Quantification of myocardial perfusion with self-gated cardiovascular magnetic resonance

Free-breathing motion-informed locally low-rank quantitative 3D myocardial perfusion imaging

Purpose

To propose respiratory motion‐informed locally low‐rank reconstruction (MI‐LLR) for robust free‐breathing single‐bolus quantitative 3D myocardial perfusion CMR imaging. Simulation and in‐vivo results are compared to locally low‐rank (LLR) and compressed sensing reconstructions (CS) for reference.

Methods

Data were acquired using a 3D Cartesian pseudo‐spiral in‐out k‐t undersampling scheme (R = 10) and reconstructed using MI‐LLR, which encompasses two stages. In the first stage, approximate displacement fields are derived from an initial LLR reconstruction to feed a motion‐compensated reference system to a second reconstruction stage, which reduces the rank of the inverse problem. For comparison, data were also reconstructed with LLR and frame‐by‐frame CS using wavelets as sparsifying transform (ℓ1‐wavelet).

Reconstruction accuracy relative to ground truth was assessed using synthetic data for realistic ranges of breathing motion, heart rates, and SNRs. In‐vivo experiments were conducted in healthy subjects at rest and during adenosine stress. Myocardial blood flow (MBF) maps were derived using a Fermi model.

Results

Improved uniformity of MBF maps with reduced local variations was achieved with MI‐LLR. For rest and stress, intra‐volunteer variation of absolute and relative MBF was lower in MI‐LLR (±0.17 mL/g/min [26%] and ±1.07 mL/g/min [33%]) versus LLR (±0.19 mL/g/min [28%] and ±1.22 mL/g/min [36%]) and versus ℓ1‐wavelet (±1.17 mL/g/min [113%] and ±6.87 mL/g/min [115%]). At rest, intra‐subject MBF variation was reduced significantly with MI‐LLR.

Conclusion

The combination of pseudo‐spiral Cartesian undersampling and dual‐stage MI‐LLR reconstruction improves free‐breathing quantitative 3D myocardial perfusion CMR imaging under rest and stress condition.

Click here for author‐reader discussions

Whole-heart, ungated, free-breathing, cardiac-phase-resolved myocardial perfusion MRI by using Continuous Radial Interleaved simultaneous Multi-slice acquisitions at sPoiled steady-state (CRIMP)

PURPOSE

To develop a whole-heart, free-breathing, non-electrocardiograph (ECG)-gated, cardiac-phase-resolved myocardial perfusion MRI framework (CRIMP; Continuous Radial Interleaved simultaneous Multi-slice acquisitions at sPoiled steady-state) and test its quantification feasibility.

METHODS

CRIMP used interleaved radial simultaneous multi-slice (SMS) slice groups to cover the whole heart in 9 or 12 short-axis slices. The sequence continuously acquired data without magnetization preparation, ECG gating or breath-holding, and captured multiple cardiac phases. Images were reconstructed by a motion-compensated patch-based locally low-rank reconstruction. Bloch simulations were performed to study the signal-to-noise ratio/contrast-to-noise ratio (SNR/CNR) for CRIMP and to study the steady-state signal under motion. Seven patients were scanned with CRIMP at stress and rest to develop the sequence. One human and two dogs were scanned at rest with a dual-bolus method to test the quantification feasibility of CRIMP. The dual-bolus scans were performed using both CRIMP and an ungated radial SMS saturation recovery (SMS-SR) sequence with injection dose = 0.075 mmol/kg to compare the sequences in terms of SNR, cardiac phase resolution and quantitative myocardial blood flow (MBF).

RESULTS

Perfusion images with multiple cardiac phases in all image slices with a temporal resolution of 72 ms/frame were obtained. Simulations and in-vivo acquisitions showed CRIMP kept the inner slices in steady-state regardless of motion. CRIMP outperformed SMS-SR in slice coverage (9 over 6), SNR (mean 20% improvement), and provided cardiac phase resolution. CRIMP and SMS-SR sequences provided comparable MBF values (rest systolic CRIMP = 0.58 ± 0.07, SMS-SR = 0.61 ± 0.16).

CONCLUSION

CRIMP allows for whole-heart, cardiac-phase-resolved myocardial perfusion images without ECG-gating or breath-holding. The sequence can provide MBF if an accurate arterial input function is obtained separately.

Respiratory motion compensation using data binning in dynamic contrast enhanced golden-angle radial MRI

