Home | About | Journals | Submit | Contact Us | Français |

**|**J Cardiovasc Magn Reson**|**v.14(1); 2012**|**PMC3519741

Formats

Article sections

- Abstract
- Background
- Methods
- Results
- Discussion
- Conclusion
- Appendix
- Competing interests
- Authors’ contributions
- References

Authors

Related links

J Cardiovasc Magn Reson. 2012; 14(1): 73.

Published online 2012 October 24. doi: 10.1186/1532-429X-14-73

PMCID: PMC3519741

Sumeda Abeykoon: gro.cmhcc@nookyeba.ademus; Michelle Sargent: gro.cmhcc@tnegras.ellehcim; Janaka P Wansapura: gro.cmhcc@arupasnaw.akanaj

Received 2012 May 16; Accepted 2012 October 11.

Copyright ©2012 Abeykoon et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article has been cited by other articles in PMC.

In the conventional approach to arterial spin labeling in the rodent heart, the relative difference in the apparent *T*_{1} relaxation times corresponding to selective and non-selective inversion is related to perfusion via a two compartment model of tissue. But accurate determination of *T*_{1} in small animal hearts is difficult and prone to errors due to long scan times and high heart rates. In this study we introduce the theoretical frame work for an alternative method (*SI-method*) based purely on the signal intensity of slice-select and non-select inversion recovery images at a single inversion time at short repetition time.

A modified Bloch equation was solved to derive perfusion as a function of signal intensity of flow sensitized segmented gradient echo acquisitions. A two compartment fast exchanging model of tissue was assumed. To test the new technique first it was implemented on a flow phantom and then it was compared with the conventional *T*_{1} method in an *in vivo* study of healthy C57BL/6 mice (n=12). Finally the *SI-method* was used in comparison to a Late Gadolinium Enhanced (LGE) method to qualitatively and quantitatively assess perfusion deficits in an ischemia-reperfusion mouse model (n=4).

The myocardial perfusion of healthy mice obtained by the *SI-method*, 5.6 ± 0.5 ml/g/min, (mean ± standard deviation) was similar (p=0.38) to that obtained by the conventional method, 5.6 ± 0.3 ml/g/min. The variance in perfusion within the left ventricle was less for the *SI-method* than that for the conventional method (p<0.0001). The mean percentage standard deviation among repeated measures was 3.6%. The LGE regions of the ischemia reperfusion model were matched with regions of hypo-perfusion in the perfusion map. The average perfusion in the hypo perfused region among all four IR mice was 1.2 ± 0.9 ml/g/min and that of the remote region was 4.4 ± 1.2 ml/g/min.

The proposed signal intensity based ASL method with a segmented acquisition scheme allows accurate high resolution perfusion mapping in small animals. It’s short scan time, high reproducibility and ease of post process makes it a robust alternative to the conventional ASL technique that relies on T_{1} measurements.

To quantify perfusion, the Arterial Spin Labeling (ASL) technique uses water protons as a freely diffusible contrast agent. The contrast due to perfusion is obtained by preparation of magnetization by selective and non-selective inversion. Slice-selective inversion causes signal enhancement due to in-flow of thermally stabilized spins whereas global inversion is insensitive to flow. In the cardiac implementation of ASL, spin labeling is generally done within the detection slice. One approach of ASL is to mathematically relate the relative difference in the apparent relaxation times corresponding to selective and non-selective inversion, to perfusion via a two compartment model of tissue [1]. In this case perfusion, *P*, is given by:

$$P=\frac{\lambda}{{T}_{1c}}\left(\frac{{T}_{1g}}{{T}_{1s}}-1\right)$$

(1)

Where, T_{ls} and T_{lg} are the apparent longitudinal relaxation times for slice select and non-select inversion respectively. T_{lc} is the longitudinal relaxation time of blood. This ASL method, hereafter referred to as the *T1-method*, has been used to quantify myocardial perfusion in humans [2,3] as well as in small animal models [4-7]. Accurate determination of T_{l} relaxation times require long scan time to allow for full relaxation of the longitudinal magnetization. Long scan times in turn leads to varied inversion times for different read outs of the k-space in segmented acquisitions in the heart, causing errors in T_{l} maps and therefore in perfusion maps.

In this study we propose an alternative ASL method, which we will call the *SI-method*, based purely on the signal intensity of slice-select and non-select inversion recovery images at a single inversion time at short repetition time. A steady state image is also acquired to normalize the receiver characteristics. A new expression for perfusion is derived for this acquisition scheme assuming a two compartment model of tissue [1]. Similar signal intensity based methods with single inversion times have been used in the past in brain and cardiac perfusion imaging [8-11]. A theoretical frame work for ASL with short repetition time has also been introduced by Pell et al. [11]. What is different in the current study is the use of short repetition time approach combined with a dedicated two-compartment theory.

