Email updates

Keep up to date with the latest news and content from Journal of Cardiovascular MR and BioMed Central.

Open Access Research

Quantification of myocardial perfusion using CMR with a radial data acquisition: comparison with a dual-bolus method

Tae Ho Kim12, Nathan A Pack1, Liyong Chen1 and Edward VR DiBella1*

Author Affiliations

1 Utah Center for Advanced Imaging Research, Department of Radiology, 729 Arapeen Dr. University of Utah, Salt Lake City, Utah 84108, USA

2 Current Address: Department of Radiation Oncology, Stanford University, Stanford, CA 94305, USA

For all author emails, please log on.

Journal of Cardiovascular Magnetic Resonance 2010, 12:45  doi:10.1186/1532-429X-12-45


The electronic version of this article is the complete one and can be found online at: http://www.jcmr-online.com/content/12/1/45


Received:11 February 2010
Accepted:23 July 2010
Published:23 July 2010

© 2010 Kim 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.

Abstract

Background

Quantitative estimates of myocardial perfusion generally require accurate measurement of the arterial input function (AIF). The saturation of signal intensity in the blood that occurs with most doses of contrast agent makes obtaining an accurate AIF challenging. This work seeks to evaluate the performance of a method that uses a radial k-space perfusion sequence and multiple saturation recovery times (SRT) to quantify myocardial perfusion with cardiovascular magnetic resonance (CMR).

Methods

Perfusion CMR was performed at 3 Tesla with a saturation recovery radial turboFLASH sequence with 72 rays. Fourteen subjects were given a low dose (0.004 mmol/kg) of dilute (1/5 concentration) contrast agent (Gd-BOPTA) and then a higher non-dilute dose of the same volume (0.02 mmol/kg). AIFs were calculated from the blood signal in three sub-images with differing effective saturation recovery times. The full and sub-images were reconstructed iteratively with a total variation constraint. The images from the full 72 ray data were processed to obtain six tissue enhancement curves in two slices of the left ventricle in each subject. A 2-compartment model was used to determine absolute flows

Results

The proposed multi-SRT method resulted in AIFs that were similar to those obtained with the dual-bolus method. Myocardial blood flow (MBF) estimates from the dual-bolus and the multi-SRT methods were related by MBFmulti-SRT = 0.85MBFdual-bolus + 0.18 (r = 0.91).

Conclusions

The multi-SRT method, which uses a radial k-space perfusion sequence, can be used to obtain an accurate AIF and thus quantify myocardial perfusion for doses of contrast agent that result in a relatively saturated AIF.

Background

Myocardial perfusion cardiovascular magnetic resonance (CMR) is a useful modality for measuring perfusion to detect coronary artery disease and myocardial ischemia [1,2]. Quantitative estimates of perfusion require that the gadolinium concentration in the arterial input function (AIF) and the myocardial tissue curves be known. This is achieved by estimating the T1 values at each time frame of the blood and tissue. If the pre-contrast T1 value of blood and tissue is also known, these T1 values can be used to calculate the gadolinium concentrations over time. Equivalently, if the change in signal intensity does not saturate (meaning a loss in linearity between the gadolinium concentration and the change in the signal from pre-contrast), the change in signal intensity time curves can be used without explicit conversion to gadolinium concentration. However, the saturation of signal intensity in the blood during the first pass that occurs with most doses of gadolinium-based contrast agents hinders quantification. Dual-bolus methods with a 1 ml first injection [3] and with dilute matched volume injections [4-6] have been used to enable quantitative perfusion CMR. With these methods, an unsaturated AIF is obtained from a low dose injection, which is then scaled and used to replace the saturated AIF from the high dose injection. The requirement of two injections can be logistically challenging and may reflect different physiological states, particularly during vasodilation.

Recently, an alternative method for estimating T1 and thus gadolinium concentration from a single injection using a radial k-space acquisition was proposed [7]. That method created four sub-images each with different effective saturation recovery times (eSRTs) from 96 ray acquisitions to estimate the AIF accurately. Extracting images with different eSRTs from a single dataset is possible with a radial acquisition if one uses a subset of the rays to perform the image reconstruction, and it is assumed that each ray of the subset contributes equally to the image [7].

However, no measures of truth were used to evaluate the method in vivo. In this work, we employ a similar multiple effective saturation recovery time (multi-SRT) approach for dynamic estimation of contrast agent concentration in the left ventricle blood pool. The purpose of this study was to study the performance of the radial k-space perfusion sequence multi-SRT method by comparison to the established dual-bolus approach.

