Modeling regional myocardial flows from residue functions of an intravascular indicator

Effects of bolus injection duration on perfusion estimates in dynamic CT and dynamic susceptibility contrast MRI

Estimates of cerebral blood flow (CBF) and tissue mean transit time (MTT) have been shown to differ between dynamic CT perfusion (CTP) and dynamic susceptibility contrast MRI (DSC-MRI). This study investigates whether these discrepancies regarding CBF and MTT between CTP and DSC-MRI can be attributed to the different injection durations of these techniques. Five subjects were scanned using CTP and DSC-MRI. Region-wise estimates of CBF, MTT, and cerebral blood volume (CBV) were derived based on oscillatory index regularized singular value decomposition. A parametric model that reproduced the shape of measured time curves and characteristics of resulting perfusion parameter estimates was developed and used to simulate data with injection durations typical for CTP and DSC-MRI for a clinically relevant set of perfusion scenarios and noise levels. In simulations, estimates of CBF/MTT showed larger negative/positive bias and increasing variability for CTP when compared to DSC-MRI, especially for high CBF levels. While noise also affected estimates, at clinically relevant levels, the injection duration effect was larger. There are several methodological differences between CTP and DSC-MRI. The results of this study suggest that the injection duration is among those that can explain differences in estimates of CBF and MTT between these bolus tracking techniques.

Quantification of balanced SSFP myocardial perfusion imaging at 1.5 T: Impact of the reference image

PURPOSE

To investigate the use of a high flip-angle (HFA) balanced SSFP (bSSFP) reference image (in comparison to conventional proton density [PD]-weighted reference images) for conversion of bSSFP myocardial perfusion images into dynamic T1 maps for improved myocardial blood flow (MBF) quantification at 1.5 T.

METHODS

The HFA-bSSFP (flip angle [FA] = 50°), PD gradient-echo (PD-GRE; FA = 5°), and PD-bSSFP (FA = 8°) reference images were acquired before a dual-sequence bSSFP perfusion acquisition. Simulations were used to study accuracy and precision of T1 and MBF quantification using the three techniques. The accuracy and precision of T1 , and the precision and intersegment variability of MBF were compared among the three techniques in 8 patients under rest conditions.

RESULTS

In simulations, HFA-bSSFP demonstrated improved T1 /MBF precision (higher T1 /MBF SD of 30%-80%/50%-100% and 30%-90%/60%-115% for PD-GRE and PD-bSSFP, respectively). Proton density-GRE and PD-bSSFP were more sensitive to effective FA than HFA-bSSFP (maximum T1 /MBF errors of 13%/43%, 20%/43%, and 1%/3%, respectively). Sensitivity of all techniques (defined as T1 /MBF errors) to native T1 , native T2 , and effective saturation efficiency were negligible (<1%/<1%), moderate (<14%/<19%), and high (<63%/<94%), respectively. In vivo, no difference in T1 accuracy was observed among HFA-bSSFP, PD-GRE, and PD-bSSFP (-9 ± 44 ms vs -28 ± 55 ms vs -22 ± 71 ms, respectively; p > .08). The HFA-bSSFP led to improved T1 /MBF precision (T1 /MBF SD: 41 ± 19 ms/0.24 ± 0.08 mL/g/min vs PD-GRE: 48 ± 20 ms/0.29 ± 0.09 mL/g/min and PD-bSSFP: 59 ± 23 ms/0.33 ± 0.11 mL/g/min; p ≤ .02) and lower MBF intersegment variability (0.14 ± 0.09 mL/g/min vs PD-GRE: 0.21 ± 0.09 mL/g/min and PD-bSSFP: 0.20 ± 0.10 mL/g/min; p ≤ .046).

CONCLUSION

We have demonstrated the feasibility of using a HFA-bSSFP reference image for MBF quantification of bSSFP perfusion imaging at 1.5 T. Results from simulations demonstrate that the HFA-bSSFP reference image results in improved precision and reduced sensitivity to effective FA compared with conventional techniques using a PD reference image. Preliminary in vivo data acquired at rest also demonstrate improved precision and intersegment variability using the HFA-bSSFP technique compared with PD techniques; however, a clinical study in patients with coronary artery disease under stress conditions is required to determine the clinical significance of this finding.

Clinical Translation of Three-Dimensional Scar, Diffusion Tensor Imaging, Four-Dimensional Flow, and Quantitative Perfusion in Cardiac MRI: A Comprehensive Review

Cardiovascular magnetic resonance (CMR) imaging is a versatile tool that has established itself as the reference method for functional assessment and tissue characterisation. CMR helps to diagnose, monitor disease course and sub-phenotype disease states. Several emerging CMR methods have the potential to offer a personalised medicine approach to treatment. CMR tissue characterisation is used to assess myocardial oedema, inflammation or thrombus in various disease conditions. CMR derived scar maps have the potential to inform ablation therapy—both in atrial and ventricular arrhythmias. Quantitative CMR is pushing boundaries with motion corrections in tissue characterisation and first-pass perfusion. Advanced tissue characterisation by imaging the myocardial fibre orientation using diffusion tensor imaging (DTI), has also demonstrated novel insights in patients with cardiomyopathies. Enhanced flow assessment using four-dimensional flow (4D flow) CMR, where time is the fourth dimension, allows quantification of transvalvular flow to a high degree of accuracy for all four-valves within the same cardiac cycle. This review discusses these emerging methods and others in detail and gives the reader a foresight of how CMR will evolve into a powerful clinical tool in offering a precision medicine approach to treatment, diagnosis, and detection of disease.

What Is the Clinical Impact of Stress CMR After the ISCHEMIA Trial?

After progressively receding for decades, cardiovascular mortality due to coronary artery disease has recently increased, and the associated healthcare costs are projected to double by 2030. While the 2019 European Society of Cardiology guidelines for chronic coronary syndromes recommend non-invasive cardiac imaging for patients with suspected coronary artery disease, the impact of non-invasive imaging strategies to guide initial coronary revascularization and improve long-term outcomes is still under debate. Recently, the ISCHEMIA trial has highlighted the fundamental role of optimized medical therapy and the lack of overall benefit of early invasive strategies at a median follow-up of 3.2 years. However, sub-group analyses excluding procedural infarctions with longer follow-ups of up to 5 years have suggested that patients undergoing revascularization had better outcomes than those receiving medical therapy alone. A recent sub-study of ISCHEMIA in patients with heart failure or reduced left ventricular ejection fraction (LVEF <45%) indicated that revascularization improved clinical outcomes compared to medical therapy alone. Furthermore, other large observational studies have suggested a favorable prognostic impact of coronary revascularization in patients with severe inducible ischemia assessed by stress cardiovascular magnetic resonance (CMR). Indeed, some data suggest that stress CMR-guided revascularization assessing the extent of the ischemia could be useful in identifying patients who would most benefit from invasive procedures such as myocardial revascularization. Interestingly, the MR-INFORM trial has recently shown that a first-line stress CMR-based non-invasive assessment was non-inferior in terms of outcomes, with a lower incidence of coronary revascularization compared to an initial invasive approach guided by fractional flow reserve in patients with stable angina. In the present review, we will discuss the current state-of-the-art data on the prognostic value of stress CMR assessment of myocardial ischemia in light of the ISCHEMIA trial results, highlighting meaningful sub-analyses, and still unanswered opportunities of this pivotal study. We will also review the available evidence for the potential clinical application of quantifying the extent of ischemia to stratify cardiovascular risk and to best guide invasive and non-invasive treatment strategies.