We adapted the two compartment model proposed by Bauer et al. [1] and derived perfusion as a function of signal intensity of slice select and non-select inversion acquisitions by solving a modified Bloch equation. To test the new technique we first implemented it on a flow phantom with two compartments with no exchange between them. Since it is difficult to accurately quantify spin exchange rates in materials, this phantom experiment allowed us to know perfusion exactly given a specific rate of flow. Second, we compared the proposed *SI-method* with the conventional *T1-method* by carrying out *in vivo* myocardial perfusion measurements in healthy mice. Finally the *SI-method* was used in comparison to a LGE method to qualitatively and quantitatively assess perfusion deficits in an ischemia-reperfusion mouse model.

The data acquisition scheme for this method consists of three scans: (i) Slice selective inversion-TI-acquisition, (ii) Non-select inversion-TI-acquisition, (iii) steady state acquisition with TR set to TI. The method as it was implemented in mice is depicted in Figure Figure1.1. A segmented k-space acquisition is assumed. A finite TR is assumed for (i) and (ii). We will derive an expression for magnetization for each of these cases in terms of perfusion in a two compartment model of tissue (Figure (Figure2)2) where we assume myocardial tissue to be consisting of homogenous intra and extra-capillary space with respect to relaxation time. The rate of change of magnetization in the intra-capillary region m_{c}(t) depends on loss of magnetization due to T_{1} relaxation, exchange of magnetic spins between intra and extra-capillary region and gain of magnetization due to in flow of blood. The magnetization in the extra-capillary region m_{e}(t) is unaffected by blood flow. The spin dynamics of the two compartments is described by the following modified Bloch equations.

$$\begin{array}{c}\frac{d{m}_{c}\left(t\right)}{dt}=\left(\frac{{m}_{c}\left(0\right)-{m}_{c}\left(t\right)}{{T}_{1c}}\right)+\frac{F{m}_{p}}{{V}_{c}}-\frac{F{m}_{c}\left(t\right)}{{V}_{c}}\hfill \\ \phantom{\rule{4.5em}{0ex}}-{K}_{c}{m}_{c}\left(t\right)+{K}_{e}{m}_{e}\left(t\right)\hfill \end{array}$$

(2)

$$\frac{d{m}_{e}\left(t\right)}{dt}=\left(\frac{{m}_{e}\left(0\right)-{m}_{e}\left(t\right)}{{T}_{1e}}\right)+{K}_{e}{m}_{c}\left(t\right)-{K}_{e}{m}_{e}\left(t\right)$$

(3)

Where T_{lx}, m_{x}(t), V_{x} and K_{x} are longitudinal relaxation time, magnetization, volume, spin exchange rate respectively with the subscript *x=c* for intra-capillary region and *x=e* for extra-capillary region. m_{p}(t) is the magnetization of the incoming blood. The total magnetization of tissue, M(t), is calculated as,

$$M\left(t\right)=\left(\frac{RBV}{\lambda}\right){m}_{c}\left(t\right)+\left(1-\frac{RBV}{\lambda}\right){m}_{e}\left(t\right)$$

(4)

where *λ* is the blood tissue water partition coefficient [4], $RBV=\frac{{V}_{c}}{{V}_{t}}$ and *V*_{T}=*V*_{c}+*V*_{e}. Equations 13 describe the time evaluation of the total magnetization of tissue and can be solved for a given input magnetization, m_{p}(t), depending on the type of spin preparation. For slice-select inversion, m_{p}(t) is the equilibrium magnetization of capillary blood hence is equal to the constant m_{c}(0) in magnitude. For non-select inversion, at steady state, ${m}_{p}\left(t\right)={m}_{c}\left(0\right)\left(1-2{e}^{-\frac{TI}{{T}_{1c}}}\right)+{m}_{c}\left(0\right){e}^{-\frac{TR}{{T}_{1c}}}$, where *TI* and *TR* are the inversion and repetition times respectively (see Appendix).

The total magnetization for each spin preparation can be written in matrix from as follows:

$$\frac{d}{dt}\left({M}_{s}\left(t\right){e}^{-Rt}\right)=\frac{P}{RBV}{e}^{-Rt}{m}_{c}\left(0\right)+{e}^{-Rt}\left(\begin{array}{c}\frac{{m}_{c}\left(0\right)}{{T}_{1c}}\\ \frac{{m}_{e}\left(0\right)}{{T}_{1e}}\end{array}\right)$$

(5)

for slice-select inversion and

$$\begin{array}{c}\frac{d}{dt}\left({M}_{g}\left(t\right){e}^{-Rt}\right)=\frac{P}{RBV}{m}_{c}\left(0\right)\left(1-2{e}^{-\left(\frac{TI}{{T}_{1c}}\right)}+{e}^{-\left(\frac{TR}{{T}_{1c}}\right)}\right){e}^{-Rt}\\ \phantom{\rule{9.3em}{0ex}}+{e}^{-Rt}\left(\begin{array}{c}\frac{{m}_{c}\left(0\right)}{{T}_{1c}}\\ \frac{{m}_{e}\left(0\right)}{{T}_{1e}}\end{array}\right)\hfill \end{array}$$

(6)

