Quantitative contrast-enhanced first-pass cardiac perfusion MRI at 3 tesla with accurate arterial input function and myocardial wall enhancement

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.

Myocardial perfusion quantification using simultaneously acquired 13 NH3 -ammonia PET and dynamic contrast-enhanced MRI in patients at rest and stress

PURPOSE

Systematic differences with respect to myocardial perfusion quantification exist between DCE-MRI and PET. Using the potential of integrated PET/MRI, this study was conceived to compare perfusion quantification on the basis of simultaneously acquired ¹³ NH₃ -ammonia PET and DCE-MRI data in patients at rest and stress.

METHODS

Twenty-nine patients were examined on a 3T PET/MRI scanner. DCE-MRI was implemented in dual-sequence design and additional T₁ mapping for signal normalization. Four different deconvolution methods including a modified version of the Fermi technique were compared against ¹³ NH₃ -ammonia results.

RESULTS

Cohort-average flow comparison yielded higher resting flows for DCE-MRI than for PET and, therefore, significantly lower DCE-MRI perfusion ratios under the common assumption of equal arterial and tissue hematocrit. Absolute flow values were strongly correlated in both slice-average (R²  = 0.82) and regional (R²  = 0.7) evaluations. Different DCE-MRI deconvolution methods yielded similar flow result with exception of an unconstrained Fermi method exhibiting outliers at high flows when compared with PET.

CONCLUSION

Thresholds for Ischemia classification may not be directly tradable between PET and MRI flow values. Differences in perfusion ratios between PET and DCE-MRI may be lifted by using stress/rest-specific hematocrit conversion. Proper physiological constraints are advised in model-constrained deconvolution.

Analysis of spatiotemporal fidelity in quantitative 3D first-pass perfusion cardiovascular magnetic resonance

BACKGROUND

Whole-heart first-pass perfusion cardiovascular magnetic resonance (CMR) relies on highly accelerated image acquisition. The influence of undersampling on myocardial blood flow (MBF) quantification has not been systematically investigated yet. In the present work, the effect of spatiotemporal scan acceleration on image reconstruction accuracy and MBF error was studied using a numerical phantom and validated in-vivo.

METHODS

Up to 10-fold scan acceleration using k-t PCA and k-t SPARSE-SENSE was simulated using the MRXCAT CMR numerical phantom framework. Image reconstruction results were compared to ground truth data in the k-f domain by means of modulation transfer function (MTF) analysis. In the x-t domain, errors pertaining to specific features of signal intensity-time curves and MBF values derived using Fermi model deconvolution were analysed. In-vivo first-pass CMR data were acquired in ten healthy volunteers using a dual-sequence approach assessing the arterial input function (AIF) and myocardial enhancement. 10x accelerated 3D k-t PCA and k-t SPARSE-SENSE were compared and related to non-accelerated 2D reference images.

RESULTS

MTF analysis revealed good recovery of data upon k-t PCA reconstruction at 10x undersampling with some attenuation of higher temporal frequencies. For 10x k-t SPARSE-SENSE the MTF was found to decrease to zero at high spatial frequencies for all temporal frequencies indicating a loss in spatial resolution. Signal intensity-time curve errors were most prominent in AIFs from 10x k-t PCA, thereby emphasizing the need for separate AIF acquisition using a dual-sequence approach. These findings were confirmed by MBF estimation based on AIFs from fully sampled and undersampled simulations. Average in-vivo MBF estimates were in good agreement between both accelerated and the fully sampled methods. Intra-volunteer MBF variation for fully sampled 2D scans was lower compared to 10x k-t PCA and k-t SPARSE-SENSE data.

CONCLUSION

Quantification of highly undersampled 3D first-pass perfusion CMR yields accurate MBF estimates provided the AIF is obtained using fully sampled or moderately undersampled scans as part of a dual-sequence approach. However, relative to fully sampled 2D perfusion imaging, intra-volunteer variation is increased using 3D approaches prompting for further developments.

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

Purpose

Perfusion analysis from first-pass contrast enhancement kinetics requires modeling tissue contrast exchange. This study presents a new approach for numerical implementation of the tissue homogeneity model, incorporating flexible distance steps along the capillary (NTH_f).

Methods

The proposed NTH_f model considers contrast exchange in fluid packets flowing along the capillary, incorporating flexible distance steps, thus allowing more efficient and stable calculations of the transit of tracer through the tissue. We prospectively studied 8 patients (62 ± 13 years old) with suspected CAD, who underwent first-pass perfusion CMR imaging at rest and stress prior to angiography. Myocardial blood flow (MBF) and myocardial perfusion reserve index (MPRI) were estimated using both the NTH_f and the conventional adiabatic approximation of the TH models. Coronary artery lesions detected at angiography were clinically assigned to one of three categories of stenosis severity (‘insignificant’, ‘mild to moderate’ and ‘severe’) and related to corresponding myocardial territories.

