Quantitative Perfusion Analysis of First-Pass Contrast Enhancement Kinetics: Application to MRI of Myocardial Perfusion in Coronary Artery Disease

Sub-segmental quantification of single (stress)-pass perfusion CMR improves the diagnostic accuracy for detection of obstructive coronary artery disease

BACKGROUND

Myocardial perfusion with cardiovascular magnetic resonance (CMR) imaging is an established diagnostic test for evaluation of myocardial ischaemia. For quantification purposes, the 16 segment American Heart Association (AHA) model poses limitations in terms of extracting relevant information on the extent/severity of ischaemia as perfusion deficits will not always fall within an individual segment, which reduces its diagnostic value, and makes an accurate assessment of outcome data or a result comparison across various studies difficult. We hypothesised that division of the myocardial segments into epi- and endocardial layers and a further circumferential subdivision, resulting in a total of 96 segments, would improve the accuracy of detecting myocardial hypoperfusion. Higher (sub-)subsegmental recording of perfusion abnormalities, which are defined relatively to the normal reference using the subsegment with the highest value, may improve the spatial encoding of myocardial blood flow, based on a single stress perfusion acquisition.

OBJECTIVE

A proof of concept comparison study of subsegmentation approaches based on transmural segments (16 AHA and 48 segments) vs. subdivision into epi- and endocardial (32) subsegments vs. further circumferential subdivision into 96 (sub-)subsegments for diagnostic accuracy against invasively defined obstructive coronary artery disease (CAD).

METHODS

Thirty patients with obstructive CAD and 20 healthy controls underwent perfusion stress CMR imaging at 3 T during maximal adenosine vasodilation and a dual bolus injection of 0.1 mmol/kg gadobutrol. Using Fermi deconvolution for blood flow estimation, (sub-)subsegmental values were expressed relative to the (sub-)subsegment with the highest flow. In addition, endo-/epicardial flow ratios were calculated based on 32 and 96 (sub-)subsegments. A receiver operating characteristics (ROC) curve analysis was performed to compare the diagnostic performance of discrimination between patients with CAD and healthy controls. Observer reproducibility was assessed using Bland-Altman approaches.

RESULTS

Subdivision into more and smaller segments revealed greater accuracy for #32, #48 and # 96 compared to the standard #16 approach (area under the curve (AUC): 0.937, 0.973 and 0.993 vs 0.820, p < 0.05). The #96-based endo-/epicardial ratio was superior to the #32 endo-/epicardial ratio (AUC 0.979, vs. 0.932, p < 0.05). Measurements for the #16 model showed marginally better reproducibility compared to #32, #48 and #96 (mean difference ± standard deviation: 2.0 ± 3.6 vs. 2.3 ± 4.0 vs 2.5 ± 4.4 vs. 4.1 ± 5.6).

CONCLUSIONS

Subsegmentation of the myocardium improves diagnostic accuracy and facilitates an objective cut-off-based description of hypoperfusion, and facilitates an objective description of hypoperfusion, including the extent and severity of myocardial ischaemia. Quantification based on a single (stress-only) pass reduces the overall amount of gadolinium contrast agent required and the length of the overall diagnostic study.

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.

State-of-the-Art Quantitative Assessment of Myocardial Ischemia by Stress Perfusion Cardiac Magnetic Resonance

Ischemic heart disease remains the foremost determinant of death and disability across the world. Quantification of the ischemia burden is currently the preferred approach to predict event risk and to trigger adequate treatment. Cardiac magnetic resonance (CMR) can be a prime protagonist in this scenario due to its synergistic features. It allows assessment of wall motility, myocardial perfusion, and tissue scar by means of late gadolinium enhancement imaging. We discuss the clinical and preclinical aspects of gadolinium-based, perfusion CMR imaging, including the relevance of high spatial resolution and 3-dimensional whole-heart coverage, among important features of this auspicious method.

Automated Quantitative Stress Perfusion in a Clinical Routine

Cardiovascular magnetic resonance (CMR) perfusion imaging is a robust noninvasive technique to evaluate ischemia in patients with coronary artery disease (CAD). Although qualitative and semiquantitative methods have shown that CMR has high accuracy in diagnosing flow-obstructing lesions in CAD, quantitative ischemic burden is an important variable used in clinical practice for treatment decisions. Quantitative CMR perfusion techniques have evolved significantly, with accuracy comparable with both PET and microsphere evaluation. Routine clinical use of these quantitative techniques has been facilitated by the introduction of automated methods that accelerate the work flow and rapidly generate pixel-based myocardial blood flow maps.