for non-select inversion, where

$$R=\left[\stackrel{-\left(\frac{1}{{T}_{1c}}+\frac{F}{{V}_{c}}+{K}_{c}\right)}{{K}_{e}}\stackrel{{K}_{c}}{-\left(\frac{1}{{T}_{1e}}+{K}_{e}\right)}\right]$$

(none1)

Equations (5) and (6) can be solved to derive the magnetization for slice-select M_{s}(t=TI) and non-select inversion M_{g}(t=TI) using the initial condition that all magnetization within the selected slice is inverted at t=0 in both preparation schemes. The term *R* can be simplified following a series of steps proposed by Bauer et al. [1,12] assuming fast exchange between intra and extra capillary regions i.e. Kc≈Ke≈∞. This assumption reduces the complexity of the model but is appropriate for myocardial tissue that has a high exchange rate compared to other tissues of the body [12]. The difference in the signal intensities between the two preparation schemes can now be defined as *ΔM*=*M*_{g}(*t*=*TI*)−*M*_{s}(*t*=*TI*). By combining the simplified expression for *R* and the solutions to equations (5) and (6) we obtain,

$$\Delta M=\frac{2P{m}_{c}\left(0\right)\left({e}^{-\left(\frac{TI}{{T}_{1c}}\right)-{e}^{-\left(\mu ,T,I\right)}}\right)}{\left(\mu -\frac{1}{{T}_{1c}}\right)\lambda}+\frac{P{m}_{c}\left(0\right){e}^{-\left(\frac{TR}{{T}_{1c}}\right)}({e}^{-\mu TI}-1)}{\mu \lambda}$$

(7)

Where, $\mu =\left(\frac{RBV}{\lambda}\right)\frac{1}{{T}_{1c}}+\left(1-\frac{RBV}{\lambda}\right)\frac{1}{{T}_{1e}}+\frac{P}{\lambda}$

To remove the dependency on receiver characteristics, we normalize ΔM by the magnetization of a steady state gradient echo acquisition (M_{s}) with no spin preparation and with *TR* set to the same *TI* as in equation (7) (Figure (Figure1).1). Since echo time of this acquisition is small (~ 1ms) we assume that *M*_{s} is not affected by perfusion. We also assume equilibrium magnetization of intra and extra-capillary regions to be similar i.e. *m*_{c}(0)≈*m*_{e}(0). Thus we obtain,

$${M}_{s}={m}_{c}\left(0\right)\left(1-{e}^{-{\mu}_{0}TI}\right)$$

(8)

Where, ${\mu}_{0}=\left(\frac{RBV}{\lambda}\right)\frac{1}{{T}_{1c}}+\left(1-\frac{RBV}{\lambda}\right)\frac{1}{{T}_{1e}}=\frac{1}{{T}_{1}}$

Here *T*_{1} is the myocardial longitudinal relaxation time without the in-flow effects. By dividing equation (7) by (8) we obtain an expression for the normalized magnetization difference, ΔM/M_{s} Then by taking the first order Taylor expansion of ΔM/M_{s} about *TI* we obtain,

$$P=\frac{\frac{\lambda \Delta M}{{M}_{s}}}{\left(1-{e}^{-\frac{TR}{{T}_{1c}}}+{e}^{-\frac{TI}{{T}_{1c}}}\right){T}_{1}}$$

(9)

Perfusion is derived for a non-exchangeable model of two compartments based on the same acquisition method as described above. The aim is to test the *SI-method* on a flow phantom. Since water exchange rates of tubing material is difficult to assess we chose material with zero permeability. This allows perfusion to be known exactly in terms of volume of fluid / volume of phantom / minute (i.e. min^{-1}), given a specific flow rate. The phantom will have two compartments, intra- and extra capillary, but with no exchange between them. Accordingly the *SI-method* is modified for this scenario by making K_{c}=0, K_{e}=0 in equations (2) and (3). The rest of the derivation is similar to as in the previous section.

The difference in the signal intensities between the two preparation schemes, ΔM, is obtained by first order Taylor expansion as in the previous section. By solving for perfusion we get,

$$P=\frac{\lambda \Delta M}{{m}_{c}\left(0\right)\left(2{e}^{-\frac{TI}{{T}_{1c}}}-{e}^{-\frac{TR}{{T}_{1c}}}\right)TI}$$

(10)

Here we do not normalize *ΔM* by the signal intensity of a steady state acquisition, instead measure m_{c}(0)*,* the equilibrium magnetization of the capillary fluid*,* using a spin echo sequence with long *TR*.

