Combined measurement of perfusion, venous oxygen saturation, and skeletal muscle T2* during reactive hyperemia in the leg
© Englund et al.; licensee BioMed Central Ltd. 2013
Received: 20 April 2013
Accepted: 30 July 2013
Published: 19 August 2013
The function of the peripheral microvascular may be interrogated by measuring perfusion, tissue oxygen concentration, or venous oxygen saturation (SvO2) recovery dynamics following induced ischemia. The purpose of this work is to develop and evaluate a magnetic resonance (MR) technique for simultaneous measurement of perfusion, SvO2, and skeletal muscle T2*.
Perfusion, Intravascular Venous Oxygen saturation, and T2* (PIVOT) is comprised of interleaved pulsed arterial spin labeling (PASL) and multi-echo gradient-recalled echo (GRE) sequences. During the PASL post-labeling delay, images are acquired with a multi-echo GRE to quantify SvO2 and T2* at a downstream slice location. Thus time-courses of perfusion, SvO2, and T2* are quantified simultaneously within a single scan. The new sequence was compared to separately measured PASL or multi-echo GRE data during reactive hyperemia in five young healthy subjects. To explore the impairment present in peripheral artery disease patients, five patients were evaluated with PIVOT.
Comparison of PIVOT-derived data to the standard techniques shows that there was no significant bias in any of the time-course-derived metrics. Preliminary data show that PAD patients exhibited alterations in perfusion, SvO2, and T2* time-courses compared to young healthy subjects.
Simultaneous quantification of perfusion, SvO2, and T2* is possible with PIVOT. Kinetics of perfusion, SvO2, and T2* during reactive hyperemia may help to provide insight into the function of the peripheral microvasculature in patients with PAD.
Peripheral artery disease (PAD), a common manifestation of atherosclerosis in the lower limbs, causes significant morbidity and mortality in the United States [1–4]. Atherosclerotic plaques tend to develop at branch points in the peripheral arteries, increasing vascular resistance and limiting blood flow in the affected arteries . Baseline blood flow to skeletal muscle is generally maintained through the recruitment of collateral arteries [5, 6], however the vasculature is unable to quickly respond to changes in metabolic demand, such as those that occur with exercise or following a period of ischemia. Analogous to cardiac stress testing, the functional integrity of the peripheral vasculature can be interrogated by measuring the dynamic response to a period of induced ischemia.
Skeletal muscle can be stressed by means of an ischemia-reperfusion paradigm, which induces reactive hyperemia. Suprasystolic pressure applied by a cuff secured around the thigh halts distal blood flow for several minutes. Oxygen extraction in the stagnant capillary blood continues until it reaches a steady state, approximately 140 seconds into the period of arterial occlusion . When pressure in the cuff is released, reactive hyperemia ensues with a surge in arterial flow resulting in an increase in perfusion  and oxygen concentration at the level of the capillary bed [7, 9]. Additionally, venous oxygen saturation (SvO2) within the large draining veins sharply decreases as the desaturated blood formerly trapped in the capillary bed enters. As the tissue oxygenation recovers, SvO2 rises and eventually surpasses its baseline value during the time that oxygen delivery exceeds the oxygen extraction rate , or during which time physiologic shunting may be occurring.
Several magnetic resonance (MR) techniques can non-invasively evaluate the hyperemic response. Arterial spin labeling (ASL) is a well-known method to investigate perfusion in many vascular territories [11–16]. Skeletal muscle perfusion is an important parameter that quantifies microvascular blood flow thereby providing information on delivery of oxygen and nutrients to tissue. It has been shown that perfusion dynamics during reactive hyperemia are altered in PAD . These findings correlate with both disease presence and severity.
A relative measure of tissue oxygenation can be obtained by measuring changes in the apparent transverse relaxation rate (T2*), also known as ‘blood oxygen-level dependent’ (BOLD) signal . T2* can be measured using a multi-echo gradient-recalled echo (GRE) sequence. BOLD imaging has been extensively applied for functional activation studies in the brain , and can also provide information regarding activation and oxygenation of many other tissues including the kidneys  and skeletal muscle [7, 8, 21–23].
Dynamic measurement of SvO2 provides information about oxygen utilization in tissue. When continuously measured throughout an ischemia-reperfusion paradigm, intravascular blood can act as an endogenous tracer, allowing SvO2 time-course kinetics to provide information on endothelial function and vascular reactivity . SvO2 can be dynamically measured using MR susceptometry . Each of the aforementioned MR methods will be discussed in detail in the Theory section.
The kinetics of the hyperemic response can provide information on microvascular integrity and endothelial function. Healthy subjects are able to rapidly respond to increases in oxygen demand, and recover back to baseline more quickly than patients with PAD [17, 25, 26]. While methods for simultaneous quantification of perfusion and T2* or other markers of tissue oxygenation have been implemented [27, 28], none have been able to investigate SvO2 time-course kinetics as well. Recovery dynamics are altered in each of these parameters in situations of impaired vascular function, as in PAD; therefore there is potentially added benefit to concurrent measurement. In this work, we developed a method to measure perfusion, SvO2, and T2* simultaneously. Such a technique allows for a full functional assessment of the peripheral vasculature during a single scan providing information on the temporal relationships between these various functional parameters.
Pulsed arterial spin labeling perfusion imaging
Pulsed arterial spin labeling (PASL) MRI is a well-established method for noninvasive perfusion imaging [13, 14, 29]. In one PASL variant developed by Raynaud, et al. , termed Saturation Inversion Recovery (SATIR), slice-selective (SS) and non-selective (NS) inversion pulses are applied for the tag and control images, respectively. Every acquisition is followed by a slice-selective saturation pulse to ensure the same initial magnetization in each dynamic image.
where MSS and MNS are the signal intensities (SI) in the image acquired after SS and NS inversion, respectively, λ is the tissue partition coefficient (0.9 mL/g), and T is the post-labeling delay (PLD), defined as the time between the inversion pulse and readout. At 3T, T1blood = 1664 ms  and T1tissue = 1420 ms . A numerical calculation based on actual values for blood and tissue T1 shows that for perfusion values up to 100 mL/min/100g, using T1tissue ≈ T1blood = 1420 ms, underestimates perfusion by less than 5%.
SATIR has been successfully applied for perfusion imaging in the leg in healthy individuals , and the sequence has been employed to simultaneously measure BOLD [28, 32], however it has not been integrated with a multi-echo GRE for venous oxygen saturation and BOLD measurements.
MR susceptometry-based dynamic oximetry
where Δχ do represents the susceptibility difference between fully oxygenated and fully deoxygenated blood (Δχ do = 4π•0.27 ppm (SI units) [35, 36]), B0 is the main magnetic field strength, and θ is the angle of the vessel with respect to B0[24, 33]. It has been shown that for small angles, the induced field outside of the vessel is approximately homogeneous and independent of the susceptibility difference between the vein and tissue . Therefore SvO2 can be measured using a field-mapping sequence, such as a multi-echo GRE. With this technique, one can quickly, directly, and noninvasively quantify intravascular SvO2 at high-temporal resolution.
Skeletal muscle BOLD
In the microvasculature, as in the large veins, the paramagnetism of deoxyhemoglobin causes inhomogeneities in the local magnetic field. Because the microvessels are so small, it is not possible to directly measure the change in phase signal as is required for susceptometry-based oximetry. However, this local field perturbation also results in intravoxel phase dispersion, which causes the MR signal to decay more rapidly . For instance, as oxygen saturation in the capillary bed decreases, the concentration of deoxygenated hemoglobin increases, resulting in an increased intravoxel phase dispersion thereby lowering the apparent transverse relaxation time, T2*. Thus changes in T2* in response to an ischemia-reperfusion paradigm can serve as a relative marker of tissue oxygenation . The BOLD signal originates not only from changes in blood oxygen level, it is also sensitive to changes in perfusion, cellular pH, vessel diameter, and vessel orientation [9, 21, 23, 38–40]. It has been postulated, however, that the BOLD signal changes primarily result from changes in the concentration of deoxyhemoglobin in the capillary bed . Quantification of T2* can be achieved by fitting signal intensity data from a multi-echo GRE to a monoexponential function. Prior studies have shown the utility of investigating dynamic skeletal muscle BOLD during exercise [41–43], ischemia , reactive hyperemia , and in disease states [26, 44].
Simultaneous measurement of perfusion, SvO2, and T2*
In PIVOT, as in SATIR, tag and control conditions for perfusion imaging are achieved using SS and NS inversion pulses, respectively. During the PLD, a keyhole  multi-echo GRE sequence acquires data at a distal slice for SvO2 and T2* quantification. The distal location was chosen to ensure that the multi-echo GRE interleave does not impact the signal from previously labeled perfusing blood. Because the NS inversion affects both the PASL and multi-echo GRE slices, only multi-echo GRE data acquired following SS inversion are analyzed, though the interleave is run every PLD to control for magnetization transfer effects. The multi-echo GRE is immediately followed by a GRE-EPI readout at isocenter to capture data for perfusion quantification.
The University of Pennsylvania’s Institutional Review Board approved all imaging procedures, and each subject provided informed consent prior to his or her participation. To evaluate PIVOT compared to the standard measurement methods, five young healthy male subjects (27 ± 2 years old) were recruited and imaged on two separate occasions (Visit 1 and Visit 2). Four ten-minute scans were run in both sessions, each scan consisting of one minute baseline, three minute arterial occlusion, and six minute recovery. PIVOT, a repeat of the PIVOT scan (PIVOT Repeat), an otherwise identical PASL-only sequence, or otherwise identical multi-echo GRE-only sequence were run in a randomized order. To ensure the PASL interleave did not impact quantification of SvO2 and T2*, dynamic SvO2 and T2* results obtained with PIVOT were compared to the multi-echo GRE-only derived SvO2 and T2* results. Similarly, to ensure the multi-echo GRE interleave did not confound perfusion, the perfusion results obtained with PIVOT were compared to the PASL-only derived perfusion data. PIVOT was repeated to provide data on the intra-session variability of these parameters.
In addition, five PAD patients (67.2 ± 6.8 years old, ankle-brachial index (ABI) = 0.61 ± 0.14, 3 male) were drawn from an ongoing study and PIVOT imaging was performed during a single ischemia-reperfusion paradigm. For experiments in PAD patients, the total scan time was 12 minutes, with 2 minutes of baseline, 5 minutes of arterial occlusion, and 6 minutes of recovery. Since repeated arterial occlusions were not performed in PAD patients, a longer ischemic duration was used to ensure a maximal hyperemic response.
PIVOT, PASL-only, and multi-echo GRE-only sequences were written in SequenceTree  and exported for use on a 3T scanner (Siemens Medical Equipment; Erlanger, Germany). Each subject was positioned with the maximum girth of the calf centered in an 8-channel transmit/receive knee coil (Invivo, Inc; Gainesville, FL). For proximal arterial occlusion, a cuff was secured around the thigh and was rapidly inflated to 75 mmHg above the systolic pressure using the Hokanson E20 AG101 Rapid Cuff Inflation System (D. E. Hokanson, Inc; Bellevue, WA). Perfusion images were acquired with a partial Fourier GRE-EPI readout with the following parameters: FOV = 250 × 250 mm2; acquired matrix = 80×50, reconstructed to 80 × 80; slice thickness = 1 cm; slice location = isocenter; TR/TE = 1 s/8.05 ms; PLD = 952 ms. The keyhole multi-echo GRE used the following parameters: FOV = 96 × 96 mm2; keyhole acquired matrix = 96 × 24, (for SvO2 data analysis, reconstructed matrix = 96 × 96 using a fully sampled reference image obtained immediately after the dynamic PIVOT or multi-echo GRE acquisition; only dynamic data were used for T2* analysis); slice thickness = 1 cm; slice location = 3 cm inferior from isocenter; TR/TE1/TE2/TE3/TE4/TE5 = 38.12/3.78/6.99/12.32/19.32/26.32 ms. Perfusion, SvO2, and T2* each were quantified with two-second temporal resolution.
Oximetry: Dynamic SvO2 images were reconstructed to a matrix size of 96 × 96 using outer k-space data from a fully sampled reference image acquired immediately after the dynamic scan . A phase difference image was generated for each of the dynamic time-points and the low spatial-frequency phase modulation was removed as described in . An ROI was prescribed in the larger of the peroneal veins, and reference tissue was selected in an ROI immediately surrounding the peroneal vein. The phase accumulation was calculated from echoes at TE1 and TE2, with ΔTE = 3.21 ms and the difference between the intravascular and extravascular phase accumulation (∆φ) was computed. SvO2 was calculated from equation 2 . Hematocrit of 0.45 was assumed for the healthy subjects, and in the PAD patients hematocrit was measured by blood draw. Washout time (time at which minimum SvO2 occurs), upslope, representing the rate of venous resaturation (maximum slope during recovery), and overshoot (peak SvO2 – baseline SvO2) were recorded (Figure 2b).
T 2 *: T2* analysis was performed on keyhole-only data since high spatial resolution is not as critical (acquired matrix = reconstructed matrix = 96 × 24). T2* was calculated by fitting a monoexponential function to magnitude SI from echoes TE2-TE5. Even though TE1 should have highest SNR, TE1 was not included in the monoexponential fitting because large switching gradients just prior to TE1 induced significant eddy current effects that would potentially confound T2* quantification. For BOLD analysis, average SI in an ROI prescribed in the soleus muscle for each of the four echoes was fitted to a mono-exponential function to determine T2* at each time-point. T2* values were normalized to the baseline average and relative T2*min, relative T2*max, and time to peak (TTPT2*) were determined (Figure 2c).
The average and standard deviation of every time-point across all subjects and both sessions was calculated for PIVOT, PASL, and multi-echo GRE time-course data. Pearson’s correlation coefficient was calculated to compare the time-courses measured with PIVOT and the standard methods.
For each key time-course-derived parameter, Wilcoxon signed-rank tests were used to assess whether statistically significant differences exist between PIVOT and the standard method. Specifically, PIVOT-derived perfusion parameters were compared to results obtained with PASL-only, and PIVOT-derived SvO2 or T2* parameters were compared to results from the multi-echo GRE-only scan. Wilcoxon signed-rank tests were also used to determine whether significant differences exist between key parameters measured with PIVOT and PIVOT Repeat. Wilcoxon signed-rank tests were used in lieu of a standard paired Student’s t-test as only five subjects were enrolled in the evaluation study and thus it cannot be assumed that the data are normally distributed. Statistical significance was set at p < 0.05.
To assess both intra-session and inter-session repeatability, the average within-subject coefficient of variation (CV) was calculated. Specifically, to calculate intra-session repeatability, for each parameter the within-subject standard deviation across PIVOT and PIVOT Repeat from Visit 1 was averaged across subjects, then divided by the between-subject mean parameter value from Visit 1. The same analysis was performed for intra-session repeatability on data acquired during Visit 2. Similarly, to calculate inter-session repeatability, the within-subject standard deviation across all PIVOT scans from Visit 1 and Visit 2 (4 measurements per subject) was averaged across all subjects and divided by the between-subject mean parameter value from both visits.
In the PAD patients, all time-course parameters described above were calculated, but no statistical analyses were performed since patients included in this preliminary study have varying disease severity and the ischemia-reperfusion paradigm was slightly different. The purpose of including PAD patient data was for proof of principle and to explore differences that exist between patients and healthy subjects.
PIVOT evaluation in young healthy subjects
Means and standard deviations (in parentheses) of key time-course metrics measured with PIVOT and standard measurement methods in five young healthy subjects on two separate occasions
Perfusion (Standard Method = PASL)
Hyperemic duration (s)
SvO 2 (Standard Method = multi-echo GRE)
Washout time (s)
T 2 * (Standard Method = multi-echo GRE)
Baseline T2* (ms)
Summary of intra-session and inter-session repeatability for all time-course-derived metrics is presented as the within-subject coefficient of variation
Visit 1 vs. 2
Hyperemic duration (s)
Washout time (s)
T 2 *
Baseline T2* (ms)
PIVOT in PAD patients
PIVOT results in individual PAD patients and for the average of all young healthy subjects
Hyperemic duration (s)
Washout time (s)
T 2 *
Baseline T2* (ms)
The repeatability assessment comparing perfusion, SvO2, and T2* metrics derived from the two successive PIVOT scans show some intra-session variability. This variability could be physiologic in nature or could be due to noise in the time-course data. Given the low average baseline signal standard deviation (perfusion: 1.6 mL/min/100g; SvO2: 1.5 %HbO2; T2*: 0.5%) it is likely that intra-session variations of key measured parameters are outweighed by physiologic variability during separate ischemia-reperfusion episodes. In all cases there were no significant differences detected between time-course parameters measured with PIVOT and PIVOT Repeat (p>0.05), suggesting that there was no training effect due to multiple periods of ischemia. This finding justifies the comparison of within-session PIVOT to PASL-only or to multi-echo GRE-only sequences, even though the data were acquired separately.
PIVOT effect on perfusion quantification
In comparing perfusion metrics measured with PIVOT to those measured with an otherwise identical PASL-only sequence, it is evident that the two measurements provide similar values of all perfusion time-course-derived metrics. Specifically, there were no significant differences between any of the key parameters measured with PIVOT or with PASL. This suggests that the multi-echo GRE interleave does not impact the quantification of perfusion.
Perfusion time-course metrics had expected results for longitudinal reproducibility. A CV of approximately 20% has been reported in previous studies [15, 49]. Perfusion varies physiologically with time of day , hydration level , and hormonal fluctuations  among other factors . Care was taken to schedule Visits 1 and 2 at the same time, and all subjects were instructed to refrain from caffeine intake and vigorous activity for 12 hours prior. Yet even within a single scan there was variability in the hyperemic response from one ischemia-reperfusion episode to another, suggesting that these physiologic factors cannot be completely controlled. The time-course data averaged over all experiments and all healthy subjects showed good correlation in the shape and magnitude of the reactive hyperemia response, reflected by the high correlation coefficients. Thus not only were the key time-course-derived metrics not significantly different, but that the overall response was highly similar.
Peak hyperemic flow measured with either PIVOT or PASL was somewhat lower than previously reported values, where PHF was measured to be 50 ± 13 mL/min/100 g using the same SATIR perfusion preparation, but with a RARE readout . PHF in the soleus in our study was lower, reaching only 37.0 ± 6.1 mL/min/100 g. This difference in perfusion could be attributed to the fact that Raynaud et al. quantified whole-leg perfusion , while individual muscle perfusion was calculated here. Additionally, care was taken to exclude vessels from the ROI, as their inclusion would increase measured perfusion. In another study by Proctor and colleagues, perfusion was measured in the calf of 64 men using strain gauge plethysmography , a technique considered to be a standard for limb perfusion measurement. Peak perfusion in the calf was found to be 35.1 ± 1.1 mL/min/100 g following a period of 10 minutes of ischemia. This work and that by Raynaud  and Proctor  report much lower PHF than a similar study by Wu et al. using continuous arterial spin labeling (CASL). In Wu’s study, PHF in the soleus muscle was found to be 116 ± 57 mL/min/100 g , however the temporal resolution in CASL was limited to 16 seconds. Using an EPI readout, the magnitude signal over the course of an ischemia reperfusion paradigm varies substantially due to the BOLD effect. If the time course is not sampled with high enough temporal resolution, the changing signal intensity due to the BOLD effect may contaminate the quantification of perfusion . The higher temporal resolution data acquired with SATIR can track the changes in signal intensity better and thus may be less susceptible to BOLD contamination, allowing more accurate quantification of perfusion.
Time to peak perfusion agreed with Raynaud’s values , however it was much shorter than the TTP reported by Wu . Again, this discrepancy between our work and that of Wu et al. could be due to the better temporal resolution of SATIR over CASL. The temporal resolution for PIVOT and SATIR was 2 seconds, while that of the CASL sequence employed in Wu’s work was 16 s. The improvement in temporal resolution was in part due to the pulsed tagging scheme used in SATIR, in which arterial tagging takes only 8 ms instead of 2 seconds as in CASL.
Hyperemic flow volume quantified the total blood delivery during hyperemia, and along with the hyperemic duration may provide more qualitative measures of hyperemia. HFV in particular may be less sensitive to time-course noise compared to PHF and TTP. PHF and TTP are determined based on a single data-point, whereas HFV is the total integrated area. The values we reported for HFV and hyperemic duration were lower and shorter than those reported by Wu, et al. , which is not surprising since our measured PHF and TTP are lower and shorter as well.
PIVOT effect on SvO2 quantification
High spatial resolution is necessary to measure the phase in the vein; therefore the keyhole multi-echo GRE data was supplemented with outer k-space data from a fully-phase-encoded reference scan that was run at the end of the dynamic acquisition. Langham et al. have previously shown that keyhole reconstruction provides accurate SvO2 results with high temporal resolution during reactive hyperemia in the femoral vein .
No significant differences were detected between time-course metrics measured with PIVOT or with the multi-echo GRE. While inter-session variability was present, it was comparable to the intra-session variability measured between PIVOT and PIVOT Repeat on the same visit. Average washout time was slightly lower than previously reported in the femoral vein of young healthy subjects (17 ± 7 s), but upslope and overshoot were in agreement with prior results . The lower washout time could be explained by the fact that the tourniquet system used in this study deflated much more quickly than that used in . The resulting decrease in the resistance to arterial flow could potentially shorten the washout time. Another potential reason that washout time differs is that we investigated SvO2 in the peroneal vein, as opposed to the more superior femoral vein. Thus not only does the peroneal vein collect from a smaller volume of muscle, but the distance between the capillary bed and peroneal vein was also smaller. Overall, key parameters measured with PIVOT agree with multi-echo GRE-derived data, suggesting that inclusion of the PASL interleave does not impact quantification of SvO2.
PIVOT effect on T2* quantification
In both PIVOT and multi-echo GRE T2* data, the signal intensity in the ROI was first averaged then fit to a mono-exponential function. The fitting of average signal intensity was nearly perfect with an overall average R2 for all T2* fits of 0.999.
During the period of arterial occlusion, relative T2* was found to decrease as the deoxyhemoglobin concentration in the capillary bed increases. Following cuff release, hyperemic arterial inflow replenishes blood in the capillary bed, bringing oxygenated arterial blood in and moving desaturated blood into the large draining veins. Relative T2* increased during reactive hyperemia due to the increase in perfusion and the decrease in deoxyhemoglobin in the capillary bed. Even though perfusion and T2* were measured at separate slice locations, they both represent changes that occur in the soleus muscle. TTPT2* was longer than perfusion TTP, suggesting that the increase in oxygen concentration at the level of the capillary bed lasted longer than the increase in microvascular flow. This finding is in agreement with a prior study investigating the combined measurement of perfusion and BOLD during reactive hyperemia . Duteil et al. suggest that because the brief period of arterial occlusion does not cause significant oxygen debt in muscle (as shown by ), the muscle’s demand for oxygen remains unchanged and thus the increase in perfusion results in a decrease of oxygen extraction, causing the BOLD signal to increase . The decrease in oxygen extraction also physiologically manifests as the SvO2 overshoot.
The measured values are consistent with previous literature reported values for baseline T2* , relative T2*min, T2*max[26, 59], and TTPT2*. Comparisons between PIVOT and multi-echo GRE key time-course parameters yielded no significant differences, and the T2* time-course data measured with PIVOT and multi-echo GRE were highly correlated (Pearson’s r = 0.99). These results indicate that T2* quantification in a downstream slice with PIVOT is not affected by the PASL interleave.
Considerations for applying PIVOT in PAD patient studies
Statistical comparisons between PAD patients and young healthy subjects were not made because the ischemic duration differed in the two cohorts. As repeated periods of arterial occlusion were used to compare PIVOT to standard measurement methods in healthy subjects, a shortened ischemic duration was used, allowing enough time for the four scans to be performed during one hour of scanning. PAD patients included in this preliminary evaluation had varying disease severity, represented by the diverse ABIs. A five-minute period of arterial occlusion was employed in PAD patients to ensure the maximal hyperemic stimulus. The PAD patient data were, however, included to explore the range of values that exist between states of health and disease. In agreement with the results of Wu et al., PAD patients exhibited a decrease in PHF and an increase in TTP . SvO2 data showed a blunted and delayed response, in agreement with Langham et al. [10, 25]. T2* data showed the characteristic alterations during the period of ischemia  and reactive hyperemia .
Quantification of multiple parameters may improve diagnosis and enhance power for detection of the response to therapeutic intervention. The traditional marker of disease severity in PAD is the ankle-brachial index (ABI), which is the ratio of systolic blood pressures measured at the level of the ankle and in the brachial artery. The ABI primarily represents occlusions and stenoses on the macrovascular level. Previous studies have shown that physiologic improvements such as increased peak walking time do not correlate with clinically significant changes in ABI [60, 61]. PIVOT provides a measure of microvascular function, thus may be more sensitive to early treatment effects. By measuring many parameters, PIVOT will provide insight into the relationship between impairments in perfusion, SvO2, and T2*.
In this study we showed that no major measurement bias in key time-course parameters was introduced by using PIVOT instead of the standard individual measurement methods. However, it should be noted that the small sample size affects the power to detect such a bias. Even though no measurement bias was detected, the precision and thus statistical power of PIVOT-derived measures was limited by physiologic variability, which cannot be completely controlled. In order to use this method to assess disease presence or monitor a treatment effect, the study must be well designed, controlling for factors that are known to affect perfusion and the hyperemic response [50–53].
The preliminary PAD patient data showed that for many of the key time-course parameters there was a relatively wide range of values that exist between states of health and disease. A recent study performed at 1.5T showed relatively poor repeatability of BOLD measurements with the exception of TTPT2* during reactive hyperemia in both healthy subjects and PAD patients . In our study a lower inter-session CV, and hence better repeatability, was found for several key T2* parameters in healthy subjects. The improvement in repeatability may be due to the higher field strength, which confers increased signal to noise ratio and greater BOLD signal contrast. A longitudinal study will be necessary to determine the repeatability of PIVOT measures in PAD patients. Additionally vascular reactivity decreases with age , therefore it will be important to compare PIVOT results obtained for PAD patients to age-matched healthy controls to assess the impact of age on PIVOT-derived measures of microvascular function.
In summary, we have introduced a quantitative MR method that measures perfusion, SvO2, and T2* simultaneously, thereby allowing a comprehensive assessment of the functional integrity of the peripheral microvasculature during a single ischemia-reperfusion paradigm. The added value of the proposed approach will require rigorous evaluation in cohorts of patients with impaired peripheral circulation in comparison to their healthy peers. In the future PIVOT could possibly serve as a means to monitor disease progression and effectiveness of intervention.
We would like to acknowledge Dr. Kathleen Propert, who provided assistance in the statistical analysis of the results. This research was supported by an award from the American Heart Association (EKE) and NIH grants R01 HL075649, R01 HL109545, and K25 HL111422.
- Criqui MH, Fronek A, Barrett-Connor E, Klauber MR, Gabriel S, Goodman D: The prevalence of peripheral arterial disease in a defined population. Circulation. 1985, 71: 510-5. 10.1161/01.CIR.71.3.510.View ArticlePubMedGoogle Scholar
- Meijer WT, Hoes AW, Rutgers D, Bots ML, Hofman A, Grobbee DE: Peripheral arterial disease in the elderly: the Rotterdam study. Arterioscler Thromb Vasc Biol. 1998, 18: 185-92. 10.1161/01.ATV.18.2.185.View ArticlePubMedGoogle Scholar
- Hirsch AT HCM, Treat-Jacobson D, Regensteiner JG, Creager MA, Olin JW, Krook SH, Hunninghake DB, Comerota AJ, Walsh ME: Peripheral arterial disease detection, awareness, and treatment in primary care. JAMA: J Am Med Assoc. 2001, 286: 1317-24. 10.1001/jama.286.11.1317.View ArticleGoogle Scholar
- Hirsch AT, Haskal ZJ, Hertzer NR, Bakal CW, Creager MA, Halperin JL, Hiratzka LF, Murphy WRC, Olin JW, Puschett JB, Rosenfield KA, Sacks D, Stanley JC, Taylor LM, White CJ, White J, White RA, Antman EM, Smith SC, Adams CD, Anderson JL, Faxon DP, Fuster V, Gibbons RJ, Hunt SA, Jacobs AK, Nishimura R, Ornato JP, Page RL, Riegel B: ACC/AHA 2005 Practice Guidelines for the management of patients with peripheral arterial disease (lower extremity, renal, mesenteric, and abdominal aortic): a collaborative report from the American Association for Vascular Surgery/Society for Vascular Surgery, Society for Cardiovascular Angiography and Interventions, Society for Vascular Medicine and Biology, Society of Interventional Radiology, and the ACC/AHA Task Force on Practice Guidelines (Writing Committee to Develop Guidelines for the Management of Patients With Peripheral Arterial Disease): endorsed by the American Association of Cardiovascular and Pulmonary Rehabilitation; National Heart, Lung, and Blood Institute; Society for Vascular Nursing; TransAtlantic Inter-Society Consensus; and Vascular Disease Foundation. Circulation. 2005, 2006: e463-654.Google Scholar
- Faxon DP, Fuster V, Libby P, Beckman JA, Hiatt WR, Thompson RW, Topper JN, Annex BH, Rundback JH, Fabunmi RP, Robertson RM, Loscalzo J, American Heart Association: Atherosclerotic vascular disease conference: writing group III: pathophysiology. Circulation. 2004, 109: 2617-25. 10.1161/01.CIR.0000128520.37674.EF.View ArticlePubMedGoogle Scholar
- Bragadeesh T, Sari I, Pascotto M, Micari A, Kaul S, Lindner JR: Detection of peripheral vascular stenosis by assessing skeletal muscle flow reserve. J Am Coll Cardiol. 2005, 45: 780-5. 10.1016/j.jacc.2004.11.045.View ArticlePubMedGoogle Scholar
- Ledermann HP, Heidecker H-G, Schulte A-C, Thalhammer C, Aschwanden M, Jaeger KA, Scheffler K, Bilecen D: Calf muscles imaged at BOLD MR: correlation with TcPO2 and flowmetry measurements during ischemia and reactive hyperemia–initial experience. Radiology. 2006, 241: 477-84. 10.1148/radiol.2412050701.View ArticlePubMedGoogle Scholar
- Toussaint J-F, Kwong KK, M'Kparu F, Weisskoff RM, LaRaia PJ, Kantor HL: Perfusion changes in human skeletal muscle during reactive hyperemia measured by echo-planar imaging. Magn Reson Med. 1996, 35: 62-9. 10.1002/mrm.1910350109.View ArticlePubMedGoogle Scholar
- Damon BM, Hornberger JL, Wadington MC, Lansdown DA, Kent-Braun JA: Dual gradient-echo MRI of post-contraction changes in skeletal muscle blood volume and oxygenation. Magn Reson Med. 2007, 57: 670-9. 10.1002/mrm.21191.PubMed CentralView ArticlePubMedGoogle Scholar
- Langham MC, Floyd TF, Mohler ER, Magland JF, Wehrli FW: Evaluation of cuff-induced ischemia in the lower extremity by magnetic resonance oximetry. J Am Coll Cardiol. 2010, 55: 598-606. 10.1016/j.jacc.2009.08.068.PubMed CentralView ArticlePubMedGoogle Scholar
- Detre J, Leigh J, Williams D: Detre, MRM 1992 Perfusion imaging. Magn Reson Med. 1992, 23: 1-9. 10.1002/mrm.1910230102.View ArticleGoogle Scholar
- Williams DS, Detre JA, Leigh JS, Koretsky AP: Magnetic resonance imaging of perfusion using spin inversion of arterial water. Proc Natl Acad Sci. 1992, 89: 212-6. 10.1073/pnas.89.1.212.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim SG: Quantification of relative cerebral blood flow change by flow-sensitive alternating inversion recovery (FAIR) technique: application to functional mapping. Magn Reson Med. 1995, 34: 293-301. 10.1002/mrm.1910340303.View ArticlePubMedGoogle Scholar
- Kim SG, Tsekos NV: Perfusion imaging by a flow-sensitive alternating inversion recovery (FAIR) technique: application to functional brain imaging. Magn Reson Med. 1997, 37: 425-35. 10.1002/mrm.1910370321.View ArticlePubMedGoogle Scholar
- Wu W-C, Wang J, Detre JA, Wehrli FW, Mohler E, Ratcliffe SJ, Floyd TF: Hyperemic flow heterogeneity within the calf, foot, and forearm measured with continuous arterial spin labeling MRI. Am J Physiol Heart Circ Physiol. 2008, 294: H2129-36. 10.1152/ajpheart.01399.2007.PubMed CentralView ArticlePubMedGoogle Scholar
- Raynaud JS, Duteil S, Vaughan JT, Hennel F, Wary C, Leroy-Willig A, Carlier PG: Determination of skeletal muscle perfusion using arterial spin labeling NMRI: validation by comparison with venous occlusion plethysmography. Magn Reson Med. 2001, 46: 305-11. 10.1002/mrm.1192.View ArticlePubMedGoogle Scholar
- Wu W-C, Mohler E, Ratcliffe SJ, Wehrli FW, Detre JA, Floyd TF: Skeletal muscle microvascular flow in progressive peripheral artery disease: assessment with continuous arterial spin-labeling perfusion magnetic resonance imaging. J Am Coll Cardiol. 2009, 53: 2372-7. 10.1016/j.jacc.2009.03.033.View ArticlePubMedGoogle Scholar
- Ogawa S, Lee T, Kay A: Brain magnetic resonance imaging with contrast dependent on blood oxygenation. Proc Natl Acad Sci. 1990, 87: 9868-72. 10.1073/pnas.87.24.9868.PubMed CentralView ArticlePubMedGoogle Scholar
- Ogawa S, Menon RS, Tank DW, Kim SG, Merkle H, Ellermann JM, Uğurbil K: Functional brain mapping by blood oxygenation level-dependent contrast magnetic resonance imaging. A comparison of signal characteristics with a biophysical model. Biophys J. 1993, 64: 803-12. 10.1016/S0006-3495(93)81441-3.PubMed CentralView ArticlePubMedGoogle Scholar
- Prasad PV, Edelman RR, Epstein FH: Noninvasive evaluation of intrarenal oxygenation with BOLD MRI. Circulation. 1996, 94: 3271-5. 10.1161/01.CIR.94.12.3271.View ArticlePubMedGoogle Scholar
- Noseworthy M, Bulte DP, Alfonsi J: BOLD magnetic resonance imaging in skeletal muscle. Semin Musculoskelet Radiol. 2003, 7: 307-15.View ArticlePubMedGoogle Scholar
- Damon BM, Gore JC: Physiological basis of muscle functional MRI: predictions using a computer model. J Appl Physiol. 2004, 98: 264-73. 10.1152/japplphysiol.00369.2004.View ArticlePubMedGoogle Scholar
- Lebon V, Brillault-Salvat C, Bloch G, Leroy-Willig A, Carlier PG: Evidence of muscle BOLD effect revealed by simultaneous interleaved gradient-echo NMRI and myoglobin NMRS during leg ischemia. Magn Reson Med. 1998, 40: 551-8. 10.1002/mrm.1910400408.View ArticlePubMedGoogle Scholar
- Fernández-Seara MA, Techawiboonwong A, Detre JA, Wehrli FW: MR susceptometry for measuring global brain oxygen extraction. Magn Reson Med. 2006, 55: 967-73. 10.1002/mrm.20892.View ArticlePubMedGoogle Scholar
- Langham MC, Englund EK, Mohler ER, Li C, Rodgers ZB, Floyd TF, Wehrli FW: Quantitative CMR markers of impaired vascular reactivity associated with age and peripheral artery disease. J Cardiov Magn Reson. 2013, 15: 17-10.1186/1532-429X-15-17.View ArticleGoogle Scholar
- Ledermann HP, Schulte A-C, Heidecker H-G, Aschwanden M, Jäger KA, Scheffler K, Steinbrich W, Bilecen D: Blood oxygenation level-dependent magnetic resonance imaging of the skeletal muscle in patients with peripheral arterial occlusive disease. Circulation. 2006, 113: 2929-35. 10.1161/CIRCULATIONAHA.105.605717.View ArticlePubMedGoogle Scholar
- Lebon V, Carlier PG, Brillault-Salvat C, Leroy-Willig A: Simultaneous measurement of perfusion and oxygenation changes using a multiple gradient-echo sequence: application to human muscle study. Magn Reson Imaging. 1998, 16: 721-9. 10.1016/S0730-725X(98)00088-5.View ArticlePubMedGoogle Scholar
- Duteil S, Wary C, Raynaud JS, Lebon V, Lesage D, Leroy-Willig A, Carlier PG: Influence of vascular filling and perfusion on BOLD contrast during reactive hyperemia in human skeletal muscle. Magn Reson Med. 2006, 55: 450-4. 10.1002/mrm.20760.View ArticlePubMedGoogle Scholar
- Wong EC, Buxton RB, Frank LR: Implementation of quantitative perfusion imaging techniques for functional brain mapping using pulsed arterial spin labeling. NMR Biomed. 1997, 10: 237-49. 10.1002/(SICI)1099-1492(199706/08)10:4/5<237::AID-NBM475>3.0.CO;2-X.View ArticlePubMedGoogle Scholar
- Lu H, Clingman C, Golay X, van Zijl PCM: Determining the longitudinal relaxation time (T1) of blood at 3.0 Tesla. Magn Reson Med. 2004, 52: 679-82. 10.1002/mrm.20178.View ArticlePubMedGoogle Scholar
- Gold GE, Han E, Stainsby J, Wright G, Brittain J, Beaulieu C: Musculoskeletal MRI at 3.0 T: relaxation times and image contrast. Am J Roentgenol. 2004, 183: 343-51. 10.2214/ajr.183.2.1830343.View ArticleGoogle Scholar
- Carlier PG, Bertoldi D, Baligand C, Wary C, Fromes Y: Muscle blood flow and oxygenation measured by NMR imaging and spectroscopy. NMR Biomed. 2006, 19: 954-67. 10.1002/nbm.1081.View ArticlePubMedGoogle Scholar
- Haacke EM, Lai S, Reichenbach JR, Kuppusamy K, Hoogenraad FG, Takeichi H, Lin W: In vivo measurement of blood oxygen saturation using magnetic resonance imaging: a direct validation of the blood oxygen level-dependent concept in functional brain imaging. Hum Brain Mapp. 1997, 5: 341-6. 10.1002/(SICI)1097-0193(1997)5:5<341::AID-HBM2>3.0.CO;2-3.View ArticlePubMedGoogle Scholar
- Langham MC, Magland JF, Epstein CL, Floyd TF, Wehrli FW: Accuracy and precision of MR blood oximetry based on the long paramagnetic cylinder approximation of large vessels. Magn Reson Med. 2009, 62: 333-40. 10.1002/mrm.21981.PubMed CentralView ArticlePubMedGoogle Scholar
- Spees WM, Yablonskiy DA, Oswood MC, Ackerman JJ: Water proton MR properties of human blood at 1.5 Tesla: magnetic susceptibility, T(1), T(2), T*(2), and non-Lorentzian signal behavior. Magn Reson Med. 2001, 45: 533-42. 10.1002/mrm.1072.View ArticlePubMedGoogle Scholar
- Jain V, Abdulmalik O, Propert KJ, Wehrli FW: Investigating the magnetic susceptibility properties of fresh human blood for noninvasive oxygen saturation quantification. Magn Reson Med. 2011, 68: 863-7.PubMed CentralView ArticlePubMedGoogle Scholar
- Li C, Langham MC, Epstein CL, Magland JF, Wu J, Gee J, Wehrli FW: Accuracy of the cylinder approximation for susceptometric measurement of intravascular oxygen saturation. Magn Reson Med. 2011, 67: 808-13.View ArticlePubMedGoogle Scholar
- Donahue KM, Van Kylen J, Guven S, El Bershawi A, Luh WM, Bandettini PA, Cox RW, Hyde JS, Kissebah AH: Simultaneous gradient‒echo/spin‒echo EPI of graded ischemia in human skeletal muscle. J Magn Reson Imaging. 1998, 8: 1106-13. 10.1002/jmri.1880080516.View ArticlePubMedGoogle Scholar
- Sanchez OA, Copenhaver EA, Elder CP, Damon BM: Absence of a significant extravascular contribution to the skeletal muscle BOLD effect at 3 T. Magn Reson Med. 2010, 64: 527-35.PubMed CentralPubMedGoogle Scholar
- Partovi S, Karimi S, Jacobi B, Schulte A-C, Aschwanden M, Zipp L, Lyo JK, Karmonik C, Müller-Eschner M, Huegli RW, Bongartz G, Bilecen D: Clinical implications of skeletal muscle blood-oxygenation-level-dependent (BOLD) MRI. Mag Reson Mater Phy. 2012, 25: 251-61. 10.1007/s10334-012-0306-y.View ArticleGoogle Scholar
- Slade JM, Towse TF, Gossain VV, Meyer RA: Peripheral microvascular response to muscle contraction is unaltered by early diabetes but decreases with age. J Appl Physiol. 2011, 111: 1361-71.PubMed CentralView ArticlePubMedGoogle Scholar
- Towse TF, Slade JM, Ambrose JA, DeLano MC, Meyer RA: Quantitative analysis of the postcontractile blood-oxygenation-level-dependent (BOLD) effect in skeletal muscle. J Appl Physiol. 2011, 111: 27-39. 10.1152/japplphysiol.01054.2009.PubMed CentralView ArticlePubMedGoogle Scholar
- Elder CP, Cook RN, Wilkens KL, Chance MA, Sanchez OA, Damon BM: A method for detecting the temporal sequence of muscle activation during cycling using MRI. J Appl Physiol. 2011, 110: 826-33. 10.1152/japplphysiol.00185.2010.PubMed CentralView ArticlePubMedGoogle Scholar
- Potthast S, Schulte A, Kos S, Aschwanden M, Bilecen D: Blood oxygenation level-dependent MRI of the skeletal muscle during ischemia in patients with peripheral arterial occlusive disease. Rofo. 2009, 181: 1157-61. 10.1055/s-0028-1109786.View ArticlePubMedGoogle Scholar
- Van Vaals JJ, Brummer ME, Thomas Dixon W, Tuithof HH, Engels H, Nelson RC, Gerety BM, Chezmar JL, Den Boer JA: “Keyhole” method for accelerating imaging of contrast agent uptake. J Magn Reson Imaging. 2005, 3: 671-5.View ArticleGoogle Scholar
- Magland J, Wehrli F: Pulse sequence programming in a dynamic visual environment. Proc Int Soc Magnet Reson Med. 2006, 3032-Google Scholar
- Wu W-C, Wang J, Detre JA, Ratcliffe SJ, Floyd TF: Transit delay and flow quantification in muscle with continuous arterial spin labeling perfusion-MRI. J Magn Reson Imaging. 2008, 28: 445-52. 10.1002/jmri.21322.PubMed CentralView ArticlePubMedGoogle Scholar
- Langham MC, Magland JF, Floyd TF, Wehrli FW: Retrospective correction for induced magnetic field inhomogeneity in measurements of large-vessel hemoglobin oxygen saturation by MR susceptometry. Magn Reson Med. 2009, 61: 626-33. 10.1002/mrm.21499.View ArticlePubMedGoogle Scholar
- Versluis B, Backes WH, van Eupen MGA, Jaspers K, Nelemans PJ, Rouwet EV, Teijink JAW, Mali WPTM, Schurink G-W, Wildberger JE, Leiner T: Magnetic resonance imaging in peripheral arterial disease: reproducibility of the assessment of morphological and functional vascular status. Invest Radiol. 2011, 46: 11-24. 10.1097/RLI.0b013e3181f2bfb8.View ArticlePubMedGoogle Scholar
- Etsuda H, Takase B, Uehata A, Kusano H, Hamabe A, Kuhara R, Akima T, Matsushima Y, Arakawa K, Satomura K, Kurita A, Ohsuzu F: Morning attenuation of endothelium-dependent, flow-mediated dilation in healthy young men: possible connection to morning peak of cardiac events?. Clin Cardiol. 1999, 22: 417-21. 10.1002/clc.4960220610.View ArticlePubMedGoogle Scholar
- Nadel ER, Fortney SM, Wenger CB: Effect of hydration state of circulatory and thermal regulations. J Appl Physiol. 1980, 49: 715-21.PubMedGoogle Scholar
- Bungum L, Kvernebo K, Oian P, Maltau JM: Laser doppler-recorded reactive hyperaemia in the forearm skin during the menstrual cycle. Brit J Obstet Gynaec. 1996, 103: 70-5. 10.1111/j.1471-0528.1996.tb09517.x.View ArticlePubMedGoogle Scholar
- Parkes LM, Rashid W, Chard DT, Tofts PS: Normal cerebral perfusion measurements using arterial spin labeling: reproducibility, stability, and age and gender effects. Magn Reson Med. 2004, 51: 736-43. 10.1002/mrm.20023.View ArticlePubMedGoogle Scholar
- Proctor DN, Le KU, Ridout SJ: Age and regional specificity of peak limb vascular conductance in men. J Appl Physiol. 2005, 98: 193-202.View ArticlePubMedGoogle Scholar
- Langham MC, Englund EK, Li C, Floyd TF, Mohler ER, Wehrli FW: Balanced tissue magnetization reduces confounding BOLD effect in post-ischemic muscle perfusion quantification. Proc Int Soc Magnet Reson Med. 2012, 2027-Google Scholar
- Langham MC, Wehrli FW: mproved temporal resolution of dynamic oximetry via keyhole acquisition for quantifying reactive hyperemia. Proc Int Soc Magnet Reson Med. 2012, 1147-Google Scholar
- Brillault-Salvat C, Giacomini E, Jouvensal L, Wary C, Bloch G, Carlier PG: Simultaneous determination of muscle perfusion and oxygenation by interleaved NMR plethysmography and deoxymyoglobin spectroscopy. NMR Biomed. 1997, 10: 315-23. 10.1002/(SICI)1099-1492(199710)10:7<315::AID-NBM489>3.0.CO;2-X.View ArticlePubMedGoogle Scholar
- Partovi S, Schulte A-C, Jacobi B, Klarhöfer M, Lumsden AB, Loebe M, Davies MG, Noon GP, Karmonik C, Zipp L, Bongartz G, Bilecen D: Blood oxygenation level-dependent (BOLD) MRI of human skeletal muscle at 1.5 and 3 T. J Magn Reson Imaging. 2012, 35: 1227-32. 10.1002/jmri.23583.View ArticlePubMedGoogle Scholar
- Schulte AC, Aschwanden M, Bilecen D: Calf muscles at blood oxygen level-dependent MR imaging: aging effects at postocclusive reactive hyperemia. Radiology. 2008, 247: 482-9. 10.1148/radiol.2472070828.View ArticlePubMedGoogle Scholar
- Mohler ER, Beebe HG, Salles-Cuhna S, Zimet R, Zhang P, Heckman J, Forbes WP: Effects of cilostazol on resting ankle pressures and exercise-induced ischemia in patients with intermittent claudication. Vasc Med. 2001, 6: 151-6. 10.1177/1358836X0100600305.View ArticlePubMedGoogle Scholar
- Murphy TP, Cutlip DE, Regensteiner JG, Mohler ER, Cohen DJ, Reynolds MR, Massaro JM, Lewis BA, Cerezo J, Oldenburg NC, Thum CC, Goldberg S, Jaff MR, Steffes MW, Comerota AJ, Ehrman J, Treat-Jacobson D, Walsh ME, Collins T, Badenhop DT, Bronas U, Hirsch AT, for the CLEVER Study Investigators: Supervised exercise versus primary stenting for claudication resulting from aortoiliac peripheral artery disease: six-month outcomes from the claudication: exercise versus endoluminal revascularization (CLEVER) study. Circulation. 2012, 125: 130-9. 10.1161/CIRCULATIONAHA.111.075770.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.