Guidelines for Cardiovascular Magnetic Resonance Imaging from the Korean Society of Cardiovascular Imaging-Part 3: Perfusion, Delayed Enhancement, and T1- and T2 Mapping

This document is the third part of the guidelines for the protocol, the interpretation and post-processing of cardiac magnetic resonance (CMR) studies. These consensus recommendations have been developed by the Consensus Committee of the Korean Society of Cardiovascular Imaging to standardize the requirements for image interpretation and post-processing of CMR. This third part of the recommendations describes tissue characterization modules, including perfusion, late gadolinium enhancement, and T1- and T2 mapping. Additionally, this document provides guidance for visual and quantitative assessment consisting of “What-to-See,” “How-To,” and common pitfalls for the analysis of each module. The Consensus Committee hopes that this document will contribute to the standardization of image interpretation and post-processing of CMR studies.

Automatic in-line quantitative myocardial perfusion mapping: Processing algorithm and implementation

PURPOSE

Quantitative myocardial perfusion mapping has advantages over qualitative assessment, including the ability to detect global flow reduction. However, it is not clinically available and remains a research tool. Building upon the previously described imaging sequence, this study presents algorithm and implementation of an automated solution for inline perfusion flow mapping with step by step performance characterization.

METHODS

Proposed workflow consists of motion correction (MOCO), arterial input function blood detection, intensity to gadolinium concentration conversion, and pixel-wise mapping. A distributed kinetics model, blood-tissue exchange model, is implemented, computing pixel-wise maps of myocardial blood flow (mL/min/g), permeability-surface-area product (mL/min/g), blood volume (mL/g), and interstitial volume (mL/g).

RESULTS

Thirty healthy subjects (11 men; 26.4 ± 10.4 years) were recruited and underwent adenosine stress perfusion cardiovascular MR. Mean MOCO quality score was 3.6 ± 0.4 for stress and 3.7 ± 0.4 for rest. Myocardial Dice similarity coefficients after MOCO were significantly improved (P < 1e-6), 0.87 ± 0.05 for stress and 0.86 ± 0.06 for rest. Arterial input function peak gadolinium concentration was 4.4 ± 1.3 mmol/L at stress and 5.2 ± 1.5 mmol/L at rest. Mean myocardial blood flow at stress and rest were 2.82 ± 0.47 mL/min/g and 0.68 ± 0.16 mL/min/g, respectively. The permeability-surface-area product was 1.32 ± 0.26 mL/min/g at stress and 1.09 ± 0.21 mL/min/g at rest (P < 1e-3). Blood volume was 12.0 ± 0.8 mL/100 g at stress and 9.7 ± 1.0 mL/100 g at rest (P < 1e-9), indicating good adenosine vasodilation response. Interstitial volume was 20.8 ± 2.5 mL/100 g at stress and 20.3 ± 2.9 mL/100 g at rest (P = 0.50).

CONCLUSIONS

An inline perfusion flow mapping workflow is proposed and demonstrated on normal volunteers. Initial evaluation demonstrates this fully automated solution for the respiratory MOCO, arterial input function left ventricle mask detection, and pixel-wise mapping, from free-breathing myocardial perfusion imaging.

Quantitative myocardial perfusion from static cardiac and dynamic arterial CT

Quantitative myocardial blood flow (MBF) estimation by dynamic contrast enhanced cardiac computed tomography (CT) requires multi-frame acquisition of contrast transit through the blood pool and myocardium to inform the arterial input and tissue response functions. Both the input and the tissue response functions for the entire myocardium are sampled with each acquisition. However, the long breath holds and frequent sampling can result in significant motion artifacts and relatively high radiation dose. To address these limitations, we propose and evaluate a new static cardiac and dynamic arterial (SCDA) quantitative MBF approach where (1) the input function is well sampled using either prediction from pre-scan timing bolus data or measured from dynamic thin slice 'bolus tracking' acquisitions, and (2) the whole-heart tissue response data is limited to one contrast enhanced CT acquisition. A perfusion model uses the dynamic arterial input function to generate a family of possible myocardial contrast enhancement curves corresponding to a range of MBF values. Combined with the timing of the single whole-heart acquisition, these curves generate a lookup table relating myocardial contrast enhancement to quantitative MBF. We tested the SCDA approach in 28 patients that underwent a full dynamic CT protocol both at rest and vasodilator stress conditions. Using measured input function plus single (enhanced CT only) or plus double (enhanced and contrast free baseline CT's) myocardial acquisitions yielded MBF estimates with root mean square (RMS) error of 1.2 ml/min/g and 0.35 ml/min/g, and radiation dose reductions of 90% and 83%, respectively. The prediction of the input function based on timing bolus data and the static acquisition had an RMS error compared to the measured input function of 26.0% which led to MBF estimation errors greater than threefold higher than using the measured input function. SCDA presents a new, simplified approach for quantitative perfusion imaging with an acquisition strategy offering substantial radiation dose and computational complexity savings over dynamic CT.

SLIC robust (SLICR) processing for fast, robust CT myocardial blood flow quantification

There are several computational methods for estimating myocardial blood flow (MBF) using CT myocardial perfusion imaging (CT-MPI). Previous work has shown that model-based deconvolution methods are more accurate and precise than model-independent methods such as singular value decomposition and max-upslope. However, iterative optimization is computationally expensive and models are sensitive to image noise, thus limiting the utility of low x-ray dose acquisitions. We propose a new processing method, SLICR, which segments the myocardium into super-voxels using a modified simple linear iterative clustering (SLIC) algorithm and quantifies MBF via a robust physiologic model (RPM). We compared SLICR against voxel-wise SVD and voxel-wise model-based deconvolution methods (RPM, single-compartment and Johnson-Wilson). We used image data from a digital CT-MPI phantom to evaluate robustness of processing methods to noise at reduced x-ray dose. We validate SLICR in a porcine model with and without partial occlusion of the LAD coronary artery with known pressure-wire fractional flow reserve. SLICR was ~50 times faster than voxel-wise RPM and other model-based methods while retaining sufficient resolution to show all clinically interesting features (e.g., a flow deficit in the endocardial wall). SLICR showed much better precision and accuracy than the other methods. For example, at simulated MBF=100 mL/min/100g and 100 mAs exposure (50% of nominal dose) in the digital simulator, MBF estimates were 101 ± 12 mL/min/100g, 160 ± 54 mL/min/100g, and 122 ± 99 mL/min/100g for SLICR, SVD, and Johnson-Wilson, respectively. SLICR even gave excellent results (103 ± 23 ml/min/100g) at 50 mAs, corresponding to 25% nominal dose.

Estimating extraction fraction and blood flow by combining first-pass myocardial perfusion and T1 mapping results

Background

Quantifying myocardial perfusion is complicated by the complexity of pharmacokinetic model being used and the reliability of perfusion parameter estimates. More complex modeling provides more information about the underlying physiology, but too many parameters in complex models introduce a new problem of reliable estimation. To overcome the problem of multiple parameters, we have developed a technique that combines knowledge from two different cardiac magnetic resonance (MR) imaging techniques: dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) and T1 mapping. Using extracellular volume (ECV) estimates from T1 mapping may allow more robust model parameter estimates.