GRASP (Golden-Angle Radial Sparse Parallel MRI) is a data acquisition and reconstruction technique that combines parallel imaging and golden-angle radial sampling. The continuously acquired free breathing Dynamic Contrast Enhanced (DCE) golden-angle radial MRI data of liver and abdomen has artifacts due to respiratory motion, resulting in low vessel-tissue contrast that makes GRASP reconstructed images less suitable for diagnosis. In this paper, DCE golden-angle radial MRI data of abdomen and liver perfusion is sorted into different motion states using the self-gating property of radial acquisition and then reconstructed using GRASP. Three methods of amplitude-based data binning namely uniform binning, adaptive binning and optimal binning are applied on the DCE golden-angle radial data to extract different motion states and a comparison is performed with the conventional GRASP reconstruction. Also, a comparison among the amplitude-based data binning techniques is performed and benefits of each of these binning techniques are discussed from a clinical perspective. The image quality assessment in terms of hepatic vessel clarity, liver edge sharpness, contrast enhancement clarity and streaking artifacts is performed by a certified radiologist. The results show that DCE golden-angle radial trajectories benefit from all the three types of amplitude-based data binning methods providing improved reconstruction results. The choice of binning technique depends upon the clinical application e.g. uniform and adaptive binning are helpful for a detailed analysis of lesion characteristic and contrast enhancement in different motion states while optimal binning can be used when clinical analysis requires a single image per contrast enhancement phase with no motion blurring artifacts.

Review of Journal of Cardiovascular Magnetic Resonance (JCMR) 2015-2016 and transition of the JCMR office to Boston

The Journal of Cardiovascular Magnetic Resonance (JCMR) is the official publication of the Society for Cardiovascular Magnetic Resonance (SCMR). In 2016, the JCMR published 93 manuscripts, including 80 research papers, 6 reviews, 5 technical notes, 1 protocol, and 1 case report. The number of manuscripts published was similar to 2015 though with a 12% increase in manuscript submissions to an all-time high of 369. This reflects a decrease in the overall acceptance rate to <25% (excluding solicited reviews). The quality of submissions to JCMR continues to be high. The 2016 JCMR Impact Factor (which is published in June 2016 by Thomson Reuters) was steady at 5.601 (vs. 5.71 for 2015; as published in June 2016), which is the second highest impact factor ever recorded for JCMR. The 2016 impact factor means that the JCMR papers that were published in 2014 and 2015 were on-average cited 5.71 times in 2016.In accordance with Open-Access publishing of Biomed Central, the JCMR articles are published on-line in the order that they are accepted with no collating of the articles into sections or special thematic issues. For this reason, over the years, the Editors have felt that it is useful to annually summarize the publications into broad areas of interest or themes, so that readers can view areas of interest in a single article in relation to each other and other recent JCMR articles. The papers are presented in broad themes with previously published JCMR papers to guide continuity of thought in the journal. In addition, I have elected to open this publication with information for the readership regarding the transition of the JCMR editorial office to the Beth Israel Deaconess Medical Center, Boston and the editorial process.Though there is an author publication charge (APC) associated with open-access to cover the publisher's expenses, this format provides a much wider distribution/availability of the author's work and greater manuscript citation. For SCMR members, there is a substantial discount in the APC. I hope that you will continue to send your high quality manuscripts to JCMR for consideration. Importantly, I also ask that you consider referencing recent JCMR publications in your submissions to the JCMR and elsewhere as these contribute to our impact factor. I also thank our dedicated Associate Editors, Guest Editors, and reviewers for their many efforts to ensure that the review process occurs in a timely and responsible manner and that the JCMR continues to be recognized as the leading publication in our field.

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.

Technical Note: Evaluation of pre-reconstruction interpolation methods for iterative reconstruction of radial k-space data

PURPOSE

To evaluate the use of three different pre-reconstruction interpolation methods to convert non-Cartesian k-space data to Cartesian samples such that iterative reconstructions can be performed more simply and more rapidly.

METHODS

Phantom as well as cardiac perfusion radial datasets were reconstructed by four different methods. Three of the methods used pre-reconstruction interpolation once followed by a fast Fourier transform (FFT) at each iteration. The methods were: bilinear interpolation of nearest-neighbor points (BINN), 3-point interpolation, and a multi-coil interpolator called GRAPPA Operator Gridding (GROG). The fourth method performed a full non-Uniform FFT (NUFFT) at each iteration. An iterative reconstruction with spatiotemporal total variation constraints was used with each method. Differences in the images were quantified and compared.

RESULTS

The GROG multicoil interpolation, the 3-point interpolation, and the NUFFT-at-each-iteration approaches produced high quality images compared to BINN, with the GROG-derived images having the fewest streaks among the three preinterpolation approaches. However, all reconstruction methods produced approximately equal results when applied to perfusion quantitation tasks. Pre-reconstruction interpolation gave approximately an 83% reduction in reconstruction time.

CONCLUSION

