Simultaneous multislice imaging for native myocardial T1 mapping: Improved spatial coverage in a single breath-hold

Motion-compensated T1 mapping in cardiovascular magnetic resonance imaging: a technical review

T 1 mapping is becoming a staple magnetic resonance imaging method for diagnosing myocardial diseases such as ischemic cardiomyopathy, hypertrophic cardiomyopathy, myocarditis, and more. Clinically, most T 1 mapping sequences acquire a single slice at a single cardiac phase across a 10 to 15-heartbeat breath-hold, with one to three slices acquired in total. This leaves opportunities for improving patient comfort and information density by acquiring data across multiple cardiac phases in free-running acquisitions and across multiple respiratory phases in free-breathing acquisitions. Scanning in the presence of cardiac and respiratory motion requires more complex motion characterization and compensation. Most clinical mapping sequences use 2D single-slice acquisitions; however newer techniques allow for motion-compensated reconstructions in three dimensions and beyond. To further address confounding factors and improve measurement accuracy, T 1 maps can be acquired jointly with other quantitative parameters such as T 2 , T 2 ∗ , fat fraction, and more. These multiparametric acquisitions allow for constrained reconstruction approaches that isolate contributions to T 1 from other motion and relaxation mechanisms. In this review, we examine the state of the literature in motion-corrected and motion-resolved T 1 mapping, with potential future directions for further technical development and clinical translation.

Improved Simultaneous Multi-slice imaging with Composition of k-space Interpolations (SMS-COOKIE) for myocardial T1 mapping

The aim of this study is to develop and evaluate a regularized Simultaneous Multi-Slice (SMS) reconstruction method for improved Cardiac Magnetic Resonance Imaging (CMR). The proposed reconstruction method, SMS with COmpOsition of k-space IntErpolations (SMS-COOKIE) combines the advantages of Iterative Self-consistent Parallel Imaging Reconstruction (SPIRiT) and split slice-Generalized Autocalibrating Partially Parallel Acquisitions (GRAPPA), while allowing regularization for further noise reduction. The proposed SMS-COOKIE was implemented with and without regularization, and validated using a Saturation Pulse-Prepared Heart rate Independent inversion REcovery (SAPPHIRE) myocardial T1 mapping sequence. The performance of the proposed reconstruction method was compared to ReadOut (RO)-SENSE-GRAPPA and split slice-GRAPPA, on both retrospectively and prospectively three-fold SMS-accelerated data with an additional two-fold in-plane acceleration. All SMS reconstruction methods yielded similar T1 values compared to single band imaging. SMS-COOKIE showed lower spatial variability in myocardial T1 with significant improvement over RO-SENSE-GRAPPA and split slice-GRAPPA (P < 10-4). The proposed method with additional locally low rank (LLR) regularization reduced the spatial variability, again with significant improvement over RO-SENSE-GRAPPA and split slice-GRAPPA (P < 10-4). In conclusion, improved reconstruction quality was achieved with the proposed SMS-COOKIE, which also provided lower spatial variability with significant improvement over split slice-GRAPPA.

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.

Signal intensity informed multi-coil encoding operator for physics-guided deep learning reconstruction of highly accelerated myocardial perfusion CMR

PURPOSE

To develop a physics-guided deep learning (PG-DL) reconstruction strategy based on a signal intensity informed multi-coil (SIIM) encoding operator for highly-accelerated simultaneous multislice (SMS) myocardial perfusion cardiac MRI (CMR).

METHODS

First-pass perfusion CMR acquires highly-accelerated images with dynamically varying signal intensity/SNR following the administration of a gadolinium-based contrast agent. Thus, using PG-DL reconstruction with a conventional multi-coil encoding operator leads to analogous signal intensity variations across different time-frames at the network output, creating difficulties in generalization for varying SNR levels. We propose to use a SIIM encoding operator to capture the signal intensity/SNR variations across time-frames in a reformulated encoding operator. This leads to a more uniform/flat contrast at the output of the PG-DL network, facilitating generalizability across time-frames. PG-DL reconstruction with the proposed SIIM encoding operator is compared to PG-DL with conventional encoding operator, split slice-GRAPPA, locally low-rank (LLR) regularized reconstruction, low-rank plus sparse (L + S) reconstruction, and regularized ROCK-SPIRiT.

RESULTS

Results on highly accelerated free-breathing first pass myocardial perfusion CMR at three-fold SMS and four-fold in-plane acceleration show that the proposed method improves upon the reconstruction methods use for comparison. Substantial noise reduction is achieved compared to split slice-GRAPPA, and aliasing artifacts reduction compared to LLR regularized reconstruction, L + S reconstruction and PG-DL with conventional encoding. Furthermore, a qualitative reader study indicated that proposed method outperformed all methods.