Methods

Simulations and human scans were performed. The myocardial perfusion scans used an ungated saturation recovery prepared TurboFLASH pulse sequence. Four short-axis (SA) slices were acquired after a single saturation pulse with a saturation recovery time of ~25 ms before the first slice. Gadoteridol was injected and ~240 frames were acquired over a minute with shallow breathing and no electrocardiograph (ECG) gating. This was followed 20±5 minutes later by an injection of regadenoson to induce hyperemia. The data were acquired using an under-sampled golden angle radial acquisition. Modified look-locker inversion recovery (MOLLI) T1 mapping was performed in 3 slices pre- and post-contrast. The pre- and post-contrast T1 maps were used for ECV estimation. Quantification of perfusion was done using a 4-parameter model with additional information about ECV supplied during model fitting. Phase contrast scans of the coronary sinus (CS) were acquired at rest and immediately after the stress perfusion acquisition to estimate global flow.

Results

Without ECV information, the 5-parameter model fails to converge to a unique solution and often gives incorrect estimates for the perfusion parameters. The myocardial blood flow (MBF) estimates during rest and stress were 0.9±0.1 and 2.3±0.6 mL/min/g, respectively. The extraction fraction estimates were 0.49±0.04 and 0.34±0.05 during rest and stress, respectively.

Conclusions

These results show that it is possible to successfully fit a dynamic perfusion model with an extraction fraction parameter by using information from T1 mapping scans. This hybrid approach is especially important when the 5-parameter model alone fails to converge on a unique solution. This work is a good example of exploiting information overlaps between various cardiac MR imaging techniques.

Introduction to Cardiovascular Magnetic Resonance: Technical Principles and Clinical Applications

Cardiovascular magnetic resonance (CMR) is a set of magnetic resonance imaging (MRI) techniques designed to assess cardiovascular morphology, ventricular function, myocardial perfusion, tissue characterization, flow quantification and coronary artery disease. Since MRI is a non-invasive tool and free of radiation, it is suitable for longitudinal monitoring of treatment effect and follow-up of disease progress. Compared to MRI of other body parts, CMR faces specific challenges from cardiac and respiratory motion. Therefore, CMR requires synchronous cardiac and respiratory gating or breath-holding techniques to overcome motion artifacts. This article will review the basic principles of MRI and introduce the CMR techniques that can be optimized for enhanced clinical assessment.

KEY WORDS

Cardiovascular MR • Coronary arteries • Flow quantification • Myocardial fibrosis • Myocardial perfusion • Myocardial scarring • Regional wall motion • Ventricular function.

Quantification of myocardial blood flow using non-electrocardiogram-triggered MRI with three-slice coverage

PURPOSE

Accurate quantification of myocardial perfusion is dependent on reliable electrocardiogram (ECG) triggering. Measuring myocardial blood flow (MBF) in patients with arrhythmias or poor ECGs is currently infeasible with MR. The purpose of this study was to demonstrate the feasibility of a non-ECG-triggered method with clinically useful three-slice ventricular coverage for measurement of MBF in healthy volunteers.

METHODS

A saturation recovery magnetization-prepared gradient recalled echo acquisition was continuously repeated during first-pass imaging. A slice-interleaved radial trajectory was employed to enable image-based retrospective triggering. The arterial input function was generated using a beat-by-beat T1 estimation method. The proposed technique was validated against a conventional ECG-triggered dual-bolus technique in 10 healthy volunteers. The technique was further demonstrated under adenosine stress in 12 healthy volunteers.

RESULTS

The proposed method produced MBF with no significant difference compared with the ECG-triggered technique. The proposed method yielded mean myocardial perfusion reserve comparable to published literature.

CONCLUSION

We have developed a non-ECG-triggered quantitative perfusion imaging method. In this preliminary study, our results demonstrate that our method yields comparable MBF compared with the conventional ECG-triggered method and that it is feasible for stress imaging.

Simulation of contrast agent transport in arteries with multilayer arterial wall: impact of arterial transmural transport on the bolus delay and dispersion

One assumption of DSC-MRI is that the injected contrast agent is kept totally intravascular and the arterial wall is impermeable to contrast agent. The assumption is unreal for such small contrast agent as Gd-DTPA can leak into the arterial wall. To investigate whether the unreal assumption is valid for the estimation of the delay and dispersion of the contrast agent bolus, we simulated flow and Gd-DTPA transport in a model with multilayer arterial wall and analyzed the bolus delay and dispersion qualified by mean vascular transit time (MVTT) and the variance of the vascular transport function. Factors that may affect Gd-DTPA transport hence the delay and dispersion were further investigated, such as integrity of endothelium and disturbed flow. The results revealed that arterial transmural transport would slightly affect MVTT and moderately increase the variance. In addition, although the integrity of endothelium can significantly affect the accumulation of contrast agent in the arterial wall, it had small effects on the bolus delay and dispersion. However, the disturbed flow would significantly increase both MVTT and the variance. In conclusion, arterial transmural transport may have a small effect on the bolus delay and dispersion when compared to the flow pattern in the artery.

Quantification of myocardial blood flow using non-ECG-triggered MR imaging

PURPOSE

MR myocardial perfusion imaging is dependent on reliable electrocardiogram (ECG) triggering for accurate measurement of myocardial blood flow (MBF). A non-ECG-triggered method for quantitative first-pass imaging may improve clinical feasibility in patients with poor ECG signal. The purpose of this study is to evaluate the feasibility of a non-ECG-triggered method for myocardial perfusion imaging in a single slice.

METHODS

The proposed non-ECG-triggered technique uses a saturation-recovery magnetization preparation and golden-angle radial acquisition for integrated arterial input function (AIF) measurement. Image based self-gating with a temporal resolution of 42.6 ms is used to generate a first-pass image series with consistent cardiac phase. The AIF is measured using beat-by-beat T1 estimation of the ventricular blood pool. The proposed technique was performed on 14 healthy volunteers and compared against a conventional ECG-triggered dual-bolus acquisition.

RESULTS

The proposed method produced MBF with no significant difference compared with ECG-triggered technique (mean of 0.63 ± 0.22 mL/min/g to 0.73 ± 0.21 mL/min/g).