Results

The mean MBF (ml/g/min) at rest/stress and MPRI were 0.80 ± 0.33/1.25 ± 0.45 and 1.68 ± 0.54 in the insignificant regions, 0.74 ± 0.21/1.09 ± 0.28 and 1.54 ± 0.46 in the mild to moderate regions, and 0.79 ± 0.28/0.63 ± 0.34 and 0.85 ± 0.48 in the severe regions, respectively. The correlation coefficients of MBFs at rest/stress and MPRI between the NTH_f and AATH models were r = 0.97/0.93 and r = 0.91, respectively.

Conclusions

The proposed NTH_f model allows efficient quantitative analysis of the transit of tracer through tissue, particularly at higher flow. Results of initial application to MRI of myocardial perfusion in CAD are encouraging.

Measurement of extracellular volume and transit time heterogeneity using contrast-enhanced myocardial perfusion MRI in patients after acute myocardial infarction

PURPOSE

To assess the ability of dynamic contrast-enhanced myocardial perfusion MRI to measure extracellular volume (ECV) and to investigate the possibility of estimating capillary transit time heterogeneity (CTH) in patients after myocardial infarction and successful revascularization.

METHODS

Twenty-four perfusion data sets were acquired on a 3 Tesla positron emission tomography (PET)/MRI scanner. Three perfusion models of different complexity were implemented in a hierarchical fashion with an Akaike information criterion being used to determine the number of fit parameters supported by the data. Results were compared sector-wise to ECV from an equilibrium T₁ mapping method (modified look-locker inversion recovery (MOLLI)).

RESULTS

ECV derived from the perfusion analysis correlated well with equilibrium measurements (R² = 0.76). Estimation of CTH was supported in 16% of sectors (mostly remote). Inclusion of a nonzero CTH parameter usually led to lower estimates of first-pass extraction and slightly higher estimates of blood volume and flow. Estimation of the capillary permeability-surface area product was feasible in 81% of sectors.

CONCLUSION

Transit time heterogeneity has a measurable effect on the kinetic analysis of myocardial perfusion MRI data, and Gd-DTPA extravasation in the myocardium is usually not flow-limited in infarct-related pathology. Measurement of myocardial ECV using perfusion imaging could provide a scan-time efficient alternative to methods based on T₁ mapping. Magn Reson Med 77:2320-2330, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Quantitative three-dimensional myocardial perfusion cardiovascular magnetic resonance with accurate two-dimensional arterial input function assessment

BACKGROUND

Quantification of myocardial perfusion from first-pass cardiovascular magnetic resonance (CMR) images at high contrast agent (CA) dose requires separate acquisition of blood pool and myocardial tissue enhancement. In this study, a dual-sequence approach interleaving 2D imaging of the arterial input function with high-resolution 3D imaging for myocardial perfusion assessment is presented and validated for low and high CA dose.

METHODS

A dual-sequence approach interleaving 2D imaging of the aortic root and 3D imaging of the whole left ventricle using highly accelerated k-t PCA was implemented. Rest perfusion imaging was performed in ten healthy volunteers after administration of a Gadolinium-based CA at low (0.025 mmol/kg b.w.) and high dose (0.1 mmol/kg b.w.). Arterial input functions extracted from the 2D and 3D images were analysed for both doses. Myocardial contrast-to-noise ratios (CNR) were compared across volunteers and doses. Variations of myocardial perfusion estimates between volunteers and across myocardial territories were studied.

RESULTS

High CA dose imaging resulted in strong non-linearity of the arterial input function in the 3D images at peak CA concentration, which was avoided when the input function was derived from the 2D images. Myocardial CNR was significantly increased at high dose compared to low dose, with a 2.6-fold mean CNR gain. Most robust myocardial blood flow estimation was achieved using the arterial input function extracted from the 2D image at high CA dose. In this case, myocardial blood flow estimates varied by 24% between volunteers and by 20% between myocardial territories when analysed on a per-volunteer basis.

CONCLUSION

Interleaving 2D imaging for arterial input function assessment enables robust quantitative 3D myocardial perfusion imaging at high CA dose.

Wideband arrhythmia-Insensitive-rapid (AIR) pulse sequence for cardiac T1 mapping without image artifacts induced by an implantable-cardioverter-defibrillator

PURPOSE