CONCLUSION

PG-DL reconstruction with the proposed SIIM encoding operator improves generalization across different time-frames /SNRs in highly accelerated perfusion CMR.

Cardiac phase-resolved late gadolinium enhancement imaging

Late gadolinium enhancement (LGE) with cardiac magnetic resonance (CMR) imaging is the clinical reference for assessment of myocardial scar and focal fibrosis. However, current LGE techniques are confined to imaging of a single cardiac phase, which hampers assessment of scar motility and does not allow cross-comparison between multiple phases. In this work, we investigate a three step approach to obtain cardiac phase-resolved LGE images: (1) Acquisition of cardiac phase-resolved imaging data with varying T₁ weighting. (2) Generation of semi-quantitative T1* maps for each cardiac phase. (3) Synthetization of LGE contrast to obtain functional LGE images. The proposed method is evaluated in phantom imaging, six healthy subjects at 3T and 20 patients at 1.5T. Phantom imaging at 3T demonstrates consistent contrast throughout the cardiac cycle with a coefficient of variation of 2.55 ± 0.42%. In-vivo results show reliable LGE contrast with thorough suppression of the myocardial tissue is healthy subjects. The contrast between blood and myocardium showed moderate variation throughout the cardiac cycle in healthy subjects (coefficient of variation 18.2 ± 3.51%). Images were acquired at 40–60 ms and 80 ms temporal resolution, at 3T and 1.5, respectively. Functional LGE images acquired in patients with myocardial scar visualized scar tissue throughout the cardiac cycle, albeit at noticeably lower imaging resolution and noise resilience than the reference technique. The proposed technique bears the promise of integrating the advantages of phase-resolved CMR with LGE imaging, but further improvements in the acquisition quality are warranted for clinical use.

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 Multi-Slice Cardiac MR Multitasking for Motion-Resolved, Non-ECG, Free-Breathing T1–T2 Mapping

The aim of this study is to simultaneously quantify T1/T2 across three slices of the left-ventricular myocardium without breath-holds or ECG monitoring, all within a 3 min scan. Radial simultaneous multi-slice (SMS) encoding, self-gating, and image reconstruction was incorporated into the cardiovascular magnetic resonance (CMR) Multitasking framework to simultaneously image three short-axis slices. A T2prep-IR FLASH sequence with two flip angles was designed and implemented to allow B1+-robust T1 and T2 mapping. The proposed Multitasking-SMS method was validated in a standardized phantom and 10 healthy volunteers, comparing T1 and T2 measurements and scan-rescan repeatability against corresponding reference methods in one layer of phantom vials and in 16 American Heart Association (AHA) myocardial segments. In phantom, Multitasking-SMS T1/T2 measurements showed substantial correlation (R² > 0.996) and excellent agreement [intraclass correlation coefficients (ICC) ≥ 0.999)] with reference measurements. In healthy volunteers, Multitasking-SMS T1/T2 maps reported similar myocardial T1/T2 values (1,215 ± 91.0/41.5 ± 6.3 ms) to the reference myocardial T1/T2 values (1,239 ± 67.5/42.7 ± 4.1 ms), with P = 0.347 and P = 0.296, respectively. Bland–Altman analyses also demonstrated good in vivo repeatability in both the multitasking and references, with segment-wise coefficients of variation of 4.7% (multitasking T1), 8.9% (multitasking T2), 2.4% [modified look-locker inversion recovery (MOLLI)], and 4.6% (T2-prep FLASH), respectively. In summary, multitasking-SMS is feasible for free-breathing, non-ECG, myocardial T1/T2 quantification in 16 AHA segments over 3 short-axis slices in 3 min. The method shows the great potential for reducing exam time for quantitative CMR without ECG or breath-holds.

Free-breathing 3D cardiac T1 mapping with transmit B1 correction at 3T

PURPOSE

To develop a cardiac T₁ mapping method for free-breathing 3D T₁ mapping of the whole heart at 3 T with transmit B₁ ( B 1 + ) correction.

METHODS

A free-breathing, electrocardiogram-gated inversion-recovery sequence with spoiled gradient-echo readout was developed and optimized for cardiac T₁ mapping at 3 T. High-frame-rate dynamic images were reconstructed from sparse (k,t)-space data acquired along a stack-of-stars trajectory using a subspace-based method for accelerated imaging. Joint T₁ and flip-angle estimation was performed in T₁ mapping to improve its robustness to B 1 + inhomogeneity. Subject-specific timing of data acquisition was used in the estimation to account for natural heart-rate variations during the imaging experiment.

RESULTS