CONCLUSION

We have developed a non-ECG-triggered perfusion imaging method with T1 based measurement of the AIF in a single slice. In this preliminary study, our results demonstrate that MBF measured using the proposed method is comparable to the conventional ECG-triggered method.

Free-breathing cardiac MR stress perfusion with real-time slice tracking

PURPOSE

To develop a free-breathing cardiac MR perfusion sequence with slice tracking for use after physical exercise.

METHODS

We propose to use a leading navigator, placed immediately before each 2D slice acquisition, for tracking the respiratory motion and updating the slice location in real-time. The proposed sequence was used to acquire CMR perfusion datasets in 12 healthy adult subjects and 8 patients. Images were compared with the conventional perfusion (i.e., without slice tracking) results from the same subjects. The location and geometry of the myocardium were quantitatively analyzed, and the perfusion signal curves were calculated from both sequences to show the efficacy of the proposed sequence.

RESULTS

The proposed sequence was significantly better compared with the conventional perfusion sequence in terms of qualitative image scores. Changes in the myocardial location and geometry decreased by 50% in the slice tracking sequence. Furthermore, the proposed sequence had signal curves that are smoother and less noisy.

CONCLUSION

The proposed sequence significantly reduces the effect of the respiratory motion on the image acquisition in both rest and stress perfusion scans.

Uncertainty in MR tracer kinetic parameters and water exchange rates estimated from T1-weighted dynamic contrast enhanced MRI

PURPOSE

The aim of this study was to assess the uncertainty in estimation of MR tracer kinetic parameters and water exchange rates in T1-weighted dynamic contrast enhanced (DCE) MRI.

METHODS

Simulated DCE-MRI data were used to assess four kinetic models; general kinetic model with a vascular compartment (GKM2), GKM2 combined with water exchange (SSM2), adiabatic approximation of the tissue homogeneity model (ATH), and ATH combined with water exchange (ATHX).

RESULTS

In GKM2 and SSM2, increase in transfer constant (K(trans)) led to underestimation of vascular volume fraction (vb), and increase in vb led to overestimation of K(trans). Such coupling between K(trans) and vb was not observed in ATH and ATHX. The precision of estimated intracellular water lifetime (τi) was substantially improved in both SSM2 and ATHX when K(trans) > 0.3 min(-1). K(trans) and vb from ATHX model had significantly smaller errors than those from ATH model (P < 0.05).

CONCLUSION

The results of this study demonstrated the feasibility of measuring τi from DCE-MRI data albeit low precision. While the inclusion of water exchange improved the accuracy of K(trans), vb, and the interstitial volume fraction estimation (ve), it lowered the precision of other kinetic model parameters within the conditions investigated in this study.

The role of capillary transit time heterogeneity in myocardial oxygenation and ischemic heart disease

Ischemic heart disease (IHD) is characterized by an imbalance between oxygen supply and demand, most frequently caused by coronary artery disease (CAD) that reduces myocardial perfusion. In some patients, IHD is ascribed to microvascular dysfunction (MVD): microcirculatory disturbances that reduce myocardial perfusion at the level of myocardial pre-arterioles and arterioles. In a minority of cases, chest pain and reductions in myocardial flow reserve may even occur in patients without any other demonstrable systemic or cardiac disease. In this topical review, we address whether these findings might be caused by impaired myocardial oxygen extraction, caused by capillary flow disturbances further downstream. Myocardial blood flow (MBF) increases approximately linearly with oxygen utilization, but efficient oxygen extraction at high MBF values is known to depend on the parallel reduction of capillary transit time heterogeneity (CTH). Consequently, changes in capillary wall morphology or blood viscosity may impair myocardial oxygen extraction by preventing capillary flow homogenization. Indeed, a recent re-analysis of oxygen transport in tissue shows that elevated CTH can reduce tissue oxygenation by causing a functional shunt of oxygenated blood through the tissue. We review the combined effects of MBF, CTH, and tissue oxygen tension on myocardial oxygen supply. We show that as CTH increases, normal vasodilator responses must be attenuated in order to reduce the degree of functional shunting and improve blood-tissue oxygen concentration gradients to allow sufficient myocardial oxygenation. Theoretically, CTH can reach levels such that increased metabolic demands cannot be met, resulting in tissue hypoxia and angina in the absence of flow-limiting CAD or MVD. We discuss these predictions in the context of MVD, myocardial infarction, and reperfusion injury.

Comparison of blood flow models and acquisitions for quantitative myocardial perfusion estimation from dynamic CT

Myocardial blood flow (MBF) can be estimated from dynamic contrast enhanced (DCE) cardiac CT acquisitions, leading to quantitative assessment of regional perfusion. The need for low radiation dose and the lack of consensus on MBF estimation methods motivates this study to refine the selection of acquisition protocols and models for CT-derived MBF. DCE cardiac CT acquisitions were simulated for a range of flow states (MBF = 0.5, 1, 2, 3 ml (min g)(-1), cardiac output = 3, 5, 8 L min(-1)). Patient kinetics were generated by a mathematical model of iodine exchange incorporating numerous physiological features including heterogenenous microvascular flow, permeability and capillary contrast gradients. CT acquisitions were simulated for multiple realizations of realistic x-ray flux levels. CT acquisitions that reduce radiation exposure were implemented by varying both temporal sampling (1, 2, and 3 s sampling intervals) and tube currents (140, 70, and 25 mAs). For all acquisitions, we compared three quantitative MBF estimation methods (two-compartment model, an axially-distributed model, and the adiabatic approximation to the tissue homogeneous model) and a qualitative slope-based method. In total, over 11 000 time attenuation curves were used to evaluate MBF estimation in multiple patient and imaging scenarios. After iodine-based beam hardening correction, the slope method consistently underestimated flow by on average 47.5% and the quantitative models provided estimates with less than 6.5% average bias and increasing variance with increasing dose reductions. The three quantitative models performed equally well, offering estimates with essentially identical root mean squared error (RMSE) for matched acquisitions. MBF estimates using the qualitative slope method were inferior in terms of bias and RMSE compared to the quantitative methods. MBF estimate error was equal at matched dose reductions for all quantitative methods and range of techniques evaluated. This suggests that there is no particular advantage between quantitative estimation methods nor to performing dose reduction via tube current reduction compared to temporal sampling reduction. These data are important for optimizing implementation of cardiac dynamic CT in clinical practice and in prospective CT MBF trials.

Uncertainty in MR tracer kinetic parameters and water exchange rates estimated from T1-weighted dynamic contrast enhanced MRI

PURPOSE

The aim of this study was to assess the uncertainty in estimation of MR tracer kinetic parameters and water exchange rates in T1-weighted dynamic contrast enhanced (DCE) MRI.

METHODS

Simulated DCE-MRI data were used to assess four kinetic models; general kinetic model with a vascular compartment (GKM2), GKM2 combined with water exchange (SSM2), adiabatic approximation of the tissue homogeneity model (ATH), and ATH combined with water exchange (ATHX).

