SMS MUSSELS: A navigator-free reconstruction for simultaneous multi-slice-accelerated multi-shot diffusion weighted imaging

Motion-compensated diffusion encoding in multi-shot human brain acquisitions: Insights using high-performance gradients

PURPOSE

To evaluate the utility of up to second-order motion-compensated diffusion encoding in multi-shot human brain acquisitions.

METHODS

Experiments were performed with high-performance gradients using three forms of diffusion encoding motion-compensated through different orders: conventional zeroth-order-compensated pulsed gradients (PG), first-order-compensated gradients (MC1), and second-order-compensated gradients (MC2). Single-shot acquisitions were conducted to correlate the order of motion compensation with resultant phase variability. Then, multi-shot acquisitions were performed at varying interleaving factors. Multi-shot images were reconstructed using three levels of shot-to-shot phase correction: no correction, channel-wise phase correction based on FID navigation, and correction based on explicit phase mapping (MUSE).

RESULTS

In single-shot acquisitions, MC2 diffusion encoding most effectively suppressed phase variability and sensitivity to brain pulsation, yielding residual variations of about 10° and of low spatial order. Consequently, multi-shot MC2 images were largely satisfactory without phase correction and consistently improved with the navigator correction, which yielded repeatable high-quality images; contrarily, PG and MC1 images were inadequately corrected using the navigator approach. With respect to MUSE reconstructions, the MC2 navigator-corrected images were in close agreement for a standard interleaving factor and considerably more reliable for higher interleaving factors, for which MUSE images were corrupted. Finally, owing to the advanced gradient hardware, the relative SNR penalty of motion-compensated diffusion sensitization was substantially more tolerable than that faced previously.

CONCLUSION

Second-order motion-compensated diffusion encoding mitigates and simplifies shot-to-shot phase variability in the human brain, rendering the multi-shot acquisition strategy an effective means to circumvent limitations of retrospective phase correction methods.

Multi-band multi-shot diffusion MRI reconstruction with joint usage of structured low-rank constraints and explicit phase mapping

PURPOSE

To develop a joint reconstruction method for multi-band multi-shot diffusion MRI.

THEORY AND METHODS

Multi-band multi-shot EPI acquisition is an effective approach for high-resolution diffusion MRI, but requires specific algorithms to correct the inter-shot phase variations. The phase correction can be done by first estimating the explicit phase map and then feeding it into the k-space signal formulation model. Alternatively, the phase information can be used indirectly as structured low-rank constraints in k-space. The 2 methods differ in reconstruction accuracy and efficiency. We aim to combine the 2 different approaches for improved image quality and reconstruction efficiency simultaneously, termed "joint usage of structured low-rank constraints and explicit phase mapping" (JULEP). The proposed JULEP reconstruction is tested on both single-band and multi-band, multi-shot diffusion data, with different resolutions and b values. The results of JULEP are compared with conventional methods with explicit phase mapping (i.e., multiplexed sensitivity-encoding [MUSE]) and structured low-rank constraints (i.e., MUSSELS), and another joint reconstruction method (i.e., network estimated artifacts for tempered reconstruction [NEATR]).

RESULTS

JULEP improves the quality of the navigator and subsequently facilitates the reconstruction of final diffusion images. Compared with all 3 other methods (MUSE, MUSSELS, and NEATR), JULEP mitigates residual structural bias and improves temporal SNRs in the final diffusion image, particularly at high multi-band factors. Compared with MUSSELS, JULEP also improves computational efficiency.

CONCLUSION

The proposed JULEP method significantly improves the image quality and reconstruction efficiency of multi-band multi-shot diffusion MRI, which can promote a broader application of high-resolution diffusion MRI.

Preliminary results of abdominal simultaneous multi-slice accelerated diffusion-weighted imaging with motion-correction in patients with cystic fibrosis and impaired compliance

OBJECTIVES

The aim of this prospective study was to compare scan time, image quality, signal-to-noise Ratio (SNR), and apparent diffusion coefficient (ADC) values of simultaneous multi-slice accelerated diffusion-weighted imaging with motion-correction (DWI SMS Moco) to standard diffusion-weighted imaging (sDWI) in free-breathing abdominal magnetic resonance imaging (MRI) in pediatric and young adult patients with cystic fibrosis (CF).