Simulations showed that accuracy and precision of T₁ mapping can be improved with joint T₁ and flip-angle estimation and optimized electrocardiogram-gated spoiled gradient echo-based inversion-recovery acquisition scheme. The phantom study showed good agreement between the T₁ maps from the proposed method and the reference method. Three-dimensional cardiac T₁ maps (40 slices) were obtained at a 1.9-mm in-plane and 4.5-mm through-plane spatial resolution from healthy subjects (n = 6) with an average imaging time of 14.2 ± 1.6 minutes (heartbeat rate: 64.2 ± 7.1 bpm), showing myocardial T₁ values comparable to those obtained from modified Look-Locker inversion recovery. The proposed method generated B 1 + maps with spatially smooth variation showing 21%-32% and 11%-15% variations across the septal-lateral and inferior-anterior regions of the myocardium in the left ventricle.

CONCLUSION

The proposed method allows free-breathing 3D T₁ mapping of the whole heart with transmit B₁ correction in a practical imaging time.

Free-breathing simultaneous T1 , T2 , and T2∗ quantification in the myocardium

PURPOSE

To implement a free-breathing sequence for simultaneous quantification of T 1 , T 2 , and T 2 ∗ for comprehensive tissue characterization of the myocardium in a single scan using a multi-gradient-echo readout with saturation and T 2 preparation pulses.

METHODS

In the proposed Saturation And T 2 -prepared Relaxometry with Navigator-gating (SATURN) technique, a series of multi-gradient-echo (GRE) images with different magnetization preparations was acquired during free breathing. A total of 35 images were acquired in 26.5 ± 14.9 seconds using multiple saturation times and T 2 preparation durations and with imaging at 5 echo times. Bloch simulations and phantom experiments were used to validate a 5-parameter fit model for accurate relaxometry. Free-breathing simultaneous T 1 , T 2 , and T 2 ∗ measurements were performed in 10 healthy volunteers and 2 patients using SATURN at 3T and quantitatively compared to conventional single-parameter methods such as SASHA for T 1 , T 2 -prepared bSSFP, and multi-GRE for T 2 ∗ .

RESULTS

Simulations confirmed accurate fitting with the 5-parameter model. Phantom measurements showed good agreement with the reference methods in the relevant range for in vivo measurements. Compared to single-parameter methods comparable accuracy was achieved. SATURN produced in vivo parameter maps that were visually comparable to single-parameter methods. No significant difference between T 1 , T 2 , and T 2 ∗ times acquired with SATURN and single-parameter methods was shown in quantitative measurements (SATURN T 1 = 1573 ± 86 ms , T 2 = 33.2 ± 3.6 ms , T 2 ∗ = 25.3 ± 6.1 ms ; conventional methods: T 1 = 1544 ± 107 ms , T 2 = 33.2 ± 3.6 ms , T 2 ∗ = 23.8 ± 5.5 ms ; P > . 2 ) CONCLUSION: SATURN enables simultaneous quantification of T 1 , T 2 , and T 2 ∗ in the myocardium for comprehensive tissue characterization with co-registered maps, in a single scan with good agreement to single-parameter methods.

Model-based reconstruction for simultaneous multi-slice T1 mapping using single-shot inversion-recovery radial FLASH

PURPOSE

To develop a single-shot multi-slice T 1 mapping method by combing simultaneous multi-slice (SMS) excitations, single-shot inversion-recovery (IR) radial fast low-angle shot (FLASH), and a nonlinear model-based reconstruction method.

METHODS

SMS excitations are combined with a single-shot IR radial FLASH sequence for data acquisition. A previously developed single-slice calibrationless model-based reconstruction is extended to SMS, formulating the estimation of parameter maps and coil sensitivities from all slices as a single nonlinear inverse problem. Joint-sparsity constraints are further applied to the parameter maps to improve T 1 precision. Validations of the proposed method are performed for a phantom and for the human brain and liver in 6 healthy adult subjects.

RESULTS

Phantom results confirm good T 1 accuracy and precision of the simultaneously acquired multi-slice T 1 maps in comparison to single-slice references. In vivo human brain studies demonstrate the better performance of SMS acquisitions compared to the conventional spoke-interleaved multi-slice acquisition using model-based reconstruction. Aside from good accuracy and precision, the results of 6 healthy subjects in both brain and abdominal studies confirm good repeatability between scan and re-scans. The proposed method can simultaneously acquire T 1 maps for 5 slices of a human brain ( 0.75 × 0.75 × 5 mm 3 ) or 3 slices of the abdomen ( 1.25 × 1.25 × 6 mm 3 ) within 4 seconds.

CONCLUSIONS

The IR SMS radial FLASH acquisition together with a nonlinear model-based reconstruction enable rapid high-resolution multi-slice T 1 mapping with good accuracy, precision, and repeatability.

Minimum TR radiofrequency-pulse design for rapid gradient echo sequences