Methods

MR of Vials

To validate the multi-SRT method in a controlled setting, we estimated the T1 of vials with known concentrations of contrast agent using the saturation recovery radial acquisition described in the human protocol below. Eighteen vials with different concentrations of Gd-BOPTA (Multihance, Bracco Diagnostics Inc.) were prepared with the range of T1 spanning 10-1000 ms, covering the range of T1 of blood after a typical contrast injection. Each vial was positioned in a water bath containing gadolinium (~0.3 mM) to reduce magnetic susceptibility artifacts in the images. The image signal in a region of interest in each vial was averaged over 10 time frames to improve SNR. Reconstruction and multi-SRT processing were performed as detailed for the human studies below, with the T1 of the different vials being estimated jointly.

For comparison, T1 was also measured using an inversion recovery pulse sequence with three inversion times. The imaging parameters for this fast spin echo sequence were TR/TE = 4000/15 ms, TI = 25/100/300 ms, 9 echoes, and resolution 2.0 × 2.0 × 8.0 mm3.

Human Subject CMR Protocol

Fourteen volunteers (8 female, 6 male, 52 ± 13 yrs, 183 ± 55 lbs) assumed to be without ischemia were imaged in this study. Informed consent from all the subjects was obtained in accordance with the University of Utah's human subject policies. All studies employed an extra "backcheck" or one-way valve on the saline side of the power injector, in addition to the standard backcheck valve on the contrast agent side of the power injector. This was to prevent any flow of contrast agent into the saline syringe, which has been observed to occur (though may be a negligible amount, depending on dose). All subjects were given a low dose (0.004 mmol/kg) of dilute (1/5 concentration) contrast agent (Gd-BOPTA, Multihance) and perfusion data acquired. Approximately 3-5 minutes later, a higher non-dilute dose (0.02 mmol/kg) was given and perfusion data again acquired. Both injections were given at 5 cc/sec and flushed with 25 ml of saline at 5 cc/sec. The studies were performed on a Siemens 3 Tesla Trio (n = 12) or a Siemens 3 Tesla Verio (n = 2) system (Siemens Medical Solutions, Erlangen, Germany), using the spine and body phased array coils. We employed a saturation recovery radial turboFLASH sequence that was very similar to the standard saturation recovery Cartesian turboFLASH sequences, just with a radial non-Cartesian readout. The radial readout acquired 72 rays over 180° in an interleaved manner (12 subsets of 6 rays each) approximately 10 msec after each saturation pulse to give the data for one slice. In order to sample more of k-space over time, a small offset of 180/(4*72) = 0.625° was used for sequential time frame acquisitions, with a period of four [7]. The acquisition parameters were TR/TE = 2.6/1.14 ms, prescribed flip angle = 14°, slice-thickness = 8 mm, number of readouts = 256, FOV = ~300 mm and ~70 time frames. Three to four slices were acquired per beat, and subjects were requested to breathe shallowly.

Reconstruction of CMR Data

An iterative total variation constrained reconstruction was used with the 72 ray dynamic datasets to obtain the tissue enhancement curves. It was assumed for the dose used here that the tissue curves were not saturated. The reconstruction method was similar to that used in [8], although in this work a forward and inverse non-uniform FFT from [9] was used at every iteration. The same method was used to reconstruct three subsets of 24 rays each, in order to create three sub-images with different effective saturation recovery times. This is not possible with conventional Cartesian sequences - that is, using a subset of Cartesian lines (phase encodes) will result in aliased images that will mostly reflect edges if the center of k-space is not included. On the other hand, for radial acquisitions, any subset of rays can be used to reconstruct an image, albeit with more streaking than if all of the rays were used. The chosen rays will dictate the contrast of the sub-image because each ray samples the center of k-space which controls image contrast. Here, the three sub-images that were reconstructed were used solely to estimate the AIF.

Analysis of CMR Perfusion Data

First, the 72 ray (full data) reconstruction was processed to obtain myocardial tissue curves and saturated AIFs. This was done with freely available custom software [10]. Manual correction of respiratory motion over the time frames was done with the software, as was selection of endocardial and epicardial borders in each short axis slice. A region of interest was also drawn on the left ventricular blood pool to provide the AIF signal. The mean tissue enhancement signals from six equiangular regions of the myocardium in each of the short axis slices were obtained from the full 72 ray images. Two relatively basal slices were selected for further processing. A uniform region of signal that comprised ~20% of the left ventricular blood pool and that did not contain papillary muscles was manually selected for the AIF.

