Quantitative myocardial perfusion in mice based on the signal intensity of flow sensitized CMR
© Abeykoon et al.; licensee BioMed Central Ltd. 2012
Received: 16 May 2012
Accepted: 11 October 2012
Published: 24 October 2012
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 T1 measurements.
Where, Tls and Tlg are the apparent longitudinal relaxation times for slice select and non-select inversion respectively. Tlc 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 Tl 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 Tl 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 . 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. . 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.  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.
SI-method for a fast exchanging two compartment model of tissue
where λ is the blood tissue water partition coefficient , 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, mp(t), depending on the type of spin preparation. For slice-select inversion, mp(t) is the equilibrium magnetization of capillary blood hence is equal to the constant mc(0) in magnitude. For non-select inversion, at steady state, , where TI and TR are the inversion and repetition times respectively (see Appendix).
SI-method for a non-exchangeable two compartment model of tissue
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 Kc=0, Ke=0 in equations (2) and (3). The rest of the derivation is similar to as in the previous section.
Here we do not normalize ΔM by the signal intensity of a steady state acquisition, instead measure mc(0), the equilibrium magnetization of the capillary fluid, using a spin echo sequence with long TR.
T1-method for a non-exchangeable two compartment model of tissue
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.  for a fast exchanging two compartment model of tissue. Here we modify it to reflect no exchange between compartments. Again we begin with setting Kc=0, Ke=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. mc(0)=me(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, mp(t), depending on the type of spin preparation. According to the normalizing scheme introduced here the input magnetization, mp(t), for slice-select inversion is zero and , for non-select inversion.
The RBV in the context of the said phantom is equal to the volume ratio Vc/VT, 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.
Two compartmented flow phantom study
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 cm2, 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 . 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 .
In vivo mouse study
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 1. A gated, Look-Locker sequence  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 cm2, 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: T1 of blood, T1c = 1800 ms, T1 of myocardial tissue in the absence of flow effects, T1 = 1400 ms, λ=0.95 ml/g . 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 .
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 error bars shown in Figure 3 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.
Comparison of in vivo results for each mouse
Reproducibility of the SI-method
Mean ± std. dev
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  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  and its derivative ASL methods by Kwong , Buxton  and Zun  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 1 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.
Published values of myocardial perfusion of healthy mice using CMR methods
 Makowski M. et al., 2010
7.3 ± 0.8 *
 Streif JU. et al., 2005
7.0 ± 0.5 *
 Nahrendorf M. et al., 2006
6.7 ± 0.3 *
 Kober F. et al., 2005
6.0 ± 1.9 +
5.6 ± 0.5 +
 Vandsburger MH. et al., 2010
5.2 ± 0.8 *
 Vandsburger MH. et al., 2007
4.3 ± 0.3 *
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 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 .
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. . 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.
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. 10.1002/mrm.1910350107.View ArticlePubMedGoogle Scholar
- 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-10.1186/1532-429X-10-53.PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1002/jmri.10386.View ArticlePubMedGoogle Scholar
- 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. 10.1002/jmri.1880080610.View ArticlePubMedGoogle Scholar
- 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. 10.1002/mrm.10676.View ArticlePubMedGoogle Scholar
- 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. 10.1002/mrm.20327.View ArticlePubMedGoogle Scholar
- 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. 10.1002/mrm.20461.View ArticlePubMedGoogle Scholar
- 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. 10.1002/mrm.1910400308.View ArticlePubMedGoogle Scholar
- 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. 10.1002/mrm.1910340613.View ArticlePubMedGoogle Scholar
- 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. 10.1002/mrm.22088.View ArticlePubMedGoogle Scholar
- 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. 10.1002/(SICI)1522-2594(199904)41:4<829::AID-MRM24>3.0.CO;2-U.View ArticlePubMedGoogle Scholar
- 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. 10.1002/mrm.1910380610.View ArticlePubMedGoogle Scholar
- 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. 10.1002/mrm.1910260104.View ArticlePubMedGoogle Scholar
- 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. 10.1002/mrm.20373.View ArticlePubMedGoogle Scholar
- Detre JA, Leigh JS, Williams DS, Koretsky AP: Perfusion imaging. Magn Reson Med. 1992, 23 (1): 37-45. 10.1002/mrm.1910230106.View ArticlePubMedGoogle Scholar
- 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. 10.1097/00005344-200009000-00017.View ArticlePubMedGoogle Scholar
- 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. 10.1002/mrm.22280.PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1002/mrm.22470.PubMed CentralView ArticlePubMedGoogle Scholar
- Nahrendorf M, Streif JU, Hiller KH, Hu K, Nordbeck P, Ritter O, Sosnovik D, Bauer L, Neubauer S, Jakob PM: 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. 10.1152/ajpheart.01038.2005.View ArticlePubMedGoogle Scholar
- 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. 10.1093/eurheartj/ehm241.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.