In vivo characterization of rodent cyclic myocardial perfusion variation at rest and during adenosine-induced stress using cine-ASL cardiovascular magnetic resonance
Journal of Cardiovascular Magnetic Resonance volume 16, Article number: 18 (2014)
Assessment of cyclic myocardial blood flow (MBF) variations can be an interesting addition to the characterization of microvascular function and its alterations. To date, totally non-invasive in vivo methods with this capability are still lacking. As an original technique, a cine arterial spin labeling (ASL) cardiovascular magnetic resonance approach is demonstrated to be able to produce dynamic MBF maps across the cardiac cycle in rats.
High-resolution MBF maps in left ventricular myocardium were computed from steady-state perfusion-dependent gradient-echo cine images produced by the cine-ASL sequence. Cyclic changes of MBF over the entire cardiac cycle in seven normal rats were analyzed quantitatively every 6ms at rest and during adenosine-induced stress.
The study showed a significant MBF increase from end-systole (ES) to end-diastole (ED) in both physiological states. Mean MBF over the cardiac cycle within the group was 5.5 ± 0.6 mL g-1 min-1 at rest (MBFMin = 4.7 ± 0.8 at ES and MBFMax = 6.5 ± 0.6 mL g-1 min-1 at ED, P = 0.0007). Mean MBF during adenosine-induced stress was 12.8 ± 0.7mL g-1 min-1 (MBFMin = 11.7±1.0 at ES and MBFMax = 14.2 ± 0.7 mL g-1 min-1 at ED, P = 0.0007). MBF percentage relative variations were significantly different with 27.2 ± 9.3% at rest and 17.8 ± 7.1% during adenosine stress (P = 0.014). The dynamic analysis also showed a time shift of peak MBF within the cardiac cycle during stress.
The cyclic change of myocardial perfusion was examined by mapping MBF with a steady-pulsed ASL approach. Dynamic MBF maps were obtained with high spatial and temporal resolution (6ms) demonstrating the feasibility of non-invasively mapping cyclic myocardial perfusion variation at rest and during adenosine stress. In a pathological context, detailed assessment of coronary responses to infused vasodilators may give valuable complementary information on microvascular functional defects in disease models.
The pulsatility of coronary flow has been studied for several decades [1, 2]: as elsewhere in the arterial circulation, it is the pressure difference between the aorta and the peripheral capillaries that provides the driving force for coronary circulation. In a paradoxical way, the heart, which is dedicated to supply blood to the organs, restricts its own blood supply during peak contraction in systole unlike most other systemic vascular beds which have a fresh oxygenated blood supply in systole. The squeezing effect of myocardial contraction causes arterial blood inflow to peak during diastole and venous blood outflow to peak during systole [3, 4]. Microvessels are abundant in the myocardium so that myocardial blood volume (MBV) and flow (MBF) vary during the cardiac cycle in response to this pulsatile change in coronary circulation and cyclic variation in myocardial tension. Diseases such as diabetes, atherosclerosis, cardiomyopathies, and arterial hypertension result in functional and morphologic microvascular alterations, which may precede clinical signs and symptoms. The mechanical interaction between coronary microcirculation, and myocardial contraction and relaxation is fundamentally important for an understanding of intramural coronary hemodynamics. The morphology and function of capillaries across myocardial wall during a cardiac cycle play a pivotal role on the mechanical control of blood perfusion to regional myocardium together with arteriolar and venular mechanics and responses . Determination of the role of systolic extravascular compression on the pattern of blood flow in the microcirculation during the cardiac cycle has been limited to measurements made in those regions of the myocardium where direct observations of these microvessels is possible [6–8]. Quantitative assessment of myocardial perfusion and characterization of cyclic myocardial perfusion variation would help better understand the mechanisms underlying myocardial contractility and may give valuable complementary information on microvascular functional defects in non-ischemic heart diseases.
Cyclic change of regional MBV in the mouse heart has already been reported using cardiovascular magnetic resonance (CMR) based on the correlation with cyclic variation in the steady-state T2* shortening effect of an ultrasmall iron oxide nanoparticle intravascular contrast agent . Another study focused on noninvasive MBV variation assessment in the canine left ventricle using displacement encoding with stimulated echoes (DENSE) technique . Nevertheless, no study has assessed MBF variation over the entire cardiac cycle in the rodent heart. In humans, however, Radjenovic et al. demonstrated the feasibility of assessing systolic and diastolic MBF changes by dynamic contrast-enhanced CMR .
Even if first-pass CMR has been proved to quantitatively assess MBF [12–15], arterial spin labeling (ASL) has become a method of choice to quantitatively and non-invasively map rodent MBF [16–24]. ASL has also shown its potential to quantify human MBF [25, 26]. Cine-ASL has recently been proposed as an original approach to assess MBF [27, 28]. This technique is based on a combination of a continuous cine CMR gradient echo readout with a steady-pulsed arterial labeling approach. This scheme provides quasi-continuous tagging of the blood feeding the coronaries while maintaining compatibility with the constraints of cardiac motion and the highly pulsatile blood flow in the ascending aorta. The steady-state approach associated with the cine readout allowed us to dynamically map MBF over the entire cardiac cycle with high temporal resolution. In a preliminary study, cyclic change of regional MBF in healthy mice was measured at rest , and a significant MBF increase of 30% from end-systole (ES) to end-diastole (ED) was shown. Based on this finding, we propose a new protocol which compares cyclic MBF changes in healthy Wistar rats between rest and adenosine-induced stress.
All experiments were conducted according to a protocol approved by the University’s animal experimentation committee. Seven healthy female Wistar rats (age 25 weeks, body weight 329 ± 8 g; Janvier Laboratories, Le Genest-Saint-Isle, France) were anesthetized in an induction chamber with 3% isoflurane. During the experiments, a 2% dose of isoflurane mixed with 1L min-1 of pure oxygen gas stream was continuously delivered through a face mask. For the present study, isoflurane was chosen as it maintains cardiac function and permits real time regulation of the depth of anesthesia whereas barbiturates such as pentobarbital are known to have depressant effects on hemodynamic parameters . Nonetheless, isoflurane has been shown to induce coronary vasodilation and MBF increase at rest compared to other anesthetics . Isoflurane concentration was regulated using a dedicated vaporizer (Ohmeda/General Electric, Milwaukee, WI, USA) so as to obtain regular breathing frequencies in the range of 70 breaths per minute. Respiration was monitored using a pressure sensor connected to an air-filled balloon positioned on the back of the rat. Body temperature was monitored using a rectal probe and maintained at 37°C using a heating blanket with hot water circulation placed around the rat. The electrocardiogram (ECG) signal was monitored by placing two subcutaneous electrodes, one in the right foreleg and one at abdominal level. The electrodes were connected to an ECG trigger unit (Rapid Biomedical, Rimpar, Germany) to record the signal and to trigger the CMR sequence. Oxygen saturation was also monitored to ensure stable physiological conditions along the experiment. A catheter was placed in the tail vein to deliver adenosine. The venous cannula had a dead volume of about 120 μL.
Cine-ASL has recently been proposed as an original ASL scheme to assess MBF in small animals. Theory of the technique and validation against another ASL approach are described elsewhere [27, 28]. Briefly, cine-ASL relies on a fast ECG-gated cine Fast Low Angle Shot (cine-FLASH) sequence repeated over several cardiac cycles for each line of k-space. At each cardiac cycle, one single gradient-echo is substituted by a spatially selective inversion pulse, labeling the arterial blood at the level of the aortic root. This labeling has to be performed just before the arterial blood enters through the myocardial tissue via the coronaries, i.e. at ES. Using this approach, a continuous cine acquisition can be performed while the blood magnetization feeding the coronaries is driven into a steady-state, leading to a perfusion-dependent stationary regime of tissue magnetization. The repetition time TR between each gradient echo was kept as short as possible to maximize temporal resolution, six milliseconds in these experiments.
As for all prospectively gated techniques, the original sequence had to be stopped before the next ECG trigger so that the observable time window within the cardiac cycle was only about two third of the cardiac cycle. Also, the waiting delay caused steady-state interruptions that might bias the perfusion signal on the short time scale observed here. Compared with the original cine-ASL version , the sequence was therefore slightly modified so as to acquire during almost two cardiac cycles (Figure 1) after each trigger pulse. We assumed that the heart rate was constant over two cardiac cycles so that ECG-triggering every two heart beats was sufficient. By appropriately adjusting the number of cardiac phases to acquire per cardiac cycle as well as the number of echoes at which the labeling pulses occur, more than one cardiac cycle coverage without interrupting the steady-pulsed labeling scheme was possible. In summary, this modification was done for two reasons: (i) the cine readout covered almost two cardiac cycles with a negligible interruption of the FLASH steady-state, and (ii) the labeling was still performed every heart beat which is mandatory to optimize labeling efficiency (see perfusion quantification section below).
All imaging experiments were performed on a Bruker Biospec Avance 4.7T/30 imager with horizontal bore (Bruker, Ettlingen Germany). Homogeneous RF excitation was achieved using a proton volume resonator (diameter 60 mm, homogeneous length 80 mm), and the animals were positioned prone on an actively decoupled surface receive coil (Rapid Biomedical, Wurzburg, Germany). Prior to the perfusion measurements, scout images were acquired to determine the short-axis plane used for the perfusion imaging sequences. A cine-FLASH sequence in four-chamber long-axis view was used to achieve precise spatial and temporal placement of the labeling pulses at ES, just before aortic valve closure.
Perfusion measurements were carried out in the short-axis plane using cine-ASL with the following parameters: imaging slice thickness: 2 mm, field of view: 40 mm, in-plane resolution: 312 × 625 μm2, TE/TR: 1.64/6 ms, flip angle α: 8° (gaussian pulses, 0.5 ms), labeling slice thickness: 6 mm (adiabatic hyperbolic secant pulses, 3.7 ms). The number of acquired cardiac phases (Nechoes) was adjusted according to the animal’s heart rate, about 30 echoes per cardiac cycle and 50 echoes per cine block. Acquisition time was approximately 12 min at 400 bpm. All animals first underwent MBF quantification at rest. Continuous adenosine infusion was accomplished using a dedicated low-flow pump injector (World Precision Instruments, Sarasota, FL, USA) and a tail vein catheter. The infusion rate was set to 280 μg kg-1 min-1, which corresponds to twice the standard clinical rate used in human clinical practice. This rate was shown in preliminary experiments to reliably ensure a maximal response to adenosine. The second cine-ASL sequence was started 5 minutes after an initial heart rate drop observed upon arrival of adenosine and infusion was maintained until completion of the sequence. The criterion for appreciation of a hyperemic state was a significant increase in MBF compared to rest condition.
Global MBF maps in left ventricular myocardium were computed from the steady-state perfusion-dependent gradient-echo cine images. As a result of the theory presented earlier , the experimentally measured relevant quantity is the stationary signal difference between control and tag scans , which is related to myocardial blood flow by the following equation:
where, β corresponds to the average labeling efficiency (β = 1 for a complete inversion and β = 0.5 for a saturation), λ = 0.95 the blood-tissue partition coefficient for water  and the apparent relaxation time measured under the influence of rapid succession of small flip angle pulses . T1 was fixed to 1.31 s as reported in a previous study . MSS is the longitudinal steady-state magnetization observed under FLASH partial saturation and can be directly calculated from M0, TR, T1 and the excitation flip angle α:
As reported in previous work , the labeling efficiency reached with this particular steady-pulsed scheme was only 0.5, although adiabatic inversion pulses were used. Considering the large extent of the inversion slice of 6mm which covers both atria, pulmonary compartments and the aorta, blood entering into the coronaries has already been inverted several times with different inversion delays before reaching the microvascular compartment. It has been shown that this multiple exposure of blood spins to inversion pulses led to reduced labeling efficiency β of 0.5, confirmed by absence of signal in the LV blood pool as also observed earlier .
To perform image analysis, an in-house developed program running in an IDL environment (ITT, Boulder, CO, USA) was used, and perfusion maps were generated by applying this model pixel-by-pixel. Transmural MBF was quantified as an average of pixel values for regions of interest drawn on every map in the series (about 50 maps per acquisition). Global myocardium was first evaluated from the MBF maps for each rat as a function of time. Since the blood signal coming from the left ventricle and from the coronaries was significantly higher than the measured signal of interest, particular attention was paid at masking this signal on each individual map. Indeed, coronary flow is known to be highly variable within the cardiac cycle , and taking blood signal into account could lead to a misinterpretation of cyclic myocardial perfusion variation.
Since the cyclic variations observed at rest and during stress were dependent on the animal’s heart rate, data were temporally normalized with respect to the individual RR interval. In order to perform a group analysis, MBF data were sorted into bins representing the mean temporal resolution after the normalization process. Finally, myocardial perfusion reserve (MPR) was computed as a ratio of stress and rest MBF for each temporal bin.
To obtain the average change of MBF over the cardiac cycle, the two highest and lowest values found at end-diastole (ED) and end-systole (ES), respectively, were averaged so as to calculate the relative cyclic variation defined as ΔMBFCycl = (MBFMax - MBFMin)/MBFMax. Measurements of systolic and diastolic MBF were done in the first cardiac cycle for each individual rat at rest and during stress before averaging within the group. To assess inter-observer variability, all measurements were performed twice by two different readers.
Data were presented either as group mean value ± SD across the entire cardiac cycle or as group mean value ± SD for each temporal bin within the cardiac cycle. All statistical processing was performed using Prism 5 software (Graph Pad, San Diego, CA, USA). Statistical differences between systolic and diastolic values of MBF and MPR, i.e. minimum and maximum, were evaluated using a paired t-test. A probability value of <0.05 was taken as statistically significant. Correlation between heart rate and relative cyclic variation ΔMBFCycl was quantified by the Pearson rank correlation coefficient rho. A paired t-test was also used to compare cyclic variation measurements between the two observers.
The Cine-ASL approach allowed for dynamic quantification of MBF across the cardiac cycle with time-frames given by the cine repetition time of 6 milliseconds. Figure 2 shows an example of a perfusion map series over the entire cardiac cycle obtained with cine-ASL from the same rat at rest and during adenosine stress. By inducing strong hyperemia with perfusion significantly increased by a factor of 2.1. In this example, mean perfusion over the entire cardiac cycle in the global myocardium increased from 6.3 ± 0.6 mL g-1 min-1 at rest to 13.7 ± 1.3 mL g-1 min-1 during stress. This result demonstrates that cine-ASL has the potential to dynamically map MBF over the cardiac cycle and to study and quantify changes in myocardial perfusion within the cardiac cycle for different physiological states. In terms of breath and heart rates, no statistical differences were found between experiments at rest and during stress (Table 1).
Capillary flow dynamics
In Figure 3A all individual MBF data sets observed in the seven normal rats studied are shown as a function of time, at rest and during adenosine stress. Figure 3B shows normalized data over the cardiac cycle. Figure 4A shows the evolution of MBF across the cardiac cycle at rest and during adenosine-induced stress. These results are plotted as a function of the normalized cardiac cycle averaged over the seven animals studied ± standard deviation within the group. Along the series starting at early-systole, perfusion gradually and significantly decreased, reaching a minimum at ES and then recovered during the diastole. This general behavior was well reproducible in all subjects. The maximal perfusion values were always found to beat ED just before the next QRS complex. One can also notice a slight shift of the maximum MBF during stress compared to rest which was found to be reached sooner within the cardiac cycle. In Figure 4B, averaged MPR within the group is reported as a function of the normalized cardiac cycle. MPR behavior was found reversed compared to MBF variation, with an initial increase during systole until a maximum at ES. This maximum value was sustained during most of the diastole. Towards the end of the cardiac cycle, MPR dropped in a sudden way reaching a minimum at ED.
For each temporal bin along the cardiac cycle, mean perfusion values were found in agreement between the two observers. No significant differences were found between the two observers, with a Pearson correlation coefficient of r = 0.88 at rest and r = 0.89 during stress. We found a mean MBF over the cardiac cycle of 5.5 ± 0.6 mL g-1 min-1 at rest (observer 1, 5.7 ± 0.6 mL g-1 min-1 and observer 2, 5.4 ± 0.6 mL g-1 min-1, P = 0.12) which was significantly different from 12.8 ± 0.7 mL g-1 min-1 obtained during stress (observer 1, 13.1 ± 1.2 mL g-1 min-1 and observer 2, 12.6 ± 0.7 mL g-1 min-1, P = 0.38). Results are presented in Figure 5A. A significant difference in myocardial perfusion was found between rest and stress conditions (P<0.0001), leading to a mean perfusion reserve of 2.4 ± 0.2 within the group.
The two lowest values and the two highest values obtained within the cardiac cycle were identified and averaged for each studied animal. The minimum and maximum MBF values found at ES and ED were defined as and respectively. Individual results are reported in Figure 5B. Resting transmural MBF extrema were significantly different with = 4.7 ± 0.8 mL g-1 min-1 and = 6.5 ± 0.6 mL g-1 min-1 (P = 0.0007). During stress, they were also significantly different with = 11.7 ± 1.0 mL g-1 min-1 and = 14.2 ± 0.7 mL g-1 min-1 (P = 0.0007). The mean relative variation ΔMBFCycl within the group (Figure 5C), was found significantly different between both conditions. We found 27.2 ± 9.3% decrease from ED to ES at rest and 17.8 ± 7.1% during adenosine stress, with paired samples t-test P = 0.014. On the other hand, MPR was minimum at ED, 1.9 ± 0.1 and maximum at ES, 2.8 ± 0.3 (P < 0.0001). Individual quantifications as well as a statistic summary are reported in Table 2 for MBF and Table 3 for MPR.
Relationship between heart rate and cyclic MBF variation amplitude
In Figure 6A, the impact of the heart rate on the amplitude of MBF variation is highlighted. This graph shows the two extreme cases found within the group in terms of heart rate, i.e. the animal with the highest heart rate (Rat 1, 370 bpm) and the one with the lowest heart rate (Rat 3, 260 bpm). For clarity, a six orders polynomial regression was used to fit the data. ΔMBFCycl (Figure 6B) was lower for Rat 1 (24% at rest and 9% during stress) compared to the Rat 3 (42% at rest and 31% during stress). We verified whether a linear group correlation existed between individual heart rate and the corresponding relative variation ΔMBFCycl for both rest and stress conditions (Figure 7). There was a significant correlation between ΔMBFCycl and heart rate during stress (rho = −0.83, P = 0.02) but not at rest (rho = −0.67, P = 0.10).
In this study we report time-resolved absolute quantification of myocardial blood flow by CMR with a new approach based on a cine-FLASH labeling and readout module. The spatial resolution (313 × 625 μm2) and the sensitivity were sufficient to measure global myocardial perfusion and hemodynamic changes within the cardiac cycle in approximately 12 min. MBF maps were obtained with a high temporal resolution of 6 ms covering almost two cardiac cycles. The tagging scheme used here consists in rapidly repeated local inversions with precise timing within the cardiac cycle. It allows maintaining the labeled blood state throughout the experiment while a cine readout can be achieved during the remaining time with a steady-state acquisition. Tissue magnetization is thus driven into a perfusion-dependent equilibrium state so that signal changes that occur during the cardiac cycle can only reflect changes in momentary tissue blood flow. Arterial spin labeling measures the integrated effect of tagged blood flowing through the vasculature and remaining in the capillary system during the mean transit time (MTT). MTT is of the order of a second in the rodent heart and therefore much longer than a cardiac cycle. The averaged blood flow across the cardiac cycle reported here corresponds to the mean blood flow commonly defined by MBF = MBV/MTT which is understood on a much longer time scale than a cardiac cycle. By using high-temporal resolution cine readouts of 6 ms, signal that is acquired within a narrow window allows capturing the time-variant changes in the microvascular status of the myocardium, which is related to the number of patent myocardial microvessels. Thus, the variations observed at different time-points within the cardiac cycle have to be distinguished from the mean blood flow and are due to cyclic variations of momentary flow.
Capillary flow dynamics
Figure 4 shows that cyclic temporal variations were measurable using the cine-ASL technique.
Left-coronary arterial blood flow is characterized by a diastolic-dominant inflow  whereas small distal venous flow exhibits a systolic-dominant pattern . Our results are consistent with the expected physiological pattern of preferential coronary diastolic filling.
A clear and complete description of coronary flow patterns across the cardiac cycle in humans has been proposed by Davies et al. , who used intracoronary wires for measuring pressure and Doppler velocity waveforms. Microvascular MBF is closely related to macrovascular coronary flow such that these results can be confronted to ours. During systole, the authors observed a dominant forward-travelling pushing wave in proximal coronary arteries that is reflected when reaching the microvascular bed in distal coronary arterioles, causing blood to move in the opposite direction. Coronary arteries are also known to have a significant capacitance to accommodate for blood supply . Also, during ventricular contraction, the compressive forces on the small vessels lying in the myocardium increase while venous outflow is in its dominant phase. One can therefore expect decreasing capillary flow during systole, which is consistent with our observations. During diastole, however, Davies and coworkers highlighted a dominant backward-travelling ‘suction’ wave initiated at the beginning of ventricular relaxation along with decreased compressive forces on capillaries and resistance of microvessels. Closure of the aortic valve in early diastole generates a late forward-travelling pushing wave which accelerates blood still further towards the myocardium. Most of the arterial inflow occurs during diastole with a velocity peak at early to mid-diastole. We observed the same pattern with a delay that corresponds to the time for the blood to reach the microvascular compartment.
We showed that momentary blood flow is non-uniform through the cardiac cycle and occurs mainly in diastole when myocardial tension is low. Therefore, the timing of data acquisition in perfusion CMR experiments may affect the measured MBF and might be considered in the evaluation of perfusion exams in both research and clinical practice. We note here that for studying temporal variations of perfusion using arterial spin labeling CMR, the underlying spin-exchange model has to be kept in mind. One-compartment models are valid for perfusion quantification in the heart , in which free diffusivity and complete exchange of the blood magnetization with the tissue during the presence of a spin in the arterial capillary compartment is assumed. However, the time scale relevant here is the delay between successive cardiac frames, which is much shorter than the time for spins to travel across the arterial capillary compartment, i.e. the time relevant for perfusion model validity. Taking into account spin exchange between perfusate and tissue, rapid MBF variations are observable if the exchange is fast on the scale of our frame rate. The significant cyclic perfusion change found in this study supports the assumption of sufficiently fast capillary/tissue spin exchange in myocardium. In the literature, there are few studies regarding changes in spin exchange rates in the context of cardiovascular pathologies. Since arterial spin labeling only uses water molecules, we can hypothesize that the fast exchange assumption should be preserved even across pathologies.
Rest versus stress
Changes in MBF across the cardiac cycle were significant in this study with, in general, a similar behavior between rest and stress experiments. Nevertheless, two major differences arose from these results. First, maximum MBF appeared sooner within the cardiac cycle at stress than at rest. This was the case not only within the averaged group results, but also for each studied animal. This observation is also reflected in the MPR variations, in which a drop was found at the end of the cardiac cycle, corresponding to the delay between both physiological states. Second, we found a mean group decrease of MBF from ED to ES (ΔMBFCycl) of 27.2 ± 9.3% at rest, which was significantly higher from ΔMBFCycl obtained during stress (17.8 ± 7.1%, P = 0.014).
In a study published by Kajiya et al. , the authors focused on the phase opposition of velocity waveforms between coronary arteries and veins. They introduced the unstressed volume (UV), which refers to a capacitance accommodating for blood supply during diastole without significant increase in UV pressure to propel venous outflow. But UV can become saturated when diastolic arterial inflow is higher than normal, such as in pharmacologic stress situations causing vasodilation. In this case, the pressure of ordinary capacitance exceeds the venous pressure causing diastolic venous outflow. Such theory might explain why the maximum MBF was found to appear sooner within the cardiac cycle during stress.
When administered intravenously, adenosine effect is mixed, endothelial dependent and non-dependent  and thus causes relaxation of smooth muscles as found inside the arterial walls. This action leads to vasodilation and increases arterial inflow along with a decreased resistance and an increased capacitance of microvessels. Here, the term vascular resistance refers to the impediment offered by the vascular bed to a unidirectional and constant blood flow. The concept of vascular impedance, however, involves the pressure-flow relations at each instantaneous timing change and refers to the impediment offered to flow at the input of a vascular bed where pulsatile flow is involved. Thus, vascular impedance takes into account not only the resistance but also dynamic effects and seems more appropriate in this context. It is defined as the ratio of arterial pressure and flow. By reaching the maximum coronary response using adenosine, microvessels are already subjected to very high blood flow but also to reduced arterial pressure , which in our experiments resulted in lower vascular impedance represented in our experiments by a lower ΔMBFCycl during stress.
Relationship between heart rate and cyclic MBF variation amplitude
There was a linear correlation when comparing heart rate and ΔMBFCycl amplitude during stress (Figures 6 and 7). However, no significant correlation was found at rest which might be explained by the lower sensitivity of ASL at rest compared to stress. This might also be due to an insufficient number of animals to reach significance. There is a dependence of ΔMBFCycl on heart rate under stress, which could be present at rest as well. A likely explanation for this behavior is that the vascular system of the heart has a capacitive property when submitted to pressure variations at its entry. Analogies between electrical circuits and the vascular system, such as the Windkessel model, have been often used for modeling approaches . In a similar way as electric current in a capacitor submitted to voltage variations, the pressure oscillations lead to flow oscillations whose amplitude depend on the frequency of the pressure variations, here given by the heart rate.
Absolute MBF and MPR quantifications
It should be noted that all MBF values reported in this study, especially at rest, cannot be considered as physiological baseline values due to the influence of isoflurane anesthesia on coronary vasodilation [31, 41]. Resting MBF within the group was 5.5 ± 0.6 mL g-1 min-1 and therefore similar to previously reported values in rats, which are in the range of 3.83−7.89 mL g-1 min-1, using CMR [18, 23, 42] or fluorescent microspheres [23, 43, 44]. MBF obtained during stress was 12.8 ± 0.7 mL g-1 min-1, and mean MPR was 2.4 ± 0.2. To our knowledge, MPR has not been reported to be more than 2.5 in rodents and our results were in accordance with earlier obtained values using CMR in comparison with fluorescent microspheres  or contrast-enhanced ultrasound [45, 46]. However, MPR in this experiment was found higher than that reported by Croteau et al. (1.4 ± 0.5) in a study using positron emission tomography and 13N-NH3 as radiotracer . In their article, the authors discussed the fact that vasodilation induced by the pharmacologic stress was not maximal in their experimental design with the commonly used clinical injection rate of 140 μg kg-1 min-1. In our first experimental trial, we also encountered difficulties to obtain a reproducible stress within a group of healthy rats using this typical infusion rate. Increasing twice the infusion rate and using a small injection line with inner diameter of 0.58 mm resulted in obtaining a maximal MPR more reliably.
Study comparison and limitations
In humans, Radjenovic et al.  found no differences between systolic and diastolic MBF at rest and a significant MBF increase in diastole during adenosine-induced stress. The authors discussed that cyclic MBF changes may be present also at rest in humans but that the precision of the first-pass gadolinium-based MRI method used in their work might not be sufficient to assess them.
MBF quantification errors may be present by various sources such as respiratory motion. No respiratory gating was done with this sequence. Likely owing to the high number of signal averages performed with cine-ASL along with randomization of breathing events with myocardial contraction, impact of respiratory motion on MBF quantification and cyclic behavior is expected to be negligible. Although the heart lengthens longitudinally from systole to diastole, through-plane motion affects both control and tag scans in the same way so that their effects should cancel when calculating the difference maps. Effect of through-plane motion on the magnetization steady-state should also be relatively small considering the rather thick (2 mm) slices. With cine-ASL, temporal placement of the labeling pulse cannot be corrected once the acquisition started. Particular attention was paid to heart rate stability so as to ensure regular labeling efficiency throughout the entire 12 minutes scan. As shown in Table 1, the heart rate variation in this group of healthy rats did not exceed 15% at rest and 20% during stress, and one can therefore expect limited impact of label pulse timing imperfections on MBF quantifications. Heart rate variations can, however, be larger in models of cardiac pathology and would necessitate individual labeling efficiency determination.
Another potential source of error, especially when assessing dynamic perfusion changes over the cardiac cycle, are the short waiting periods before each ECG trigger that lead to interruptions of the longitudinal magnetization steady-state. Again, these interruptions affect tag and control scans in the same way, but in addition, the acquisition was extended over two cardiac cycles instead of one, such that the second cycle was not affected by the steady-state interruption. Figures 3 and 4 indeed prove good concordance between the dynamics of the initial phases of the first and the second observed cycle. Finally, it cannot be excluded that tissue and/or blood T1 as well as T2* are submitted to changes throughout the cardiac cycle. T1 might in addition be altered in a pathological context when tissue structural changes play a role. As a limitation common to many ASL techniques, cine-ASL would indeed require the acquisition of a T1 map to be fully quantitative. In terms of cyclic changes, the influence of T1 variations across the cardiac cycle can, however, be expected to remain small and would only linearly affect the resulting MBF values. Given the short echo time used in the cine-ASL sequence, the influence of T2* on the signal is also expected to remain small.
To our knowledge, this study presents the first attempt to quantify and monitor cyclic MBF variation in vivo in rats using a new and fully non-invasive arterial spin labeling CMR protocol. These experiments were carried out in a group of healthy rats with a focus on cyclic changes at rest and during adenosine-induced stress. Our observations clearly show that, in rodents, blood flow is phase-dependent and differs significantly between systole and diastole, for both rest and stress conditions. We discussed rat myocardial perfusion changes in response to adenosine as an infused vasodilator in detail and provided complementary information on how cyclic perfusion changes are affected by the vasodilator. Also, maximum absolute blood flow occurs earlier under stress conditions than under resting conditions. This might become relevant when calculating perfusion reserve based on two measurements with equal trigger delays. Considering the nature of coronary hemodynamics, cyclic MBF changes may reveal new physiologic information because they are a function of coronary flow, myocardial contraction and microvascular condition. This technique may be useful in a pathologic context to study microvascular defects in non-ischemic heart disease models.
Myocardial blood volume
Myocardial blood flow
Mean transit time
Arterial spin labeling
Cardiovascular magnetic resonance
Fast low-angle shot
Myocardial perfusion reserve
Relative variation of myocardial blood flow across the cardiac cycle.
Gregg DE, Sabiston DC: Effect of cardiac contraction on coronary blood flow. Circulation. 1957, 15: 14-20. 10.1161/01.CIR.15.1.14.
Hoffman JI, Spaan JA: Pressure-flow relations in coronary circulation. Physiol Rev. 1990, 70: 331-390.
Morgenstern C, Holjes U, Arnold G, Lochner W: The influence of coronary pressure and coronary flow on intracoronary blood volume and geometry of the left ventricle. Pflugers Arch. 1973, 340: 101-111. 10.1007/BF00588169.
Spaan JA, Breuls NP, Laird JD: Diastolic-systolic coronary flow differences are caused by intramyocardial pump action in the anesthetized dog. Circ Res. 1981, 49: 584-593. 10.1161/01.RES.49.3.584.
Klassen GA, Barclay KD, Wong R, Paton B, Wong AY: Red cell flux during the cardiac cycle in the rabbit myocardial microcirculation. Cardiovasc Res. 1997, 34: 504-514. 10.1016/S0008-6363(97)00062-X.
Ashikawa K, Kanatsuka H, Suzuki T, Takishima T: Phasic blood flow velocity pattern in epimyocardial microvessels in the beating canine left ventricle. Circ Res. 1986, 59: 704-711. 10.1161/01.RES.59.6.704.
Kajiya F, Tsujioka K, Ogasawara Y, Hiramatsu O, Wada Y, Goto M, Yanaka M: Analysis of the characteristics of the flow velocity waveforms in left atrial small arteries and veins in the dog. Circ Res. 1989, 65: 1172-1181. 10.1161/01.RES.65.5.1172.
Yada T, Hiramatsu O, Kimura A, Goto M, Ogasawara Y, Tsujioka K, Yamamori S, Ohno K, Hosaka H, Kajiya F: In vivo observation of subendocardial microvessels of the beating porcine heart using a needle-probe videomicroscope with a CCD camera. Circ Res. 1993, 72: 939-946. 10.1161/01.RES.72.5.939.
Wu EX, Tang H, Wong KK, Wang J: Mapping cyclic change of regional myocardial blood volume using steady-state susceptibility effect of iron oxide nanoparticles. J Magn Reson Imaging. 2004, 19: 50-58. 10.1002/jmri.10426.
Rodriguez I, Ennis DB, Wen H: Noninvasive measurement of myocardial tissue volume change during systolic contraction and diastolic relaxation in the canine left ventricle. Magn Reson Med. 2006, 55: 484-490. 10.1002/mrm.20786.
Radjenovic A, Biglands JD, Larghat A, Ridgway JP, Ball SG, Greenwood JP, Jerosch-Herold M, Plein S: Estimates of systolic and diastolic myocardial blood flow by dynamic contrast-enhanced MRI. Magn Reson Med. 2010, 64: 1696-1703. 10.1002/mrm.22538.
Coolen BF, Moonen RP, Paulis LE, Geelen T, Nicolay K, Strijkers GJ: Mouse myocardial first-pass perfusion MR imaging. Magn Reson Med. 2010, 64: 1658-1663. 10.1002/mrm.22588.
Makowski M, Jansen C, Webb I, Chiribiri A, Nagel E, Botnar R, Kozerke S, Plein S: First-pass contrast-enhanced myocardial perfusion MRI in mice on a 3-T clinical MR scanner. Magn Reson Med. 2010, 64: 1592-1598. 10.1002/mrm.22470.
Jogiya R, Makowski M, Phinikaridou A, Patel AS, Jansen C, Zarinabad N, Chiribiri A, Botnar R, Nagel E, Kozerke S, Plein S: Hyperemic stress myocardial perfusion cardiovascular magnetic resonance in mice at 3 Tesla: initial experience and validation against microspheres. J Cardiovasc Magn Reson. 2013, 15 (1): 62-10.1186/1532-429X-15-62.
van Nierop BJ, Coolen BF, Dijk WJ, Hendriks AD, de Graaf L, Nicolay K, Strijkers GJ: Quantitative first-pass perfusion MRI of the mouse myocardium. Magn Reson Med. 2013, 69: 1735-1744. 10.1002/mrm.24424.
Belle V, Kahler E, Waller C, Rommel E, Voll S, Hiller KH, Bauer WR, Haase A: In vivo quantitative mapping of cardiac perfusion in rats using a noninvasive MR spin-labeling method. J Magn Reson Imaging. 1998, 8: 1240-1245. 10.1002/jmri.1880080610.
Waller C, Hiller KH, Voll S, Haase A, Ertl G, Bauer WR: Myocardial perfusion imaging using a non-contrast agent MR imaging technique. Int J Cardiovasc Imaging. 2001, 17: 123-132. 10.1023/A:1010603217785.
Kober F, Iltis I, Izquierdo M, Desrois M, Ibarrola D, Cozzone PJ, Bernard M: High-resolution myocardial perfusion mapping in small animals in vivo by spin-labeling gradient-echo imaging. Magn Reson Med. 2004, 51: 62-67. 10.1002/mrm.10676.
Kober F, Iltis I, Cozzone PJ, Bernard M: Myocardial blood flow mapping in mice using high-resolution spin labeling magnetic resonance imaging: influence of ketamine/xylazine and isoflurane anesthesia. Magn Reson Med. 2005, 53: 601-606. 10.1002/mrm.20373.
Streif JU, Nahrendorf M, Hiller KH, Waller C, Wiesmann F, Rommel E, Haase A, Bauer WR: In vivo assessment of absolute perfusion and intracapillary blood volume in the murine myocardium by spin labeling magnetic resonance imaging. Magn Reson Med. 2005, 53: 584-592. 10.1002/mrm.20327.
Nahrendorf M, Streif JU, Hiller KH, Hu K, Nordbeck P, Ritter O, Sosnovik D, Bauer L, Neubauer S, Jakob PM, Ertl G, Spindler M, Bauer WR: Multimodal functional cardiac MRI in creatine kinase-deficient mice reveals subtle abnormalities in myocardial perfusion and mechanics. Am J Physiol Heart Circ Physiol. 2006, 290: H2516-2521. 10.1152/ajpheart.01038.2005.
Vandsburger MH, Janiczek RL, Xu Y, French BA, Meyer CH, Kramer CM, Epstein FH: Improved arterial spin labeling after myocardial infarction in mice using cardiac and respiratory gated look-locker imaging with fuzzy C-means clustering. Magn Reson Med. 2010, 63: 648-657. 10.1002/mrm.22280.
Jacquier A, Kober F, Bun S, Giorgi R, Cozzone PJ, Bernard M: Quantification of myocardial blood flow and flow reserve in rats using arterial spin labeling MRI: comparison with a fluorescent microsphere technique. NMR Biomed. 2011, 24 (9): 1047-1053. 10.1002/nbm.1645.
Campbell-Washburn AE, Price AN, Wells JA, Thomas DL, Ordidge RJ, Lythgoe MF: Cardiac arterial spin labeling using segmented ECG-gated Look-Locker FAIR: variability and repeatability in preclinical studies. Magn Reson Med. 2013, 69 (1): 238-247. 10.1002/mrm.24243.
Northrup BE, McCommis KS, Zhang HS, Ray S, Woodard PK, Gropler RJ, Zheng J: Resting myocardial perfusion quantification with CMR arterial spin labeling at 1.5 T and 3.0 T. J Cardiovasc Magn Reson. 2008, 10: 53-10.1186/1532-429X-10-53.
Zun ZH, Varadarajan P, Pai RG, Wong EC, Nayak KS: Arterial spin labeled CMR detects clinically relevant increase in myocardial blood flow with vasodilation. Jacc-Cardiovasc Imag. 2011, 4: 1253-1261.
Troalen T, Capron T, Cozzone PJ, Bernard M, Kober F: Cine-ASL: A steady-pulsed arterial spin labeling method for myocardial perfusion mapping in mice. Part I. Experimental study. Magn Reson Med. 2013, 70: 1389-1398. 10.1002/mrm.24565.
Capron T, Troalen T, Cozzone PJ, Bernard M, Kober F: Cine-ASL: A steady-pulsed arterial spin labeling method for myocardial perfusion mapping in mice. Part II. Theoretical model and sensitivity optimization. Magn Reson Med. 2013, 70: 1399-1408. 10.1002/mrm.24588.
Troalen T, Capron T, Cozzone PJ, Bernard M, Kober F: Cyclic Variation of Myocardial Blood Flow Assessed with cine-ASL. Abstracts of ESMRMB (European Society for Magnetic Resonance in Medicine and Biology) 2012, the 29th Annual Scientific Meeting. Lisbon, Portugal. October 4–6, 2012. MAGMA. 2012, 25 (Suppl 1): 94-
Takuma S, Suehiro K, Cardinale C, Hozumi T, Yano H, Shimizu J, Mullis-Jansson S, Sciacca R, Wang J, Burkhoff D, Di Tullio MR, Homma S, et al: Anesthetic inhibition in ischemic and non ischemic murine heart: comparison with conscious echocardiographic approach. Am J Physiol Heart Circ Physiol. 2001, 280: H2364-2370.
Iltis I, Kober F, Dalmasso C, Lan C, Cozzone PJ, Bernard M: In vivo assessment of myocardial blood flow in rat heart using magnetic resonance imaging: effect of anesthesia. J Magn Reson Imaging. 2005, 22: 242-247. 10.1002/jmri.20352.
Waller C, Kahler E, Hiller KH, Hu K, Nahrendorf M, Voll S, Haase A, Ertl G, Bauer WR: Myocardial perfusion and intracapillary blood volume in rats at rest and with coronary dilatation: MR imaging in vivo with use of a spin-labeling technique. Radiology. 2000, 215: 189-197. 10.1148/radiology.215.1.r00ap07189.
Deichmann R, Haase A: Quantification of T1 Values by Snapshot-Flash Nmr Imaging. J Magn Reson. 1992, 96: 608-612.
Davies JE, Whinnett ZI, Francis DP, Manisty CH, Aguado-Sierra J, Willson K, Foale RA, Malik IS, Hughes AD, Parker KH, Mayet J: Evidence of a dominant backward-propagating “suction” wave responsible for diastolic coronary filling in humans, attenuated in left ventricular hypertrophy. Circulation. 2006, 113: 1768-1778. 10.1161/CIRCULATIONAHA.105.603050.
Chilian WM, Eastham CL, Layne SM, Marcus ML: Small vessel phenomena in the coronary microcirculation: phasic intramyocardial perfusion and coronary microvascular dynamics. Prog Cardiovasc Dis. 1988, 31: 17-38. 10.1016/0033-0620(88)90009-6.
Kajiya M, Hirota M, Inai Y, Kiyooka T, Morimoto T, Iwasaki T, Endo K, Mohri S, Shimizu J, Yada T, Ogasawara Y, Naruse K, Ohe T, Kajiya F: Impaired NO-mediated vasodilation with increased superoxide but robust EDHF function in right ventricular arterial microvessels of pulmonary hypertensive rats. Am J Physiol Heart Circ Physiol. 2007, 292: H2737-2744. 10.1152/ajpheart.00548.2006.
Kajiya F, Yada T, Hiramatsu O, Ogasawara Y, Inai Y, Kajiya M: Coronary microcirculation in the beating heart. Med Biol Eng Comput. 2008, 46: 411-419. 10.1007/s11517-008-0335-x.
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-249. 10.1002/(SICI)1099-1492(199706/08)10:4/5<237::AID-NBM475>3.0.CO;2-X.
Sato A, Terata K, Miura H, Toyama K, Loberiza FR, Hatoum OA, Saito T, Sakuma I, Gutterman DD: Mechanism of vasodilation to adenosine in coronary arterioles from patients with heart disease. Am J Physiol Heart Circ Physiol. 2005, 288: H1633-1640.
Kissling G, Blickle B, Pascht U: Modified heart-lung preparation for the evaluation of systolic and diastolic coronary flow in rats. Am J Physiol Heart Circ Physiol. 2000, 278: H277-284.
Crystal GJ: Vasomotor effects of isoflurane in the coronary circulation. Anesthesiology. 1996, 84: 1516-1518.
Iltis I, Kober F, Desrois M, Dalmasso C, Lan C, Portha B, Cozzone PJ, Bernard M: Defective myocardial blood flow and altered function of the left ventricle in type 2 diabetic rats: a noninvasive in vivo study using perfusion and cine magnetic resonance imaging. Invest Radiol. 2005, 40: 19-26.
Wicker P, Tarazi RC: Coronary blood flow measurements with left atrial injection of microspheres in conscious rats. Cardiovasc Res. 1982, 16: 580-586. 10.1093/cvr/16.10.580.
Debaene B, Goldfarb G, Braillon A, Jolis P, Lebrec D: Effects of ketamine, halothane, enflurane, and isoflurane on systemic and splanchnic hemodynamics in normovolemic and hypovolemic cirrhotic rats. Anesthesiology. 1990, 73: 118-124. 10.1097/00000542-199007000-00017.
Oshita A, Ohmori K, Yu Y, Kondo I, Takeuchi H, Takagi Y, Wada Y, Yukiiri K, Mizushige K, Kohno M: Myocardial blood flow measurements in rats with simple pulsing contrast echocardiography. Ultrasound Med Biol. 2002, 28: 459-466. 10.1016/S0301-5629(02)00481-7.
Raher MJ, Thibault H, Poh KK, Liu R, Halpern EF, Derumeaux G, Ichinose F, Zapol WM, Bloch KD, Picard MH, Scherrer-Crosbie M: In vivo characterization of murine myocardial perfusion with myocardial contrast echocardiography: validation and application in nitric oxide synthase 3 deficient mice. Circulation. 2007, 116: 1250-1257. 10.1161/CIRCULATIONAHA.107.707737.
Croteau E, Benard F, Bentourkia M, Rousseau J, Paquette M, Lecomte R: Quantitative myocardial perfusion and coronary reserve in rats with 13N-ammonia and small animal PET: impact of anesthesia and pharmacologic stress agents. J Nucl Med. 2004, 45: 1924-1930.
Grant support: this work was supported by the French ‘Agence Nationale de la Recherche’, ANR Blanc 2008 ‘QASAREM’ and by the Centre National de la Recherche Scientifique “CNRS-UMR 7339”.
TT was a recipient of a PhD fellowship “CIFRE” from Siemens Healthcare, France.
TT carried out the sequence design, conceived the protocol, performed the acquisition and analysis of data, and drafted the manuscript. TC created the theoretical model for perfusion quantification and participated in the revision of the manuscript. MB helped in the interpretation of data. FK participated in the study design, managed coordination and participated in the revision of the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Troalen, T., Capron, T., Bernard, M. et al. In vivo characterization of rodent cyclic myocardial perfusion variation at rest and during adenosine-induced stress using cine-ASL cardiovascular magnetic resonance. J Cardiovasc Magn Reson 16, 18 (2014). https://doi.org/10.1186/1532-429X-16-18