RESULTS

In GKM2 and SSM2, increase in transfer constant (K(trans)) led to underestimation of vascular volume fraction (vb), and increase in vb led to overestimation of K(trans). Such coupling between K(trans) and vb was not observed in ATH and ATHX. The precision of estimated intracellular water lifetime (τi) was substantially improved in both SSM2 and ATHX when K(trans) > 0.3 min(-1). K(trans) and vb from ATHX model had significantly smaller errors than those from ATH model (P < 0.05).

CONCLUSION

The results of this study demonstrated the feasibility of measuring τi from DCE-MRI data albeit low precision. While the inclusion of water exchange improved the accuracy of K(trans), vb, and the interstitial volume fraction estimation (ve), it lowered the precision of other kinetic model parameters within the conditions investigated in this study.

Computational fluid dynamics simulations of contrast agent bolus dispersion in a coronary bifurcation: impact on MRI-based quantification of myocardial perfusion

Contrast-enhanced first-pass magnetic resonance imaging (MRI) in combination with a tracer kinetic model, for example, MMID4, can be used to determine myocardial blood flow (MBF) and myocardial perfusion reserve (MPR). Typically, the arterial input function (AIF) required for this methodology is estimated from the left ventricle (LV). Dispersion of the contrast agent bolus might occur between the LV and the myocardial tissue. Negligence of bolus dispersion could cause an error in MBF determination. The aim of this study was to investigate the influence of bolus dispersion in a simplified coronary bifurcation geometry including one healthy and one stenotic branch on the quantification of MBF and MPR. Computational fluid dynamics (CFD) simulations were combined with MMID4. Different inlet boundary conditions describing pulsatile and constant flows for rest and hyperemia and differing outflow conditions have been investigated. In the bifurcation region, the increase of the dispersion was smaller than inside the straight vessels. A systematic underestimation of MBF values up to -16.1% for pulsatile flow and an overestimation of MPR up to 7.5% were found. It was shown that, under the conditions considered in this study, bolus dispersion can significantly influence the results of quantitative myocardial MR-perfusion measurements.

Development of a universal dual-bolus injection scheme for the quantitative assessment of myocardial perfusion cardiovascular magnetic resonance

BACKGROUND

The dual-bolus protocol enables accurate quantification of myocardial blood flow (MBF) by first-pass perfusion cardiovascular magnetic resonance (CMR). However, despite the advantages and increasing demand for the dual-bolus method for accurate quantification of MBF, thus far, it has not been widely used in the field of quantitative perfusion CMR. The main reasons for this are that the setup for the dual-bolus method is complex and requires a state-of-the-art injector and there is also a lack of post processing software. As a solution to one of these problems, we have devised a universal dual-bolus injection scheme for use in a clinical setting. The purpose of this study is to show the setup and feasibility of the universal dual-bolus injection scheme.

METHODS

The universal dual-bolus injection scheme was tested using multiple combinations of different contrast agents, contrast agent dose, power injectors, perfusion sequences, and CMR scanners. This included 3 different contrast agents (Gd-DO3A-butrol, Gd-DTPA and Gd-DOTA), 4 different doses (0.025 mmol/kg, 0.05 mmol/kg, 0.075 mmol/kg and 0.1 mmol/kg), 2 different types of injectors (with and without "pause" function), 5 different sequences (turbo field echo (TFE), balanced TFE, k-space and time (k-t) accelerated TFE, k-t accelerated balanced TFE, turbo fast low-angle shot) and 3 different CMR scanners from 2 different manufacturers. The relation between the time width of dilute contrast agent bolus curve and cardiac output was obtained to determine the optimal predefined pause duration between dilute and neat contrast agent injection.

RESULTS

161 dual-bolus perfusion scans were performed. Three non-injector-related technical errors were observed (1.9%). No injector-related errors were observed. The dual-bolus scheme worked well in all the combinations of parameters if the optimal predefined pause was used. Linear regression analysis showed that the optimal duration for the predefined pause is 25s to separate the dilute and neat contrast agent bolus curves if 0.1 mmol/kg dose of Gd-DO3A-butrol is used.

CONCLUSION

The universal dual-bolus injection scheme does not require sophisticated double-head power injector function and is a feasible technique to obtain reasonable arterial input function curves for absolute MBF quantification.

Comparison of myocardial perfusion estimates from dynamic contrast-enhanced magnetic resonance imaging with four quantitative analysis methods

Dynamic contrast-enhanced MRI has been used to quantify myocardial perfusion in recent years. Published results have varied widely, possibly depending on the method used to analyze the dynamic perfusion data. Here, four quantitative analysis methods (two-compartment modeling, Fermi function modeling, model-independent analysis, and Patlak plot analysis) were implemented and compared for quantifying myocardial perfusion. Dynamic contrast-enhanced MRI data were acquired in 20 human subjects at rest with low-dose (0.019 +/- 0.005 mmol/kg) bolus injections of gadolinium. Fourteen of these subjects were also imaged at adenosine stress (0.021 +/- 0.005 mmol/kg). Aggregate rest perfusion estimates were not significantly different between all four analysis methods. At stress, perfusion estimates were not significantly different between two-compartment modeling, model-independent analysis, and Patlak plot analysis. Stress estimates from the Fermi model were significantly higher (approximately 20%) than the other three methods. Myocardial perfusion reserve values were not significantly different between all four methods. Model-independent analysis resulted in the lowest model curve-fit errors. When more than just the first pass of data was analyzed, perfusion estimates from two-compartment modeling and model-independent analysis did not change significantly, unlike results from Fermi function modeling.

Quantification of myocardial perfusion by cardiovascular magnetic resonance

The potential of contrast-enhanced cardiovascular magnetic resonance (CMR) for a quantitative assessment of myocardial perfusion has been explored for more than a decade now, with encouraging results from comparisons with accepted "gold standards", such as microspheres used in the physiology laboratory. This has generated an increasing interest in the requirements and methodological approaches for the non-invasive quantification of myocardial blood flow by CMR. This review provides a synopsis of the current status of the field, and introduces the reader to the technical aspects of perfusion quantification by CMR. The field has reached a stage, where quantification of myocardial perfusion is no longer a claim exclusive to nuclear imaging techniques. CMR may in fact offer important advantages like the absence of ionizing radiation, high spatial resolution, and an unmatched versatility to combine the interrogation of the perfusion status with a comprehensive tissue characterization. Further progress will depend on successful dissemination of the techniques for perfusion quantification among the CMR community.

Nonparametric Residue Analysis of Dynamic PET Data With Application to Cerebral FDG Studies in Normals

