Feasibility of multiple-view myocardial perfusion MRI using radial simultaneous multi-slice acquisitions

Quantitative myocardial perfusion with a hybrid 2D simultaneous multi-slice sequence

PURPOSE

To evaluate a novel 2D simultaneous multi-slice (SMS) myocardial perfusion acquisition and compare directly to a published quantitative 3D stack-of-stars (SoS) acquisition.

METHODS

A hybrid saturation recovery radial 2D SMS sequence following a single saturation was created for the quantification of myocardial blood flow (MBF). This sequence acquired three slices simultaneously and generated an arterial input function (AIF) using the first 24 rays. Validation was done in a novel way by alternating heartbeats between the hybrid 2D SMS and the 3D SoS acquisitions. Initial studies were done to study the effects of using only every other beat for the 2D SMS in two subjects, and for the 3D SoS in four subjects. The proposed alternating acquisitions were then performed in ten dog studies at rest, four dog studies at adenosine stress, and two human resting studies. Quantitative MBF analysis was performed for 2D SMS and 3D SoS separately, using a compartment model.

RESULTS

Acquiring every-other-beat data resulted in 6 ± 5% ("ideal") and 11 ± 8% ("practical") perfusion changes for both 2D SMS and 3D SoS methods. For alternating acquisitions, 2D SMS and 3D SoS quantitative perfusion values were comparable for both the twelve rest studies (2D SMS: 0.69 ± 0.16 vs 3D: 0.69 ± 0.15 ml/g/min, p = 0.55) and the four stress studies (2D SMS: 1.28 ± 0.22 vs 3D: 1.30 ± 0.24 ml/g/min, p = 0.61).

CONCLUSION

Every-other-beat acquisition changed estimated perfusion values relatively little for both sequences. The quantitative hybrid radial 2D SMS myocardial first-pass perfusion imaging sequence gave results similar to 3D perfusion when compared directly with an alternating beat acquisition.

Contrast-optimal simultaneous multi-slice bSSFP cine cardiac imaging at 0.55 T

PURPOSE

To determine if contemporary 0.55 T MRI supports the use of contrast-optimal flip angles (FA) for simultaneous multi-slice (SMS) balanced SSFP (bSSFP) cardiac function assessment, which is impractical at conventional field strengths because of excessive SAR and/or banding artifacts.

METHODS

Blipped-CAIPI bSSFP was combined with spiral sampling for ventricular function assessment at 0.55 T. Cine movies with single band and SMS factors of 2 and 3 (SMS 2 and 3), and FA ranging from 60° to 160°, were acquired in seven healthy volunteers. Left ventricular blood and myocardial signal intensity (SI) normalized by background noise and blood-myocardium contrast were measured and compared across acquisition settings.

RESULTS

Myocardial SI was slightly higher in single band than in SMS and decreased with an increasing FA. Blood SI increased as the FA increased for single band, and increment was small for FA ≥120°. Blood SI for SMS 2 and 3 increased with an increasing FA up to ∼100°. Blood-myocardium contrast increased with an increasing FA for single band, peaked at FA = 160° (systole: 28.43, diastole: 29.15), attributed mainly to reduced myocardial SI when FA ≥120°. For SMS 2, contrast peaked at 120° (systole: 21.43, diastole: 19.85). For SMS 3, contrast peaked at 120° in systole (16.62) and 100° in diastole (19.04).

CONCLUSIONS

Contemporary 0.55 T MR scanners equipped with high-performance gradient systems allow the use of contrast-optimal FA for SMS accelerated bSSFP cine examinations without compromising image quality. The contrast-optimal FA was found to be 140° to 160° for single band and 100° to 120° for SMS 2 and 3.

Simultaneous multi-slice T1 mapping using MOLLI with blipped CAIPIRINHA bSSFP

BACKGROUND

This study evaluates the possibility for replacing conventional 3 slices, 3 breath-holds MOLLI cardiac T1 mapping with single breath-hold 3 simultaneous multi-slice (SMS3) T1 mapping using blipped-CAIPIRINHA SMS-bSSFP MOLLI sequence. As a major drawback, SMS-bSSFP presents unique artefacts arising from side-lobe slice excitations that are explained by imperfect RF modulation rendering and bSSFP low flip angle enhancement. Amplitude-only RF modulation (AM) is proposed to reduce these artefacts in SMS-MOLLI compared to conventional Wong multi-band RF modulation (WM).

MATERIALS AND METHODS