PURPOSE

A framework to design radiofrequency (RF) pulses specifically to minimize the TR of gradient echo sequences is presented, subject to hardware and physiological constraints.

METHODS

Single-band and multiband (MB) RF pulses can be reduced in duration using variable-rate selective excitation (VERSE) VERSE for a range of flip angles; however, minimum-duration pulses do not guarantee minimum TR because these can lead to a high specific absorption rate (SAR). The optimal RF pulse is found by meeting spatial encoding, peripheral nerve stimulation (PNS) and SAR constraints. A TR reduction for a range of designs is achieved and an application of this in an MB cardiac balanced steady-state free-precession (bSSFP) experiment is presented. Gradient imperfections and their imaging effects are also considered.

RESULTS

Sequence TR with low-time bandwidth product (TBP) pulses, as used in bSSFP, was reduced up to 14%, and the TR when using high TBP pulses, as used in slab-selective imaging, was reduced by up to 72%. A breath-hold cardiac exam was reduced by 46% using both MB and the TR-optimal framework. The importance of RF-based correction of gradient imperfections is demonstrated. PNS was not a practical limitation.

CONCLUSION

The TR-optimal framework designs RF pulses for a range of pulse parameters, specifically to minimize sequence TR.

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.

Cardiac cine magnetic resonance fingerprinting for combined ejection fraction, T1 and T2 quantification

This study introduces a technique called cine magnetic resonance fingerprinting (cine-MRF) for simultaneous T₁ , T₂ and ejection fraction (EF) quantification. Data acquired with a free-running MRF sequence are retrospectively sorted into different cardiac phases using an external electrocardiogram (ECG) signal. A low-rank reconstruction with a finite difference sparsity constraint along the cardiac motion dimension yields images resolved by cardiac phase. To improve SNR and precision in the parameter maps, these images are nonrigidly registered to the same phase and matched to a dictionary to generate T₁ and T₂ maps. Cine images for computing left ventricular volumes and EF are also derived from the same data. Cine-MRF was tested in simulations using a numerical relaxation phantom. Phantom and in vivo scans of 19 subjects were performed at 3 T during a 10.9 seconds breath-hold with an in-plane resolution of 1.6 x 1.6 mm² and 24 cardiac phases. Left ventricular EF values obtained with cine-MRF agreed with the conventional cine images (mean bias -1.0%). Average myocardial T₁ times in diastole/systole were 1398/1391 ms with cine-MRF, 1394/1378 ms with ECG-triggered cardiac MRF (cMRF) and 1234/1212 ms with MOLLI; and T₂ values were 30.7/30.3 ms with cine-MRF, 32.6/32.9 ms with ECG-triggered cMRF and 37.6/41.0 ms with T₂ -prepared FLASH. Cine-MRF and ECG-triggered cMRF relaxation times were in good agreement. Cine-MRF T₁ values were significantly longer than MOLLI, and cine-MRF T₂ values were significantly shorter than T₂ -prepared FLASH. In summary, cine-MRF can potentially streamline cardiac MRI exams by combining left ventricle functional assessment and T₁ -T₂ mapping into one time-efficient acquisition.

Autocalibrated cardiac tissue phase mapping with multiband imaging and k-t acceleration

PURPOSE

To develop an autocalibrated multiband (MB) CAIPIRINHA acquisition scheme with in-plane k-t acceleration enabling multislice three-directional tissue phase mapping in one breath-hold.

METHODS

A k-t undersampling scheme was integrated into a time-resolved electrocardiographic-triggered autocalibrated MB gradient-echo sequence. The sequence was used to acquire data on 4 healthy volunteers with MB factors of two (MB2) and three (MB3), which were reconstructed using a joint reconstruction algorithm that tackles both k-t and MB acceleration. Forward simulations of the imaging process were used to tune the reconstruction model hyperparameters. Direct comparisons between MB and single-band tissue phase-mapping measurements were performed.

RESULTS

Simulations showed that the velocities could be accurately reproduced with MB2 k-t (average ± twice the SD of the RMS error of 0.08 ± 0.22 cm/s and velocity peak reduction of 1.03% ± 6.47% compared with fully sampled velocities), whereas acceptable results were obtained with MB3 k-t (RMS error of 0.13 ± 0.58 cm/s and peak reduction of 2.21% ± 13.45%). When applied to tissue phase-mapping data, the proposed technique allowed three-directional velocity encoding to be simultaneously acquired at two/three slices in a single breath-hold of 18 heartbeats. No statistically significant differences were detected between MB2/MB3 k-t and single-band k-t motion traces averaged over the myocardium. Regional differences were found, however, when using the American Heart Association model for segmentation.

CONCLUSION