Multi-SRT Processing to obtain the AIF

The (saturated) AIFs from the three sub-images were created from the ROI chosen on the full data reconstruction. These three curves were then processed to obtain a single non-saturated AIF. Each of the sub-images has a different eSRT (43 ms for the first sub-image, 107 ms for the second sub-image, and 172 ms for the third sub-image), and can be expressed by the equation in [7] in terms of the readout flip angle α, a proton density term M, and T1. By simultaneously fitting the signal intensity from each sub-image curve to this signal equation, T1 was estimated for each time point. Only a single α and M were estimated for each time curve [7]. The T1 values at each time point were converted to concentration using the standard relation 1/T1-1/T1(0) = ß[Gd], with relaxivity ß = 5.5 mM-1s-1 and T1(0) defined as the pre-contrast T1 value. This equation states that changes in (inverse) T1 from pre-contrast are linearly related to changes in gadolinium contrast agent concentration ([Gd]). T1(0) was obtained by fitting the low (1/5 concentration) dose blood pool multi-SRT curves. The average of the first 7 pre-contrast time points was used to find T1(0).

The multi-SRT AIF concentration curve was then scaled so that its average value in the final four time points matched the average of the four final frames of the blood signal from the full 72 ray image [7]. This procedure assumed that in the final frames there was no signal saturation of the blood signal in the full image. Note that this scaling step makes the choice of the relaxivity ß irrelevant in the conversion to concentration equation given above, and was done to match to the scale of the tissue curves.

Dual-bolus Processing to obtain the AIF

For the dual-bolus method, the AIF for each subject was obtained from a manually drawn region of the left ventricular blood pool in the low dose 72 ray images. The AIF was upscaled by five to replace the volume matched high dose 72 ray AIF blood signal.

Fitting for MBF

For perfusion quantification, the AIFs obtained from the multi-SRT method and the tissue enhancement curves from the 72 ray higher dose acquisition were input to a 2-compartment model [11] using the software mentioned previously [10] in order to obtain Ktrans, a measure of absolute MBF. The fitting procedure was repeated using the dual-bolus AIF.

Statistical Analysis

Paired comparisons of MBF estimated using the two different AIFs were tested using the Students t-test. Pearson's correlation was calculated along with linear regression analysis, and Bland-Altman analysis was performed.

Results

Vials

T1 estimates from the 18 vials are presented in Figure 1. The multi-SRT method gave T1 values similar to the reference standard inversion recovery sequence, T1est = 1.02 T1ref+3.26 msec.

thumbnailFigure 1. The correlation between estimated T1 values in 18 vials using the multi-SRT method and the reference standard method (inversion recovery with three inversion recovery times).

Human Studies

Three sub-images and the corresponding full 72 ray image from a single time frame of a typical perfusion dataset are presented in Figure 2. Although the 72 ray combined image shows better discrimination of the myocardium and the left ventricular blood pool, each of the sub-images is of sufficient quality to obtain the multi-SRT blood enhancement signals. An example of the time curves that are obtained from each of the sub-images is shown in Figure 3. The curve from the shortest eSRT has the lowest signal and the signal increases with eSRT. Fits to the curves from the different eSRT images are also shown.

thumbnailFigure 2. Three sub-images with different eSRTs and the corresponding 72 ray combined image from a single time frame of a perfusion dataset are shown.

thumbnailFigure 3. The measured blood enhancement (saturated AIF) curves corresponding to the three sub-images with different eSRTs are shown, along with their fits obtained using the proposed multi-SRT measurement method.

In Figure 4, the corrected AIF using the proposed multi-SRT estimate method is shown, along with the saturated AIF from the "high" dose (0.02 mmol/kg) injection, and the upscaled low dose AIF from the dual-bolus method. The peak of the measured AIF from the 0.02 mmol/kg scan is approximately 40% below the upscaled low dose AIF. The proposed multi-SRT estimate method gave AIFs that were similar to those obtained with the dual-bolus method. The same type of comparison is shown for all 14 subjects in Figure 5. The other parameters estimated along with T1 and concentration were T1(0) = 2049 ± 434 ms, flip angle α = 8.5° ± 2.2°, and M= 812 ± 434.