The aim is to compare the results of the *T1-method* with that of the *SI-method* using the flow phantom described above. The *T1-method* was originally developed by Bauer et al. [4] for a fast exchanging two compartment model of tissue. Here we modify it to reflect no exchange between compartments. Again we begin with setting K_{c}=0, K_{e}=0 in equations (2) and (3). In the *T1-method* we are interested in the apparent longitudinal relaxation times and not the magnitudes of magnetization. To make the derivation of relaxation times easier we normalize the magnetization of intra and extra-capillary regions at the time of inversion to one. The equilibrium magnetization on the other hand is set to zero i.e. m_{c}(0)=m_{e}(0)=0 in equations (2) and (3). The total magnetization of tissue, M(t), is given by equation (4) as before. Modified Equations 24 can be solved for a given input magnetization, m_{p}(t), depending on the type of spin preparation. According to the normalizing scheme introduced here the input magnetization, m_{p}(t)*,* for slice-select inversion is zero and ${m}_{p}\left(t\right)={e}^{-\frac{t}{{T}_{1c}}}$, for non-select inversion.

By taking into account the mean relaxation time approximation that *T*_{1}=∫_{0}^{∞}M(*t*)*dt*, [4,13], we obtain the following expressions for apparent relaxation times

$${T}_{1s}=\frac{RBV}{\lambda}\left(\frac{RBV{T}_{1c}}{P{T}_{1c}+RBV}\right)+\left(1-\frac{RBV}{\lambda}\right){T}_{1e}$$

(11)

for slice-select inversion, and

$${T}_{1g}=\frac{RBV}{\lambda}{T}_{1c}+\left(1-\frac{RBV}{\lambda}\right){T}_{1e}$$

(12)

for non-select inversion. From equations (11) and (12) we derive an expression for perfusion in a two compartment model of tissue with no exchange.

$$P=\frac{1}{{T}_{1c}\left\{\frac{1}{\lambda}\left(\frac{{T}_{1c}}{{T}_{1g}-{T}_{1s}}\right)-\frac{1}{RBV}\right\}}$$

(13)

The *RBV* in the context of the said phantom is equal to the volume ratio V_{c}/V_{T}, where *V*_{c} is the volume of capillary tube and *V*_{T} is the total volume of the phantom as measured from a cross section of the phantom.