An autocalibrated MB k-t acquisition/reconstruction framework is presented that allows three-directional velocity encoding of the myocardial velocities at multiple slices in one breath-hold.

Faster 3D saturation-recovery based myocardial T1 mapping using a reduced number of saturation points and denoising

PURPOSE

To accelerate the acquisition of free-breathing 3D saturation-recovery-based (SASHA) myocardial T1 mapping by acquiring fewer saturation points in combination with a post-processing 3D denoising technique to maintain high accuracy and precision.

METHODS

3D SASHA T1 mapping acquires nine T1-weighted images along the saturation recovery curve, resulting in long acquisition times. In this work, we propose to accelerate conventional cardiac T1 mapping by reducing the number of saturation points. High T1 accuracy and low standard deviation (as a surrogate for precision) is maintained by applying a 3D denoising technique to the T1-weighted images prior to pixel-wise T1 fitting. The proposed approach was evaluated on a T1 phantom and 20 healthy subjects, by varying the number of T1-weighted images acquired between three and nine, both prospectively and retrospectively. Following the results from the healthy subjects, three patients with suspected cardiovascular disease were acquired using five T1-weighted images. T1 accuracy and precision was determined for all the acquisitions before and after denoising.

RESULTS

In the T1 phantom, no statistical difference was found in terms of accuracy and precision for the different number of T1-weighted images before or after denoising (P = 0.99 and P = 0.99 for accuracy, P = 0.64 and P = 0.42 for precision, respectively). In vivo, both prospectively and retrospectively, the precision improved considerably with the number of T1-weighted images employed before denoising (P<0.05) but was independent on the number of T1-weighted images after denoising.

CONCLUSION

We demonstrate the feasibility of accelerating 3D SASHA T1 mapping by reducing the number of acquired T1-weighted images in combination with an efficient 3D denoising, without affecting accuracy and precision of T1 values.

2D multislice MP2RAGE sequence for fast T1 mapping at 7 T: Application to mouse imaging and MR thermometry

PURPOSE

To develop a 2D radial multislice MP2RAGE sequence for fast and reliable T₁ mapping at 7 T in mice and for MR thermometry.

METHODS

The 2D-MP2RAGE sequence was performed with the following parameters: TI₁ -TI₂ -MP2RAGE_TR = 1000-3000-9000 ms. The multiple dead times within the sequence were used for interleaved multislice acquisition, enabling one to acquire six slices in 9 seconds. The excitation pulse shape, inversion selectivity, and interslice gap were optimized. In vitro comparison with the inversion-recovery sequence was performed. The T₁ variations with temperature were measured on tubes with T₁ ranging from 800 ms to 2000 ms. The sequence was used to acquire T₁ maps continuously during 30 minutes on the brain and abdomen of healthy mice.

RESULTS

A three-lobe cardinal sine excitation pulse, combined with an inversion slice thickness and an interslice gap of respectively 150% and 50% of the imaging slice thickness, led to a SD and bias of the T₁ measurements below 1% and 2%, respectively. A linear dependence of T₁ with temperature was measured between 10°C and 60°C. In vivo, less than 1% variation was measured between successive T₁ maps in the mouse brain. In the abdomen, no obvious in-plane motion artifacts were observed but respiratory motion in the slice dimension led to 6% T₁ underestimation.

CONCLUSION

The multislice MP2RAGE sequence could be used for fast whole-body T₁ mapping and MR thermometry. Its reconstruction method would enable on-the-fly reconstruction.

Free-Breathing Three-Dimensional T1 Mapping of the Heart Using Subspace-Based Data Acquisition and Image Reconstruction

Mapping the longitudinal relaxation time constant (T1) of the myocardium using Magnetic Resonance Imaging (MRI) is an emerging technique for quantitative assessment of the morphology and viability of the myocardium. However, three-dimensional (3D) T1 mapping of the heart is challenging due to the high dimensionality of the signal and the presence of cardiac and respiratory motions. We propose a subspace-based method for free-breathing 3D T1 mapping of the heart without respiratory gating. The image function is represented as a high-order partially separable (PS) function to explore the inherent spatiotemporal correlations of the underlying signal. A special data acquisition scheme enabled by the high-order PS model is used for sparse sampling of the (k,t)-space, where complementary sparse datasets are acquired, each covering only a small portion of the (k,t)-space to characterize a single subspace (spatial or temporal). High-resolution dynamic MR images are reconstructed from the highly undersampled (k,t)-space using low-rank tensor and sparsity constraints. We demonstrate the feasibility of our proposed method using in vivo data obtained from healthy subjects on a 3T MR scanner. The proposed method can enable new clinical applications of T1 mapping in cardiac MR.

Improved Regularized Reconstruction for Simultaneous Multi-Slice Cardiac MRI T 1 Mapping

Myocardial T ₁ mapping is a quantitative MRI technique that has found great clinical utility in the detection of various heart disease. These acquisitions typically require three breath-holds, leading to long scan durations and patient discomfort. Simultaneous multi-slice (SMS) imaging has been shown to reduce the scan time of myocardial T ₁ mapping to a single breath-hold without sacrificing coverage, albeit at reduced precision. In this work, we propose a new reconstruction strategy for SMS imaging that combines the advantages of two different k-space interpolation strategies, while allowing for regularization, in order to improve the precision of accelerated mycordial T ₁ mapping.

Non-Cartesian slice-GRAPPA and slice-SPIRiT reconstruction methods for multiband spiral cardiac MRI

PURPOSE

Spiral MRI has advantages for cardiac imaging, and multiband (MB) spiral MRI of the heart shows promise. However, current reconstruction methods for MB spiral imaging have limitations. We sought to develop improved reconstruction methods for MB spiral cardiac MRI.

METHODS

Two reconstruction methods were developed. The first is non-Cartesian slice-GRAPPA (NCSG), which uses phase demodulation and gridding operations before application of a Cartesian slice-separating kernel. The second method, slice-SPIRiT, formulates the reconstruction as a minimization problem that enforces in-plane coil consistency and consistency with the acquired MB data, and uses through-plane coil sensitivity information in the iterative solution. These methods were compared with conjugate-gradient SENSE in phantoms and volunteers. Temporal alternation of CAIPIRINHA (controlled aliasing in parallel imaging results in higher acceleration) phase and the use of a temporal filter were also investigated.

RESULTS

Phantom experiments with 3 simultaneous slices (MB = 3) showed that mean artifact power was highest for conjugate-gradient SENSE, lower for NCSG, and lowest for slice-SPIRiT. For volunteer cine imaging (MB = 3, N = 5), the artifact power was 0.182 ± 0.037, 0.148 ± 0.036, and 0.139 ± 0.034 for conjugate-gradient SENSE, NCSG, and slice-SPIRiT, respectively (P < .05, analysis of variance). Temporal alternation of CAIPIRINHA reduced artifacts for both NCSG and slice-SPIRiT.

CONCLUSION

The NCSG and slice-SPIRiT methods provide more accurate reconstructions for MB spiral cine imaging compared with conjugate-gradient SENSE. These methods hold promise for non-Cartesian MB imaging.

Model-based myocardial T1 mapping with sparsity constraints using single-shot inversion-recovery radial FLASH cardiovascular magnetic resonance

BACKGROUND

This study develops a model-based myocardial T1 mapping technique with sparsity constraints which employs a single-shot inversion-recovery (IR) radial fast low angle shot (FLASH) cardiovascular magnetic resonance (CMR) acquisition. The method should offer high resolution, accuracy, precision and reproducibility.

METHODS

The proposed reconstruction estimates myocardial parameter maps directly from undersampled k-space which is continuously measured by IR radial FLASH with a 4 s breathhold and retrospectively sorted based on a cardiac trigger signal. Joint sparsity constraints are imposed on the parameter maps to further improve T1 precision. Validations involved studies of an experimental phantom and 8 healthy adult subjects.

RESULTS

In comparison to an IR spin-echo reference method, phantom experiments with T1 values ranging from 300 to 1500 ms revealed good accuracy and precision at simulated heart rates between 40 and 100 bpm. In vivo T1 maps achieved better precision and qualitatively better preservation of image features for the proposed method than a real-time CMR approach followed by pixelwise fitting. Apart from good inter-observer reproducibility (0.6% of the mean), in vivo results confirmed good intra-subject reproducibility (1.05% of the mean for intra-scan and 1.17, 1.51% of the means for the two inter-scans, respectively) of the proposed method.

CONCLUSION

Model-based reconstructions with sparsity constraints allow for single-shot myocardial T1 maps with high spatial resolution, accuracy, precision and reproducibility within a 4 s breathhold. Clinical trials are warranted.

MR susceptibility contrast imaging using a 2D simultaneous multi-slice gradient-echo sequence at 7T

Purpose

To develop a 7T simultaneous multi-slice (SMS) 2D gradient-echo sequence for susceptibility contrast imaging, and to compare its quality to 3D imaging.

Methods

A frequency modulated and phase cycled RF pulse was designed to simultaneously excite multiple slices in multi-echo 2D gradient-echo imaging. The imaging parameters were chosen to generate images with susceptibility contrast, including T₂*-weighted magnitude/phase images, susceptibility-weighted images and quantitative susceptibility/R₂* maps. To compare their image quality with 3D gradient-echo imaging, both 2D and 3D imaging were performed on 11 healthy volunteers and 4 patients with multiple sclerosis (MS). The signal to noise ratio (SNR) in gray and white matter and their contrast to noise ratio (CNR) was simulated for the 2D and 3D magnitude images using parameters from the imaging. The experimental SNRs and CNRs were measured in gray/white matter and deep gray matter structures on magnitude, phase, R₂* and QSM images from volunteers and the visibility of MS lesions on these images from patients was visually rated. All SNRs and CNRs were compared between the 2D and 3D imaging using a paired t-test.

Results

Although the 3D magnitude images still had significantly higher SNRs (by 13.0~17.6%), the 2D magnitude and QSM images generated significantly higher gray/white matter or globus pallidus/putamen contrast (by 13.3~87.5%) and significantly higher MS lesion contrast (by 5.9~17.3%).

Conclusion

2D SMS gradient-echo imaging can serve as an alternative to often used 3D imaging to obtain susceptibility-contrast-weighted images, with an advantage of providing better image contrast and MS lesion sensitivity.

Whole-heart spiral simultaneous multi-slice first-pass myocardial perfusion imaging

PURPOSE

To develop and evaluate a simultaneous multislice (SMS) spiral perfusion pulse sequence with whole-heart coverage.

METHODS

An orthogonal set of phase cycling angles following a Hadamard pattern was incorporated into a golden-angle (GA) variable density spiral perfusion sequence to perform SMS imaging at different multiband (MB) factors. Images were reconstructed using an SMS extension of L1-SPIRiT that we have termed SMS-L1-SPIRiT. The proposed sequence was evaluated in 40 subjects (10 each for MB factors of 1, 2, 3, and 4). Images were blindly graded by 2 cardiologists on a 5-point scale (5, excellent). To quantitatively evaluate the reconstruction performance against images acquired without SMS, the MB =1 data were used to retrospectively simulate data acquired at MB factors of 2 to 4.

RESULTS

Analysis of the SMS point-spread function for the desired slice showed that the proposed sampling strategy significantly canceled the main-lobe energy of the other slices and has low side-lobe energy resulting in an incoherent temporal aliasing pattern when rotated by the GA. Retrospective experiments demonstrated the SMS-L1-SPIRiT method removed aliasing from the interfering slices and showed excellent agreement with the ground-truth MB =1 images. Clinical evaluation demonstrated high-quality perfusion images with average image-quality scores of 4.3 ± 0.5 (MB =2), 4.2 ± 0.5 (MB =3), and 4.4 ± 0.4 (MB =4) with no significant quality difference in image quality between MB factors (P = 0.38).

CONCLUSION

SMS spiral perfusion at MB factors 2, 3, and 4 produces high-quality perfusion images with whole-heart coverage in a clinical setting with high sampling efficiency.

Fast Interleaved Multislice T1 Mapping: Model-Based Reconstruction of Single-Shot Inversion-Recovery Radial FLASH

Purpose

To develop a high-speed multislice T1 mapping method based on a single-shot inversion-recovery (IR) radial FLASH acquisition and a regularized model-based reconstruction.

Methods

Multislice radial k-space data are continuously acquired after a single nonselective inversion pulse using a golden-angle sampling scheme in a spoke-interleaved manner with optimized flip angles. Parameter maps and coil sensitivities of each slice are estimated directly from highly undersampled radial k-space data using a model-based nonlinear inverse reconstruction in conjunction with joint sparsity constraints. The performance of the method has been validated using a numerical and experimental T1 phantom as well as demonstrated for studies of the human brain and liver at 3T.

Results

The proposed method allows for 7 simultaneous T1 maps of the brain at 0.5 × 0.5 × 4 mm³ resolution within a single IR experiment of 4 s duration. Phantom studies confirm similar accuracy and precision as obtained for a single-slice acquisition. For abdominal applications, the proposed method yields three simultaneous T1 maps at 1.25 × 1.25 × 6 mm³ resolution within a 4 s breath hold.

Conclusion

Rapid, robust, accurate, and precise multislice T1 mapping may be achieved by combining the advantages of a model-based nonlinear inverse reconstruction, radial sampling, parallel imaging, and compressed sensing.

Time efficient whole-brain coverage with MR Fingerprinting using slice-interleaved echo-planar-imaging