MATERIAL AND METHODS

16 patients (7 male and 9 female, 12-41 years old) with CF were examined prospectively in a single-center from November 2020 to March 2021 on a 1.5 Tesla clinical MR scanner. The characteristics of overall image quality and delimitability of mesenteric lymph nodes were evaluated using a 5-point Likert scale by two experienced pediatric radiologists independently from each other. Quantitative parameters with SNR and ADC values were assessed in 8 different locations and compared using a Wilcoxon signed-rank test.

RESULTS

The acquisition time for DWI SMS Moco was 32% shorter than for sDWI. Regarding quality comparison, overall image quality and delimitability of mesenteric lymph nodes were significant higher in DWI SMS Moco (p ≤ 0.05 for both readers). The readers preferred DWI SMS Moco to sDWI in all cases (16/16). Mean SNR values from DWI SMS Moco and sDWI were similar in 7 from 8 locations. The ADC values showed no significant difference between DWI SMS Moco and sDWI in any of the evaluated locations (p > 0.05).

CONCLUSIONS

The DWI SMS Moco improves overall image quality and delimitability of mesenteric lymph nodes compared to sDWI with similar SNR and ADC values and a distinguished reduction of scan time in free-breathing by one third. We conclude that MRI with DWI SMS Moco could be helpful in monitoring the effect of the high-efficiency modulator (HEM) therapy in cystic fibrosis (CF) patients homozygous or heterozygous for F508del in the abdomen.

Multi-band- and in-plane-accelerated diffusion MRI enabled by model-based deep learning in q-space and its extension to learning in the spherical harmonic domain

PURPOSE

To propose a new method for the recovery of combined in-plane- and multi-band (MB)-accelerated diffusion MRI data.

METHODS

Combining MB acceleration with in-plane acceleration is crucial to improve the time efficiency of high (angular and spatial) resolution diffusion scans. However, as the MB factor and in-plane acceleration increase, the reconstruction becomes challenging due to the heavy aliasing. The new reconstruction utilizes an additional q-space prior to constrain the recovery, which is derived from the previously proposed qModeL framework. Specifically, the qModeL prior provides a pre-learned representation of the diffusion signal space to which the measured data belongs. We show that the pre-learned q-space prior along with a model-based iterative reconstruction that accommodate multi-band unaliasing, can efficiently reconstruct the in-plane- and MB-accelerated data. The power of joint reconstruction is maximally utilized by using an incoherent under-sampling pattern in the k-q domain. We tested the proposed method on single- and multi-shell data, with high/low angular resolution, high/low spatial resolution, healthy/abnormal tissues, and 3T/7T field strengths. Furthermore, the learning is extended to the spherical harmonic basis, to provide a rotational invariant learning framework.

RESULTS

The qModeL joint reconstruction is shown to simultaneously unalias and jointly recover DWIs with reasonable accuracy in all the cases studied. The reconstruction error from 18-fold accelerated multi-shell datasets was <3%. The microstructural maps derived from the accelerated acquisitions also exhibit reasonable accuracy for both healthy and abnormal tissues. The deep learning (DL)-enabled reconstructions are comparable to those derived using traditional methods.

CONCLUSION

qModeL enables the joint recovery of combined in-plane- and MB-accelerated dMRI utilizing DL.

qModeL: A plug-and-play model-based reconstruction for highly accelerated multi-shot diffusion MRI using learned priors

PURPOSE

To introduce a joint reconstruction method for highly undersampled multi-shot diffusion weighted (msDW) scans.

METHODS