thumbnailFigure 4. A corrected AIF using the proposed multi-SRT method (black), along with the saturated AIF from the higher dose (0.02 mmol/kg) injection (blue) and the upscaled low dose AIF is shown (red).

thumbnailFigure 5. Comparison of the corrected AIFs obtained using the multi-SRT method (red) and the low dose (0.004 mmol/kg) upscaled AIF used for the dual-bolus method (green), along with the saturated AIF from the high dose (0.02 mmol/kg) injection (blue), in all 14 subjects in the study. Peak height of the saturated AIF varies from approximately 30% reduction of the dual-bolus or multi-SRT AIF, to 70% reductions.

Figure 6 shows the correlation between MBF values obtained using the multi-SRT method and the dual-bolus method. The MBF results from the proposed method were similar to the dual-bolus method: MBFmulti-SRT = 0.85MBFdual-bolus + 0.18 (r = 0.91). The aggregate MBF values estimated from the dual-bolus and the multi-SRT methods for all 14 subjects were 0.66 ± 0.24 and 0.75 ± 0.23 ml/min/g, respectively. These means were significantly different, p < 0.05. Figure 7 shows a Bland-Altman plot of the MBF estimates from the multi-SRT method and the dual-bolus method. For these 14 subjects, there was a mean overestimation in the multi-SRT MBF estimates of 0.08 ± 0.1 ml/min/g compared to the dual bolus results.

thumbnailFigure 6. The correlation between MBF estimates using the proposed multi-SRT measurement method and dual-bolus imaging.

thumbnailFigure 7. A Bland-Altman plot showing the mean MBF estimates and the difference between MBF estimates from the multi-SRT method and dual-bolus imaging.

Discussion

We report on a multi-SRT quantitative perfusion CMR method using a radial k-space acquisition with 72 rays. The method uses an iterative reconstruction with a total variation constraint, which has been shown to produce image quality similar or better than Cartesian acquisitions when only 24 rays per slice and multiple slices are acquired after a single saturation pulse and the first slice is not considered [8]. In this work, the reconstructed images used all 72 rays to obtain tissue enhancement curves, and used three subset images of 24 rays each to estimate the non-saturated AIF. There were relatively small signal changes from the second sub-image (eSRT = 107 ms) to the third sub-image (eSRT = 172 ms) suggesting that the third sub-image may not be essential for the approach. However, all three sub-images were used to estimate T1 in order to minimize the effects of undersampling and noise in the sub-images.

Other methods have been proposed to handle the saturated AIF problem. Low doses [12-17] can be used but give low signal. Analytical methods convert the measured signal intensities into gadolinium concentrations [18,19], but can be sensitive to noise when the signal-concentration relationship is not close to linear. Alternatively, two groups have developed custom sequences to measure an additional slice at low resolution and short SRT so that signal in that slice will be more linear with concentration and that slice can be used solely for the estimation of the AIF [20,21].

The proposed multi-SRT method was shown here to give accurate T1 estimates in vials, when compared to T1 values using an inversion recovery pulse sequence, (Figure 1). Ref. [7] found a similar relationship, although a 96 ray acquisition and four subsets were used, and a 3 D gradient echo acquisition was used as the reference standard.

When estimating gadolinium concentrations of the AIF in this work, we did not assume the prescribed flip angle was accurate, and instead estimated it as an unknown parameter in the multi-SRT curve-fitting process. Both spatial excitation (flip angle) variations and an imperfect slice profile likely played a role in the findings of a net flip angle for the blood region of 8.5 ± 2.2 degrees, instead of the 14° that was designated in the radial perfusion sequence. This nearly 40% reduction in flip angle is similar to that found previously in the heart at 3T [22]. A 23-48% variation in flip angle across the heart has also been reported at 3T [23]. Flip angle variation over the region of interest may have led to some variation in the multi-SRT gadolinium concentration estimates of the AIF. Changes in the AIF would alter the MBF estimates of all tissue curves of that subject.

T1 values prior to contrast injection were also estimated in this work. For reference, Lu et al. reported the T1 of ex-vivo bovine blood at 3T as 1664 ± 14 msec [24]. Studies in ex-vivo human blood at 37°C and 3T found 1932 ± 85 msec [25]. In our studies, the aggregate pre-contrast T1 estimates in the left ventricular blood pool were 2049 ± 434 ms, using the radial perfusion sequence. The mean of this result is in reasonable agreement; the relatively high standard deviation of the T1 estimates may be due to inter-subject variability.