Image quality suffers little from using a pre-reconstruction interpolation approach compared to the more accurate NUFFT-based approach. GROG-based pre-reconstruction interpolation appears to offer the best compromise by using multicoil information to perform the interpolation to Cartesian sample points prior to image reconstruction. Speed gains depend on the implementation and relatively standard optimizations on a MATLAB platform result in preinterpolation speedups of ~ 6 compared to using NUFFT at every iteration, reducing the reconstruction time from around 42 min to 7 min.

Rapid rest/stress regadenoson ungated perfusion CMR for detection of coronary artery disease in patients with atrial fibrillation

Cardiovascular magnetic resonance (CMR) perfusion has been established as a useful imaging modality for the detection of coronary artery disease (CAD). However, there are several limitations when applying standard, ECG-gated stress/rest perfusion CMR to patients with atrial fibrillation (AF). In this study we investigate an approach with no ECG gating and a rapid rest/stress perfusion protocol to determine its accuracy for detection of CAD in patients with AF. 26 patients with AF underwent a rapid rest/regadenoson stress CMR perfusion imaging protocol, and all patients had X-ray coronary angiography. An ungated radial myocardial perfusion sequence was used. Imaging protocol included: rest perfusion image acquisition, followed nearly immediately by administration of regadenoson to induce hyperemia, 60 s wait, and stress image acquisition. CMR perfusion images were interpreted by three blinded readers as normal or abnormal. Diagnostic accuracy was evaluated by comparison to X-ray angiography. 21 of the CMR rest/stress perfusion scans were negative, and 5 were positive by angiography criteria. Majority results of the ungated datasets from all of the readers showed a sensitivity, specificity and accuracy of 80, 100 and 96%, respectively, for detection of CAD. An ungated, rapid rest/stress regadenoson perfusion CMR protocol appears to be useful for the diagnosis of obstructive CAD in patients with AF.

Split Bregman multicoil accelerated reconstruction technique: A new framework for rapid reconstruction of cardiac perfusion MRI

PURPOSE

Rapid reconstruction of undersampled multicoil MRI data with iterative constrained reconstruction method is a challenge. The authors sought to develop a new substitution based variable splitting algorithm for faster reconstruction of multicoil cardiac perfusion MRI data.

METHODS

The new method, split Bregman multicoil accelerated reconstruction technique (SMART), uses a combination of split Bregman based variable splitting and iterative reweighting techniques to achieve fast convergence. Total variation constraints are used along the spatial and temporal dimensions. The method is tested on nine ECG-gated dog perfusion datasets, acquired with a 30-ray golden ratio radial sampling pattern and ten ungated human perfusion datasets, acquired with a 24-ray golden ratio radial sampling pattern. Image quality and reconstruction speed are evaluated and compared to a gradient descent (GD) implementation and to multicoil k-t SLR, a reconstruction technique that uses a combination of sparsity and low rank constraints.

RESULTS

Comparisons based on blur metric and visual inspection showed that SMART images had lower blur and better texture as compared to the GD implementation. On average, the GD based images had an ∼18% higher blur metric as compared to SMART images. Reconstruction of dynamic contrast enhanced (DCE) cardiac perfusion images using the SMART method was ∼6 times faster than standard gradient descent methods. k-t SLR and SMART produced images with comparable image quality, though SMART was ∼6.8 times faster than k-t SLR.

CONCLUSIONS

The SMART method is a promising approach to reconstruct good quality multicoil images from undersampled DCE cardiac perfusion data rapidly.

Review of Journal of Cardiovascular Magnetic Resonance 2015

There were 116 articles published in the Journal of Cardiovascular Magnetic Resonance (JCMR) in 2015, which is a 14 % increase on the 102 articles published in 2014. The quality of the submissions continues to increase. The 2015 JCMR Impact Factor (which is published in June 2016) rose to 5.75 from 4.72 for 2014 (as published in June 2015), which is the highest impact factor ever recorded for JCMR. The 2015 impact factor means that the JCMR papers that were published in 2013 and 2014 were cited on average 5.75 times in 2015. The impact factor undergoes natural variation according to citation rates of papers in the 2 years following publication, and is significantly influenced by highly cited papers such as official reports. However, the progress of the journal's impact over the last 5 years has been impressive. Our acceptance rate is <25 % and has been falling because the number of articles being submitted has been increasing. In accordance with Open-Access publishing, the JCMR articles go on-line as they are accepted with no collating of the articles into sections or special thematic issues. For this reason, the Editors have felt that it is useful once per calendar year to summarize the papers for the readership into broad areas of interest or theme, so that areas of interest can be reviewed in a single article in relation to each other and other recent JCMR articles. The papers are presented in broad themes and set in context with related literature and previously published JCMR papers to guide continuity of thought in the journal. We hope that you find the open-access system increases wider reading and citation of your papers, and that you will continue to send your quality papers to JCMR for publication.