Kinetic analysis is used to extract metabolic information from dynamic positron emission tomography (PET) uptake data. The theory of indicator dilutions, developed in the seminal work of Meier and Zierler (1954), provides a probabilistic framework for representation of PET tracer uptake data in terms of a convolution between an arterial input function and a tissue residue. The residue is a scaled survival function associated with tracer residence in the tissue. Nonparametric inference for the residue, a deconvolution problem, provides a novel approach to kinetic analysis-critically one that is not reliant on specific compartmental modeling assumptions. A practical computational technique based on regularized cubic B-spline approximation of the residence time distribution is proposed. Nonparametric residue analysis allows formal statistical evaluation of specific parametric models to be considered. This analysis needs to properly account for the increased flexibility of the nonparametric estimator. The methodology is illustrated using data from a series of cerebral studies with PET and fluorodeoxyglucose (FDG) in normal subjects. Comparisons are made between key functionals of the residue, tracer flux, flow, etc., resulting from a parametric (the standard two-compartment of Phelps et al. 1979) and a nonparametric analysis. Strong statistical evidence against the compartment model is found. Primarily these differences relate to the representation of the early temporal structure of the tracer residence-largely a function of the vascular supply network. There are convincing physiological arguments against the representations implied by the compartmental approach but this is the first time that a rigorous statistical confirmation using PET data has been reported. The compartmental analysis produces suspect values for flow but, notably, the impact on the metabolic flux, though statistically significant, is limited to deviations on the order of 3%-4%. The general advantage of the nonparametric residue analysis is the ability to provide a valid kinetic quantitation in the context of studies where there may be heterogeneity or other uncertainty about the accuracy of a compartmental model approximation of the tissue residue.

Quantitative analysis of dynamic contrast-enhanced MR images based on Bayesian P-splines

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) is an important tool for detecting subtle kinetic changes in cancerous tissue. Quantitative analysis of DCE-MRI typically involves the convolution of an arterial input function (AIF) with a nonlinear pharmacokinetic model of the contrast agent concentration. Parameters of the kinetic model are biologically meaningful, but the optimization of the nonlinear model has significant computational issues. In practice, convergence of the optimization algorithm is not guaranteed and the accuracy of the model fitting may be compromised. To overcome these problems, this paper proposes a semi-parametric penalized spline smoothing approach, where the AIF is convolved with a set of B-splines to produce a design matrix using locally adaptive smoothing parameters based on Bayesian penalized spline models (P-splines). It has been shown that kinetic parameter estimation can be obtained from the resulting deconvolved response function, which also includes the onset of contrast enhancement. Detailed validation of the method, both with simulated and in vivo data, is provided.

Attenuation resilient AIF estimation based on hierarchical Bayesian modelling for first pass myocardial perfusion MRI

Non-linear attenuation of the Arterial Input Function (AIF) is a major problem in first-pass MR perfusion imaging due to the high concentration of the contrast agent in the blood pool. This paper presents a technique to reconstruct the true AIF using signal intensities in the myocardium and the attenuated AIF based on a Hierarchical Bayesian Model (HBM). With the proposed method, both the AIF and the response function are modeled as smoothed functions by using Bayesian penalty splines (P-Splines). The derived AIF is then used to estimate the impulse response of the myocardium based on deconvolution analysis. The proposed technique is validated both with simulated data using the MMID4 model and ten in vivo data sets for estimating myocardial perfusion reserve rates. The results demonstrate the ability of the proposed technique in accurately reconstructing the desired AIF for myocardial perfusion quantification. The method does not involve any MRI pulse sequence modification, and thus is expected to have wider clinical impact.

Renal function measurements from MR renography and a simplified multicompartmental model

The purpose of this study was to determine the accuracy and sources of error in estimating single-kidney glomerular filtration rate (GFR) derived from low-dose gadolinium-enhanced T1-weighted MR renography. To analyze imaging data, MR signal intensity curves were converted to concentration vs. time curves, and a three-compartment, six-parameter model of the vascular-nephron system was used to analyze measured aortic, cortical, and medullary enhancement curves. Reliability of the parameter estimates was evaluated by sensitivity analysis and by Monte Carlo analyses of model solutions to which random noise had been added. The dominant sensitivity of the medullary enhancement curve to GFR 1–4 min after tracer injection was supported by a low coefficient of variation in model-fit GFR values (4%) when measured data were subjected to 5% noise. These analyses also showed the minimal effects of bolus dispersion in the aorta on parameter reliability. Single-kidney GFR from MR renography analyzed by the three-compartment model (4.0–71.4 ml/min) agreed well with reference measurements from 99mTc-DTPA clearance and scintigraphy ( r = 0.84, P < 0.001). Bland-Altman analysis showed an average difference of 11.9 ml/min (95% confidence interval = 5.8–17.9 ml/min) between model and reference values. We conclude that a nephron-based multicompartmental model can be used to derive clinically useful estimates of single-kidney GFR from low-dose MR renography.

The delay of contrast arrival in magnetic resonance first-pass perfusion imaging: a novel non-invasive parameter detecting collateral-dependent myocardium

AIM

To establish the regional delay of contrast arrival in magnetic resonance perfusion imaging (MRPI) for the detection of collateral-dependent myocardium in patients with coronary artery disease.

DESIGN AND SETTING

Observational study, case series; single centre, university hospital.

PATIENTS

30 patients with coronary artery disease and collateral-dependent myocardium and 17 healthy volunteers.

METHODS

Resting and hyperaemic (adenosine) MRPI was used to determine the delay time (Deltat(d)) of contrast arrival between the left ventricle and collateral-dependent or antegradely perfused myocardium, and myocardial perfusion (MP, ml/min/g).

RESULTS

In healthy volunteers, mean (SD) Deltat(d) at rest and during hyperaemia were 0.8 (0.4) and 0.3 (0.3) s, and MP was 1.14 (0.21) and 4.23 (1.12) ml/min/g. In patients Deltat(d) in antegradely perfused vs collateral-dependent myocardium was 0.9 (0.7) vs 1.7 (1.0) s at rest (p<0.001), and 0.4 (0.3) vs 1.1 (0.6) s (p<0.001) during hyperaemia. MP was 1.12 (0.11) and 0.98 (0.28) ml/min/g (p = NS) at rest and 2.46 (0.85) vs 1.86 (0.91) ml/min/g (p<0.01) during hyperaemia. Receiver operating characteristics analysis showed the best sensitivity and specificity of 90% and 83% for hyperaemic Deltat(d) of >0.6 s (area under the curve (AUC) = 0.89) to detect collateral-dependent myocardium, while resting Deltat(d) (AUC = 0.77) and perfusion (AUC = 0.69 at rest or 0.70 during hyperaemia) were less accurate.

CONCLUSIONS

MRPI-derived hyperaemic delay of contrast arrival detects collateral-dependent myocardium with high sensitivity and specificity. Perfusion was less sensitive, emphasising the clinical role of Deltat(d) in non-invasive detection of collateral-dependent myocardium.

High-resolution dynamic imaging of contrast agent uptake in a beating heart

