- Open Access
Evaluation of skeletal muscle microvascular perfusion of lower extremities by cardiovascular magnetic resonance arterial spin labeling, blood oxygenation level-dependent, and intravoxel incoherent motion techniques
Journal of Cardiovascular Magnetic Resonance volume 20, Article number: 18 (2018)
Noninvasive cardiovascular magnetic resonance (CMR) techniques including arterial spin labeling (ASL), blood oxygenation level-dependent (BOLD), and intravoxel incoherent motion (IVIM), are capable of measuring tissue perfusion-related parameters. We sought to evaluate and compare these three CMR techniques in characterizing skeletal muscle perfusion in lower extremities and to investigate their abilities to diagnose and assess the severity of peripheral arterial disease (PAD).
Fifteen healthy young subjects, 14 patients with PAD, and 10 age-matched healthy old subjects underwent ASL, BOLD, and IVIM CMR perfusion imaging. Healthy young and healthy old participants were subjected to a cuff-induced ischemia experiment with pressures of 20 mmHg and 40 mmHg above systolic pressure during imaging. Perfusion-related metrics, including blood flow, T2* relaxation time, perfusion fraction f, diffusion coefficient D, and pseudodiffusion coefficient D*, were measured in the anterior, lateral, soleus, and gastrocnemius muscle groups. Friedman, Mann-Whitney, Wilcoxon signed rank, and Spearman rank correlation tests were used for statistical analysis.
In cases of significant differences determined by the Friedman test (P < 0.05), blood flow, T2*, and D values gradually decreased, while f values showed a tendency to increase in healthy subjects under cuff compression. No significant correlations were found among the ASL, BOLD, and IVIM parameters (all P > 0.05). Blood flow and T2* values showed significant positive correlations with transcutaneous oxygen pressure measurements (ρ = 0.465 and 0.522, respectively; both P ≤ 0.001), while f values showed a significant negative correlation in healthy young subjects (ρ = − 0.351; P = 0.018). T2* was independent of age in every muscle group. T2* values were significantly decreased in PAD patients compared with healthy old subjects and severe PAD patients compared with mild-to-moderate PAD patients (all P < 0.0125). Significant correlations were found between T2* and ankle–brachial index values in all muscle groups in PAD patients (ρ = 0.644–0.837; all P < 0.0125). Other imaging parameters failed to show benefits towards the diagnosis and disease severity evaluation of PAD.
ASL, BOLD, and IVIM provide complementary information regarding tissue perfusion. Compared with ASL and IVIM, BOLD may be a more reliable technique for assessing PAD in the resting state and could thus be applied together with angiography in clinical studies as a tool to comprehensively assess microvascular and macrovascular properties in PAD patients.
Peripheral arterial disease (PAD) is a highly prevalent and severe atherosclerotic condition characterized by progressive peripheral arterial development of lower extremeity stenosis/occlusions . Patients affected with PAD suffer from reduced quality of life, and more importantly, increased risk of cardiovascular and cerebrovascular events . Therefore, a noninvasive and objective method is desirable for diagnostic, prognostic, and therapeutic purposes, such as early detection of physiological function changes, clinical risk stratification for predicting myocardial infarction or stroke, and intervention planning for symptomatic patients.
Noninvasive testing of flow-limiting stenosis typically includes measurement of the ankle–brachial index (ABI), the ratio of ankle systolic blood pressure to arm systolic pressure . PAD is considered to be present when the ABI is ≤0.90 and severe when the ABI is ≤0.50 . However, ABI has low sensitivity for PAD diagnosis and may not be necessarily associated with symptom relief after interventions [4, 5]. Transcutaneous oxygen pressure (TcPO2) measurement is an additional method used to indirectly assess the degree of ischemia in ischemic skeletal muscle by measuring tissue oxygenation . The use of TcPO2 measurement is limited because it is confined to the skin and thus does not accurately reflect muscle perfusion .
Medical imaging has emerged as an important tool in the diagnosis and management of PAD. Imaging modalities, including computer tomography (CT) angiography, cardiovascular magnetic resonance (CMR) angiography, and digital subtraction angiography, are commonly used to assess abnormal blood vessels and blood flow to the lower extremities. These techniques fail to provide information regarding skeletal muscle microvascular perfusion in the affected extremity . However, as PAD extends beyond the large-vessels, blood flow impairment leads to microvascular dysfunction. Precise assessment of skeletal muscle perfusion would facilitate the comprehensive evaluation of PAD and could be combined with conventional angiography to reveal both functional and anatomical characteristics.
Several CMR techniques can noninvasively measure microvascular perfusion using endogenous tracers, including arterial spin labeling (ASL), blood oxygenation level-dependent (BOLD), and intravoxel incoherent motion (IVIM). ASL magnetically tags arterial blood using radiofrequency pulses, and the perfusion contrast is given by the signal difference between the tagged image and the nontagged control image obtained without net magnetization perturbation in arterial blood . BOLD uses the paramagnetic effect of deoxygenated hemoglobin as an intrinsic contrast agent, which decreases the T2* relaxation signal . IVIM is a variant of conventional diffusion-weighted imaging by separating the effect of blood flow in the randomly oriented capillary network from that of thermally driven water molecular diffusion . ASL, BOLD, and IVIM have been successfully applied to measure skeletal muscle perfusion in previous studies [12,13,14]. Since perfusion-related metrics derived from these different CMR techniques are based on completely distinct mechanisms, they may depict different aspects of muscle perfusion properties. ASL is more related to the function of blood delivery to target tissues, BOLD is related to tissue oxygenation, and IVIM is related to pseudodiffusion within capillary beds. We hypothesize that multi-parametric CMR techniques, including ASL, BOLD, and IVIM, could provide complementary information regarding perfusion in skeletal muscles and would represent various alteration patterns in the presence of perfusion deficits.
Hence, this study aimed to 1) test the feasibility of using ASL, BOLD, and IVIM to measure perfusion changes in the lower extremities of healthy subjects under different external compression statuses; 2) validate the associations between ASL-, BOLD-, and IVIM-derived parameters and TcPO2 measurements; 3) evaluate the effects of age on imaging parameter measurements of ASL, BOLD, and IVIM in healthy subjects at rest; 4) use ASL, BOLD, and IVIM to compare perfusion in affected and contralateral (asymptomatic) lower extremities in PAD patients at rest; 5) compare the capabilities of resting-state ASL, BOLD, and IVIM in detecting perfusion differences between PAD patients and age-matched healthy subjects and between mild-to-moderate and severe PAD patients; and 6) investigate the associations between ASL-, BOLD-, and IVIM-derived parameters and ABI in PAD patients.
The local institutional review board approved the prospective study, and written consent was obtained from all subjects prior to participation. Technical support for imaging sequence optimization was provided by a GE Healthcare employee (JQ). Authors not associated with GE Healthcare had full control of the data and information submitted for publication.
Between February 2016 and October 2017, three groups of subjects were enrolled: 1) healthy young subjects (n = 15); 2) PAD patients (n = 14); and 3) age-matched healthy old subjects (n = 10). None of the healthy subjects showed clinical evidence of PAD, cardiac insufficiency, hypoxic pulmonary diseases, or lower extremity venous disorders. Each of these subjects had a normal peripheral pulse status and an ABI > 0.90, and they were all non-smokers. The healthy old group was age-matched to the PAD group to eliminate the confounding effect of age on CMR perfusion parameter measurements. Patients with PAD were recruited from the department of vascular surgery with symptoms of intermittent claudication, rest pain, or critical limb ischemia and an ABI ≤ 0.90. The PAD group was then stratified into two disease severity subgroups based on ABI: 1) the mild-to-moderate group corresponded to an ABI of 0.51 to 0.90 (n = 7); and 2) the severe group corresponded to an ABI ≤ 0.50 (n = 7). Details for each group are listed in Table 1.
In preparation for the scans, all subjects were asked to refrain from alcohol, caffeine, and vigorous exercise 12 h before imaging. The healthy young and healthy old subjects were subjected to a cuff-induced arterial occlusion experiment during the CMR imaging scans to test the sensitivities of ASL, BOLD, and IVIM to pressure variations. Ischemia via arterial occlusion was induced in the lower extremity by a sphygmomanometer cuff tied around the middle of one thigh. The contralateral lower extremity without intervention was imaged simultaneously as the control side. ASL, BOLD, and IVIM were conducted four times in the following order: baseline (Pre), cuff compression with a pressure of 20 mmHg above systolic pressure (Cuff-20), cuff compression with a pressure of 40 mmHg above systolic pressure (Cuff-40), and recovery period (Post). Two different cuff pressures were used to verify whether different air pressures could modulate the perfusion signal intensity and thus provoke different degrees of ischemia. The pressure of 40 mmHg above systolic pressure was chosen in accordance with the previously reported range of 30–50 mmHg above systolic pressure, which is recommended for provoking complete and reproducible ischemia . During the Cuff-20 and Cuff-40 sessions, the cuff was kept inflated until all scans were completed (6 min 20 s). The last three sessions were performed at 30-min intervals to avoid the interference from the preceding session (Fig. 1). To test the interscan reproducibility, 5 healthy subjects were subjected to a second scan within 1 week. Patients with PAD received baseline examinations only at rest. All subjects were asked to lie still in the supine position for approximately 15 min before the onset of CMR imaging to ensure that their legs were at heart height , and remain still during the entire examination.
All the CMR measurements were carried out on a 3-T CMR system (HDxt, General Electric Healthcare, Waukesha, Wisconsin, USA) with an eight-channel cardiac coil. The subjects were investigated in the supine position. Prior to acquisition of the CMR perfusion images, axial three-dimensional (3D) spoiled gradient-recalled echo T1-weighted images (repetition time ms/echo time ms 4.1/1.5; flip angle = 12°; matrix, 320 × 320; FOV, 32 × 32 mm2; slice thickness = 5 mm; number of slices = 12) were acquired for use as anatomical landmarks. A FOV of 32 × 32 mm2 could cover both legs without leaving anything outside, thus prohibiting wrapping artifacts . Pseudocontinuous ASL was performed using an interleaved 3D stack of spiral fast spin-echo sequence with background suppression. Each spiral arm included 512 sampling points in k-space and a total of 6 spiral arms were acquired. Background suppression was achieved via 5 inversion pulses placed 0, 1465, 2100, 2600, and 2880 ms after the labeling start point to suppress a broad range of T1 values . Other ASL parameters were as follows: repetition time ms/echo time ms = 4316/9.4; bandwidth = 62.5 kHz; FOV = 32 × 32 mm2; slice thickness = 5 mm; number of slices = 12; number of averages = 2; post-labeling delay time = 1525 ms. A post-labeling delay time of 1525 ms before imaging was used to allow the blood to perfuse all muscle groups, which was similar to that used in a previous study . BOLD was performed using a multi-echo gradient-recalled echo sequence implementing the following parameters: repetition time ms/echo times ms 875/(2.5, 6.7, 10.9, 15.1, 19.3, 23.5, 27.7, 31.9, 36.1, 40.3, 44.5, 48.7, 52.9, 57.1, 61.3, and 65.5); matrix = 256 × 128; FOV = 32 × 32 mm2; slice thickness/gap = 5/0 mm; number of slices = 12. IVIM imaging was performed using a single shot spin-echo echo-planar imaging sequence at 9 b-values (0, 20, 50, 100, 150, 200, 300, 500, and 800 s/mm2) in three orthogonal gradient directions, with the following parameters: repetition time ms/echo time ms = 2800/70; matrix = 192 × 192; FOV = 32 × 32 mm2; slice thickness/gap = 5/0 mm; number of slices = 12; number of averages = 3; parallel imaging factor = two. A standard monopolar Stejskal-Tanner diffusion encoding scheme was applied with diffusion gradient pulse duration of 16 ms. The acquisition times for ASL, BOLD, and IVIM were 2 min 18 s, 1 min 59 s, and 2 min 3 s, respectively. All imaging was performed in the axial plane at the level of the middle calf.
ASL, BOLD, and IVIM data were post-processed on a pixel-by-pixel basis on a workstation (Advantage Workstation 4.5; General Electric Healthcare) to obtain corresponding parametric maps of blood flow, T2* relaxation time, perfusion fraction f, diffusion coefficient D, and pseudodiffusion coefficient D*.
Pseudocontinuous ASL perfusion was calculated using a one-compartment model for blood after subtracting the tagged images from the nontagged control images. Blood flow was quantified using the following equation:
where PLD is the post-labeling delay (1525 ms), λ is the tissue partition coefficient (0.9 ml/g) , T1, blood is the longitudinal relaxation time of blood (1600 ms) , α is the labeling efficiency (0.80 × 0.75) (label pseudocontinuous ASL × background suppression) , τ is the labeling duration (1450 ms), SI PD is the proton density reference without labeling or background saturation, and SI control and SI label are the control and tagged signals, respectively. The arterial transit time was ignored in the quantification.
T2* relaxation time was calculated from the multi-echo T2* gradient-recalled echo data using the least-square fit of monoexponential decay , according to the following equation: S(TE) = STE0·exp.(−TE/T2*), where TE is the gradient echo time and S(TE) and STE0 are the measured signal intensities for TE ≠ 0 and TE = 0, respectively.
The multi-b-value diffusion-weighted images were analyzed using the IVIM model according to the following equation: S(b) = Sb0·[f·exp.(−bD*) + (1-f) ·exp.(−bD)], where S(b) is the measured signal intensity obtained with a nonzero b-value and Sb0 is the measured signal intensity for b = 0. With this equation, perfusion fraction f together with diffusion coefficient D and pseudodiffusion coefficient D* were calculated using a nonlinear biexponential fit based on the Levenberg-Marquardt technique .
For BOLD and IVIM analysis, goodness of fit was assessed using R2 = 1 − ESS/TSS, where ESS is the sum of squared errors between the data points and the fitting curve and TSS is the sum of squared differences between the data points and the mean value of all data points. Additional computations were performed to assess the signal-to-noise ratios (SNRs) of the images obtained with TE of 65.5 ms for BOLD and b-value of 800 s/mm2 for IVIM. Noisy images were excluded from curve fitting.
Regions of interest (ROIs) were manually drawn on T1-weighted images around the 4 muscle groups (anterior, lateral, soleus, and gastrocnemius) at the largest cross-sectional area of the calf on both the experimental and control sides (Fig. 1). Attention was given to exclude areas influenced by bones and large vessels, and the inter-osseous muscle was not investigated because it contains a relatively large number of vessels . ROIs were then copied and pasted into the corresponding functional perfusion maps. Average values within ROIs were recorded. Independent analysis of perfusion maps (from 7 randomly selected healthy subjects) by 2 radiologists blinded to the clinical outcomes was conducted to evaluate the interreader reproducibility. In addition, the repeat scans of 5 healthy subjects were analyzed by the same blinded radiologist to test the interscan reproducibility. Normalized values were obtained by dividing each imaging parameter value obtained under the experimental statuses (Cuff-20/Cuff-40/Post) by the baseline measurement (Pre).
TcPO2 measurements were acquired in all healthy young subjects the day after the CMR examinations using a TcPO2 monitoring system (Periflux System 5000; Perimed, Jarfalla, Sweden) in an air-conditioned room maintained at 22 °C. The cuff-induced ischemia paradigm followed the same process as that described for the CMR experiments. One author (LZ, 12 years of experience in TcPO2 measurement) placed the electrode at the same spot at which the CMR measurements were taken (i.e., the medial upper third of the calf adjacent to the gastrocnemius muscle at the maximal calf diameter) .
Statistical analyses were performed using SPSS version 20 (International Business Machines, Armonk, New York, USA), OriginPro 2016 (OriginLab Corp., Northampton, Massachusetts, USA), and Prism 5 (GraphPad Software Inc., La Jolla, California, USA). P-values less than 0.05 were considered to indicate statistical significance.
All data were expressed as the median (range) owing to non-normal data distributions. To assess interreader and interscan reproducibility, the intraclass correlation coefficient (ICC) was calculated. ICC values less than 0.40 indicated poor reproducibility, those ranging from 0.40 to 0.75 indicated fair to good reproducibility, and those higher than 0.75 indicated excellent reproducibility.
To determine whether differences existed in each imaging parameter value under different statuses (Pre/Cuff-20/Cuff-40/Post) in the experimental and control lower extremities, the Friedman test was used. In cases of statistical significance, further pairwise comparisons with the Dunn test were performed. The Wilcoxon signed rank test was used to compare imaging parameter measurements between left and right lower extremities in all healthy subjects under different statuses. The non-parametric Spearman rank correlation test was performed to assess correlations between imaging parameters derived from different methods, as well as between imaging parameters in gastrocnemius and TcPO2 measurements. Strength of correlation based on the Spearman rank correlation coefficient (ρ) was interpreted as follows: 0.00 to 0.20, very weak to negligible correlation; 0.21 to 0.40, weak correlation; 0.41 to 0.70, moderate correlation; 0.71 to 0.90, strong correlation; and 0.91 to 1.00, very strong correlation . The effects of age on the baseline imaging parameter measurements were investigated by comparing the values between the healthy young and healthy old groups with the Mann-Whitney U test. The Wilcoxon signed rank test was used for comparisons of imaging parameter measurements between the left and right lower extremities in PAD patients. Then, imaging parameters were compared between PAD patients and healthy old subjects, as well as between mild-to-moderate and severe PAD patients using the Mann-Whitney U test. Finally, correlations between imaging parameters and ABI in PAD patients were assessed using the Spearman rank correlation test. Bonferroni correction for multiple comparisons for the number of muscle groups was applied when necessary.
The sample size of patients included in this study was estimated using one-sided calculations with α of 0.05 and a power of 80% to detect an absolute T2* decrease of 4 ms (compared with normal) with a standard deviation of 2 ms based on the results of previous studies [24, 25]. Assuming a 20% dropout rate, it was determined that 10 participants were required.
All subjects successfully completed the CMR examinations (Table 1). Cuff compression of the thigh was well tolerated. Quantitative image analysis was conducted for each participant. Representative source images and parametric maps are illustrated in Additional file 1: Figure S1. The median R2 values for BOLD and IVIM fittings were 0.90 (range, 0.80–0.98) and 0.98 (range, 0.92–1.00), respectively. The median SNR values on the TE = 65.5 ms and b = 800 s/mm2 images were 19.2 (range, 14.4–23.5) and 17.2 (range, 14.2–32.0), respectively. The overall interreader reproducibility was fair to excellent, with ICC values of 0.83 for blood flow, 0.93 for T2*, 0.77 for f, 0.92 for D, and 0.64 for D*. The overall interscan reproducibility was fair to excellent, with ICC values of 0.73 for blood flow, 0.85 for T2*, 0.67 for f, 0.84 for D, and 0.55 for D*.
Functional imaging parameter variation under the cuff compression paradigm
Changes in quantitative imaging parameters under the cuff compression paradigm in the healthy young and healthy old groups are illustrated in Figs. 2 and 3 for the experimental lower extremity and in Additional file 2: Figure S2 and Additional file 3: Figure S3 for the control side.
In the healthy young group on the experimental side, only blood flow and T2* values showed significant differences among the 4 statuses in all muscle groups (all P < 0.01 for blood flow, and all P < 0.001 for T2*). f and D values showed significant differences in the lateral, soleus, and gastrocnemius muscle groups, and D* values were significantly different in the soleus and gastrocnemius muscle groups (all P < 0.05). Under the Cuff-20/Cuff-40 compression statuses, T2* values were significantly lower than the baseline measurements in all muscle groups (all P < 0.05), and blood flow, f, D, and D* values showed marked differences in only some muscle groups. No significant differences for any parameter in either muscle group were observed between the Post and Pre statuses (all P > 0.05) (Fig. 2).
In the healthy old group on the experimental side, only T2*, f, and D values showed significant differences among the 4 statuses in all muscle groups (all P < 0.001 for T2*, and all P < 0.01 for f and D). D* values showed significant differences in the soleus and gastrocnemius muscle groups (both P < 0.01). Blood flow values did not significantly differ in any muscle group (all P > 0.05). Under the Cuff-20/Cuff-40 compression statuses, significant differences were found for T2* values compared with those at baseline in all muscle groups except for the anterior muscle group under the 20-mmHg status (all P < 0.05). Furthermore, f, D, and D* values were significantly different in several muscle groups under these statuses (all P < 0.05). Additionally, no significant differences between Pre and Post measurements were observed for any of the parameters (all P > 0.05) (Fig. 3).
In both the healthy young and healthy old groups on the experimental side, in the cases of significant differences, blood flow, T2*, and D values gradually decreased under the Cuff-20 and Cuff-40 compression statuses, while f values showed a tendency to increase. During the recovery period, all parameters nearly returned to normal (Figs. 2 and 3).
In both the healthy young and healthy old groups on the control side, among all parameters, only T2* values in some muscle groups showed significant differences among the 4 statuses (P < 0.05) (Additional file 2: Figure S2 and Additional file 3: Figure S3).
Comparison of imaging parameters between the experimental and control sides
Results of Wilcoxon signed rank test comparing functional imaging parameters between the experimental and control sides under different statuses are illustrated in Table 2. No significant differences in any parameter were observed between the left and right lower extremities in healthy subjects at rest (all P > 0.0125). Most muscle groups exhibited significant differences in perfusion-related parameter values especially when the cuff compression pressure was increased to 40 mmHg above systolic pressure.
Correlation between functional imaging parameters derived from ASL, BOLD, and IVIM
No significant correlations between functional imaging parameters derived from the different methods were observed (all P > 0.05) (Additional file 4: Table S1).
Correlation between functional imaging parameters and TcPO2 measurements
The normalized CMR imaging parameters blood flow, T2*, and f were all correlated with normalized TcPO2 measurements (Fig. 4). Blood flow and T2* showed significant moderate correlations with TcPO2 measurements (ρ = 0.465 (P = 0.001) and ρ = 0.522 (P < 0.001), respectively). A significant negative correlation was observed between f and TcPO2 measurements (P = 0.018), although the correlation was weak (ρ = − 0.351). No significant correlation was found for D (P = 0.054) or D* values (P = 0.340).
Effects of age on resting-state functional imaging parameter measurements in healthy subjects
Within the muscle groups studied, blood flow for all muscle groups, f for the soleus and gastrocnemius groups, and D and D* for the gastrocnemius group demonstrated significant differences between the healthy young and healthy old groups at rest (all P < 0.0125). Only T2* was found to be independent of age in every muscle group (all P > 0.0125) (Table 3).
Comparison of resting-state functional imaging parameters between the affected and contralateral sides in PAD patients
In PAD patients, T2* was markedly reduced on the affected side compared with that on the contralateral side; significance was reached in the anterior muscle group (P = 0.005), and no significance was reached in the lateral, soleus, or gastrocnemius muscle groups (all P > 0.0125) (Additional file 5: Table S2).
Using ASL, BOLD, and IVIM to measure perfusion in PAD patients
Results of Mann-Whitney U tests comparing functional perfusion parameters between PAD patients and age-matched healthy subjects and between mild-to-moderate and severe PAD patients at rest are illustrated in Table 4. Among all parameters, T2* showed the best performance for the discrimination, with significantly reduced values observed in PAD patients compared with age-matched healthy old subjects and severe PAD patients compared with mild-to-moderate patients in all muscle groups (all P < 0.0125). No significant differences in other parameters in any muscle group were observed except D* values between PAD patients and healthy old subjects in the lateral muscle group.
Spearman rank correlation analysis showed that in PAD patients, T2* was significantly correlated with ABI in the anterior (ρ = 0.837; P < 0.001), lateral (ρ = 0.820; P < 0.001), soleus (ρ = 0.785; P = 0.001), and gastrocnemius (ρ = 0.644; P = 0.012) muscle groups (Fig. 5), whereas no significant correlation was observed for the other parameters (all P > 0.05).
Noninvasive monitoring of skeletal muscle perfusion in the lower extremities is critical for PAD patient management as perfusion can provide insight into microvascular function and endothelial integrity . Advanced CMR techniques, including ASL, BOLD, and IVIM, have been utilized in the assessment of skeletal muscle perfusion [12,13,14]; however, to our knowledge, few studies have been performed to directly compare these techniques in the same subject cohort [5, 26]. Furthermore, there is limited literature on comparisons of these CMR parameters with routinely used parameters in clinical practice such as TcPO2 and ABI [5, 7, 13]. The results of our study suggested that 1) ASL, BOLD, and IVIM could respond to cuff-induced ischemia in healthy subjects—that is, when the difference reached a significant level, ASL-derived blood flow values, BOLD-derived T2* values, and IVIM-derived D values tended to decrease with increasing external pressure while IVIM-derived f values tended to increase under cuff compression; 2) blood flow, T2*, and f values were all correlated with TcPO2 measurements; 3) ASL of all muscle groups and IVIM of the gastrocnemius group were influenced by age; only BOLD was independent of age in every muscle group; 4) BOLD could detect perfusion differences between the affected and contralateral lower extremities in PAD patients; 5) BOLD could separate PAD patients from healthy old subjects and PAD patients with different severities; and 6) BOLD-derived T2* was correlated with ABI in PAD patients.
Similar to the results in published literature [19, 27], baseline blood flow in skeletal muscle as measured by ASL in the lower extremity was mostly near or less than 20 ml/100 g/min in our study. For BOLD CMR, we observed a baseline mean T2* value of approximately 25 ms in healthy subjects, which lies within the previously reported range of 22–27 ms [13, 25, 28, 29]. IVIM imaging of skeletal muscle in the lower extremity has rarely been studied. In other parts of the body, IVIM-derived f, D, and D* values were reported to be 3%, 1.45 × 10− 3 mm2/s, and 28.5 × 10− 3 mm2/s, respectively, in the forearm muscle at rest by Filli et al. , and 6.6%, 1.45 × 10− 3 mm2/s, and 11.7× 10− 3 mm2/s, respectively, in the shoulder muscle by Nguyen et al. . These findings are consistent with our observations (e.g., 5.9%, 1.52 × 10− 3 mm2/s, and 18.6 × 10− 3 mm2/s, respectively, in the soleus muscle in healthy young subjects).
Under cuff compression conditions, negative ASL and BOLD contrasts in healthy young subjects developed due to ischemic insult, which agreed with previous studies [7, 12, 13, 25, 31]. ASL is capable of measuring blood flow through muscle tissue microvasculature given that ASL and radiolabeled microsphere measurements in rat leg muscle have shown good correlation for perfusion . Cuff compression interrupted both arterial inflow and venous outflow simultaneously, thus provoking reduced blood flow obtained by ASL. Although the exact source of the BOLD signal in skeletal muscle is not yet fully understood, it is generally accepted that the signal is primarily associated with capillary and blood oxygenation state . Lebon et al. also found that the T2* signal in muscle rapidly decreased during ischemia and attributed this change to early hemoglobin desaturation . This finding is logical given that the BOLD signal changed almost synchronously with hemoglobin desaturation but preceded myoglobin desaturation , as the dissociation constant of hemoglobin is more than 10 times higher than that of myoglobin. IVIM-derived D values also showed a decreasing trend under cuff-induced ischemic conditions. Local ischemia leads to decreased diffusivity of water molecules within muscles, and this decreased D was mostly likely attributed to this physical effect. Conversely, the perfusion fraction f derived from IVIM showed a tendency to increase in the case of arterial occlusion. Given that f reflects the signal fraction of capillary blood flow in entire water molecule diffusion pool within each voxel , it can be hypothesized that obstruction of venous reflux is probably responsible for the altered f values. In addition, it has been suggested that decrease in venous oxygen saturation may release relaxing factors than can cause microvascular dilation , which may also increase the f value. In healthy old subjects, changing trends for the imaging parameters were similar to those observed in the healthy young group except for ASL. ASL is limited by the intrinsic low SNR in skeletal muscle, wherein the ASL signal represents only 0.5%–1% of the raw image intensity . This notion may account for the lack of statistical significance of ASL measurements, especially given that the healthy old sample number was small (n = 10).
Interestingly, we also observed aberrant T2* signal changes on the control sides of healthy subjects during the cuff compression experiment on the other lower extremity, in accordance with the findings of a previous study . Yeung et al. suggested that this might be because of the high sensitivity of BOLD CMR imaging to local magnetic field disturbances caused by magnetic susceptibility effects, which may be induced by oxygen in the air at high pressure during cuff inflation .
Ledermann et al. reported that BOLD CMR imaging correlated with TcPO2 measurements in healthy volunteers during muscle ischemia , which was consistent with our results. However, the correlation observed by Ledermann et al. was stronger (correlation coefficient, 0.96) than that in our study (correlation coefficient, 0.522), primarily because Ledermann et al. averaged signal intensities across all volunteers for statistics, not for individuals. In addition, ASL-derived blood flow and IVIM-derived f values were also found to correlate with TcPO2 measurements in our study, with correlation coefficients of 0.465 and − 0.351, respectively. The correlation between BOLD-derived T2* and TcPO2 measurements was stronger than those between the other parameters, which may be attributed to the fact that both BOLD CMR imaging and TcPO2 examination are directly associated with the oxygenation state at the microvascular level.
Age effects on parameter measurements varied among the different sequences and muscles. In our study, no significant differences between the healthy young and healthy old groups were observed for the baseline T2* value, indicating that age did not appear to affect BOLD in healthy subjects at rest. This finding was in accordance with other observations . In contrast, ASL and IVIM were more easily influenced by the age factor especially in the gastrocnemius muscle group. The gastrocnemius is a fast-twitch type muscle, whereas the soleus belongs to the slow-twitch type. Degenerative processes of muscle fibers have been demonstrated to differ with fiber type, and the fast-twitch muscle is more prone to aging and fatigue [36, 37].
The ABI is a measure providing objective data for diagnosing PAD. When applying these imaging techniques on PAD patients, we found that BOLD was capable of detecting perfusion deficits at rest better than ASL and IVIM. Lower T2* values were related to the presence of PAD and more disease severity stratified by ABI. As discussed earlier, BOLD effect is generally assumed to reflect blood oxygenation state influenced by the ratio of oxygenated to deoxygenated hemoglobin, which is determined by the balance between oxygen supply and consumption . In PAD, arterial blood flow in the lower extremities is limited, leading to reduced oxygenated hemoglobin. Moreover, the impaired vascular function causes a longer contact time between blood and myocytes, leading to more efficient deoxygenation of hemoglobin . These two effects both contribute to a reduced T2* value. In a previous study by Englund et al., BOLD-derived metrics under the ischemia-reperfusion paradigm were also found to be correlated with ABI, suggestive of disease severity-dependent impairment of vascular response in PAD patients . Our results suggest that even at rest, vascular function at the tissue level could be indicative of disease progression. However, unlike CMR imaging, ABI cannot be obtained from all patients especially in patients with critical limb ischemia.
In the cuff-induced arterial occlusion experiment on healthy subjects, varying degrees of ischemia were factitiously induced by different pressures above the systolic pressure, and CMR imaging techniques were able to detect these changes. However, in healthy old subjects and PAD patients with differing degrees of ischemia, only BOLD was effective for this discrimination at rest. One possible explanation for this finding could be that the degree of ischemic insult in PAD patients was less severe than that induced by cuff occlusion. Additionally, collateral arteries developed in skeletal muscles in PAD patients would compensate for the perfusion deficit. BOLD is more sensitive to these less dramatic changes, which may be ASL- or IVIM-insensitive.
Numerous previous studies have used ASL or BOLD to monitor dynamic perfusion changes in skeletal muscles at rest and during ischemia and hyperemia, which allows the measurement of key parameters, such as peak hyperemic value (PHV) and time-to-peak (TTP). In our study, we did not measure continuous temporal changes of ASL, BOLD, and IVIM in our subjects mainly because it was not technically feasible to perform these three sequences sequentially at a high temporal resolution. Nevertheless, our protocol can be regarded as a simplified approach to the dynamic scanning method; a similar approach was also used in a previous study . Although PHV and TTP proved useful for the assessment of PAD in most published studies, conflicting results also exist . In addition, the reproducibility of the data is another concern. Versluis et al. investigated the reproducibility of BOLD-derived PHV and TTP values in healthy subjects and PAD patients . The reproducibility was unsatisfying with a coefficient of variation up to 26.7% and an ICC value as low as 0.59 . Moreover, due to massive pain or risk of worsening the clinical condition caused by cuff compression, patients with critical limb ischemia, ulceration, necrosis, or gangrene should be considered with caution or even excluded from the study . Compared with the cuff compression paradigm, the resting-state imaging scheme is simple to perform, less time-consuming, and more acceptable to PAD patients . In our study, we wanted to investigate and were most concerned with whether the baseline measurements of these techniques had value in assessing PAD.
The present study has several limitations. First, the number of study subjects was relatively small. Second, no gold standard for blood flow, oxygenation, or microvascular perfusion in skeletal muscle could be established in our subjects. For example, mixed venous oxygen saturation measurements would be informative regarding the confirmation of BOLD results. Nevertheless, our results revealed significant relationships with TcPO2 measurements, which are commonly used in clinical routines. However, the use of TcPO2 measurement is limited since it is confined to the skin microvasculature and thus fails to directly analyze the skeletal muscle . Therefore, once further validated, noninvasive CMR techniques might be used to gain direct information regarding skeletal muscle perfusion in lower extremities. Third, the cuff compression paradigm was used in the study instead of exercise. Exercise is more physiologically and clinically relevant. A previous study showed that muscle perfusion at peak exercise was correlated with 6-min walk distance in lower extremities . Compared with exercise, cuff compression has improved test-retest reproducibility  and less motion artifacts . Moreover, cuff compression allow the assessment of muscles which respond less to commonly used ankle flexion exercise . Lopez et al. suggested the use of cuff compression in a more general study population and the use of exercise in specific PAD therapies in claudicants based on its physiologic and clinical relevance . Fourth, skeletal muscle energetics were not investigated in the current study. 31P CMR spectroscopy is a useful tool to noninvasively probe skeletal muscle energetics, including adenosine triphosphate and creatine phosphate metabolism , which would help to better understand the perfusion results in our study. Further studies are warranted. Finally, CMR angiography was not performed in the current study.
In conclusion, the present study shows that multi-parametric CMR techniques including ASL, BOLD, and IVIM provide useful and complementary information regarding tissue perfusion in the lower extremities of healthy subjects. Perfusion-related metrics derived from these techniques correlate with TcPO2 measurements. In PAD patients, BOLD is a more reliable imaging technique for the detection and stratification of alterations in microvascular function at rest compared with ASL and IVIM.
Arterial spin labeling
Blood oxygenation level-dependent
Cardiovascular magnetic resonance
Sum of squared errors
Field of view
Intraclass correlation coefficient
Intravoxel incoherent motion
Peripheral arterial disease
Peak hyperemic value
Region of interest
Transcutaneous oxygen pressure
Sum of squared differences
Hirsch AT, Criqui MH, Treat-Jacobson D, Regensteiner JG, Creager MA, Olin JW, et al. Peripheral arterial disease detection, awareness, and treatment in primary care. JAMA. 2001;286:1317–24.
Criqui MH, Langer RD, Fronek A, Feigelson HS, Klauber MR, McCann TJ, et al. Mortality over a period of 10 years in patients with peripheral arterial disease. New Engl J Med. 1992;326:381–6.
Kullo IJ, Rooke TW. Peripheral artery disease. N Engl J Med. 2016;374:861–71.
Stein R, Hriljac I, Halperin JL, Gustavson SM, Teodorescu V, Olin JW. Limitation of the resting ankle-brachial index in symptomatic patients with peripheral arterial disease. Vasc Med. 2006;11:29–33.
Englund EK, Langham MC, Ratcliffe SJ, Fanning MJ, Wehrli FW, Mohler ER, et al. Multiparametric assessment of vascular function in peripheral artery disease dynamic measurement of skeletal muscle perfusion, blood-oxygen-level dependent signal, and venous oxygen saturation. Circ Cardiovasc Imaging. 2015;8:e002673.
Bunt T, Holloway GA. TcPO2 as an accurate predictor of therapy in limb salvage. Ann Vasc Surg. 1996;10:224–7.
Ledermann HP, Heidecker H-G, Schulte A-C, Thalhammer C, Aschwanden M, Jaeger KA, et al. Calf muscles imaged at BOLD MR: correlation with TcPo2 and Flowmetry measurements during ischemia and reactive hyperemia—initial experience. Radiology. 2006;241:477–84.
Bajwa A, Wesolowski R, Patel A, Saha P, Ludwinski F, Smith A, et al. Assessment of tissue perfusion in the lower limb current methods and techniques under development. Circ Cardiovasc Imaging. 2014;7:836–43.
Detre JA, Leigh JS, Williams DS, Koretsky AP. Perfusion imaging. Magn Reson Med. 1992;23:37–45.
Ogawa S, Lee T-M, Kay AR, Tank DW. Brain magnetic resonance imaging with contrast dependent on blood oxygenation. Proc Natl Acad Sci U S A. 1990;87:9868–72.
Le Bihan D, Breton E, Lallemand D, Aubin M, Vignaud J, Laval-Jeantet M. Separation of diffusion and perfusion in intravoxel incoherent motion MR imaging. Radiology. 1988;168:497–505.
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.
Bajwa A, Wesolowski R, Patel A, Saha P, Ludwinski F, Ikram M, et al. Blood oxygenation level-dependent CMR-derived measures in critical limb ischemia and changes with revascularization. J Am Coll Cardiol. 2016;67:420–31.
Filli L, Boss A, Wurnig MC, Kenkel D, Andreisek G, Guggenberger R. Dynamic intravoxel incoherent motion imaging of skeletal muscle at rest and after exercise. NMR Biomed. 2015;28:240–6.
Jacobi B, Bongartz G, Partovi S, Schulte AC, Aschwanden M, Lumsden AB, et al. Skeletal muscle BOLD MRI: from underlying physiological concepts to its usefulness in clinical conditions. J Magn Reson Imaging. 2012;35:1253–65.
Duteil S, Wary C, Raynaud JS, Lebon V, Lesage D, Leroy-Willig A, et al. Influence of vascular filling and perfusion on BOLD contrast during reactive hyperemia in human skeletal muscle. Magn Reson Med. 2006;55:450–4.
Loughran T, Higgins DM, McCallum M, Coombs A, Straub V, Hollingsworth KG. Improving highly accelerated fat fraction measurements for clinical trials in muscular dystrophy: origin and quantitative effect of R2* changes. Radiology. 2015;275:570–8.
Mutsaerts HJ, Steketee RM, Heijtel DF, Kuijer JP, van Osch MJ, Majoie CB, et al. Inter-vendor reproducibility of pseudo-continuous arterial spin labeling at 3 tesla. PLoS One. 2014;9:e104108.
Grozinger G, Pohmann R, Schick F, Grosse U, Syha R, Brechtel K, et al. Perfusion measurements of the calf in patients with peripheral arterial occlusive disease before and after percutaneous transluminal angioplasty using MR arterial spin labeling. J Magn Reson Imaging. 2014;40:980–7.
Ledermann H-P, Schulte A-C, Heidecker H-G, Aschwanden M, Jäger KA, Scheffler K, et al. Blood oxygenation level–dependent magnetic resonance imaging of the skeletal muscle in patients with peripheral arterial occlusive disease. Circulation. 2006;113:2929–35.
Cao M, Suo S, Han X, Jin K, Sun Y, Wang Y, et al. Application of a simplified method for estimating perfusion derived from diffusion-weighted MR imaging in glioma grading. Front Aging Neurosci. 2018;9:432.
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.
Heye S, Maleux G, Oyen RH, Claes K, Kuypers DR. Occupational radiation dose: percutaneous interventional procedures on hemodialysis arteriovenous fistulas and grafts. Radiology. 2012;264:278–84.
Zuo CS, Sung YH, Simonson DC, Habecker E, Wang J, Haws C, et al. Reduced T2* values in soleus muscle of patients with type 2 diabetes mellitus. PLoS One. 2012;7:e49337.
Englund EK, Langham MC, Li C, Rodgers ZB, Floyd TF, Mohler ER, et al. Combined measurement of perfusion, venous oxygen saturation, and skeletal muscle T2* during reactive hyperemia in the leg. J Cardiovasc Magn Reson. 2013;15:70.
Andreisek G, White LM, Sussman MS, Langer DL, Patel C, Su JW, et al. T2*-weighted and arterial spin labeling MRI of calf muscles in healthy volunteers and patients with chronic exertional compartment syndrome: preliminary experience. AJR Am J Roentgenol. 2009;193:W327–33.
Wu WC, 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.
Kos S, Klarhofer M, Aschwanden M, Scheffler K, Jacob AL, Bilecen D. Simultaneous dynamic blood oxygen level-dependent magnetic resonance imaging of foot and calf muscles: aging effects at ischemia and postocclusive hyperemia in healthy volunteers. Investig Radiol. 2009;44:741–7.
Partovi S, Schulte AC, Jacobi B, Klarhofer M, Lumsden AB, Loebe M, et al. Blood oxygenation level-dependent (BOLD) MRI of human skeletal muscle at 1.5 and 3 T. J Magn Reson Imaging. 2012;35:1227–32.
Nguyen A, Ledoux JB, Omoumi P, Becce F, Forget J, Federau C. Application of intravoxel incoherent motion perfusion imaging to shoulder muscles after a lift-off test of varying duration. NMR Biomed. 2016;29:66–73.
Yu G, Floyd TF, Durduran T, Zhou C, Wang J, Detre JA, et al. Validation of diffuse correlation spectroscopy for muscle blood flow with concurrent arterial spin labeled perfusion MRI. Opt Express. 2007;15:1064–75.
Pohmann R, Künnecke B, Fingerle J, Kienlin MV. Fast perfusion measurements in rat skeletal muscle at rest and during exercise with single-voxel FAIR (flow-sensitive alternating inversion recovery). Magn Reson Med. 2006;55:108–15.
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.
Clifford PS, Hellsten Y. Vasodilatory mechanisms in contracting skeletal muscle. J Appl Physiol. 2004;97:393–403.
Yeung DK, Griffith JF, Li AF, Ma HT, Yuan J. Air pressure-induced susceptibility changes in vascular reactivity studies using BOLD MRI. J Magn Reson Imaging. 2013;38:976–80.
Ishihara A, Naitoh H, Katsuta S. Effects of ageing on the total number of muscle fibers and motoneurons of the tibialis anterior and soleus muscles in the rat. Brain Res. 1987;435:355–8.
Noseworthy MD, Bulte DP, Alfonsi J. BOLD magnetic resonance imaging of skeletal muscle. Semin Musculoskelet Radiol. 2003;7:307–15.
Wang C, Zhang R, Zhang X, Wang H, Zhao K, Jin L, et al. Noninvasive measurement of lower extremity muscle oxygen extraction fraction under cuff compression paradigm. J Magn Reson Imaging. 2016;43:1148–58.
Versluis B, Backes WH, van Eupen MG, Jaspers K, Nelemans PJ, Rouwet EV, et al. Magnetic resonance imaging in peripheral arterial disease: reproducibility of the assessment of morphological and functional vascular status. Investig Radiol. 2011;46:11–24.
Partovi S, Karimi S, Jacobi B, Schulte AC, Aschwanden M, Zipp L, et al. Clinical implications of skeletal muscle blood-oxygenation-level-dependent (BOLD) MRI. MAGMA. 2012;25:251–61.
Anderson JD, Epstein FH, Meyer CH, Hagspiel KD, Wang H, Berr SS, et al. Multifactorial determinants of functional capacity in peripheral arterial disease: uncoupling of calf muscle perfusion and metabolism. J Am Coll Cardiol. 2009;54:628–35.
Lopez D, Pollak AW, Meyer CH, Epstein FH, Zhao L, Pesch AJ, et al. Arterial spin labeling perfusion cardiovascular magnetic resonance of the calf in peripheral arterial disease: cuff occlusion hyperemia vs exercise. J Cardiovasc Magn Reson. 2015;17:23.
Weiss K, Schar M, Panjrath GS, Zhang Y, Sharma K, Bottomley PA, et al. Fatigability, exercise intolerance, and abnormal skeletal muscle energetics in heart failure. Circ Heart Fail. 2017;10:e004129.
This work was supported by National Natural Science Foundation of China (grants 81501458, 81271638, and 81571630) and the Shanghai Pujiang Program (grant 15PJ1405200).
Availability of data and materials
The data that support the findings of this study are available from corresponding author upon reasonable request.
Ethics approval and consent to participate
The protocol was approved by the Institutional Review Board at Renji Hospital, School of Medicine, Shanghai Jiao Tong University, and all participants provided written informed consent.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Example of source images and processed maps for arterial spin labeling (ASL), blood-oxygen level dependent (BOLD), and intravoxel incoherent motion (IVIM) cardiovascular magnetic resonance (CMR). (PPTX 11134 kb)
Figure S2. Graphs depicting serial measurements of imaging parameters from ASL, BOLD, and IVIM for the anterior, lateral, soleus, and gastrocnemius muscle groups in the control side in healthy young subjects. (PPTX 727 kb)
Figure S3. Graphs depicting serial measurements of imaging parameters from ASL, BOLD, and IVIM for the anterior, lateral, soleus, and gastrocnemius muscle groups in the control side in healthy old subjects. (PPTX 654 kb)
Table S1. Spearman rank correlation coefficients for ASL, BOLD, and IVIM imaging parameters. (DOCX 16 kb)
Table S2. Comparison of functional imaging parameters between affected and contralateral lower extremities in PAD patients. (DOCX 18 kb)
About this article
Cite this article
Suo, S., Zhang, L., Tang, H. et al. Evaluation of skeletal muscle microvascular perfusion of lower extremities by cardiovascular magnetic resonance arterial spin labeling, blood oxygenation level-dependent, and intravoxel incoherent motion techniques. J Cardiovasc Magn Reson 20, 18 (2018). https://doi.org/10.1186/s12968-018-0441-3
- Cardiovascular magnetic resonance
- Blood flow
- Diffusion-weighted imaging