As shown in Figs. 4-5, when no signal intensity correction is performed, the peak of the measured AIF from the 0.02 mmol/kg injection scan is saturated by approximately 30-70%. This is likely a larger effect than found with 0.02 mmol/kg scans using Gd-DTPA or Gd-DTPA-BMA (Magnevist or Omniscan), which have lower relaxivity [26]. Although done at 1.5T, others have reported that Gd-BOPTA saturates at significantly lower concentrations than Gd-DTPA [27].

The proposed multi-SRT method can correct for signal saturation in the AIF using a single injection of contrast agent that overcomes the complexity of multi-injection protocols. In addition, by using a single injection, there is no physiological change between the corrected AIF and the measured AIF that could be confounding with some dual-bolus imaging protocols. For example, in this study there were sometimes variations in a subject's heart rate (and likely the hemodynamics of blood and contrast agent flow in the left ventricle) during the two separate injections, indicating some degree of physiological change. With a single injection of contrast agent using the proposed method, these changes are not an issue. This is a larger concern in stress studies where it can be more difficult to have vasodilation kept constant, particularly if breath-holds also change perfusion to some degree.

The new multi-SRT method resulted in AIFs that were similar to those obtained with the dual-bolus method, and gave MBF estimates that on the average differed by less than 11%. This difference is small and likely not critical for clinical applications of the method.

The MBF values are slightly lower than that reported in some CMR studies and with other modalities such as PET. Model-independent and Fermi model methods also result in very similar perfusion estimates [28]. The relatively low perfusion values may be due in part from coil sensitivity- preliminary results with using proton density measurements to compensate for coil sensitivity variations tended to scale the blood pool region differently than the tissue, typically increasing MBF ~15%. As well, flow effects from arterial blood that has experienced fewer alpha pulses than the stationary tissue would tend to decrease MBF estimates due to increasing the AIF. This could be modeled with the equation used here if it could be determined how many fewer alpha pulses affected the blood pool. Adjustment for rate-pressure product [29] left mean MBF values virtually unchanged, with higher standard deviations: 0.67 ± 0.32 for the dual-bolus method, and 0.75 ± 0.33 for the multi-SRT method. The rate-pressure product adjusted MBFs did have a closer relationship: MBFmulti-SRT = 0.97MBFdualbolus + 0.1, r = 0.95.

Another source of error with the multi-SRT method are T2* effects during the first pass of contrast agent through the left ventricular blood pool [30]. In our study, T2* effects were assumed to be negligible for the moderate doses used. More study is needed to determine how well the method works with larger doses, and if T2* effects at the peak of the bolus can be accurately estimated [7].

Conclusions

The multi-SRT estimation method using an undersampled radial k-space perfusion sequence accurately quantifies myocardial perfusion in the presence of 30-70% peak signal decrease due to saturation of the AIF. Unlike the dual-bolus imaging method, the multi-SRT method requires only a single contrast agent injection, which can greatly simplify myocardial perfusion studies, especially during stress imaging.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

TK performed image acquisition, image reconstruction, developed portions of the analysis, and created the initial draft of the manuscript. NAP participated in image acquisition, reconstruction, and analysis. LC performed key pulse sequence modifications and reconstruction design. EVRD conceived of the study, was involved with all aspects and edited the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We appreciate the assistance of Henry Buswell with data collection. Funding was received from Funding NIH R01 EB00177. Portions of this work have appeared at Society Cardiovascular Magnetic Resonance, 2010.