Multi-shot EPI methods enable higher spatial resolution for diffusion MRI, but at the expense of long scan-time. Highly accelerated msDW scans are needed to enable their utilization in advanced microstructure studies, which require high q-space coverage. Previously, joint k-q undersampling methods coupled with compressed sensing were shown to enable very high acceleration factors. However, the reconstruction of this data using sparsity priors is challenging and is not suited for multi-shell data. We propose a new reconstruction that recovers images from the combined k-q data jointly. The proposed qModeL reconstruction brings together the advantages of model-based iterative reconstruction and machine learning, extending the idea of plug-and-play algorithms. Specifically, qModeL works by prelearning the signal manifold corresponding to the diffusion measurement space using deep learning. The prelearned manifold prior is incorporated into a model-based reconstruction to provide a voxel-wise regularization along the q-dimension during the joint recovery. Notably, the learning does not require in vivo training data and is derived exclusively from biophysical modeling. Additionally, a plug-and-play total variation denoising provides regularization along the spatial dimension. The proposed framework is tested on k-q undersampled single-shell and multi-shell msDW acquisition at various acceleration factors.

RESULTS

The qModeL joint reconstruction is shown to recover DWIs from 8-fold accelerated msDW acquisitions with error less than 5% for both single-shell and multi-shell data. Advanced microstructural analysis performed using the undersampled reconstruction also report reasonable accuracy.

CONCLUSION

qModeL enables the joint recovery of highly accelerated multi-shot dMRI utilizing learning-based priors. The bio-physically driven approach enables the use of accelerated multi-shot imaging for multi-shell sampling and advanced microstructure studies.

Segmented simultaneous multi-slice diffusion-weighted imaging with navigated 3D rigid motion correction

PURPOSE

To improve the robustness of diffusion-weighted imaging (DWI) data acquired with segmented simultaneous multi-slice (SMS) echo-planar imaging (EPI) against in-plane and through-plane rigid motion.

THEORY AND METHODS

The proposed algorithm incorporates a 3D rigid motion correction and wavelet denoising into the image reconstruction of segmented SMS-EPI diffusion data. Low-resolution navigators are used to estimate shot-specific diffusion phase corruptions and 3D rigid motion parameters through SMS-to-volume registration. The shot-wise rigid motion and phase parameters are integrated into a SENSE-based full-volume reconstruction for each diffusion direction. The algorithm is compared to a navigated SMS reconstruction without gross motion correction in simulations and in vivo studies with four-fold interleaved 3-SMS diffusion tensor acquisitions.

RESULTS

Simulations demonstrate high fidelity was achieved in the SMS-to-volume registration, with submillimeter registration errors and improved image reconstruction quality. In vivo experiments validate successful artifact reduction in 3D motion-compromised in vivo scans with a temporal motion resolution of approximately 0.3 s.

CONCLUSION

This work demonstrates the feasibility of retrospective 3D rigid motion correction from shot navigators for segmented SMS DWI.

Magnetic resonance tractography of the lumbosacral plexus: Step-by-step

ABSTRACT

MR tractography of the lumbosacral plexus (LSP) is challenging due to the difficulty of acquiring high quality data and accurately estimating the neuronal tracts. We proposed an algorithm for an accurate visualization and assessment of the major LSP bundles using the segmentation of the cauda equina as seed points for the initial starting area for the fiber tracking algorithm.Twenty-six healthy volunteers underwent MRI examinations on a 3T MR scanner using the phased array coils with optimized measurement protocols for diffusion-weighted images and coronal T2 weighted 3D short-term inversion recovery sampling perfection with application optimized contrast using varying flip angle evaluation sequences used for LSP fiber reconstruction and MR neurography (MRN).The fiber bundles reconstruction was optimized in terms of eliminating the muscle fibers contamination using the segmentation of cauda equina, the effects of the normalized quantitative anisotropy (NQA) and angular threshold on reconstruction of the LSP. In this study, the NQA parameter has been used for fiber tracking instead of fractional anisotropy (FA) and the regions of interest positioning was precisely adjusted bilaterally and symmetrically in each individual subject.The diffusion data were processed in individual L3-S2 nerve fibers using the generalized Q-sampling imaging algorithm. Data (mean FA, mean diffusivity, axial diffusivity and radial diffusivity, and normalized quantitative anisotropy) were statistically analyzed using the linear mixed-effects model. The MR neurography was performed in MedINRIA and post-processed using the maximum intensity projection method to demonstrate LSP tracts in multiple planes.FA values significantly decreased towards the sacral region (P < .001); by contrast, mean diffusivity, axial diffusivity, radial diffusivity and NQA values significantly increased towards the sacral region (P < .001).Fiber tractography of the LSP was feasible in all examined subjects and closely corresponded with the nerves visible in the maximum intensity projection images of MR neurography. Usage of NQA instead of FA in the proposed algorithm enabled better separation of muscle and nerve fibers.The presented algorithm yields a high quality reconstruction of the LSP bundles that may be helpful both in research and clinical practice.

Joint calibrationless reconstruction of highly undersampled multicontrast MR datasets using a low-rank Hankel tensor completion framework

PURPOSE

To jointly reconstruct highly undersampled multicontrast two-dimensional (2D) datasets through a low-rank Hankel tensor completion framework.

METHODS

A multicontrast Hankel tensor completion (MC-HTC) framework is proposed to exploit the shareable information in multicontrast datasets with respect to their highly correlated image structure, common spatial support, and shared coil sensitivity for joint reconstruction. This is achieved by first organizing multicontrast k-space datasets into a single block-wise Hankel tensor. Subsequent low-rank tensor approximation via higher-order singular value decomposition (HOSVD) uses the image structural correlation by considering different contrasts as virtual channels. Meanwhile, the HOSVD imposes common spatial support and shared coil sensitivity by treating data from different contrasts as from additional k-space kernels. The missing k-space data are then recovered by iteratively performing such low-rank approximation and enforcing data consistency. This joint reconstruction framework was evaluated using multicontrast multichannel 2D human brain datasets (T₁ -weighted, T₂ -weighted, fluid-attenuated inversion recovery, and T₁ -weighted-inversion recovery) of identical image geometry with random and uniform undersampling schemes.

RESULTS

The proposed method offered high acceleration, exhibiting significantly less residual errors when compared with both single-contrast SAKE (simultaneous autocalibrating and k-space estimation) and multicontrast J-LORAKS (joint parallel-imaging-low-rank matrix modeling of local k-space neighborhoods) low-rank reconstruction. Furthermore, the MC-HTC framework was applied uniquely to Cartesian uniform undersampling by incorporating a novel complementary k-space sampling strategy where the phase-encoding direction among different contrasts is orthogonally alternated.

CONCLUSION

The proposed MC-HTC approach presents an effective tensor completion framework to jointly reconstruct highly undersampled multicontrast 2D datasets without coil-sensitivity calibration.

Improved MUSSELS reconstruction for high-resolution multi-shot diffusion weighted imaging

PURPOSE

MUSSELS is a one-step iterative reconstruction method for multishot diffusion weighted (msDW) imaging. The current work presents an efficient implementation, termed IRLS MUSSELS, that enables faster reconstruction to enhance its utility for high-resolution diffusion MRI studies.

METHODS

The recently proposed MUSSELS reconstruction belongs to a new class of parallel imaging-based methods that recover artifact-free DWIs from msDW data without needing phase compensation. The reconstruction is achieved via structured low-rank matrix completion algorithms, which are computationally demanding due to the large size of the Hankel matrices and their associated computations involving singular value decompositions. Because of this, computational demands of the MUSSELS reconstruction scales as the matrix size and the number of shots increases, which hinders its practical utility for high-resolution applications. In this work, we derive a computationally efficient MUSSELS formulation by modifying the iterative reweighted least squares (IRLS) method that were proposed earlier to solve such problems. Using whole-brain in vivo data, we show the utility of the IRLS MUSSELS for routine high-resolution studies with reduced computational burden.

RESULTS

IRLS MUSSELS provides about five times faster reconstruction for matrix sizes 192 × 192 and 256 × 256 compared to the earlier MUSSELS implementation. The widely employed conjugate symmetry priors can also be incorporated into IRLS MUSSELS to reduce blurring of the partial Fourier acquisitions, without incurring much computational burden.

CONCLUSIONS

The proposed method is observed to be computationally efficient to enable routine high-resolution studies. The computational complexity matches the traditional msDWI reconstruction methods and provides improved reconstruction results with the additional constraints.