Phantoms and ten healthy volunteers were imaged at 1.5 T using a modified blipped-CAIPIRINHA SMS-bSSFP MOLLI sequence with 3 simultaneous slices. WM-SMS3 and AM-SMS3 were compared to conventional single-slice (SMS1) MOLLI. First, SNR degradation and T1 accuracy were measured in phantoms. Second, artefacts from side-lobe excitations were evaluated in a phantom designed to reproduce fat presence near the heart. Third, the occurrence of these artefacts was observed in volunteers, and their impact on T1 quantification was compared between WM-SMS3 and AM-SMS3 with conventional MOLLI as a reference.

RESULTS

In the phantom, larger slice gaps and slice thicknesses yielded higher SNR. There was no significant difference of T1 values between conventional MOLLI and SMS3-MOLLI (both WM and AM). Positive banding artefacts were identified from fat neighbouring the targeted FOV due to side-lobe excitations from WM and the unique bSSFP signal profile. AM RF pulses reduced these artefacts by 38%. In healthy volunteers, AM-SMS3-MOLLI showed similar artefact reduction compared to WM-SMS3-MOLLI (3 ± 2 vs 5 ± 3 corrupted LV segments out of 16). In-vivo native T1 values obtained from conventional MOLLI and AM-SMS3-MOLLI were equivalent in LV myocardium (SMS1-T1 = 935.5 ± 36.1 ms; AM-SMS3-T1 = 933.8 ± 50.2 ms; P = 0.436) and LV blood pool (SMS1-T1 = 1475.4 ± 35.9 ms; AM-SMS3-T1 = 1452.5 ± 70.3 ms; P = 0.515). Identically, no differences were found between SMS1 and SMS3 postcontrast T1 values in the myocardium (SMS1-T1 = 556.0 ± 19.7 ms; SMS3-T1 = 521.3 ± 28.1 ms; P = 0.626) and the blood (SMS1-T1 = 478 ± 65.1 ms; AM-SMS3-T1 = 447.8 ± 81.5; P = 0.085).

CONCLUSIONS

Compared to WM RF modulation, AM SMS-bSSFP MOLLI was able to reduce side-lobe artefacts considerably, providing promising results to image the three levels of the heart in a single breath hold. However, few artefacts remained even using AM-SMS-bSSFP due to residual RF imperfections. The proposed blipped-CAIPIRINHA MOLLI T1 mapping sequence provides accurate in vivo T1 quantification in line with those obtained with a single slice acquisition.

A Slice-Low-Rank Plus Sparse (slice-L + S) Reconstruction Method for k-t Undersampled Multiband First-Pass Myocardial Perfusion MRI

Purpose

The synergistic use of k‐t undersampling and multiband (MB) imaging has the potential to provide extended slice coverage and high spatial resolution for first‐pass perfusion MRI. The low‐rank plus sparse (L + S) model has shown excellent performance for accelerating single‐band (SB) perfusion MRI.

Methods

A MB data consistency method employing ESPIRiT maps and through‐plane coil information was developed. This data consistency method was combined with the temporal L + S constraint to form the slice‐L + S method. Slice‐L + S was compared to SB L + S and the sequential operations of split slice‐GRAPPA and SB L + S (seq‐SG‐L + S) using synthetic data formed from multislice SB images. Prospectively k‐t undersampled MB data were also acquired and reconstructed using seq‐SG‐L + S and slice‐L + S.

Results

Using synthetic data with total acceleration rates of 6–12, slice‐L + S outperformed SB L + S and seq‐SG‐L + S (N = 7 subjects) with respect to normalized RMSE and the structural similarity index (P < 0.05 for both). For the specific case with MB factor = 3 and rate 3 undersampling, or for SB imaging with rate 9 undersampling (N = 7 subjects), the normalized RMSE values were 0.037 ± 0.007, 0.042 ± 0.005, and 0.031 ± 0.004; and the structural similarity index values were 0.88 ± 0.03, 0.85 ± 0.03, and 0.89 ± 0.02 for SB L + S, seq‐SG‐L + S, and slice‐L + S, respectively (P < 0.05 for both). For prospectively undersampled MB data, slice‐L + S provided better image quality than seq‐SG‐L + S for rate 6 (N = 7) and rate 9 acceleration (N = 7) as scored by blinded experts.

Conclusion

Slice‐L + S outperformed SB‐L + S and seq‐SG‐L + S and provides 9 slice coverage of the left ventricle with a spatial resolution of 1.5 mm × 1.5 mm with good image quality.

Simultaneous multislice steady-state free precession myocardial perfusion with full left ventricular coverage and high resolution at 1.5 T

Purpose

To implement and evaluate a simultaneous multi‐slice balanced SSFP (SMS‐bSSFP) perfusion sequence and compressed sensing reconstruction for cardiac MR perfusion imaging with full left ventricular (LV) coverage (nine slices/heartbeat) and high spatial resolution (1.4 × 1.4 mm²) at 1.5T.

Methods

A preliminary study was performed to evaluate the performance of blipped controlled aliasing in parallel imaging (CAIPI) and RF‐CAIPI with gradient‐controlled local Larmor adjustment (GC‐LOLA) in the presence of fat. A nine‐slice SMS‐bSSFP sequence using RF‐CAIPI with GC‐LOLA with high spatial resolution (1.4 × 1.4 mm²) and a conventional three‐slice sequence with conventional spatial resolution (1.9 × 1.9 mm²) were then acquired in 10 patients under rest conditions. Qualitative assessment was performed to assess image quality and perceived signal‐to‐noise ratio (SNR) on a 4‐point scale (0: poor image quality/low SNR; 3: excellent image quality/high SNR), and the number of myocardial segments with diagnostic image quality was recorded. Quantitative measurements of myocardial sharpness and upslope index were performed.

Results

Fat signal leakage was significantly higher for blipped CAIPI than for RF‐CAIPI with GC‐LOLA (7.9% vs. 1.2%, p = 0.010). All 10 SMS‐bSSFP perfusion datasets resulted in 16/16 diagnostic myocardial segments. There were no significant differences between the SMS and conventional acquisitions in terms of image quality (2.6 ± 0.6 vs. 2.7 ± 0.2, p = 0.8) or perceived SNR (2.8 ± 0.3 vs. 2.7 ± 0.3, p = 0.3). Inter‐reader variability was good for both image quality (ICC = 0.84) and perceived SNR (ICC = 0.70). Myocardial sharpness was improved using the SMS sequence compared to the conventional sequence (0.37 ± 0.08 vs 0.32 ± 0.05, p < 0.001). There was no significant difference between measurements of upslope index for the SMS and conventional sequences (0.11 ± 0.04 vs. 0.11 ± 0.03, p = 0.84).

Conclusion

SMS‐bSSFP with multiband factor 3 and compressed sensing reconstruction enables cardiac MR perfusion imaging with three‐fold increased spatial coverage and improved myocardial sharpness compared to a conventional sequence, without compromising perceived SNR, image quality, upslope index or number of diagnostic segments.

Deep learning for radial SMS myocardial perfusion reconstruction using the 3D residual booster U-net

PURPOSE

To develop an end-to-end deep learning solution for quickly reconstructing radial simultaneous multi-slice (SMS) myocardial perfusion datasets with comparable quality to the pixel tracking spatiotemporal constrained reconstruction (PT-STCR) method.

METHODS

Dynamic contrast enhanced (DCE) radial SMS myocardial perfusion data were obtained from 20 subjects who were scanned at rest and/or stress with or without ECG gating using a saturation recovery radial CAIPI turboFLASH sequence. Input to the networks consisted of complex coil combined images reconstructed using the inverse Fourier transform of undersampled radial SMS k-space data. Ground truth images were reconstructed using the PT-STCR pipeline. The performance of the residual booster 3D U-Net was tested by comparing it to state-of-the-art network architectures including MoDL, CRNN-MRI, and other U-Net variants.

RESULTS

Results demonstrate significant improvements in speed requiring approximately 8 seconds to reconstruct one radial SMS dataset which is approximately 200 times faster than the PT-STCR method. Images reconstructed with the residual booster 3D U-Net retain quality of ground truth PT-STCR images (0.963 SSIM/40.238 PSNR/0.147 NRMSE). The residual booster 3D U-Net has superior performance compared to existing network architectures in terms of image quality, temporal dynamics, and reconstruction time.

CONCLUSION

Residual and booster learning combined with the 3D U-Net architecture was shown to be an effective network for reconstructing high-quality images from undersampled radial SMS datasets while bypassing the reconstruction time of the PT-STCR method.

All-systolic first-pass myocardial rest perfusion at a long saturation time using simultaneous multi-slice imaging and compressed sensing acceleration

PURPOSE

To enable all-systolic first-pass rest myocardial perfusion with long saturation times. To investigate the change in perfusion contrast and dark rim artefacts through simulations and surrogate measurements.

METHODS

Simulations were employed to investigate optimal saturation time for myocardium-perfusion defect contrast and blood-to-myocardium signal ratios. Two saturation recovery blocks with long/short saturation times (LTS/STS) were employed to image 3 slices at end-systole and diastole. Simultaneous multi-slice balanced steady state free precession imaging and compressed sensing acceleration were combined. The sequence was compared to a 3 slice-by-slice clinical protocol in 10 patients. Quantitative assessment of myocardium-peak pre contrast and blood-to-myocardium signal ratios, as well as qualitative assessment of perceived SNR, image quality, blurring, and dark rim artefacts, were performed.

RESULTS

Simulations showed that with a bolus of 0.075 mmol/kg, a LTS of 240-470 ms led to a relative increase in myocardium-perfusion defect contrast of 34% ± 9%-28% ± 27% than a STS = 120 ms, while reducing blood-to-myocardium signal ratio by 18% ± 10%-32% ± 14% at peak myocardium. With a bolus of 0.05 mmol/kg, LTS was 320-570 ms with an increase in myocardium-perfusion defect contrast of 63% ± 13%-62% ± 29%. Across patients, LTS led to an average increase in myocardium-peak pre contrast of 59% (P < .001) at peak myocardium and a lower blood-to-myocardium signal ratio of 47% (P < .001) and 15% (P < .001) at peak blood/myocardium. LTS had improved motion robustness (P = .002), image quality (P < .001), and decreased dark rim artefacts (P = .008) than the clinical protocol.

CONCLUSION

All-systolic rest perfusion can be achieved by combining simultaneous multi-slice and compressed sensing acceleration, enabling 3-slice cardiac coverage with reduced motion and dark rim artefacts. Numerical simulations indicate that myocardium-perfusion defect contrast increases at LTS.

Aliasing artifact reduction in spiral real-time MRI

PURPOSE

To mitigate a common artifact in spiral real-time MRI, caused by aliasing of signal outside the desired FOV. This artifact frequently occurs in midsagittal speech real-time MRI.

METHODS

Simulations were performed to determine the likely origin of the artifact. Two methods to mitigate the artifact are proposed. The first approach, denoted as "large FOV" (LF), keeps an FOV that is large enough to include the artifact signal source during reconstruction. The second approach, denoted as "estimation-subtraction" (ES), estimates the artifact signal source before subtracting a synthetic signal representing that source in multicoil k-space raw data. Twenty-five midsagittal speech-production real-time MRI data sets were used to evaluate both of the proposed methods. Reconstructions without and with corrections were evaluated by two expert readers using a 5-level Likert scale assessing artifact severity. Reconstruction time was also compared.

RESULTS

The origin of the artifact was found to be a combination of gradient nonlinearity and imperfect anti-aliasing in spiral sampling. The LF and ES methods were both able to substantially reduce the artifact, with an averaged qualitative score improvement of 1.25 and 1.35 Likert levels for LF correction and ES correction, respectively. Average reconstruction time without correction, with LF correction, and with ES correction were 160.69 ± 1.56, 526.43 ± 5.17, and 171.47 ± 1.71 ms/frame.

CONCLUSION

Both proposed methods were able to reduce the spiral aliasing artifacts, with the ES-reduction method being more effective and more time efficient.

Improved simultaneous multislice cardiac MRI using readout concatenated k-space SPIRiT (ROCK-SPIRiT)

PURPOSE

To develop and evaluate a simultaneous multislice (SMS) reconstruction technique that provides noise reduction and leakage blocking for highly accelerated cardiac MRI.

METHODS

ReadOut Concatenated k-space SPIRiT (ROCK-SPIRiT) uses the concept of readout concatenation in image domain to represent SMS encoding, and performs coil self-consistency as in SPIRiT-type reconstruction in an extended k-space, while allowing regularization for further denoising. The proposed method is implemented with and without regularization, and validated on retrospectively SMS-accelerated cine imaging with three-fold SMS and two-fold in-plane acceleration. ROCK-SPIRiT is compared with two leakage-blocking SMS reconstruction methods: readout-SENSE-GRAPPA and split slice-GRAPPA. Further evaluation and comparisons are performed using prospectively SMS-accelerated cine imaging.

RESULTS

Results on retrospectively three-fold SMS and two-fold in-plane accelerated cine imaging show that ROCK-SPIRiT without regularization significantly improves on existing methods in terms of PSNR (readout-SENSE-GRAPPA: 33.5 ± 3.2, split slice-GRAPPA: 34.1 ± 3.8, ROCK-SPIRiT: 35.0 ± 3.3) and SSIM (readout-SENSE-GRAPPA: 84.4 ± 8.9, split slice-GRAPPA: 85.0 ± 8.9, ROCK-SPIRiT: 88.2 ± 6.6 [in percentage]). Regularized ROCK-SPIRiT significantly outperforms all methods, as characterized by these quantitative metrics (PSNR: 37.6 ± 3.8, SSIM: 94.2 ± 4.1 [in percentage]). The prospectively five-fold SMS and two-fold in-plane accelerated data show that ROCK-SPIRiT and regularized ROCK-SPIRiT have visually improved image quality compared with existing methods.

CONCLUSION

The proposed ROCK-SPIRiT technique reduces noise and interslice leakage in accelerated SMS cardiac cine MRI, improving on existing methods both quantitatively and qualitatively.

Combined simultaneous multislice bSSFP and compressed sensing for first-pass myocardial perfusion at 1.5 T with high spatial resolution and coverage

PURPOSE

To implement and evaluate a pseudorandom undersampling scheme for combined simultaneous multislice (SMS) balanced SSFP (bSSFP) and compressed-sensing (CS) reconstruction to enable myocardial perfusion imaging with high spatial resolution and coverage at 1.5 T.

METHODS

A prospective pseudorandom undersampling scheme that is compatible with SMS-bSSFP phase-cycling requirements and CS was developed. The SMS-bSSFP CS with pseudorandom and linear undersampling schemes were compared in a phantom. A high-resolution (1.4 × 1.4 mm2 ) six-slice SMS-bSSFP CS perfusion sequence was compared with a conventional (1.9 × 1.9 mm2 ) three-slice sequence in 10 patients. Qualitative assessment of image quality, perceived SNR, and number of diagnostic segments and quantitative measurements of sharpness, upslope index, and contrast ratio were performed.

RESULTS

In phantom experiments, pseudorandom undersampling resulted in residual artifact (RMS error) reduction by a factor of 7 compared with linear undersampling. In vivo, the proposed sequence demonstrated higher perceived SNR (2.9 ± 0.3 vs. 2.2 ± 0.6, P = .04), improved sharpness (0.35 ± 0.03 vs. 0.32 ± 0.05, P = .01), and a higher number of diagnostic segments (100% vs. 94%, P = .03) compared with the conventional sequence. There were no significant differences between the sequences in terms of image quality (2.5 ± 0.4 vs. 2.8 ± 0.2, P = .08), upslope index (0.11 ± 0.02 vs. 0.10 ± 0.01, P = .3), or contrast ratio (3.28 ± 0.35 vs. 3.36 ± 0.43, P = .7).

CONCLUSION

A pseudorandom k-space undersampling compatible with SMS-bSSFP and CS reconstruction has been developed and enables cardiac MR perfusion imaging with increased spatial resolution and myocardial coverage, increased number of diagnostic segments and perceived SNR, and no difference in image quality, upslope index, and contrast ratio.

Non-iterative image reconstruction from sparse magnetic resonance imaging radial data without priors

The state-of-the-art approaches for image reconstruction using under-sampled k-space data are compressed sensing based. They are iterative algorithms that optimize objective functions with spatial and/or temporal constraints. This paper proposes a non-iterative algorithm to estimate the un-measured data and then to reconstruct the image with the efficient filtered backprojection algorithm. The feasibility of the proposed method is demonstrated with a patient magnetic resonance imaging study. The proposed method is also compared with the state-of-the-art iterative compressed-sensing image reconstruction method using the total-variation optimization norm.

Non-invasive imaging techniques to assess myocardial perfusion

INTRODUCTION

Coronary artery disease affecting myocardial perfusion continues to be the leading cause of cardiovascular morbidity and mortality worldwide. While invasive evaluation based on coronary angiography and flow measurements still is considered the reference technique for assessing myocardial perfusion, technological evolutions in noninvasive imaging techniques resulted in them playing a more important role in assessing myocardial perfusion influencing therapeutic decision-making and prognostication.

AREAS COVERED

Different imaging modalities are used to evaluate coronary perfusion, including echocardiography, coronary computerized tomography scan, magnetic resonance imaging, and nuclear myocardial perfusion imaging. Through a combination of different techniques, it is possible to describe coronary artery anatomy and the diameter of the epicardial vessels but more recently also of the coronary microcirculation. Quantification of myocardial perfusion is feasible both at baseline and during pharmacological or physiological stress.

EXPERT OPINION

The objective of this state-of-the-art paper is to review the recent evolutions in imaging methods to estimate myocardial perfusion and to discuss the diagnostic strengths and limitations of the different techniques. The new ultrasound technologies and the hybrid approaches seem to be the future is these fields.

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.