Influence of the k-space trajectory on the dynamic T1-weighted signal in quantitative first-pass cardiac perfusion MRI at 3T

Accelerated Wideband Myocardial Perfusion Pulse Sequence with Compressed Sensing Reconstruction for Myocardial Blood Flow Quantification in Patients with a Cardiac Implantable Electronic Device

Purpose

To develop an accelerated wideband cardiac perfusion pulse sequence and test whether it can produce diagnostically acceptable image quality and whether it can be used to reliably quantify myocardial blood flow (MBF) in patients with a cardiac implantable electronic device (CIED).

Materials and Methods

A fivefold-accelerated wideband perfusion pulse sequence was developed using compressed sensing to sample one arterial input function plane and three myocardial perfusion (MP) planes per heartbeat in patients with a CIED with heart rates as high as 102 beats per minute. Resting perfusion scans were performed in 10 patients with a CIED and in 10 patients with no device as a control group. Two clinical readers compared the resulting images and retrospective images of the 10 patients with a CIED, which were obtained by using a previously described twofold-accelerated wideband perfusion pulse sequence with temporal generalized autocalibrating partially parallel acquisition. Summed visual score (SVS) was defined as the sum of conspicuity, artifact, and noise scores individually ranging from 1 (worst) to 5 (best). Resting MBF in the remote zones was quantified using Fermi deconvolution.

Results

Median SVS was significantly different (P < .05) between the prospective and retrospective CIED groups (13 vs nine) and between the nondevice group and the retrospective CIED group (13.5 vs nine); all median SVSs were nine or greater (clinically acceptable cut point). The median resting MBF in remote zones was not significantly different (P = .27) between patients with a CIED (1.1 mL/min/g; median left ventricular ejection fraction [LVEF], 52.5%) and patients with no device (1.3 mL/min/g; median LVEF, 64.0%). Mean MBF values were consistent with those (mean resting MBF range, 1.0-1.2 mL/min/g) reported by two prior state-of-the-art cardiac perfusion MRI studies.

Conclusion

The proposed scan yielded diagnostically acceptable image quality and enabled reliable quantification of MBF with three MP planes per heartbeat in patients with a CIED with heart rates as high as 102 beats per minute. Supplemental material is available for this article. © RSNA, 2020.

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.

Comparison of centric and reverse-centric trajectories for highly accelerated three-dimensional saturation recovery cardiac perfusion imaging

PURPOSE

Highly undersampled three-dimensional (3D) saturation-recovery sequences are affected by k-space trajectory since the magnetization does not reach steady state during the acquisition and the slab excitation profile yields different flip angles in different slices. This study compares centric and reverse-centric 3D cardiac perfusion imaging.

METHODS

An undersampled (98 phase encodes) 3D ECG-gated saturation-recovery sequence that alternates centric and reverse-centric acquisitions each time frame was used to image phantoms and in vivo subjects. Flip angle variation across the slices was measured, and contrast with each trajectory was analyzed via Bloch simulation.

RESULTS

Significant variations in flip angle were observed across slices, leading to larger signal variation across slices for the centric acquisition. In simulation, severe transient artifacts were observed when using the centric trajectory with higher flip angles, placing practical limits on the maximum flip angle used. The reverse-centric trajectory provided less contrast, but was more robust to flip angle variations.

CONCLUSION

Both of the k-space trajectories can provide reasonable image quality. The centric trajectory can have higher CNR, but is more sensitive to flip angle variation. The reverse-centric trajectory is more robust to flip angle variation.

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.

Myocardial perfusion MRI with an undersampled 3D stack-of-stars sequence

PURPOSE

To determine the feasibility of three-dimensional (3D) hybrid radial (stack-of-stars) MRI with spatiotemporal total variation (TV) constrained reconstruction for dynamic contrast enhanced myocardial perfusion imaging.

METHODS

An ECG-triggered saturation recovery turboFLASH sequence with undersampled stack-of-stars sampling with spatiotemporal TV constrained reconstruction was developed for dynamic contrast enhanced myocardial perfusion imaging. Simulations were performed to study the dependence of the approach to steady state on flip angle and saturation recovery time for this stack-of-stars acquisition. Phantom studies were used to show the effect of the flip angle selection and imperfect spoiling on image qualities. Studies were done in three humans to test the feasibility of the approach for myocardial perfusion imaging.

RESULTS

The simulation and phantom studies showed that imperfect spoiling and magnetization changes during the readout were a function of flip angle and nonoptimized selection of flip angle could degrade the images. Low flip angle acquisitions in the human subjects result in images with good quality similar to multislice radial 2D images.

CONCLUSIONS

3D stack-of-stars sampling with spatiotemporal TV constrained reconstruction provides a promising alternative for myocardial perfusion imaging.

Myocardial perfusion acquisition without magnetization preparation or gating

Current myocardial perfusion MRI acquisitions are performed with a saturation recovery sequence, in large part to minimize sensitivity to arrhythmia. A new approach is proposed here where the images are acquired ungated at steady state without use of a saturation pulse. The data are acquired continuously and reach steady state after the first few images. A confluence of advances has made this new paradigm of an ungated steady-state acquisition possible-very rapid undersampled readouts with new reconstruction technologies permit enough measurements that continuous acquisition becomes a feasible approach. Gating can be applied retrospectively from a logged electrocardiogram (ECG) or with self-gating methods. In this work, simulations and measurements in a concentration phantom are used to demonstrate that similar contrast and signal can be obtained with the standard saturation recovery and the proposed spoiled gradient echo (SPGR) acquisition. Specifically, for a flip angle of 14° and a saturation recovery time of 80 ms, similar signals are acquired over a range of T(1) s that reflect realistic myocardial tissue concentrations. Preliminary results in one subject are presented to show the potential of this new approach. The method may allow for cine cardiac perfusion and more signal-to-noise ratio-efficient acquisitions.

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