A phantom was prepared by placing 12 loops of capillary tube (Cole Parmer # EW 06492–02) with inner diameter 0.51 mm, inside a 3cc syringe. One end of the capillary tube was fitted to a syringe infusion pump (Model R-99E, Razel Scientific Instruments, St. Albans, VT). The space between capillary tubes was filed with a mixture of Gadolinium and Poly Vinyl Alcohol solution. Diluted Gd solution was used as the capillary liquid. Intra and extra-capillary liquids were prepared to have approximately similar longitudinal relaxation times of blood and tissue at 7 Tesla (*T*_{1-intra} = 1896 ms, *T*_{1-extra} = 1344 ms). Perfusion was calculated using equation (10) for the *SI-method* and equation (13) for the *T1-methods* for the following flow velocities 0.03, 0.07, 0.09, 0.1, 0.18, 0.3, 0.35 cm s^{-1}. The average signal intensity of a region covering the phantom was considered for calculation of *ΔM* and for curve fitting of *T*_{1}.

Cardiovascular magnetic resonance (CMR) was performed on a 7T Bruker Biospec system using an inversion recovery FLASH sequence consisting of a sinc3 excitation pulse and a hyperbolic secant pulse (sech) for inversion. For the *T1-method*, slice select inversion and non-select inversion prepared images were acquired at a constant repetition time, varying the inversion time, *TI*, in twelve steps between 50 and 5500 ms. *TR/TE* = 6000 ms/ 2.6 ms field of view = 2.5×2.5 cm^{2}, flip angle = 30°, slice thickness = 4 mm, a matrix size = 256×64. For the *SI-method*, TR = 1000 ms and *TI* = 100 ms, flip angle = 90°, All other imaging parameters were set similar to that of the *T1-method*. Imaging was performed perpendicular to the direction of flow. For both methods, the average signal intensity of a region of interest covering a cross section of the phantom was considered for calculations. *T*_{1s} and *T*_{1g} were calculated by fitting to a mono exponential function of the form $A{e}^{-\frac{TI}{{T}_{1}}}+B$. Spin density of intra (ρ_{i}) and extra-capillary (ρ_{e}) liquids were obtained by a spin echo pulse sequence with *TR*=20 sec and with all other imaging parameters including receiver gains fixed. Perfusion was calculated using equation (10) for *SI-method* and equation (13) for *T1-method* taking *λ*=(*ρ*_{i}+*ρ*_{e})/*ρ*_{i} and *m*_{c}(0)=*ρ*_{i}.

The study was conducted under a protocol approved by the Institutional Animal Care and Use Committee. To compare the proposed *SI-method* with the conventional *T1-method*, we measured myocardial perfusion of healthy ten weeks old C57BL/6 mice (n=12) using both methods. Mice were anesthetized with continuous inhaled Isoflurane (2% by volume) administered via a nose cone. Constant body temperature of 37°C was maintained using a thermocouple/heater system. Both SI and T1 methods were performed in the same session. To test the reproducibility of the *SI-method* three perfusion measurements were made on the same mouse. Each measurement was performed on a separate session (often on a different day). Three C57BL/6 mice were used for this study. Finally we used the *SI-method* to quantify perfusion in a mouse model of ischemia-reperfusion (IR). Four C57BL/6 mice were anesthetized with isoflurane and intubated through the mouth and ventilated. The left coronary artery was ligated with a silk ligature following thoracotomy. The animals were maintained in the ligated state for 30 minutes, after which the externalized silk was pulled to release the constriction and hence allowing full reperfusion. IR mice were scanned within 24 hours of surgery. The *in vivo* study was performed on a 7T Bruker Biospec system using a custom made solenoid coil having a single turn consisting of a wide conductive sheet (diameter =3.0 cm, length =3.5 cm). The sensitivity of the coil covered approximately 80%- 90% of the mouse’s body. In IR mice, in addition to myocardial perfusion measurement, LGE experiment was performed 30 minutes after injecting a bolus of Gd-DTPA (0.3-0.6 mmol/kg) intraperitonealy. *T1* weighted cine imaging was performed in the short axis using a segmented FLASH sequence to highlight enhancement. Slice thickness=1.0 mm, matrix size=256x256, in-plane resolution=117x117 μm2. TE/TR=3/5.2ms, flip angle=30º. Image acquisition was prospectively ECG gated using pediatric ECG probes attached to the paws. A pneumatic pillow was used for respiratory gating (SA Instruments, Stony Brooks, NY).

The acquisition schemes for *T1* and *SI methods* are shown in Figure Figure1.1. A gated, Look-Locker sequence [5] was implemented for the *T1-method*. A 5ms hyperbolic secant pulse was used for inversion. For the *T1-method*, for each inversion, the relaxation curve was sampled at twelve different inversion times at every third cardiac cycle i.e. *TI*=3×*n*×*RR*,*n*=1,.12 with *R*≈50×*RR*, that is approximately 6 seconds. For the *SI-method* the same sequence was used with data acquired at a single *TI* time of 10×RR with TR=16×RR, that is approximately 2 seconds. In both methods the acquisition block comprised a single excitation pulse followed by one line of k space. For the *SI-method* an additional steady state gradient echo image was acquired without the inversion pulse with TR=10×RR. All images were acquired in the short axis plane at end diastole. The following common imaging parameters were used: TE = 1.73 ms, field of view = 2.5×2.5 cm^{2}, slice thickness = 2.0 mm, matrix size = 128×64, flip angle = 90° and 10° were used for SI and T1 method respectively. For slice select inversion, the slice thickness of the inversion pulse was set to twice that of the excitation slice to eliminate effects from imperfect edges of both the inversion and the imaging pulse profiles

For the *in vivo* study we assumed myocardial tissue to consist of fast exchanging two compartments. Therefore we used equation (9) for the *SI-method*, and equation (1) for the *T1-method* for calculating perfusion. Following constant values were used: T_{1} of blood, T_{1c} = 1800 ms, T_{1} of myocardial tissue in the absence of flow effects, T_{1} = 1400 ms, λ=0.95 ml/g [14]. Perfusion was calculated on a pixel-by-pixel basis. The mean value was calculated from an ROI drawn within the left ventricle wall. Pixel values less than −5 ml/g/min and greater than 15 ml/g/min were set to zero and also threshold limit was set to obtain a noise free perfusion map. The apparent *T*_{1} values obtained by the Look Locker acquisition was corrected for saturation effects as previously described [5].

Results are expressed as mean ± standard deviation. Statistical differences were assessed with the unpaired 2-tailed Student’s t test for two experimental groups. A nonparametric test was applied when the data were not normally distributed. A 2-tailed P value of less than 0.05 was considered statistically significant.

The results of the phantom experiments are shown in Figure Figure3.3. In this case the perfusion is expressed as the flow rate per unit volume of phantom i.e. in (ml/min)/ (ml). Perfusion measured by both methods agrees well with that calculated using the flow rates in the capillary tubes. The deviation from the actual perfusion is greater for higher flow rates for both methods, the maximum being 16%. The Bland Altman plots show that the difference between the executed perfusion and the measured stay within the 95% confidence level in both cases.

The error bars shown in Figure Figure33 are derived by propagating errors of *P* in equations (10) and (13) for *SI* and *T1 methods*. Larger errors in the T1-method reflect greater uncertainty of *T*_{1s} and *T*_{1g} measurements.

Typical left ventricular perfusion maps of a normal mouse, by *SI* and *T1-method* are given in Figure Figure4.4. The average *T*_{1s} and *T*_{1g} for normal mice in our study was 1293 ± 80 ms and 1542 ± 88 ms respectively. The Table Table11 compares the perfusion values obtained by the two methods for each mouse. Perfusion was calculated from region of interest covering the entire left ventricle wall. The myocardial perfusion of healthy mice obtained by the *SI-method* is 5.6±0.5 ml/g/min, (mean ± standard deviation). This is not significantly different (p>0.05) to that obtained by the T1 method, 5.6±0.3 ml/g/min. The standard deviation given for each mouse in Table Table11 reflects the variation of *P* within the left ventricle for each mouse. The mean value of this standard deviation is significantly higher for the *T1-method* (2.3 ml/g/min) compared to that of the *SI-method* (1.6 ml/g/min), p<0.0001. As seen in Table Table22 the mean percentage standard deviation among repeated measures is less than 5% for the *SI-method*.

Figure Figure55 shows an example of LGE and perfusion in an IR mouse in the short axis plane. The enhancement in the lateral wall and the anterior septum are matched with regions of hypo-perfusion in the perfusion map. The average perfusion in the hypoperfused region among all four IR mice was 1.2 ± 0.9 ml/g/min and that of the remote region was 4.4 ± 1.2 ml/g/min.

In this study, a new ASL method to quantify myocardial perfusion is introduced based on the signal intensity of flow sensitized CMR. Following Bauer’s [1] frame work, myocardium is modeled as a two compartment with respect to *T*_{1}, which takes in to consideration that *T*_{1} of capillary blood is different to that of tissue. The imaging method consists of segmented, single gradient echo acquisitions prepared by slice select and non-select inversion pulses with finite repetition times followed by a gradient echo acquisition with no spin preparation. A simple expression for perfusion was derived by solving a modified Bloch equation for this model.

An important feature of our tissue model is that it consists of two compartments. In contrast, the Detre’s frame work [15] and its derivative ASL methods by Kwong [9], Buxton [8] and Zun [10] considers the tissue to be consists of a single compartment of blood. However, because *T*_{1} of blood is higher than that of tissue, taking tissue to be a single compartment of blood could lead to overestimation of perfusion. This is problematic if absolute values of perfusion are desired. On the other hand, we can arrive at the same expression for *P* as did Kwong and others by reducing our model to a single compartment model. To achieve this we take the two compartment non-exchange version of the *SI-method* (Equation 10) which is similar to one compartment of blood as far as *ΔM* is concerned. Then, by making *TR*→*∞* in equation 10 we obtain the same expression for *P* as in equation 5 of reference 9. Because of the two compartment model, the *SI-method* requires the knowledge of *T*_{1} of tissue with no in-flow effects (*T*_{1}) (Equation (9)) in addition to the *T*_{1} of blood (*T*_{1c}) whereas Detre’s method only requires the knowledge of *T*_{1c}. In this paper we have presented *P* calculated using literature values of *T*_{1} and *T*_{1c}. However, in our case these relaxation times could also be calculated from data acquired for the *T1-method*. Hence we calculated *T*_{1} and *T*_{1c} for each mouse in Table Table11 and used it to derive *P* using the *SI-method* (data not shown). The average *P* calculated using individually assessed *T*_{1} and *T*_{1c} was 5.7 ± 0.4 ml/g/min. Since this is not significantly different to the average *P* calculated using literature values (=5.6 ± 0.5 ml/g/min) it appears that the variability of perfusion among mice is greater than the uncertainty posed by *T*_{1} and *T*_{1c} in the *SI-method*.

The other important feature of the *SI-method* is that its k-space acquisition scheme is segmented as opposed to a single shot acquisition. Also, in deriving *P* no assumption was made about the length of the repetition time. This allows phase encoding steps to be acquired in multiple inversions with shorter TR times than it is required for full relaxation. Unlike in humans, the cardiac cycle in small animals is too short (~ 150 ms under anesthesia) to acquire a full set of k-space data in a single inversion. Therefore, by allowing for short repetition times, the *SI-method* considerably shortens the total acquisition time for high resolution experiments. Depending on the heart rate, the total scan time for the *T1-method* ranged from 6–8 minutes while that for the *SI-method* was 2–4 minutes. In human studies, the single shot acquisition limits the in-plane resolution since all k-space data must be acquired in a short period of time in diastole. In this context the segmented acquisition scheme of the *SI-method* and its theoretical model allows to achieve high resolution perfusion maps if the temporal noise could be mitigated.

The phantom and the *in vivo* study results show that the *P* values measured with the *SI-method* is comparable with that measured with the conventional *T1-method*. With the flow phantom we were able to verify the accuracy of the *SI-method* albeit with no spin exchange effects. In order to apply the *SI-method* to the flow phantom, the theoretical model was modified to reflect spin non-exchangeability. In this respect, the phantom is an exact physical replica of the modified theoretical model. Therefore the accuracy of the predicted outcome carried no uncertainty due to the design of the model. Hence the phantom experiment was an ideal method to test the relationship of *P* to flow sensitized MR signal as predicted by physics. However, in reality it is difficult to ascertain how well the two compartment model mimic myocardial tissue.

The fluorescent microspheres or fluospheres, have recently become a popular alternative to radioactive microspheres for quantitating myocardial perfusion in rodent experimental studies. The myocardial perfusion of healthy C57BL/6 mice measured using fluorescent microspheres (5.7 ± 0.3 ml/g/min**)** by Richer et al. [16] closely matches that of the *SI-method* (5.6 ± 0.4 ml/g/min) and the *T1-method* (5.6 ± 0.3 ml/g/min) in our study. As seen in Table Table3,3, our values are also within the range of published values [6,14,17-20], for myocardial perfusion in mice using standard ASL methods. The mean standard deviation of repeated measures by the *SI-method* is 3.6%. This shows that the perfusion quantified using the proposed method is highly reproducible. Interestingly we noted that the variation of *P* within the myocardium was significantly higher for the *T1-method* compared to that of the *SI-method* in the *in vivo* study. A probable cause is the variation in beat to beat time difference in the cardiac cycle resulting errors in the inversion times. This problem is exuberated by the long scan time of the *T1-method*. However, recently introduced post processing techniques could be used to mitigate this error [17].

In general, regions of LGE matched with regions of hypo-perfusion in the ischemia reperfusion study. But in some cases, low perfusion regions extended to the mid-myocardium beyond the boundary of the enhancement (arrows, Figure Figure3).3). It is possible that contrast material was not retained in these areas since it was surrounded by normally perfused tissue. It should also be noted that the perfusion maps have a lower resolution than that of the LGE images. Therefore these low perfusion extensions represents some volume averaging.

Although the slice thickness of the inversion pulse was twice that of the excitation pulse we believe it was much less in the actual implementation since the RF pulse shapes were significantly different. The goal was to compensate for miss alignment of the inversion and excitation slice profiles due to motion. However the thicker inversion slice underestimates perfusion. The shorter scan time of the *SI-method* somewhat reduces the errors due to motion thereby allowing for better alignment of the two profiles. This error can be vastly mitigated in human studies where single shot acquisitions can be used. Recently published work by Zun et al. demonstrates the feasibility of such a signal based ASL method with a single shot acquisition in humans [10].

The aim here was to introduce a new theoretical frame work for signal based ASL but we admit that much work remains to be done to improve the imaging technique. The quality of the inversion is problematic for high fields in general where B1 field inhomogeneity is much worse; thus remains a source of error in this application. Transit time effects and coil inflow time effects are associated with improper inversion pulses. These issues were not addressed in this study. We refer the readers to a rigorous analysis of these effects in the context of short TR by Pell et al. [11]. Compared to brain studies physiological noise is another major factor affecting the quality of cardiac perfusion maps. Some tasks for future consideration are: reducing physiological noise, overcoming mis-registration of differently prepared images, achieving a better slice profile for inversion, correcting for varied TI times of different k-space lines. With these improvements in place, the SI-method could potentially be used to generate high resolution pixel-by-pixel perfusion maps in human hearts.

This study demonstrates that signal intensity based ASL method is a robust alternative to the conventional *T1-method*. The proposed *SI-method* with a segmented acquisition scheme allows faster high resolution perfusion mapping in small animals. Compared to other signal intensity based ASL methods, the two compartmented model used in the present study makes it biologically more accurate.

During repeated inversion pulses at intervals of TR the longitudinal magnetization relaxes as,

$${M}_{z}\left(TI\right)={M}_{z}\left(0\right){e}^{-\frac{TI}{{T}_{1}}}+{M}_{0}\left(1-{e}^{-\frac{TI}{{T}_{1}}}\right)$$

(14)

The 90° excitation pulse saturates the M_{z} and subsequently recovers during the time TR-TI to

$${M}_{n}^{-}={M}_{0}\left(1-{e}^{-\frac{TR-TI}{{T}_{1}}}\right)$$

(15)

just before the nth inversion pulse. Therefore just after the nth inversion pulse the magnetization is given by;

$${M}_{n}^{+}=-{M}_{n}^{-}$$

(16)

Combining (1) (2) and (3) we get,

$${M}_{z}\left(TI\right)=-{M}_{0}\left(1-{e}^{-\frac{TR-TI}{{T}_{1}}}\right){e}^{-\frac{TI}{{T}_{1}}}+{M}_{0}\left(1-{e}^{-\frac{TI}{{T}_{1}}}\right)$$

(17)

and

$${M}_{Z}\left(TI\right)={M}_{0}\left(1-2{e}^{-\frac{TI}{T1}}+{e}^{-\frac{TR}{{T}_{1}}}\right)$$

(18)

The authors declare that they have no competing interests.

SA: Contributed to the design, mathematical derivations, acquisition and analysis of data, and drafting the manuscript; MS: Performed ischemia reperfusion surgery. JPW: Contributed to the conception and design, supervised acquisition, analysis and interpretation of data, involved in drafting the manuscript. All authors read and approved the final manuscript.

This work was supported by the NIH Heart, Lung and Blood Institute grant K25HL102244.

- Bauer WR, Hiller KH, Roder F, Rommel E, Ertl G, Haase A. Magnetization exchange in capillaries by microcirculation affects diffusion-controlled spin-relaxation: a model which describes the effect of perfusion on relaxation enhancement by intravascular contrast agents. Magn Reson Med. 1996;35(1):43–55. doi: 10.1002/mrm.1910350107. [PubMed] [Cross Ref]
- Northrup BE, McCommis KS, Zhang H, 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(1):53. doi: 10.1186/1532-429X-10-53. [PMC free article] [PubMed] [Cross Ref]
- Wacker CM, Fidler F, Dueren C, Hirn S, Jakob PM, Ertl G, Haase A, Bauer WR. Quantitative assessment of myocardial perfusion with a spin-labeling technique: preliminary results in patients with coronary artery disease. J Magn Reson Imaging. 2003;18(5):555–560. doi: 10.1002/jmri.10386. [PubMed] [Cross Ref]
- 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(6):1240–1245. doi: 10.1002/jmri.1880080610. [PubMed] [Cross Ref]
- 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(1):62–67. doi: 10.1002/mrm.10676. [PubMed] [Cross Ref]
- 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(3):584–592. doi: 10.1002/mrm.20327. [PubMed] [Cross Ref]
- Zhang H, Shea SM, Park V, Li D, Woodard PK, Gropler RJ, Zheng J. Accurate myocardial T1 measurements: toward quantification of myocardial blood flow with arterial spin labeling. Magn Reson Med. 2005;53(5):1135–1142. doi: 10.1002/mrm.20461. [PubMed] [Cross Ref]
- Buxton RB, Frank LR, Wong EC, Siewert B, Warach S, Edelman RR. A general kinetic model for quantitative perfusion imaging with arterial spin labeling. Magn Reson Med. 1998;40(3):383–396. doi: 10.1002/mrm.1910400308. [PubMed] [Cross Ref]
- Kwong KK, Chesler DA, Weisskoff RM, Donahue KM, Davis TL, Ostergaard L, Campbell TA, Rosen BR. MR perfusion studies with T1-weighted echo planar imaging. Magn Reson Med. 1995;34(6):878–887. doi: 10.1002/mrm.1910340613. [PubMed] [Cross Ref]
- Zun Z, Wong EC, Nayak KS. Assessment of myocardial blood flow (MBF) in humans using arterial spin labeling (ASL): feasibility and noise analysis. Magn Reson Med. 2009;62(4):975–983. doi: 10.1002/mrm.22088. [PubMed] [Cross Ref]
- Pell GS, Thomas DL, Lythgoe MF, Calamante F, Howseman AM, Gadian DG, Ordidge RJ. Implementation of quantitative FAIR perfusion imaging with a short repetition time in time-course studies. Magn Reson Med. 1999;41(4):829–840. doi: 10.1002/(SICI)1522-2594(199904)41:4<829::AID-MRM24>3.0.CO;2-U. [PubMed] [Cross Ref]
- Bauer WR, Roder F, Hiller KH, Han H, Frohlich S, Rommel E, Haase A, Ertl G. The effect of perfusion on T1 after slice-selective spin inversion in the isolated cardioplegic rat heart: measurement of a lower bound of intracapillary-extravascular water proton exchange rate. Magn Reson Med. 1997;38(6):917–923. doi: 10.1002/mrm.1910380610. [PubMed] [Cross Ref]
- Bauer WR, Schulten K. Theory of contrast agents in magnetic resonance imaging: coupling of spin relaxation and transport. Magn Reson Med. 1992;26(1):16–39. doi: 10.1002/mrm.1910260104. [PubMed] [Cross Ref]
- 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(3):601–606. doi: 10.1002/mrm.20373. [PubMed] [Cross Ref]
- Detre JA, Leigh JS, Williams DS, Koretsky AP. Perfusion imaging. Magn Reson Med. 1992;23(1):37–45. doi: 10.1002/mrm.1910230106. [PubMed] [Cross Ref]
- Richer C, Domergue V, Gervais M, Bruneval P, Giudicelli JF. Fluospheres for cardiovascular phenotyping genetically modified mice. J Cardiovasc Pharmacol. 2000;36(3):396–404. doi: 10.1097/00005344-200009000-00017. [PubMed] [Cross Ref]
- 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(3):648–657. doi: 10.1002/mrm.22280. [PMC free article] [PubMed] [Cross Ref]
- 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(6):1592–1598. doi: 10.1002/mrm.22470. [PMC free article] [PubMed] [Cross Ref]
- Nahrendorf M, Streif JU, Hiller KH, Hu K, Nordbeck P, Ritter O, Sosnovik D, Bauer L, Neubauer S, Jakob PM. et al. 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(6):H2516–H2521. doi: 10.1152/ajpheart.01038.2005. [PubMed] [Cross Ref]
- Vandsburger MH, French BA, Helm PA, Roy RJ, Kramer CM, Young AA, Epstein FH. Multi-parameter in vivo cardiac magnetic resonance imaging demonstrates normal perfusion reserve despite severely attenuated beta-adrenergic functional response in neuronal nitric oxide synthase knockout mice. Eur Heart J. 2007;28(22):2792–2798. doi: 10.1093/eurheartj/ehm241. [PubMed] [Cross Ref]

Articles from Journal of Cardiovascular Magnetic Resonance are provided here courtesy of **BioMed Central**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |