Improving parallel imaging by jointly reconstructing multi-contrast data

PURPOSE

To develop parallel imaging techniques that simultaneously exploit coil sensitivity encoding, image phase prior information, similarities across multiple images, and complementary k-space sampling for highly accelerated data acquisition.

METHODS

We introduce joint virtual coil (JVC)-generalized autocalibrating partially parallel acquisitions (GRAPPA) to jointly reconstruct data acquired with different contrast preparations, and show its application in 2D, 3D, and simultaneous multi-slice (SMS) acquisitions. We extend the joint parallel imaging concept to exploit limited support and smooth phase constraints through Joint (J-) LORAKS formulation. J-LORAKS allows joint parallel imaging from limited autocalibration signal region, as well as permitting partial Fourier sampling and calibrationless reconstruction.

RESULTS

We demonstrate highly accelerated 2D balanced steady-state free precession with phase cycling, SMS multi-echo spin echo, 3D multi-echo magnetization-prepared rapid gradient echo, and multi-echo gradient recalled echo acquisitions in vivo. Compared to conventional GRAPPA, proposed joint acquisition/reconstruction techniques provide more than 2-fold reduction in reconstruction error.

CONCLUSION

JVC-GRAPPA takes advantage of additional spatial encoding from phase information and image similarity, and employs different sampling patterns across acquisitions. J-LORAKS achieves a more parsimonious low-rank representation of local k-space by considering multiple images as additional coils. Both approaches provide dramatic improvement in artifact and noise mitigation over conventional single-contrast parallel imaging reconstruction. Magn Reson Med 80:619-632, 2018. © 2018 International Society for Magnetic Resonance in Medicine.

High-resolution whole-brain diffusion MRI at 7T using radiofrequency parallel transmission

PURPOSE

Investigating the utility of RF parallel transmission (pTx) for Human Connectome Project (HCP)-style whole-brain diffusion MRI (dMRI) data at 7 Tesla (7T).

METHODS

Healthy subjects were scanned in pTx and single-transmit (1Tx) modes. Multiband (MB), single-spoke pTx pulses were designed to image sagittal slices. HCP-style dMRI data (i.e., 1.05-mm resolutions, MB2, b-values = 1000/2000 s/mm² , 286 images and 40-min scan) and data with higher accelerations (MB3 and MB4) were acquired with pTx.

RESULTS

pTx significantly improved flip-angle detected signal uniformity across the brain, yielding ∼19% increase in temporal SNR (tSNR) averaged over the brain relative to 1Tx. This allowed significantly enhanced estimation of multiple fiber orientations (with ∼21% decrease in dispersion) in HCP-style 7T dMRI datasets. Additionally, pTx pulses achieved substantially lower power deposition, permitting higher accelerations, enabling collection of the same data in 2/3 and 1/2 the scan time or of more data in the same scan time.

CONCLUSION

pTx provides a solution to two major limitations for slice-accelerated high-resolution whole-brain dMRI at 7T; it improves flip-angle uniformity, and enables higher slice acceleration relative to current state-of-the-art. As such, pTx provides significant advantages for rapid acquisition of high-quality, high-resolution truly whole-brain dMRI data.

Functional Sensitivity of 2D Simultaneous Multi-Slice Echo-Planar Imaging: Effects of Acceleration on g-factor and Physiological Noise

Accelerated data acquisition with simultaneous multi-slice (SMS) imaging for functional MRI studies leads to interacting and opposing effects that influence the sensitivity to blood oxygen level-dependent (BOLD) signal changes. Image signal to noise ratio (SNR) is decreased with higher SMS acceleration factors and shorter repetition times (TR) due to g-factor noise penalties and saturation of longitudinal magnetization. However, the lower image SNR is counteracted by greater statistical power from more samples per unit time and a higher temporal Nyquist frequency that allows for better removal of spurious non-BOLD high frequency signal content. This study investigated the dependence of the BOLD sensitivity on these main driving factors and their interaction, and provides a framework for evaluating optimal acceleration of SMS-EPI sequences. functional magnetic resonance imaging (fMRI) data from a scenes/objects visualization task was acquired in 10 healthy volunteers at a standard neuroscience resolution of 3 mm on a 3T MRI scanner. SMS factors 1, 2, 4, and 8 were used, spanning TRs of 2800 ms to 350 ms. Two data processing methods were used to equalize the number of samples over the SMS factors. BOLD sensitivity was assessed using g-factors maps, temporal SNR (tSNR), and t-score metrics. tSNR results show a dependence on SMS factor that is highly non-uniform over the brain, with outcomes driven by g-factor noise amplification and the presence of high frequency noise. The t-score metrics also show a high degree of spatial dependence: the lower g-factor noise area of V1 shows significant improvements at higher SMS factors; the moderate-level g-factor noise area of the parahippocampal place area shows only a trend of improvement; and the high g-factor noise area of the ventral-medial pre-frontal cortex shows a trend of declining t-scores at higher SMS factors. This spatial variability suggests that the optimal SMS factor for fMRI studies is region dependent. For task fMRI studies done with similar parameters as were used here (3T scanner, 32-channel RF head coil, whole brain coverage at 3 mm isotropic resolution), we recommend SMS accelerations of 4x (conservative) to 8x (aggressive) for most studies and a more conservative acceleration of 2x for studies interested in anterior midline regions.

Optimization of simultaneous multislice EPI for concurrent functional perfusion and BOLD signal measurements at 7T

Purpose

To overcome limitations of previous ultra‐high‐field arterial spin labeling (ASL) techniques concerning temporal resolution and brain coverage by utilizing the simultaneous multi‐slice (SMS) approach.

Methods

An optimized, flow‐alternating inversion recovery quantitative imaging of perfusion using a single subtraction II scheme was developed that tackles the challenges of 7 tesla (T) ASL. The implementation of tailored labeling radiofrequency pulses reduced the effect of transmit field ( B1+) inhomogeneities. The proposed approach utilizes an SMS echo‐planar imaging (EPI) readout to efficiently achieve large brain coverage.

Results

A pulsed ASL (PASL) technique with large brain coverage is described and optimized that can be applied at temporal resolutions below 2.5 s, similar to those achievable at 1.5 and 3T magnetic field strength. The influences of within‐ and through‐slice acceleration factors and reconstruction parameters on perfusion and blood‐oxygenation‐level‐dependent (BOLD)‐signal image and temporal signal‐to‐noise ratio (SNR) are presented. The proposed approach yielded twice the brain coverage as compared to conventional PASL at 7T, without notable loss in image quality.

Conclusion

The presented SMS EPI PASL at 7T overcomes current limitations in SNR, temporal resolution, and spatial coverage for functional perfusion and BOLD signal as well as baseline perfusion measurements. Magn Reson Med 78:121–129, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution‐NonCommercial‐NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non‐commercial and no modifications or adaptations are made.

High spatial resolution spiral first-pass myocardial perfusion imaging with whole-heart coverage at 3 T

PURPOSE

To develop and evaluate a high spatial resolution (1.25 × 1.25 mm² ) spiral first-pass myocardial perfusion imaging technique with whole-heart coverage at 3T, to better assess transmural differences in perfusion between the endocardium and epicardium, to quantify the myocardial ischemic burden, and to improve the detection of obstructive coronary artery disease.

METHODS

Whole-heart high-resolution spiral perfusion pulse sequences and corresponding motion-compensated reconstruction techniques for both interleaved single-slice (SS) and simultaneous multi-slice (SMS) acquisition with or without outer-volume suppression (OVS) were developed. The proposed techniques were evaluated in 34 healthy volunteers and 8 patients (55 data sets). SS and SMS images were reconstructed using motion-compensated L1-SPIRiT and SMS-Slice-L1-SPIRiT, respectively. Images were blindly graded by 2 experienced cardiologists on a 5-point scale (5, excellent; 1, poor).

RESULTS

High-quality perfusion imaging was achieved for both SS and SMS acquisitions with or without OVS. The SS technique without OVS had the highest scores (4.5 [4, 5]), which were greater than scores for SS with OVS (3.5 [3.25, 3.75], P < .05), MB = 2 without OVS (3.75 [3.25, 4], P < .05), and MB = 2 with OVS (3.75 [2.75, 4], P < .05), but significantly higher than those for MB = 3 without OVS (4 [4, 4], P = .95). SMS image quality was improved using SMS-Slice-L1-SPIRiT as compared to SMS-L1-SPIRiT (P < .05 for both reviewers).

CONCLUSION

We demonstrated the successful implementation of whole-heart spiral perfusion imaging with high resolution at 3T. Good image quality was achieved, and the SS without OVS showed the best image quality. Evaluation in patients with expected ischemic heart disease is warranted.

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.

A SEmi-Adiabatic matched-phase spin echo (SEAMS) PINS pulse-pair for B1 -insensitive simultaneous multislice imaging

PURPOSE

Simultaneous multislice (SMS) imaging is a powerful technique that can reduce image acquisition time for anatomical, functional, and diffusion weighted magnetic resonance imaging. At higher magnetic fields, such as 7 Tesla, increased radiofrequency (RF) field inhomogeneity, power deposition, and changes in relaxation parameters make SMS spin echo imaging challenging. We designed an adiabatic 180° Power Independent of Number of Slices (PINS) pulse and a matched-phase 90° PINS pulse to generate a SEmi-Adiabatic Matched-phase Spin echo (SEAMS) PINS sequence to address these issues.

METHODS

We used the adiabatic Shinnar Le-Roux (SLR) algorithm to generate a 180° pulse. The SLR polynomials for the 180° pulse were then used to create a matched-phase 90° pulse. The pulses were sub-sampled to produce a SEAMS PINS pulse-pair and the performance of this pulse-pair was validated in phantoms and in vivo.

RESULTS

Simulations as well as phantom and in vivo results, demonstrate multislice capability and improved B1 -insensitivity of the SEAMS PINS pulse-pair when operating at RF amplitudes of up to 40% above adiabatic threshold.

CONCLUSION

The SEAMS PINS approach presented here achieves multislice spin echo profiles with improved B1 -insensitivity when compared with a conventional spin echo.

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.

High-Resolution Simultaneous Multi-Slice Accelerated Turbo Spin-Echo Musculoskeletal Imaging: A Head-to-Head Comparison With Routine Turbo Spin-Echo Imaging

Background/Aim: The turbo spin-echo (TSE) sequence is widely used for musculoskeletal (MSK) imaging; however, its acquisition speed is limited and can be easily affected by motion artifacts. We aimed to evaluate whether the use of a simultaneous multi-slice TSE (SMS-TSE) sequence can accelerate MSK imaging while maintaining image quality when compared with the routine TSE sequence.

Methods: We prospectively enrolled 71 patients [mean age, 37.43 ± 12.56 (range, 20–67) years], including 37 men and 34 women, to undergo TSE and SMS sequences. The total scanning times for the wrist, ankle and knee joint with routine sequence were 14.92, 13.97, and 13.48 min, respectively. For the SMS-TSE sequence, they were 7.52, 7.20, and 6.87 min. Quantitative parameters, including the signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR), were measured. Three experienced MSK imaging radiologists qualitatively evaluated the image quality of bone texture, cartilage, tendons, ligament, meniscus, and artifact using a 5-point evaluation system, and the diagnostic performance of the SMS-TSE sequences was evaluated.

Results: Compared with the routine TSE sequences, the scanning time was lower by 49.60, 48.46, and 49.04% using SMS-TSE sequences for the wrist, ankle, and knee joints, respectively. For the SNR comparison, the SMS-TSE sequences were significantly higher than the routine TSE sequence for wrist (except for Axial-T2WI-FS), ankle, and knee joint MR imaging (all p < 0.05), but no statistical significance was obtained for the CNR measurement (all p > 0.05, except for Sag-PDWI-FS in ankle joint). For the wrist joint, the diagnostic sensitivity, specificity, and accuracy were 88.24, 100, and 92%. For the ankle joint, they were 100, 75, and 93.33%. For the knee joint, they were 87.50, 85.71, and 87.10%.

Conclusion: The use of the SMS-TSE sequence in the wrist, ankle, and knee joints can significantly reduce the scanning time and show similar image quality when compared with the routine TSE sequence.

Simultaneous multi-slice accelerated 4D-MRI for radiotherapy guidance

4D-MRI is becoming increasingly important for daily guidance of thoracic and abdominal radiotherapy. This study exploits the simultaneous multi-slice (SMS) technique to accelerate the acquisition of a balanced turbo field echo (bTFE) and a turbo spin echo (TSE) coronal 4D-MRI sequence performed on 1.5 T MRI scanners. SMS single-shot bTFE and TSE sequences were developed to acquire a stack of 52 coronal 2D images over 30 dynamics. Simultaneously excited slices were separated by half the field of view. Slices intersecting with the liver-lung interface were used as navigator slices. For each navigator slice location, an end-exhale dynamic was automatically identified, and used to derive the self-sorting signal by rigidly registering the remaining dynamics. Navigator slices were sorted into 10 amplitude bins, and the temporal relationship of simultaneously excited slices was used to generate sorted 4D-MRIs for 12 healthy volunteers. The self-sorting signal was validated using anin vivopeak-to-peak motion analysis. The smoothness of the liver-lung interface was quantified by comparing to sagittal cine images acquired directly after the SMS-4D-MRI sequence. To ensure compatibility with the MR-linac radiotherapy workflow, the 4D-MRIs were transformed into 3D mid-position (MidP) images using deformable image registration. Consistency of the deformable vector fields was quantified in terms of the distance discordance metric (DDM) in the body. The SMS-4D-TSE sequence was additionally acquired for 3 lung cancer patients to investigate tumor visibility. SMS-4D-MRI acquisition and processing took approximately 7 min. 4D-MRI reconstruction was possible for 26 out of 27 acquired datasets. Missing data in the sorted 4D-MRIs varied from 4%-26% for the volunteers and varied from 8%-24% for the patients. Peak-to-peak (SD) amplitudes analysis agreed within 1.8 (1.1) mm and 0.9 (0.4) mm between the sorted 4D-MRIs and the self-sorting signals of the volunteers and patients, respectively. Liver-lung interface smoothness was found to be in the range of 0.6-3.1 mm for volunteers. The percentage of DDM values smaller than 2 mm was in the range of 85%-89% and 86%-92% for the volunteers and patients, respectively. Lung tumors were clearly visibility in the SMS-4D-TSE images and MidP images. Two fast SMS-accelerated 4D-MRI sequences were developed resulting in T₂/T₁or T₂weighted contrast. The SMS-4D-MRIs and derived 3D MidP-MRIs yielded anatomically plausible images and good tumor visibility. SMS-4D-MRI is therefore a strong candidate to be used for treatment simulation and daily guidance of thoracic and abdominal MR-guided radiotherapy.

High temporal resolution motion estimation using a self-navigated simultaneous multi-slice echo planar imaging acquisition

BACKGROUND

Subject motion is known to produce spurious covariance among time-series in functional connectivity that has been reported to induce distance-dependent spurious correlations.

PURPOSE

To present a feasibility study for applying the extended Kalman filter (EKF) framework for high temporal resolution motion correction of resting state functional MRI (rs-fMRI) series using each simultaneous multi-slice (SMS) echo planar imaging (EPI) shot as its own navigator.

STUDY TYPE

Prospective feasibility study.

POPULATION/SUBJECTS

Three human volunteers.

FIELD STRENGTH/SEQUENCE

3T GE DISCOVERY MR750 scanner using a 32-channel head coil. Simultaneous multi-slice rs-fMRI sequence with repetition time (TR)/echo time (TE) = 800/30 ms, and SMS factor 6.

ASSESSMENT

Motion estimates were computed using two techniques: a conventional rigid-body volume-wise registration; and a high-temporal resolution motion estimation rigid-body approach. The reference image was resampled using the estimates obtained from both approaches and the difference between these predicted volumes and the original moving series was summarized using the normalized mean squared error (NMSE).

STATISTICAL TESTS

Direct comparison of NMSE values.

RESULTS

High-temporal motion estimation was always superior to volume-wise motion estimation for the sample presented. For staged continuous rotations, the NMSE using high-temporal resolution motion estimates ranged between [0.130, 0.150] for the first volunteer (in-plane rotations), between [0.060, 0.068] for the second volunteer (in-plane rotations), and between [0.063, 0.080] for the third volunteer (through-plane rotations). These values went up to [0.384, 0.464]; [0.136, 0.179]; and [0.080, 0.096], respectively, when using volume-wise motion estimates.

DATA CONCLUSION

Accurate high-temporal rigid-body motion estimates can be obtained for rs-fMRI taking advantage of simultaneous multi-slice EPI sub-TR shots.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2018.

A simultaneous multi-slice T2 mapping framework based on overlapping-echo detachment planar imaging and deep learning reconstruction

PURPOSE

Quantitative MRI (qMRI) is of great importance to clinical medicine and scientific research. However, most qMRI techniques are time-consuming and sensitive to motion, especially when a large 3D volume is imaged. To accelerate the acquisition, a framework is proposed to realize reliable simultaneous multi-slice T₂ mapping.

METHODS

The simultaneous multi-slice T₂ mapping framework is based on overlapping-echo detachment (OLED) planar imaging (dubbed SMS-OLED). Multi-slice overlapping-echo signals were generated by multiple excitation pulses together with echo-shifting gradients. The signals were excited and acquired with a single-channel coil. U-Net was used to reconstruct T₂ maps from the acquired overlapping-echo image.

RESULTS

Single-shot double-slice and two-shot triple-slice SMS-OLED scan schemes were designed according to the framework for evaluation. Simulations, water phantom, and in vivo rat brain experiments were carried out. Overlapping-echo signals were acquired, and T₂ maps were reconstructed and compared with references. The results demonstrate the superior performance of our method.

CONCLUSION

Two slices of T₂ maps can be obtained in a single shot within hundreds of milliseconds. Higher quality multi-slice T₂ maps can be obtained via multiple shots. SMS-OLED provides a lower specific absorption rate scheme compared with conventional SMS methods with a coil with only a single receiver channel. The new method is of potential in dynamic qMRI and functional qMRI where temporal resolution is vital.

Motion-robust diffusion compartment imaging using simultaneous multi-slice acquisition

PURPOSE

To achieve motion-robust diffusion compartment imaging (DCI) in near continuously moving subjects based on simultaneous multi-slice, diffusion-weighted brain MRI.

METHODS

Simultaneous multi-slice (SMS) acquisition enables fast and dense sampling of k- and q-space. We propose to achieve motion-robust DCI via slice-level motion correction by exploiting the rigid coupling between simultaneously acquired slices. This coupling provides 3D coverage of the anatomy that substantially constraints the slice-to-volume alignment problem. This is incorporated into an explicit model of motion dynamics that handles continuous and large subject motion in robust DCI reconstruction.

RESULTS

We applied the proposed technique, called Motion Tracking based on Simultanous Multislice Registration (MT-SMR) to multi b-value SMS diffusion-weighted brain MRI of healthy volunteers and motion-corrupted scans of 20 pediatric subjects. Quantitative and qualitative evaluation based on fractional anisotropy in unidirectional fiber regions, and DCI in crossing-fiber regions show robust reconstruction in the presence of motion.

CONCLUSION

The proposed approach has the potential to extend routine use of SMS DCI in very challenging populations, such as young children, newborns, and non-cooperative patients.

RF peak power reduction in CAIPIRINHA excitation by interslice phase optimization

The purpose of this work was to show that the overall peak power of RF pulses for CAIPIRINHA excitation can be substantially reduced by applying interslice phase relaxation. The optimal phases are scan dependent and can be quickly calculated by the proposed method. The multi-band RF pulse design is implemented as the minimization of a linear objective function with quadratic constraints. The interslice phase is considered to be a variable for optimization. In the case of a phase cycling scheme (CAIPIRINHA), the peak power is considered over all pulses. The computation time (about 1 s) is compatible with online RF pulse design. It is shown that the optimal interslice phases depend on the CAIPIRINHA scheme used and that RF peak power is reduced when the CAIPIRINHA phase cycling is taken into account in the optimization. The proposed method is extremely fast and results in RF pulses with low peak power for CAIPIRINHA excitation. The MATLAB implementation is given in the appendix; it allows for online determination of scan-dependent phase parameters. Furthermore, the method can be easily extended to pTx shimming systems in the context of multi-slice excitations, and this possibility is included in the software.

Split-slice training and hyperparameter tuning of RAKI networks for simultaneous multi-slice reconstruction

PURPOSE

Simultaneous multi-slice acquisitions are essential for modern neuroimaging research, enabling high temporal resolution functional and high-resolution q-space sampling diffusion acquisitions. Recently, deep learning reconstruction techniques have been introduced for unaliasing these accelerated acquisitions, and robust artificial-neural-networks for k-space interpolation (RAKI) have shown promising capabilities. This study systematically examines the impacts of hyperparameter selections for RAKI networks, and introduces a novel technique for training data generation which is analogous to the split-slice formalism used in slice-GRAPPA.

METHODS

RAKI networks were developed with variable hyperparameters and with and without split-slice training data generation. Each network was trained and applied to five different datasets including acquisitions harmonized with Human Connectome Project lifespan protocol. Unaliasing performance was assessed through L1 errors computed between unaliased and calibration frequency-space data.

RESULTS

Split-slice training significantly improved network performance in nearly all hyperparameter configurations. Best unaliasing results were achieved with three layer RAKI networks using at least 64 convolutional filters with receptive fields of 7 voxels, 128 single-voxel filters in the penultimate RAKI layer, batch normalization, and no training dropout with the split-slice augmented training dataset. Networks trained without the split-slice technique showed symptoms of network over-fitting.

CONCLUSIONS

Split-slice training for simultaneous multi-slice RAKI networks positively impacts network performance. Hyperparameter tuning of such reconstruction networks can lead to further improvements in unaliasing performance.

Virtual slice concept for improved simultaneous multi-slice MRI employing an extended leakage constraint

PURPOSE

To develop a novel, simultaneous multi-slice (SMS) reconstruction that extends an inter-slice leakage constraint to intra-slice aliasing with a virtual slice concept for artifact reduction.

METHODS

Inter-slice leakage constraint has been used for SMS reconstruction that mitigates leakage artifacts from the adjacent slices. In this work, the leakage constraint is extended to more general framework that includes SMS and parallel MRI as special cases by viewing intra-slice aliasing artifacts from undersampling as virtual slices while imposing data fidelity to ensure the measurement consistency. In this way, the reconstruction makes it feasible to directly estimate the individual slices from the undersampled SMS acquisition as a one-step method. The performance of the extended method is evaluated with data acquired using 2D GRE and EPI sequences.

RESULTS

Compared to a two-step method that performs slice unaliasing followed by inplane unaliasing, the proposed one-step method reduces aliasing artifacts by employing the extended leakage constraint while lowering the noise amplification by improving the conditioning for the inverse problem.

CONCLUSIONS

The proposed one-step method takes advantage of virtual slices as additional encoding power for improved image quality. We successfully demonstrated that the proposed one-step method minimizes a trade-off between aliasing artifacts and amplified noises over the two-step method.

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.

Volumetric assessment of spontaneous mechanical activities by simultaneous multi-slice MRI techniques with correlation to muscle fiber orientation

The purpose of this work was assessment of volumetric characteristics of spontaneous mechanical activities in musculature (SMAMs) by diffusion-weighted simultaneous multi-slice (DW-SMS) imaging and spatial correlation to anatomical structure, as revealed by fusion to fiber tractographic information derived from diffusion-tensor imaging (DTI). The feasibility of using DW-SMS to image spontaneous events in human musculature was assessed by phantom measurements. Series of DW-SMS images and DTI datasets were recorded from the resting calf of three human subjects. Simultaneously recorded SMAMs in multiple slices were analyzed regarding spatial extension by the Kolmogorov-Smirnov test. Direct correlation of spatial distribution of SMAMs and fiber orientation was investigated by mapping of muscle fibers to multi-slice SMAM datasets. The DW-SMS strategy allows simultaneous assessment of SMAMs in several slices of resting skeletal musculature, since 73.9% of SMAM-affected volumes have shown SMAMs in multiple DW-SMS slices. Spatial extension of SMAMs was highly correlated over different simultaneously recorded DW-SMS slices, and affected areas followed the orientation of muscle fibers with a connectivity ratio up to 57.18 ± 14.80% based on event count and connectivity count maps. In 89.2% of all SMAM-affected datasets muscle fiber connectivity was shown in at least two adjacent slices. Direct correlation between SMAMs in human lower leg musculature and underlying anatomical structure was revealed by high muscle fiber connectivity (89.2%). SMAMs have shown a wide distribution along the longitudinal muscle direction (73.9% in multiple DW-SMS slices) with direct involvement of muscle fibers. Correlation between SMAMs in multiple DW-SMS slices and crossing muscular fiber tracts provides evidence that SMAMs result from physiological processes in musculature. Fusion of DW-SMS with DTI facilitates non-invasive studies of muscle fiber involvement in SMAMs in resting muscle.

Open-source myocardial T1 mapping with simultaneous multi-slice acceleration: Combining an auto-calibrated blipped-bSSFP readout with VERSE-MB pulses

PURPOSE

Enabling fast and accessible myocardial T₁ mapping is crucial for extending its clinical application. We introduce Open-MOLLI-SMS combining simultaneous multi-slice (SMS) with auto-calibration and variable-rate selective excitation (VERSE)-multiband pulses to obtain all slices in a fast single-shot T₁ mapping sequence.

METHODS