References

  1. Al-Saadi N, Nagel E, Gross M, Schnackenburg B, Paetsch I, Klein C, Fleck E: Improvement of myocardial perfusion reserve early after coronary intervention: assessment with cardiac magnetic resonance imaging.

    J Am Coll Cardiol 2000, 36:1557-1564. PubMed Abstract | Publisher Full Text OpenURL

  2. Sandstede JJW, Lipke C, Beer M, Harre K, Pabst T, Kenn W, Neubauer S, Hahn D: Analysis of first-pass and delayed contrast-enhancement patterns of dysfunctional myocardium on MR imaging: use in the prediction of myocardial viability.

    AJR 2000, 174:1737-1740. PubMed Abstract | Publisher Full Text OpenURL

  3. Kostler H, Ritter C, Lipp M, Beer M, Hahn D, Sandstede J: Prebolus quantitative MR heart perfusion imaging.

    Magn Reson Med 2004, 52:296-299. PubMed Abstract | Publisher Full Text OpenURL

  4. Christian TF, Rettmann DW, Aletras AH, Liao SL, Taylor JL, Balaban RS, Arai AE: Absolute myocardial perfusion in canines measured by using dual-bolus first-pass MR imaging.

    Radiology 2004, 232:677-684. PubMed Abstract | Publisher Full Text OpenURL

  5. Hsu LY, Rhoads KL, Holly JE, Kellman P, Aletras AH, Arai AE: Quantitative myocardial perfusion analysis with a dual-bolus contrast-enhanced first-pass MRI technique in humans.

    J Magn Reson Imaging 2006, 23:315-322. PubMed Abstract | Publisher Full Text OpenURL

  6. Christian TF, Aletras AH, Arai AE: Estimation of absolute myocardial blood flow during first-pass MR perfusion imaging using a dual-bolus injection technique: comparison to single-bolus injection method.

    J Magn Reson Imaging 2008, 27:1271-1277. PubMed Abstract | Publisher Full Text OpenURL

  7. Kholmovski EG, DiBella EVR: Perfusion MRI with Radial Acquisition for Arterial Input Function Assessment.

    Magn Reson Med 2007, 57:821-827. PubMed Abstract | Publisher Full Text OpenURL

  8. Adluru G, McGann C, Speier P, Kholmovski EG, Shaaban A, Dibella EV: Acquisition and reconstruction of undersampled radial data for myocardial perfusion magnetic resonance imaging.

    J Magn Reson Imaging 2009, 29:466-473. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Fessler JA: Matlab Tomography Toolbox.

    Book Matlab Tomography Toolbox (Editor ed.^eds.). City 2004. OpenURL

  10. Pack N, Vijayakumar S, Kim TH, McGann C, DiBella E: A Semi-Automatic Software Package for Analysis of Dynamic Contrast-Enhanced MRI Myocardial Perfusion Studies.

    Computers in Cardiology; Park City, Utah 2009. OpenURL

  11. Tofts PS, Brix G, Buckley DL, Evelhoch JL, Henderson E, Knopp MV, Larsson HBW, Lee T-Y, Mayr NA, Parker GJM, Port RE, Taylor J, Weisskoff RM: Estimating kinetic parameters from dynamic contrast-enhanced T1-weighted MRI of a diffusible tracer: standardized quantities and symbols.

    J Magn Reson Imaging 1999, 10:223-232. PubMed Abstract | Publisher Full Text OpenURL

  12. Wilke N, Jerosch-Herold M, Wang Y, Huang Y, Christensen BV, Stillman AE, Ugurbil K, McDonald K, Wilson RF: Myocardial Perfusion Reserve: Assessment with Multisection, Quantitative, First-Pass MR Imaging.

    Radiology 1997, 204:373-384. PubMed Abstract | Publisher Full Text OpenURL

  13. Jerosch-Herold M, Wilke N, Stillman AE: Magnetic resonance quantification of the myocardial perfusion reserve with a Fermi function model for constrained deconvolution.

    Med Phys 1998, 25:73-84. PubMed Abstract | Publisher Full Text OpenURL

  14. Muehling OM, Wilke NM, Panse P, Jerosch-Herold M, Wilson BV, Wilson RF, Miller LW: Reduced myocardial perfusion reserve and transmural perfusion gradient in heart transplant arteriopathy assessed by magnetic resonance imaging.

    J Am Coll Cardiol 2003, 42:1054-1060. PubMed Abstract | Publisher Full Text OpenURL

  15. Muehling OM, Jerosch-Herold M, Panse P, Zenovich A, Wilson BV, Wilson RF, Wilke N: Regional heterogeneity of myocardial perfusion in healthy human myocardium: assessment with magnetic resonance perfusion imaging.

    J Cardiovasc Magn Reson 2004, 6:499-507. PubMed Abstract | Publisher Full Text OpenURL

  16. Jerosch-Herold M, Swingen C, Seethamraju RT: Myocardial blood flow quantification with MRI by model-independent deconvolution.

    Med Phys 2002, 29:886-897. PubMed Abstract | Publisher Full Text OpenURL

  17. Selvanayagam JB, Jerosch-Herold M, Porto I, Sheridan D, Cheng AS, Petersen SE, Searle N, Channon KM, Banning AP, Neubauer S: Resting myocardial blood flow is impaired in hibernating myocardium: a magnetic resonance study of quantitative perfusion assessment.

    Circulation 2005, 112:3289-3296. PubMed Abstract | Publisher Full Text OpenURL

  18. Cernicanu A, Axel L: Theory-based signal calibration with single-point T1 measurements for first-pass quantitative perfusion MRI studies.

    Acad Radiol 2006, 13:686-693. PubMed Abstract | Publisher Full Text OpenURL

  19. Hsu LY, Kellman P, Arai AE: Nonlinear myocardial signal intensity correction improves quantification of contrast-enhanced first-pass MR perfusion in humans.

    J Magn Reson Imaging 2008, 27:793-801. PubMed Abstract | Publisher Full Text OpenURL

  20. Gatehouse PD, Elkington AG, Ablitt NA, Yang GZ, Pennell DJ, Firmin DN: Accurate assessment of the arterial input function during high-dose myocardial perfusion cardiovascular magnetic resonance.

    J Magn Reson Imaging 2004, 20:39-45. PubMed Abstract | Publisher Full Text OpenURL

  21. Kim D, Axel L: Multislice, dual-imaging sequence for increasing the dynamic range of the contrast-enhanced blood signal and CNR of myocardial enhancement at 3T.

    Journal of Magnetic Resonance Imaging 2006, 23:81-86.

    Erratum in: J Magn Reson Imaging. 2007, 25(4):878

    PubMed Abstract | Publisher Full Text OpenURL

  22. Sung K, Nayak KS: B1+ compensation in 3T cardiac imaging using short 2DRF pulses.

    Magn Reson Med 2008, 59:441-446. PubMed Abstract | Publisher Full Text OpenURL

  23. Sung K, Nayak KS: Measurement and characterization of RF nonuniformity over the heart at 3T using body coil transmission.

    J Magn Reson Imaging 2008, 27:643-648. PubMed Abstract | Publisher Full Text OpenURL

  24. Lu H, Clingman C, Golay X, van Zijl PC: Determining the longitudinal relaxation time (T1) of blood at 3.0 Tesla.

    Magn Reson Med 2004, 52:679-682. PubMed Abstract | Publisher Full Text OpenURL

  25. Stanisz GJ, Odrobina EE, Pun J, Escaravage M, Graham SJ, Bronskill MJ, Henkelman RM: T1, T2 relaxation and magnetization transfer in tissue at 3T.

    Magn Reson Med 2005, 54:507-512. PubMed Abstract | Publisher Full Text OpenURL

  26. Pintaske J, Martirosian P, Graf H, Erb G, Lodemann KP, Claussen CD, Schick F: Relaxivity of Gadopentetate Dimeglumine (Magnevist), Gadobutrol (Gadovist), and Gadobenate Dimeglumine (MultiHance) in human blood plasma at 0.2, 1.5, and 3 Tesla.

    Invest Radiol 2006, 41:213-221. PubMed Abstract | Publisher Full Text OpenURL

  27. Kostler H, Ritter C, Lipp M, Beer M, Hahn D, Sandstede J: Comparison of different contrast agents and doses for quantitative MR myocardial perfusion imaging.

    J Magn Reson Imaging 2008, 28:382-389. PubMed Abstract | Publisher Full Text OpenURL

  28. Pack N, DiBella EVR: Comparison of myocardial blood flow estimates from dynamic contrast-enhanced magnetic resonance imaging with four quantitative analysis methods.

    Magn Reson Med 2010, 64:125-137. PubMed Abstract | Publisher Full Text OpenURL

  29. Czernin J, Muller P, Chan S, Brunken RC, Porenta G, Krivokapich J, Chen K, Chan A, Phelps ME, Schelbert HR: Influence of age and hemodynamics on myocardial blood flow and flow reserve.

    Circulation 1993, 88:62-69. PubMed Abstract | Publisher Full Text OpenURL

  30. deBazelaire C, Rofsky NM, Duhamel G, Zhang J, Michaelson MD, George D, Alsop DC: Combined T2* and T1 measurements for improved perfusion and permeability studies in high field using dynamic contrast enhancement.

    Eur Radiol 2006, 16:2083-2091. PubMed Abstract | Publisher Full Text OpenURL