PURPOSE

To develop, and validate in vivo, a robust quantitative first-pass perfusion cardiovascular MR (CMR) method with accurate arterial input function (AIF) and myocardial wall enhancement.

MATERIALS AND METHODS

A saturation-recovery (SR) pulse sequence was modified to sequentially acquire multiple slices after a single nonselective saturation pulse at 3 Tesla. In each heartbeat, an AIF image is acquired in the aortic root with a short time delay (TD) (50 ms), followed by the acquisition of myocardial images with longer TD values (∼150-400 ms). Longitudinal relaxation rates (R(1) = 1/T(1)) were calculated using an ideal saturation recovery equation based on the Bloch equation, and corresponding gadolinium contrast concentrations were calculated assuming fast water exchange condition. The proposed method was validated against a reference multi-point SR method by comparing their respective R(1) measurements in the blood and left ventricular myocardium, before and at multiple time-points following contrast injections, in 7 volunteers.

RESULTS

R(1) measurements with the proposed method and reference multi-point method were strongly correlated (r > 0.88, P < 10(-5)) and in good agreement (mean difference ±1.96 standard deviation 0.131 ± 0.317/0.018 ± 0.140 s(-1) for blood/myocardium, respectively).

CONCLUSION

The proposed quantitative first-pass perfusion CMR method measured accurate R(1) values for quantification of AIF and myocardial wall contrast agent concentrations in 3 cardiac short-axis slices, in a total acquisition time of 523 ms per heartbeat.

Combination of compressed sensing and parallel imaging for highly accelerated first-pass cardiac perfusion MRI

First-pass cardiac perfusion MRI is a natural candidate for compressed sensing acceleration since its representation in the combined temporal Fourier and spatial domain is sparse and the required incoherence can be effectively accomplished by k-t random undersampling. However, the required number of samples in practice (three to five times the number of sparse coefficients) limits the acceleration for compressed sensing alone. Parallel imaging may also be used to accelerate cardiac perfusion MRI, with acceleration factors ultimately limited by noise amplification. In this work, compressed sensing and parallel imaging are combined by merging the k-t SPARSE technique with sensitivity encoding (SENSE) reconstruction to substantially increase the acceleration rate for perfusion imaging. We also present a new theoretical framework for understanding the combination of k-t SPARSE with SENSE based on distributed compressed sensing theory. This framework, which identifies parallel imaging as a distributed multisensor implementation of compressed sensing, enables an estimate of feasible acceleration for the combined approach. We demonstrate feasibility of 8-fold acceleration in vivo with whole-heart coverage and high spatial and temporal resolution using standard coil arrays. The method is relatively insensitive to respiratory motion artifacts and presents similar temporal fidelity and image quality when compared to Generalized autocalibrating partially parallel acquisitions (GRAPPA) with 2-fold acceleration.

In vivo quantification of contrast agent concentration using the induced magnetic field for time-resolved arterial input function measurement with MRI

For pharmacokinetic modeling of tissue physiology, there is great interest in measuring the arterial input function (AIF) from dynamic contrast-enhanced (DCE) magnetic resonance imaging (MRI) using paramagnetic contrast agents. Due to relaxation effects, the measured signal is a nonlinear function of the injected contrast agent concentration and depends on sequence parameters, system calibration, and time-of-flight effects, making it difficult to accurately measure the AIF during the first pass. Paramagnetic contrast agents also affect susceptibility and modify the magnetic field in proportion to their concentration. This information is contained in the MR signal phase which is discarded in a typical image reconstruction. However, quantifying AIF through contrast agent susceptibility induced phase changes is made difficult by the fact that the induced magnetic field is nonlocal and depends upon the contrast agent spatial distribution and thus on organ and vessel shapes. In this article, the contrast agent susceptibility was quantified through inversion of magnetic field shifts using a piece-wise constant model. Its feasibility is demonstrated by a determination of the AIF from the susceptibility-induced field changes of an intravenous bolus. After in vitro validation, a time-resolved two-dimensional (2D) gradient echo scan, triggered to diastole, was performed in vivo on the aortic arch during a bolus injection of 0.1 mmol/kg Gd-DTPA. An approximate geometrical model of the aortic arch constructed from the magnitude images was used to calculate the spatial variation of the field associated with the bolus. In 14 subjects, Gd concentration curves were measured dynamically (one measurement per heart beat) and indirectly validated by independent 2D cine phase contrast flow rate measurements. Flow rate measurements using indicator conservation with this novel quantitative susceptibility imaging technique were found to be in good agreement with those obtained from the cine phase contrast measurements in all subjects. Contrary to techniques that rely on intensity, the accuracy of this signal phase based method is insensitive to factors influencing signal intensity such as flip angle, coil sensitivity, relaxation changes, and time-of-flight effects extending the range of pulse sequences and contrast doses for which quantitative DCE-MRI can be applied.