Open-MOLLI-SMS was developed by integrating SMS with the open-source method Open-MOLLI previously implemented in Pulseq. Three methods were integrated for Open-MOLLI-SMS: (1) auto-calibration blip patterns to ensure consistency between the data and coil information; (2) a blipped-balanced SSFP (bSSFP) readout to induce controlled aliasing in parallel imaging shifts without disturbing the bSSFP frequency response; and (3) a VERSE-multiband pulse for minimizing the achievable TR and the specific absortion rate (SAR) impact of SMS. Two (SMS2) or three (SMS3) slices were excited simultaneously and encoded with an in-plane acceleration factor of 2. Experiments were performed in the International Society for Magnetic Resonance in Medicine/National Institute of Standards and Technology phantom and five healthy volunteers.

RESULTS

Phantom results show accurate T₁ estimates for reference values between 400 to 2200 ms. Artifacts were visible for Open-MOLLI-SMS3 but not replicated in vivo. In vivo Open-MOLLI-SMS (T₁ ^SMS2  = 993 ± 10 ms; T₁ ^SMS3  = 1031 ± 17 ms) provided similar values to mean T₁ single-band Open-MOLLI estimates (T₁ ^Open-MOLLI  = 1005 ± 47 ms). Open-MOLLI-SMS2 provided the closest estimates to the reference.

CONCLUSION

This proof-of-principle implementation study demonstrates the feasibility of speeding up T₁ -mapping acquisitions and increasing coverage by combining auto-calibration strategies with a blipped-bSFFP readout and VERSE multiband RF excitation pulses. The proposed methodology was built on the Open-MOLLI mapping sequence, which provides a fast means for prototyping and enables open-source sharing of the method.

Abdominal Diffusion-Weighted MRI With Simultaneous Multi-Slice Acquisition: Agreement and Reproducibility of Apparent Diffusion Coefficients Measurements

BACKGROUND

Simultaneous multi-slice diffusion-weighted imaging (SMS-DWI) can shorten acquisition time in abdominal imaging.

PURPOSE

To investigate the agreement and reproducibility of apparent diffusion coefficient (ADC) from abdominal SMS-DWI acquired with different vendors and different breathing schemes.

STUDY TYPE

Prospective.

SUBJECTS

Twenty volunteers and 10 patients.

FIELD STRENGTH/SEQUENCE

3.0 T, SMS-DWI with a diffusion-weighted echo-planar imaging sequence.

ASSESSMENT

SMS-DWI was acquired using breath-hold and free-breathing techniques in scanners from two vendors, yielding four scans in each participant. Average ADC values were measured in the liver, pancreas, spleen, and both kidneys. Non-normalized ADC and ADCs normalized to the spleen were compared between vendors and breathing schemes.

STATISTICAL TESTS

Paired t-test or Wilcoxon signed rank test; intraclass correlation coefficient (ICC); Bland-Altman method; coefficient of variation (CV) analysis; significance level: P < 0.05.

RESULTS

Non-normalized ADCs from the four SMS-DWI scans did not differ significantly in the spleen (P = 0.262, 0.330, 0.166, 0.122), right kidney (P = 0.167, 0.538, 0.957, 0.086), and left kidney (P = 0.182, 0.281, 0.504, 0.405), but there were significant differences in the liver and pancreas. For normalized ADCs, there were no significant differences in the liver (P = 0.315, 0.915, 0.198, 0.799), spleen (P = 0.815, 0.689, 0.347, 0.423), pancreas (P = 0.165, 0.336, 0.304, 0.584), right kidney (P = 0.165, 0.336, 0.304, 0.584), and left kidney (P = 0.496, 0.304, 0.443, 0.371). Inter-reader agreements of non-normalized ADCs were good to excellent (ICCs ranged from 0.861 to 0.983), and agreement and reproducibility were good to excellent depending on anatomic location (CVs ranged from 3.55% to 13.98%). Overall CVs for abdominal ADCs from the four scans were 6.25%, 7.62%, 7.08, and 7.60%.

DATA CONCLUSION

The normalized ADCs from abdominal SMS-DWI may be comparable between different vendors and breathing schemes, showing good agreement and reproducibility. ADC changes above approximately 8% may potentially be considered as a reliable quantitative biomarker to assess disease or treatment-related changes.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 2.

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.

3D motion strategy for online volumetric thermometry using simultaneous multi-slice EPI at 1.5T: an evaluation study

PURPOSE

In presence of respiratory motion, temperature mapping is altered by in-plane and through-plane displacements between successive acquisitions together with periodic phase variations. Fast 2D Echo Planar Imaging (EPI) sequence can accommodate intra-scan motion, but limited volume coverage and inter-scan motion remain a challenge during free-breathing acquisition since position offsets can arise between the different slices.