Robust dynamic myocardial perfusion CT deconvolution for accurate residue function estimation via adaptive-weighted tensor total variation regularization: a preclinical study

Dynamic myocardial perfusion computed tomography (MPCT) is a promising technique for quick diagnosis and risk stratification of coronary artery disease. However, one major drawback of dynamic MPCT imaging is the heavy radiation dose to patients due to its dynamic image acquisition protocol. In this work, to address this issue, we present a robust dynamic MPCT deconvolution algorithm via adaptive-weighted tensor total variation (AwTTV) regularization for accurate residue function estimation with low-mA s data acquisitions. For simplicity, the presented method is termed 'MPD-AwTTV'. More specifically, the gains of the AwTTV regularization over the original tensor total variation regularization are from the anisotropic edge property of the sequential MPCT images. To minimize the associative objective function we propose an efficient iterative optimization strategy with fast convergence rate in the framework of an iterative shrinkage/thresholding algorithm. We validate and evaluate the presented algorithm using both digital XCAT phantom and preclinical porcine data. The preliminary experimental results have demonstrated that the presented MPD-AwTTV deconvolution algorithm can achieve remarkable gains in noise-induced artifact suppression, edge detail preservation, and accurate flow-scaled residue function and MPHM estimation as compared with the other existing deconvolution algorithms in digital phantom studies, and similar gains can be obtained in the porcine data experiment.

Accelerated dynamic MRI using patch regularization for implicit motion compensation

PURPOSE

To introduce a fast algorithm for motion-compensated accelerated dynamic MRI.

METHODS

An efficient patch smoothness regularization scheme, which implicitly compensates for inter-frame motion, is introduced to recover dynamic MRI data from highly undersampled measurements. The regularization prior is a sum of distances between each rectangular patch in the dataset with other patches in the dataset using a saturating distance metric. Unlike current motion estimation and motion compensation (ME-MC) methods, the proposed scheme does not require reference frames or complex motion models. The proposed algorithm, which alternates between inter-patch shrinkage step and conjugate gradient algorithm, is considerably more computationally efficient than ME-MC methods. The reconstructions obtained using the proposed algorithm is compared against state-of-the-art methods.

RESULTS

The proposed method is observed to yield reconstructions with minimal spatiotemporal blurring and motion artifacts. In comparison to the existing state-of-the-art ME-MC methods, PRICE provides comparable or even better image quality with faster reconstruction times (approximately nine times faster).

CONCLUSION

The presented scheme enables computationally efficient and effective motion-compensated reconstruction in a variety of applications with large inter-frame motion and contrast changes. This algorithm could be seen as an alternative over the current state-of-the-art ME-MC schemes that are computationally expensive. Magn Reson Med 77:1238-1248, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Interstudy repeatability of self-gated quantitative myocardial perfusion MRI

PURPOSE

To evaluate the interstudy repeatability of multislice quantitative cardiovascular magnetic resonance myocardial blood flow (MBF), myocardial perfusion reserve (MPR), and extracellular volume (ECV). A unique saturation recovery self-gated acquisition was used for the perfusion scans.

MATERIALS AND METHODS

An ungated golden angle radial turboFLASH pulse sequence was used to scan 10 subjects on two separate days on a 3T scanner. A single saturation pulse was followed by a set of four slices. Rest and hyperemia scans were acquired during free breathing. The images were reconstructed using an iterative algorithm with spatiotemporal constraints. The ungated images were retrospectively binned (self-gated) into near-systole and near-diastole. Deformable registration was performed to adjust for respiratory and residual cardiac motion, and the data were fit with a Fermi model to estimate the interstudy repeatability of quantitative self-gated MBF and MPR.

RESULTS

The coefficient of variation (CoV) of the territorial MPR using the self-gated near-systole data was 18.6%. The self-gated near-diastole data gave less good CoV of MPR, equal to 46.2%. For MBFs, and using smaller (segmental) regions, the CoVs were 20.1% and 22.7% for the estimation of myocardial blood flow at stress and rest, respectively, using the self-gated near-systole data. The self-gated near-diastole data gave CoV = 48.6% and 44.9% for stress and rest.

CONCLUSION

The self-gated free-breathing technique for quantification of myocardial blood flow showed good repeatability for near-systole, with results comparable to published studies on interstudy repeatability of quantitative myocardial perfusion MRI using ECG-gating and breath-holds. Self-gated near-diastole data results were less repeatable. J. Magn. Reson. Imaging 2016;43:1369-1378.