Magnetic resonance fingerprinting (MRF) is a promising method for fast simultaneous quantification of multiple tissue parameters. The objective of this study is to improve the coverage of MRF based on echo-planar imaging (MRF-EPI) by using a slice-interleaved acquisition scheme. For this, the MRF-EPI is modified to acquire several slices in a randomized interleaved manner, increasing the effective repetition time of the spoiled gradient echo readout acquisition in each slice. Per-slice matching of the signal-trace to a precomputed dictionary allows the generation of T1 and T2* maps with integrated B1+ correction. Subsequent compensation for the coil sensitivity profile and normalization to the cerebrospinal fluid additionally allows for quantitative proton density (PD) mapping. Numerical simulations are performed to optimize the number of interleaved slices. Quantification accuracy is validated in phantom scans and feasibility is demonstrated in-vivo. Numerical simulations suggest the acquisition of four slices as a trade-off between quantification precision and scan-time. Phantom results indicate good agreement with reference measurements (Difference T1: -2.4 ± 1.1%, T2*: -0.5 ± 2.5%, PD: -0.5 ± 7.2%). In-vivo whole-brain coverage of T1, T2* and PD with 32 slices was acquired within 3:36 minutes, resulting in parameter maps of high visual quality and comparable performance with single-slice MRF-EPI at 4-fold scan-time reduction.

Saturation-Recovery Myocardial T1-Mapping during Systole: Accurate and Robust Quantification in the Presence of Arrhythmia

Myocardial T₁-mapping, a cardiac magnetic resonance imaging technique, facilitates a quantitative measure of fibrosis which is linked to numerous cardiovascular symptoms. To overcome the problems of common techniques, including lack of accuracy and robustness against partial-voluming and heart-rate variability, we introduce a systolic saturation-recovery T₁-mapping method. The Saturation-Pulse Prepared Heart-rate independent Inversion-Recovery (SAPPHIRE) T₁-mapping method was modified to enable imaging during systole. Phantom measurements were used to evaluate the insensitivity of systolic T₁-mapping towards heart-rate variability. In-vivo feasibility and accuracy were demonstrated in ten healthy volunteers with native and post-contrast T₁-mappping during systole and diastole. To show benefits in the presence of RR-variability, six arrhythmic patients underwent native T₁-mapping. Resulting systolic SAPPHIRE T₁-values showed no dependence on arrhythmia in phantom (CoV < 1%). In-vivo, significantly lower T₁ (1563 ± 56 ms, precision: 84.8 ms) and ECV-values (0.20 ± 0.03) than during diastole (T₁ = 1580 ± 62 ms, p = 0.0124; precision: 60.2 ms, p = 0.03; ECV = 0.21 ± 0.03, p = 0.0098) were measured, with a strong correlation of systolic and diastolic T₁ (r = 0.89). In patients, mis-triggering-induced motion caused significant imaging artifacts in diastolic T₁-maps, whereas systolic T₁-maps displayed resilience to arrythmia. In conclusion, the proposed method enables saturation-recovery T₁-mapping during systole, providing increased robustness against partial-voluming compared to diastolic imaging, for the benefit of T₁-measurements in arrhythmic patients.

Temporally resolved parametric assessment of Z-magnetization recovery (TOPAZ): Dynamic myocardial T1 mapping using a cine steady-state look-locker approach

PURPOSE

To develop and evaluate a cardiac phase-resolved myocardial T₁ mapping sequence.

METHODS

The proposed method for temporally resolved parametric assessment of Z-magnetization recovery (TOPAZ) is based on contiguous fast low-angle shot imaging readout after magnetization inversion from the pulsed steady state. Thereby, segmented k-space data are acquired over multiple heartbeats, before reaching steady state. This results in sampling of the inversion-recovery curve for each heart phase at multiple points separated by an R-R interval. Joint T₁ and B1+ estimation is performed for reconstruction of cardiac phase-resolved T₁ and B1+ maps. Sequence parameters are optimized using numerical simulations. Phantom and in vivo imaging are performed to compare the proposed sequence to a spin-echo reference and saturation pulse prepared heart rate-independent inversion-recovery (SAPPHIRE) T₁ mapping sequence in terms of accuracy and precision.

RESULTS

In phantom, TOPAZ T₁ values with integrated B1+ correction are in good agreement with spin-echo T₁ values (normalized root mean square error = 4.2%) and consistent across the cardiac cycle (coefficient of variation = 1.4 ± 0.78%) and different heart rates (coefficient of variation = 1.2 ± 1.9%). In vivo imaging shows no significant difference in TOPAZ T₁ times between the cardiac phases (analysis of variance: P = 0.14, coefficient of variation = 3.2 ± 0.8%), but underestimation compared with SAPPHIRE (T₁ time ± precision: 1431 ± 56 ms versus 1569 ± 65 ms). In vivo precision is comparable to SAPPHIRE T₁ mapping until middiastole (P > 0.07), but deteriorates in the later phases.

CONCLUSIONS

The proposed sequence allows cardiac phase-resolved T₁ mapping with integrated B1+ assessment at a temporal resolution of 40 ms. Magn Reson Med 79:2087-2100, 2018. © 2017 International Society for Magnetic Resonance in Medicine.