Dynamic contrast-enhanced (DCE) MR cardiac imaging has been recognized as a unique and powerful tool for assessing both cardiac functions and physiological conditions of the heart tissues (e.g., tissue rejection following heart transplantation). However, because of cardiac motion and the limited data acquisition speed of existing MRI techniques, it has been very difficult to acquire dynamic images of very high spatiotemporal resolution. This paper proposes a new generalized series (GS) based imaging technique to overcome this challenging problem. Specifically, the proposed technique collects two data sets: a) a sequence of highresolution reference images over several cardiac cycles using a gated cine acquisition scheme before the injection of a contrast agent (or a molecular probe), and b) a sequence of reduced data sets with very high frame rate during the transient wash-in/wash-out stage of the contrast agent. A GS model is then used to combine these two data sets to reconstruct a high-resolution image sequence, capturing both the cardiac motions and dynamic signal changes due to the interaction of the contrast agent with the cardiac tissues. The proposed technique has been validated using both simulated and experimental data, which show that high-resolution dynamic images can be acquired with as few as 8 encodings (in contrast to 256 encodings required in the traditional Fourier transform-based methods). The technique provides a very effective tool for physiological imaging of the beating heart with molecular probes.

Myocardial perfusion imaging by cardiac magnetic resonance

Cardiovascular magnetic resonance (CMR) has been shown to provide high quality data on cardiac and valvular function, perfusion, viability, blood flow, and potentially, on cardiac metabolism as well. Several of these CMR applications (eg, function and viability assessment) matured during the past years and are now established components of a cardiac workup. Perfusion-CMR is close to this status and is already a major contributor to cardiac examinations in a growing number of expert centers. Large multicenter perfusion-CMR trials comparing the diagnostic performance of CMR with other techniques were recently reported yielding areas under the receiver-operator-characteristics curve as a high as 0.85 for coronary artery disease detection (MR-IMPACT). Anticipating a growing role for perfusion-CMR in cardiology in the near future, this article discusses the principles of perfusion-CMR and its integration into the workup of patient with coronary artery disease (CAD). In addition to a functional study, this integration is mainly composed of a perfusion-CMR part, followed by a viability assessment by late enhancement CMR techniques. The principal characteristics of these CMR techniques are compared with those of single photon emission computed tomography (SPECT) and positron emission tomography (PET). After introduction into principles and techniques of perfusion-CMR, some open questions in perfusion-CMR and challenges for the future are addressed. Finally, newer CMR applications are shortly mentioned utilizing hyperpolarized carbon-13 compounds in experimental models for quantification of myocardial perfusion and for real-time assessment of metabolic pathways in postischemic myocardium.

A comparison of myocardial perfusion and rejection in cardiac transplant patients

INTRODUCTION

Although histological evaluation of the cardiac tissue is the current gold standard for evaluation of rejection, we hypothesized that cardiac perfusion MRI is a safe non-invasive method that correlates tissue blood flow changes with biopsy proven rejection in the cardiac transplant patient.

MATERIALS AND METHODS

In a retrospective study from 1984-2001, 83 patients underwent 135 MR Gd-DTPA imaging studies. In 8 patients (9%), biopsies graded 2 or higher (by ISHLT criteria) provided evidence of rejection. Patients were age and sex matched to 11 non-rejected controls for imaging analysis. Time-signal intensity curves generated for a mid-ventricle LV short axis slice during rest and adenosine stress allowed determination of myocardial blood flow (MBF, ml/min/gm). ROC curve analysis by SPSS allowed estimation of sensitivity and specificity.

RESULTS

At rest, there was no difference in MBF between patients with prior rejection vs. those without (1.18 +/- 0.26 vs. 1.16 +/- 0.29). At stress there was a decrease in MBF for patients with prior rejection episodes (3.27 +/- 0.74) compared to no rejection (3.60 +/- 0.72), P = 0.067). The area under the ROC curve was 0.82, with specificity and sensitivity of 75% and 81%, respectively.

CONCLUSION

This study suggests that perfusion MR imaging can be used in assessing the cardiac transplant patient for rejection related microvascular changes. The high specificity and sensitivity recorded from the ROC curve illustrates the potential utility of this diagnostic test for future studies.

MRI of Myocardial Perfusion

An overwhelming number of myocardial perfusion studies are done by nuclear isotope imaging. Magnetic resonance imaging during the first pass of an injected, contrast bolus has some significant advantages for detection of blood flow deficits, namely higher spatial resolution, absence of ionizing radiation, and speed of the test. Previous clinical studies have demonstrated that excellent sensitivity and specificity can be achieved with MR myocardial perfusion imaging for detecting coronary artery disease, and assessment of patients with acute chest pain. Furthermore, an absolute quantification of myocardial blood flow is feasible, as was demonstrated by comparison of MR perfusion imaging, to measurements with isotope labeled microspheres in experimental models. An integrated assessment of perfusion, function, and viability, is thus feasible by MRI to answer important clinical challenges such as the identification of stunned or hibernating, but viable myocardium.

Study of an adaptive bolus chasing CT angiography

To improve imaging quality and to reduce contrast dose and radiation exposure, an adaptive bolus chasing CT angiography was proposed so that the bolus peak position and the imaging aperture can be synchronized. The performance of the proposed adaptive bolus chasing CT angiography was experimentally evaluated based on the actual bolus dynamics. The experimental results show that the controlled table position and the bolus peak position were highly consistent. The results clearly demonstrate that the proposed adaptive bolus chasing CT angiography that synchronizes the bolus peak position with the imaging aperture by a simple adaptive system is computationally and clinically feasible. Similar techniques may also be applied to conventional angiography to improve imaging quality and to reduce contrast dose and/or radiation exposure.

Multiscale modeling of cardiac cellular energetics

Multiscale modeling is essential to integrating knowledge of human physiology starting from genomics, molecular biology, and the environment through the levels of cells, tissues, and organs all the way to integrated systems behavior. The lowest levels concern biophysical and biochemical events. The higher levels of organization in tissues, organs, and organism are complex, representing the dynamically varying behavior of billions of cells interacting together. Models integrating cellular events into tissue and organ behavior are forced to resort to simplifications to minimize computational complexity, thus reducing the model's ability to respond correctly to dynamic changes in external conditions. Adjustments at protein and gene regulatory levels shortchange the simplified higher-level representations. Our cell primitive is composed of a set of subcellular modules, each defining an intracellular function (action potential, tricarboxylic acid cycle, oxidative phosphorylation, glycolysis, calcium cycling, contraction, etc.), composing what we call the "eternal cell," which assumes that there is neither proteolysis nor protein synthesis. Within the modules are elements describing each particular component (i.e., enzymatic reactions of assorted types, transporters, ionic channels, binding sites, etc.). Cell subregions are stirred tanks, linked by diffusional or transporter-mediated exchange. The modeling uses ordinary differential equations rather than stochastic or partial differential equations. This basic model is regarded as a primitive upon which to build models encompassing gene regulation, signaling, and long-term adaptations in structure and function. During simulation, simpler forms of the model are used, when possible, to reduce computation. However, when this results in error, the more complex and detailed modules and elements need to be employed to improve model realism. The processes of error recognition and of mapping between different levels of model form complexity are challenging but are essential for successful modeling of large-scale systems in reasonable time. Currently there is to this end no established methodology from computational sciences.