METHOD

To address this limitation, we evaluated a 2D simultaneous multi-slice EPI sequence with multiband (MB) acceleration during radiofrequency ablation on a mobile gel and in the liver of a volunteer (no heating). The sequence was evaluated in terms of resulting inter-scan motion, temperature uncertainty and elevation, potential false-positive heating and repeatability. Lastly, to account for potential through-plane motion, a 3D motion compensation pipeline was implemented and evaluated.

RESULTS

In-plane motion was compensated whatever the MB factor and temperature distribution was found in agreement during both the heating and cooling periods. No obvious false-positive temperature was observed under the conditions being investigated. Repeatability of measurements results in a 95% uncertainty below 2 °C for MB1 and MB2. Uncertainty up to 4.5 °C was reported with MB3 together with the presence of aliasing artifacts. Lastly, fast simultaneous multi-slice EPI combined with 3D motion compensation reduce residual out-of-plane motion.

CONCLUSION

Volumetric temperature imaging (12 slices/700 ms) could be performed with 2 °C accuracy or less, and offer tradeoffs in acquisition time or volume coverage. Such a strategy is expected to increase procedure safety by monitoring large volumes more rapidly for MR-guided thermotherapy on mobile organs.

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.

Accelerating joint relaxation-diffusion MRI by integrating time division multiplexing and simultaneous multi-slice (TDM-SMS) strategies

PURPOSE

To accelerate the acquisition of relaxation-diffusion imaging by integrating time-division multiplexing (TDM) with simultaneous multi-slice (SMS) for EPI and evaluate imaging quality and diffusion measures.

METHODS

The time-division multiplexing (TDM) technique and SMS method were integrated to achieve a high slice-acceleration (e.g., 6×) factor for acquiring relaxation-diffusion MRI. Two variants of the sequence, referred to as TDM3e-SMS and TDM2s-SMS, were developed to simultaneously acquire slice groups with three distinct TEs and two slice groups with the same TE, respectively. Both sequences were evaluated on a 3T scanner with in vivo human brains and compared with standard single-band (SB) -EPI and SMS-EPI using diffusion measures and tractography results.

RESULTS

Experimental results showed that the TDM3e-SMS sequence with total slice acceleration of 6 (multiplexing factor (MP) = 3 × multi-band factor (MB) = 2) provided similar image intensity and microstructure measures compared to standard SMS-EPI with MB = 2, and yielded less bias in intensity compared to standard SMS-EPI with MB = 4. The three sequences showed a similar positive correlation between TE and mean kurtosis (MK) and a negative correlation between TE and mean diffusivity (MD) in white matter. Multi-fiber tractography also shows consistency of results in TE-dependent measures between different sequences. The TDM2s-SMS sequence (MP = 2, MB = 2) also provided imaging measures similar to standard SMS-EPI sequences (MB = 2) for single-TE diffusion imaging.

CONCLUSIONS

The TDM-SMS sequence can provide additional 2× to 3× acceleration to SMS without degrading imaging quality. With the significant reduction in scan time, TDM-SMS makes joint relaxation-diffusion MRI a feasible technique in neuroimaging research to investigate new markers of brain disorders.

Novel Technique to Measure Pulse Wave Velocity in Brain Vessels Using a Fast Simultaneous Multi-Slice Excitation Magnetic Resonance Sequence

In this study, we proposed a novel pulse wave velocity (PWV) technique to determine cerebrovascular stiffness using a 3-tesla magnetic resonance imaging (MRI) to overcome the various shortcomings of existing PWV techniques for cerebral-artery PWV, such as long scan times and complicated procedures. The technique was developed by combining a simultaneous multi-slice (SMS) excitation pulse sequence with keyhole acquisition and reconstruction (SMS-K). The SMS-K technique for cerebral-artery PWV was evaluated using phantom and human experiments. In the results, common and internal carotid arteries (CCA and ICA) were acquired simultaneously in an image with a high temporal resolution-of 48 ms for one measurement. Vascular signals at 500 time points acquired within 30 s could generate pulse waveforms of CCA and ICA with 26 heartbeats, allowing for the detection of PWV changes over time. The results demonstrated that the SMS-K technique could provide more PWV information with a simple procedure within a short period of time. The procedural convenience and advantages of PWV measurements will make it more appropriate for clinical applications.