To develop and evaluate a wideband arrhythmia-insensitive-rapid (AIR) pulse sequence for cardiac T1 mapping without image artifacts induced by implantable-cardioverter-defibrillator (ICD).

METHODS

We developed a wideband AIR pulse sequence by incorporating a saturation pulse with wide frequency bandwidth (8.9 kHz) to achieve uniform T1 weighting in the heart with ICD. We tested the performance of original and "wideband" AIR cardiac T1 mapping pulse sequences in phantom and human experiments at 1.5 Tesla.

RESULTS

In five phantoms representing native myocardium and blood and postcontrast blood/tissue T1 values, compared with the control T1 values measured with an inversion-recovery pulse sequence without ICD, T1 values measured with original AIR with ICD were considerably lower (absolute percent error > 29%), whereas T1 values measured with wideband AIR with ICD were similar (absolute percent error < 5%). Similarly, in 11 human subjects, compared with the control T1 values measured with original AIR without ICD, T1 measured with original AIR with ICD was significantly lower (absolute percent error > 10.1%), whereas T1 measured with wideband AIR with ICD was similar (absolute percent error < 2.0%).

CONCLUSION

This study demonstrates the feasibility of a wideband pulse sequence for cardiac T1 mapping without significant image artifacts induced by ICD.

Three-dimensional dynamic contrast enhanced imaging of the carotid artery with direct arterial input function measurement

PURPOSE

Kinetic analysis using dynamic contrast enhanced MRI to assess neovascularization of carotid plaque requires images with high spatial and temporal resolution. This work demonstrates a new three-dimensional (3D) dynamic contrast enhanced imaging sequence, which directly measures the arterial input function with high temporal resolution yet maintains the high spatial resolution required to identify areas of increased adventitial neovascularity.

THEORY AND METHODS

The sequence consists of multiple rapid acquisitions of a saturation prepared dynamic 3D gradient recalled echo (GRE) sequence temporally interleaved with multiple acquisitions of a 2D slice. The saturation recovery time was adjusted to maintain signal linearity with the very different contrast agent concentrations in the 2D slice and 3D volume. The K(trans) maps were obtained from the 3D dynamic contrast measurements while the 2D slice was used to obtain the arterial input function. Calibration and dynamic studies are presented.

RESULTS

For contrast agent concentrations up to 5 mM, a saturation recovery time for the 2D slice of 20 ms resulted in less than a 10% deviation from the desired linear response of signal intensity with contrast agent concentration. The corresponding saturation recovery time of 83 ms for the 3D volume maintained less than a 10% deviation from the linear response up to contrast agent concentrations of 2 mM while a contrast agent concentration of 5 mM had almost a 30% deviation. There was a significant improvement in signal attenuation (9 ± 3% versus 23 ± 5% at 40 cm/s) when flow compensation was added to the slice select gradients. For patient studies, volume transfer and plasma fraction maps were calculated with data from the proposed sequence.

CONCLUSION

This work demonstrated a novel sequence for 3D dynamic contrast enhanced imaging with a simultaneously acquired 2D slice that directly measures the arterial input function with high temporal resolution. Acquisition parameters can be adjusted to accommodate the full range of contrast agent concentration values to be encountered and the kinetic parameters obtained were consistent with expected values.

Is it possible to do without the study of myocardial perfusion in 2013?

The analysis of myocardial perfusion is a key step in the cardiac MRI examination. In routine work, this exploration carried out at rest is based on the qualitative first pass study of gadolinium with an ECG-triggered saturation recovery bFFE sequence. In view of recent knowledge, the analysis of the myocardial perfusion under vasodilator stress may be carried out by scintigraphy or MRI, the latter benefiting from the absence of exposure to ionizing rays and a lower cost. Besides coronary disease, the perfusion sequence provides a rich semiology to compare with the clinics and the data from other sequences. Arterial Spin Labeling (ASL) is an alternative technique used in the animal to quantify myocardial perfusion.

SHILO, a novel dual imaging approach for simultaneous HI-/LOw temporal (Low-/Hi-spatial) resolution imaging for vascular dynamic contrast enhanced cardiovascular magnetic resonance: numerical simulations and feasibility in the carotid arteries

BACKGROUND

Dynamic contrast enhanced (DCE) cardiovascular magnetic resonance (CMR) is increasingly used to quantify microvessels and permeability in atherosclerosis. Accurate quantification depends on reliable sampling of both vessel wall (VW) uptake and contrast agent dynamic in the blood plasma (the so called arterial input function, AIF). This poses specific challenges in terms of spatial/temporal resolution and matched dynamic MR signal range, which are suboptimal in current vascular DCE-CMR protocols. In this study we describe a novel dual-imaging approach, which allows acquiring simultaneously AIF and VW images using different spatial/temporal resolution and optimizes imaging parameters for the two compartments. We refer to this new acquisition as SHILO, Simultaneous HI-/LOw-temporal (low-/hi-spatial) resolution DCE-imaging.