Time delay for arrival of MR contrast agent in collateral-dependent myocardium

An analysis of the kinetics of myocardial contrast enhancement is an important component of myocardial perfusion studies. The contrast enhancement can be modeled by a linear time-invariant system, and the myocardial impulse response, calculated by deconvolution of the measured tissue response with an arterial input, gives a direct estimate of myocardial blood flow. In this paper, we analyze the effects of delays in the contrast enhancement, that occur in collateral-dependent myocardium, where the tracer reaches the tissue region only through branches from other coronary arteries that form natural bypass vessels. We investigate how the delayed arrival of tracer alters the myocardial impulse response. Model-independent deconvolution is applied to determine the lag between arterial input and tissue enhancement. Experimental data in a porcine model of collateral development indicate that the delayed arrival of an injected tracer, measured at rest, is a useful marker to identify collateral-dependent myocardium, and predict its flow capacitance.

Breast cancer: regional blood flow and blood volume measured with magnetic susceptibility-based MR imaging--initial results

The purpose of this study was to quantify microcirculation in breast neoplasms with magnetic susceptibility-based contrast material-enhanced magnetic resonance imaging. With this imaging method for invasive cancers, the mean values of the ratios of tumor to normal blood flow and blood volume were significantly higher (P <.002) than those for benign or normal tissue. The method allows independent measurement of regional blood flow and blood volume in breast cancers.

The mechanical and metabolic basis of myocardial blood flow heterogeneity

Precise measurements of regional myocardial blood flow heterogeneity had to be developed before one could seek causation for the heterogeneity. Deposition techniques (particles or molecular microspheres) are the most precise, but imaging techniques have begun to provide high enough resolution to allow in vivo studies. Assigning causation has been difficult. There is no apparent association with the regional concentrations of energy-related enzymes or substrates, but these are measures of status, not of metabolism. There is statistical correlation between flow and regional substrate uptake and utilization. Attribution of regional flow variation to vascular anatomy or to vasomotor control appears not to be causative on a long-term basis. The closest relationships appear to be with mechanical function, but one cannot say for sure whether this is related to ATP hydrolysis at the crossbridge or associated metabolic reactions such as calcium uptake by the sarcoplasmic reticulum.

Direct comparison of an intravascular and an extracellular contrast agent for quantification of myocardial perfusion. Cardiac MRI Group

A direct comparison of extracellular and intravascular contrast agents for the assessment of myocardial perfusion was carried out in a porcine model (N = 5) with a flow-limiting occluder on the left anterior descending coronary artery. Rapid imaging during the first pass of an extracellular or intravascular contrast agent with a saturation-recovery-prepared TurboFLASH sequence showed comparable peak contrast-to-noise enhancements in myocardial tissue regions with flows averaging 1.1 +/- 0.2 at baseline to 4.8 +/- 0.6 ml/min/g during hyperemia. The coefficient of variation between the MR estimates of blood flow with Gadomer-17 and the microsphere blood flow measurements was 11 +/- 11%, while the corresponding co-efficient of variation for blood flow estimates with the extracellular CA was 23 +/- 11%. Blood volume differences between rest and hyperemia observed with the intravascular tracer were significant (Vvasc(rest) = 0.078 +/- 0.013 ml/g, versus Vvasc(hyperemia) = 0.102 +/- 0.019 ml/g; p < 0.05). The effects of water exchange were minimized through the choice of pulse sequence parameters to provide blood volume estimates consistent with the changes expected between rest and hyperemia. This study represents the first application of multiple indicators in first pass imaging studies for the assessment of myocardial perfusion. The use of an intravascular instead of an extracellular contrast agent allows a reduction of the degrees of freedom for modeling tissue residue curves and results in improved accuracy of blood flow estimates.

Assessing myocardial perfusion in coronary artery disease with magnetic resonance first-pass imaging

MRFP perfusion imaging can now be used clinically on most MR scanner systems (1.0 to 1.5 T). The current experimental data demonstrate that MRFP imaging allows the quantitative assessment of myocardial blood flow changes and accurate measurements of collateral flow, including changes in the collateral dependent zones. Certain protocols, however, as outlined here have to be followed to obtain all the possible diagnostic information. Based on the current data on MRFP imaging, it is realistic to anticipate that MRFP imaging in combination with cine or tagging MR imaging will provide clinicians with better methods to distinguish stunned and hibernating, from nonviable myocardium and obtain better outcome data. Dedicated MR scanners are now being designed to meet the needs for MR imaging of patients with coronary artery disease. These scanners, small in size and with better patient access, make placement near the coronary care unit or catheterization laboratory feasible. This is a major step toward enhancing the utility of this new technique by providing the necessary infrastructure for scanning large numbers of patients. The main obstacle to wider use of these new diagnostic tools to assess perfusion is the lack of a large clinical database because there have not yet been major multicenter trials. With the development of novel intravascular contrast agents, however, larger trials are planned that should provide the clinical data mandatory for full integration of MRFP imaging into clinical practice. In particular, the development of dedicated and user-friendly perfusion analysis software will create the means to evaluate MR perfusion data accurately in large patient populations. These studies need to be conducted in a collaborative fashion by cardiologists, heart surgeons, and radiologists to be fully accepted by health care providers in an increasingly cost-averse and competitive health care environment.

MR first pass imaging: quantitative assessment of transmural perfusion and collateral flow

Recent advances with fast switching gradient coils, and the optimization of magnetic resonance techniques for multislice imaging have made it possible to apply models of contrast agent transit for the quantification of myocardial perfusion, and determination of the transmural distribution of blood flow. This article summarizes some of these recent developments and presents examples of quantitative, multi-slice myocardial perfusion imaging studies in patients and animal models. Multi-slice, true first pass imaging, with high temporal resolution, and T1-weighted, arrhythmia insensitive contrast enhancement is used for the quantification of perfusion changes accompanying mild to severe ischemia. The first pass imaging technique and the modeling approach are sufficiently robust for fitting of tissue residue curves corresponding to a wide, physiologically realistic range of myocardial blood flows. In animals this was validated by comparison to blood flow measurements with radiolabeled microspheres as gold standard. It is demonstrated that with the proposed modeling approach one can determine the myocardial perfusion reserve from two consecutive MR first pass measurements under resting and hyperemic conditions. In patients with microvascular dysfunction the MR studies show for the first time that the myocardial perfusion reserve correlates with Doppler flow measurements (linear regression with slope of 1.02 +/- 0.09; r = 0.80). Since perfusion limitations usually begin in the subendocardium as coronary flow is gradually reduced, first pass imaging with the prerequisitie spatial and temporal resolution allows early detection of a mild coronary stenosis.