Hierarchical segmentation using equivalence test (HiSET): Application to DCE image sequences

Dynamical contrast enhanced (DCE) imaging allows non invasive access to tissue micro-vascularization. It appears as a promising tool to build imaging biomarkers for diagnostic, prognosis or anti-angiogenesis treatment monitoring of cancer. However, quantitative analysis of DCE image sequences suffers from low signal to noise ratio (SNR). SNR may be improved by averaging functional information in a large region of interest when it is functionally homogeneous. We propose a novel method for automatic segmentation of DCE image sequences into functionally homogeneous regions, called DCE-HiSET. Using an observation model which depends on one parameter a and is justified a posteriori, DCE-HiSET is a hierarchical clustering algorithm. It uses the p-value of a multiple equivalence test as dissimilarity measure and consists of two steps. The first exploits the spatial neighborhood structure to reduce complexity and takes advantage of the regularity of anatomical features, while the second recovers (spatially) disconnected homogeneous structures at a larger (global) scale. Given a minimal expected homogeneity discrepancy for the multiple equivalence test, both steps stop automatically by controlling the Type I error. This provides an adaptive choice for the number of clusters. Assuming that the DCE image sequence is functionally piecewise constant with signals on each piece sufficiently separated, we prove that DCE-HiSET will retrieve the exact partition with high probability as soon as the number of images in the sequence is large enough. The minimal expected homogeneity discrepancy appears as the tuning parameter controlling the size of the segmentation. DCE-HiSET has been implemented in C++ for 2D and 3D image sequences with competitive speed.

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.

Myocardial perfusion cardiovascular magnetic resonance: optimized dual sequence and reconstruction for quantification

BACKGROUND

Quantification of myocardial blood flow requires knowledge of the amount of contrast agent in the myocardial tissue and the arterial input function (AIF) driving the delivery of this contrast agent. Accurate quantification is challenged by the lack of linearity between the measured signal and contrast agent concentration. This work characterizes sources of non-linearity and presents a systematic approach to accurate measurements of contrast agent concentration in both blood and myocardium.

METHODS

A dual sequence approach with separate pulse sequences for AIF and myocardial tissue allowed separate optimization of parameters for blood and myocardium. A systems approach to the overall design was taken to achieve linearity between signal and contrast agent concentration. Conversion of signal intensity values to contrast agent concentration was achieved through a combination of surface coil sensitivity correction, Bloch simulation based look-up table correction, and in the case of the AIF measurement, correction of T2* losses. Validation of signal correction was performed in phantoms, and values for peak AIF concentration and myocardial flow are provided for 29 normal subjects for rest and adenosine stress.

RESULTS

For phantoms, the measured fits were within 5% for both AIF and myocardium. In healthy volunteers the peak [Gd] was 3.5 ± 1.2 for stress and 4.4 ± 1.2 mmol/L for rest. The T2* in the left ventricle blood pool at peak AIF was approximately 10 ms. The peak-to-valley ratio was 5.6 for the raw signal intensities without correction, and was 8.3 for the look-up-table (LUT) corrected AIF which represents approximately 48% correction. Without T2* correction the myocardial blood flow estimates are overestimated by approximately 10%. The signal-to-noise ratio of the myocardial signal at peak enhancement (1.5 T) was 17.7 ± 6.6 at stress and the peak [Gd] was 0.49 ± 0.15 mmol/L. The estimated perfusion flow was 3.9 ± 0.38 and 1.03 ± 0.19 ml/min/g using the BTEX model and 3.4 ± 0.39 and 0.95 ± 0.16 using a Fermi model, for stress and rest, respectively.

CONCLUSIONS

A dual sequence for myocardial perfusion cardiovascular magnetic resonance and AIF measurement has been optimized for quantification of myocardial blood flow. A validation in phantoms was performed to confirm that the signal conversion to gadolinium concentration was linear. The proposed sequence was integrated with a fully automatic in-line solution for pixel-wise mapping of myocardial blood flow and evaluated in adenosine stress and rest studies on N = 29 normal healthy subjects. Reliable perfusion mapping was demonstrated and produced estimates with low variability.