METHODS

In SHILO, the acquisition of low spatial resolution single-shot AIF images is interleaved with segments of higher spatial resolution images of the VW. This allows sampling the AIF and VW with different spatial/temporal resolution and acquisition parameters, at independent spatial locations. We show the adequacy of this temporal sampling scheme by using numerical simulations. Following, we validate the MR signal of SHILO against a standard 2D spoiled gradient recalled echo (SPGR) acquisition with in vitro and in vivo experiments. Finally, we show feasibility of using SHILO imaging in subjects with carotid atherosclerosis.

RESULTS

Our simulations confirmed the superiority of the SHILO temporal sampling scheme over conventional strategies that sample AIF and tissue curves at the same time resolution. Both the median relative errors and standard deviation of absolute parameter values were lower for the SHILO than for conventional sampling schemes. We showed equivalency of the SHILO signal and conventional 2D SPGR imaging, using both in vitro phantom experiments (R2 =0.99) and in vivo acquisitions (R2 =0.95). Finally, we showed feasibility of using the newly developed SHILO sequence to acquire DCE-CMR data in subjects with carotid atherosclerosis to calculate plaque perfusion indices.

CONCLUSIONS

We successfully demonstrate the feasibility of using the newly developed SHILO dual-imaging technique for simultaneous AIF and VW imaging in DCE-CMR of atherosclerosis. Our initial results are promising and warrant further investigation of this technique in wider studies measuring kinetic parameters of plaque neovascularization with validation against gold standard techniques.

Absolute quantification of myocardial blood flow with constrained estimation of the arterial input function

PURPOSE

To evaluate the performance of the constrained alternating minimization with model (CAMM) method for estimating the input function from the myocardial tissue curves.

MATERIALS AND METHODS

Myocardial perfusion imaging was performed on seven canine models of coronary artery disease in 15 imaging sessions. In each session, stress was induced with intravenous infusion of adenosine and a variable occluder created coronary artery stenosis. A dual bolus protocol was used for each acquisition, and input functions were then estimated using the CAMM method with data acquired from the high dose scan following each imaging session. For each acquisition, myocardial blood flow was measured by injected microspheres.

RESULTS

The dual bolus and CAMM-derived flows were not significantly different (P = 0.18), and the correlation between the two methods was high (r = 0.97). The correlation between the dual bolus and CAMM methods and microsphere measurements was lower than that for the two MR methods (r = 0.53; r = 0.43, respectively).

CONCLUSION

The CAMM method presented here shows promise in estimating myocardial blood flow in patients with coronary artery disease at stress with a single injection and without any specialized acquisitions. Further work is needed to validate the approach in a clinical setting.

Arrhythmia insensitive rapid cardiac T1 mapping pulse sequence

PURPOSE

To develop an arrhythmia-insensitive rapid (AIR) cardiac T1 mapping pulse sequence for quantification of diffuse fibrosis.

METHODS

An arrhythmia-insensitive cardiac T1 mapping pulse sequence was developed based on saturation recovery T1 weighting, which is inherently insensitive to heart rate and rhythm, and two single-shot balanced steady-state free precession image acquisitions with centric k-space ordering, where T1 calculation is inherently insensitive to T2 effects. Its performance against conventional cardiac T1 mapping based on inversion recovery (i.e., MOLLI) is compared. Phantom experiments (T1 ranging from 535 to 2123 ms) were performed with heart rate and rhythm simulated at 60 and 120 beats per minute (bpm) and arrhythmia using an external triggering device. Ten human subjects and 17 large animals were scanned precontrast and 5, 10, and 15 min after contrast agent administration.

RESULTS

Compared with the reference T1 mapping, AIR yielded lower normalized root-mean-square error than MOLLI (8% vs. 3%, respectively, at 60 bpm, 28% vs. 3%, respectively, at 120 bpm, and 22% vs. 3%, respectively, at arrhythmia). In vivo studies showed that T1 measurements made by MOLLI and AIR were strongly correlated (r = 0.99) but in poor agreement (mean difference = 161.8 ms, upper and lower 95% limits of agreements = 347.5 ms and -24.0 ms).

CONCLUSION

Our AIR pulse sequence may be clinically useful for assessment of diffuse myocardial fibrosis in patients.