Analysis of phase error effects in multishot diffusion-prepared turbo spin echo imaging

Distortion-free water-fat separated diffusion-weighted imaging using spatiotemporal joint reconstruction

PURPOSE

Diffusion-weighted imaging (DWI) suffers from geometric distortion and chemical shift artifacts due to the commonly used Echo Planar Imaging (EPI) trajectory. Even with fat suppression in DWI, severe B0 and B1 variations can result in residual fat, which becomes both a source of image artifacts and a confounding factor in diffusion-weighted contrast in distinguishing benign and malignant tissues. This work presents a method for acquiring distortion-free diffusion-weighted images using spatiotemporal acquisition and joint reconstruction. Water-fat separation is performed by chemical-shift encoding.

METHODS

Spatiotemporal acquisition is employed to obtain distortion-free images at a series of echo times. Chemical-shift encoding is used for water-fat separation. Reconstruction and separation are performed jointly in the spat-spectral domain. To address the shot-to-shot motion-induced phase in DWI, an Fast Spin Echo (FSE)-based phase navigator is incorporated into the sequence to obtain distortion-free phase information. The proposed method was validated in phantoms and in vivo for the brain, head and neck, and breast.

RESULTS

The proposed method enables the acquisition of distortion-free diffusion-weighted images in the presence of B0 field inhomogenieties commonly observed in the body. Water and fat components are separated with no obvious spectral leakage artifacts. The estimated Apparent Diffusion Coefficient (ADC) is comparable to that of multishot DW-EPI.

CONCLUSION

Distortion-free, water-fat separated diffusion-weighted images in body can be obtained through the utilization of spatiotemporal acquisition and joint reconstruction methods.

Distortionless, free-breathing, and respiratory resolved 3D diffusion weighted imaging of the abdomen

PURPOSE

Abdominal imaging is frequently performed with breath holds or respiratory triggering to reduce the effects of respiratory motion. Diffusion weighted sequences provide a useful clinical contrast but have prolonged scan times due to low signal-to-noise ratio (SNR), and cannot be completed in a single breath hold. Echo-planar imaging (EPI) is the most commonly used trajectory for diffusion weighted imaging but it is susceptible to off-resonance artifacts. A respiratory resolved, three-dimensional (3D) diffusion prepared sequence that obtains distortionless diffusion weighted images during free-breathing is presented. Techniques to address the myriad of challenges including: 3D shot-to-shot phase correction, respiratory binning, diffusion encoding during free-breathing, and robustness to off-resonance are described.

METHODS

A twice-refocused, M1-nulled diffusion preparation was combined with an RF-spoiled gradient echo readout and respiratory resolved reconstruction to obtain free-breathing diffusion weighted images in the abdomen. Cartesian sampling permits a sampling density that enables 3D shot-to-shot phase navigation and reduction of transient fat artifacts. Theoretical properties of a region-based shot rejection are described. The region-based shot rejection method was evaluated with free-breathing (normal and exaggerated breathing), and respiratory triggering. The proposed sequence was compared in vivo with multishot DW-EPI.

RESULTS

The proposed sequence exhibits no evident distortion in vivo when compared to multishot DW-EPI, robustness to B0 and B1 field inhomogeneities, and robustness to motion from different respiratory patterns.

CONCLUSION

Acquisition of distortionless, diffusion weighted images is feasible during free-breathing with a b-value of 500 s/mm2, scan time of 6 min, and a clinically viable reconstruction time.

Self-calibrated subspace reconstruction for multidimensional MR fingerprinting for simultaneous relaxation and diffusion quantification

PURPOSE

To propose a new reconstruction method for multidimensional MR fingerprinting (mdMRF) to address shading artifacts caused by physiological motion-induced measurement errors without navigating or gating.

METHODS

The proposed method comprises two procedures: self-calibration and subspace reconstruction. The first procedure (self-calibration) applies temporally local matrix completion to reconstruct low-resolution images from a subset of under-sampled data extracted from the k-space center. The second procedure (subspace reconstruction) utilizes temporally global subspace reconstruction with pre-estimated temporal subspace from low-resolution images to reconstruct aliasing-free, high-resolution, and time-resolved images. After reconstruction, a customized outlier detection algorithm was employed to automatically detect and remove images corrupted by measurement errors. Feasibility, robustness, and scan efficiency were evaluated through in vivo human brain imaging experiments.

RESULTS

The proposed method successfully reconstructed aliasing-free, high-resolution, and time-resolved images, where the measurement errors were accurately represented. The corrupted images were automatically and robustly detected and removed. Artifact-free T1, T2, and ADC maps were generated simultaneously. The proposed reconstruction method demonstrated robustness across different scanners, parameter settings, and subjects. A high scan efficiency of less than 20 s per slice has been achieved.

CONCLUSION

The proposed reconstruction method can effectively alleviate shading artifacts caused by physiological motion-induced measurement errors. It enables simultaneous and artifact-free quantification of T1, T2, and ADC using mdMRF scans without prospective gating, with robustness and high scan efficiency.

A joint linear reconstruction for multishot diffusion weighted non-Carr-Purcell-Meiboom-Gill fast spin echo with full signal

PURPOSE

Diffusion weighted Fast Spin Echo (DW-FSE) is a promising approach for distortionless DW imaging that is robust to system imperfections such as eddy currents and off-resonance. Due to non-Carr-Purcell-Meiboom-Gill (CPMG) magnetization, most DW-FSE sequences discard a large fraction of the signal (2 - 2 × $$ \sqrt{2}-2\times $$ ), reducing signal-to-noise ratio (SNR) efficiency compared to DW-EPI. The full FSE signal can be preserved by quadratically incrementing the transmit phase of the refocusing pulses, but this method of resolving non-CPMG magnetization has only been applied to single-shot DW-FSE due to challenges associated with image reconstruction. We present a joint linear reconstruction for multishot quadratic phase increment data that addresses these challenges and corrects ghosting from both shot-to-shot phase and intrashot signal oscillations. Multishot imaging reduces T2 blur and joint reconstruction of shots improves g-factor performance. A thorough analysis on the condition number of the proposed linear system is described.

METHODS

A joint multishot reconstruction is derived from the non-CPMG signal model. Multishot quadratic phase increment DW-FSE was tested in a standardized diffusion phantom and compared to single-shot DW-FSE and DW-EPI in vivo in the brain, cervical spine, and prostate. The pseudo multiple replica technique was applied to generate g-factor and SNR maps.

RESULTS

The proposed joint shot reconstruction eliminates ghosting from shot-to-shot phase and intrashot oscillations. g-factor performance is improved compared to previously proposed reconstructions, permitting efficient multishot imaging. apparent diffusion coefficient estimates in phantom experiments and in vivo are comparable to those obtained with conventional methods.

CONCLUSION

Multi-shot non-CPMG DW-FSE data with full signal can be jointly reconstructed using a linear model.

MR Fingerprinting with b-Tensor Encoding for Simultaneous Quantification of Relaxation and Diffusion in a Single Scan

PURPOSE

Although both relaxation and diffusion imaging are sensitive to tissue microstructure, studies have reported limited sensitivity and robustness of using relaxation or conventional diffusion alone to characterize tissue microstructure. Recently, it has been shown that tensor-valued diffusion encoding and joint relaxation-diffusion quantification enable more reliable quantification of compartment-specific microstructural properties. However, scan times to acquire such data can be prohibitive. Here, we aim to simultaneously quantify relaxation and diffusion using MR fingerprinting (MRF) and b-tensor encoding in a clinically feasible time.

METHODS

We developed multidimensional MRF scans (mdMRF) with linear and spherical b-tensor encoding (LTE and STE) to simultaneously quantify T1, T2, and ADC maps from a single scan. The image quality, accuracy, and scan efficiency were compared between the mdMRF using LTE and STE. Moreover, we investigated the robustness of different sequence designs to signal errors and their impact on the maps.

RESULTS

T1 and T2 maps derived from the mdMRF scans have consistently high image quality, while ADC maps are sensitive to different sequence designs. Notably, the fast imaging steady state precession (FISP)-based mdMRF scan with peripheral pulse gating provides the best ADC maps that are free of image distortion and shading artifacts.

CONCLUSION

We demonstrated the feasibility of quantifying T1, T2, and ADC maps simultaneously from a single mdMRF scan in around 24 s/slice. The map quality and quantitative values are consistent with the reference scans.

Volumetric and multispectral DWI near metallic implants using a non-linear phase Carr-Purcell-Meiboom-Gill diffusion preparation

PURPOSE

DWI near metal implants has not been widely explored due to substantial challenges associated with through-slice and in-plane distortions, the increased encoding requirement of different spectral bins, and limited SNR. There is no widely adopted clinical protocol for DWI near metal since the commonly used EPI trajectory fails completely due to distortion from extreme off-resonance ranging from 2 to 20 kHz. We present a sequence that achieves DWI near metal with moderate b-values (400-500 s/mm² ) and volumetric coverage in clinically feasible scan times.

THEORY AND METHODS

Multispectral excitation with Cartesian sampling, view angle tilting, and kz phase encoding reduce in-plane and through-plane off-resonance artifacts, and Carr-Purcell-Meiboom-Gill (CPMG) spin-echo refocusing trains counteract T2* effects. The effect of random phase on the refocusing train is eliminated using a stimulated echo diffusion preparation. Root-flipped Shinnar-Le Roux refocusing pulses permits preparation of a high spectral bandwidth, which improves imaging times by reducing the number of excitations required to cover the desired spectral range. B₁ sensitivity is reduced by using an excitation that satisfies the CPMG condition in the preparation. A method for ADC quantification insensitive to background gradients is presented.

RESULTS

Non-linear phase refocusing pulses reduces the peak B₁ by 46% which allows RF bandwidth to be doubled. Simulations and phantom experiments show that a non-linear phase CPMG pulse pair reduces B₁ sensitivity. Application in vivo demonstrates complementary contrast to conventional multispectral acquisitions and improved visualization compared to DW-EPI.

CONCLUSION

Volumetric and multispectral DW imaging near metal can be achieved with a 3D encoded sequence.

Lipid droplet-size mapping in human adipose tissue using a clinical 3T system

PURPOSE

To develop a methodology for probing lipid droplet sizes with a clinical system based on a diffusion-weighted stimulated echo-prepared turbo spin-echo sequence and to validate the methodology in water-fat emulsions and show its applicability in ex vivo adipose-tissue samples.

METHODS

A diffusion-weighted stimulated echo-prepared preparation was combined with a single-shot turbo spin-echo readout for measurements at different b-values and diffusion times. The droplet size was estimated with an analytical expression, and three fitting approaches were compared: magnitude-based spatial averaging with voxel-wise residual minimization, complex-based spatial averaging with voxel-wise residual minimization, and complex-based spatial averaging with neighborhood-regularized residual minimization. Simulations were performed to characterize the fitting residual landscape and the approaches' noise performance. The applicability was assessed in oil-in-water emulsions in comparison with laser deflection and in ten human white adipose tissue samples in comparison with histology.

RESULTS

The fitting residual landscape showed a minimum valley with increasing extent as the droplet size increased. In phantoms, a very good agreement of the mean droplet size was observed between the diffusion-weighted MRI-based and the laser deflection measurements, showing the best performance with complex-based spatial averaging with neighborhood-regularized residual minimization processing (R² /P: 0.971/0.014). In the human adipose-tissue samples, complex-based spatial averaging with neighborhood-regularized residual minimization processing showed a significant correlation (R² /P: 0.531/0.017) compared with histology.

CONCLUSION

The proposed acquisition and parameter-estimation methodology was able to probe restricted diffusion effects in lipid droplets. The methodology was validated using phantoms, and its feasibility in measuring an apparent lipid droplet size was demonstrated ex vivo in white adipose tissue.

Diffusion-prepared fast spin echo for artifact-free spinal cord imaging

PURPOSE

Diffusion MRI provides unique contrast important for the detection and examination of pathophysiology after acute neurologic insults, including spinal cord injury. Diffusion weighted imaging of the rodent spinal cord has typically been evaluated with axial EPI readout. However, Diffusion weighted imaging is prone to motion artifacts, whereas EPI is prone to susceptibility artifacts. In the context of acute spinal cord injury, diffusion filtering has previously been shown to improve detection of injury by minimizing the confounding effects of edema. We propose a diffusion-preparation module combined with a rapid acquisition with relaxation enhancement readout to minimize artifacts for sagittal imaging.

METHODS

Sprague-Dawley rats with cervical contusion spinal cord injury were scanned at 9.4 Tesla. The sequence optimization included the evaluation of motion-compensated encoding diffusion gradients, gating strategy, and different spinal cord-specific diffusion-weighting schemes.

RESULTS

A diffusion-prepared rapid acquisition with relaxation enhancement achieved high-quality images free from susceptibility artifacts with both second-order motion-compensated encoding and gating necessary for reduction of motion artifacts. Axial diffusivity obtained from the filtered diffusion-encoding scheme had greater lesion-to-healthy tissue contrast (52%) compared to the similar metric from DTI (25%).

CONCLUSION

This work demonstrated the feasibility of high-quality diffusion sagittal imaging in the rodent cervical cord with diffusion-prepared relaxation enhancement. The sequence and results are expected to improve injury detection and evaluation in acute spinal cord injury.

Diffusion-prepared 3D gradient spin-echo sequence for improved oscillating gradient diffusion MRI

PURPOSE

Oscillating gradient (OG) enables the access of short diffusion times for time-dependent diffusion MRI (dMRI); however, it poses several technical challenges for clinical use. This study proposes a 3D oscillating gradient-prepared gradient spin-echo (OGprep-GRASE) sequence to improve SNR and shorten acquisition time for OG dMRI on clinical scanners.

METHODS

The 3D OGprep-GRASE sequence consisted of global saturation, diffusion encoding, fat saturation, and GRASE readout modules. Multiplexed sensitivity-encoding reconstruction was used to correct the phase errors between multiple shots. We compared the scan time and SNR of the proposed sequence and the conventional 2D-EPI sequence for OG dMRI at 30-90-mm slice coverage. We also examined the time-dependent diffusivity changes with OG dMRI acquired at frequencies of 50 Hz and 25 Hz and pulsed-gradient dMRI at diffusion times of 30 ms and 60 ms.

RESULTS

The OGprep-GRASE sequence reduced the scan time by a factor of 1.38, and increased the SNR by 1.74-2.27 times compared with 2D EPI for relatively thick slice coverage (60-90 mm). The SNR gain led to improved diffusion-tensor reconstruction in the multishot protocols. Image distortion in 2D-EPI images was also reduced in GRASE images. Diffusivity measurements from the pulsed-gradient dMRI and OG dMRI showed clear diffusion-time dependency in the white matter and gray matter of the human brain, using both the GRASE and EPI sequences.

CONCLUSION

The 3D OGprep-GRASE sequence improved scan time and SNR and reduced image distortion compared with the 2D multislice acquisition for OG dMRI on a 3T clinical system, which may facilitate the clinical translation of time-dependent dMRI.

Diffusion tensor imaging and tractography for preoperative assessment of benign peripheral nerve sheath tumors

PURPOSE

To evaluate the diagnostic value of fiber tractography and diffusivity analysis generated from 3D diffusion-weighted (DW) sequences for preoperative assessment of benign peripheral nerve sheath tumors.

METHOD

MR imaging at 3 T was performed in 22 patients (mean age 41.9 ± 17.1y, 13 women) with histologically confirmed schwannomas (N = 18) and histologically confirmed neurofibromas (N = 11), including a 3D DW turbo spin echo sequence with fat suppression. Diffusion tensor parameters were computed and fiber tracks were determined. Evaluation was performed by two radiologists and one orthopedic surgeon blinded for final diagnosis. Mean diffusivity was computed to allow further assessment of tumor microstructure. Preoperative fascicle visualization was graded, fascicles were categorized regarding anatomical location and amount of fascicles surrounding the tumor. The agreement of imaging findings with intraoperative findings was assessed.

RESULTS

On 78.3 % of the DTI images, the fascicle visualization was rated as good or very good. Tractography differences were observed in schwannomas and neurofibromas, showing schwannomas to be significantly more often located eccentrically to the nerve (94.8 %) than neurofibromas (0 %, P < 0.01). Fascicles were significantly more often continuous (87.5 %) in schwannomas, while in neurofibromas, none of the tracks was graded to be continuous (0 %, P = 0.014). A substantial agreement between fiber tracking and surgical anatomy was found regarding the fascicle courses surrounding the tumor (κ = 0.78). Mean diffusivity of schwannomas (1.5 ± 0.2 × 10^-3 mm²/s) was significantly lower than in neurofibromas (1.8 ± 0.2 × 10^-3 mm²/s; P < 0.001). The Youden index showed an optimal cutoff at 1.7 × 10^-3 mm²/s (sensitivity, 0.91; specificity, 0.78; J = 0.69).

CONCLUSIONS

Preoperative diffusion tensor imaging allowed to accurately differentiate between schwannomas and neurofibromas and to describe their location in relation to the nerve fascicles for preoperative planning.

Three-dimensional simultaneous brain T1 , T2 , and ADC mapping with MR Multitasking

PURPOSE

To develop a simultaneous T₁ , T₂ , and ADC mapping method that provides co-registered, distortion-free images and enables multiparametric quantification of 3D brain coverage in a clinically feasible scan time with the MR Multitasking framework.

METHODS

The T₁ /T₂ /diffusion weighting was generated by a series of T₂ preparations and diffusion preparations. The underlying multidimensional image containing 3 spatial dimensions, 1 T₁ weighting dimension, 1 T₂ -preparation duration dimension, 1 b-value dimension, and 1 diffusion direction dimension was modeled as a 5-way low-rank tensor. A separate real-time low-rank model incorporating time-resolved phase correction was also used to compensate for both inter-shot and intra-shot phase inconsistency induced by physiological motion. The proposed method was validated on both phantom and 16 healthy subjects. The quantification of T₁ /T₂ /ADC was evaluated for each case. Three post-surgery brain tumor patients were scanned for demonstration of clinical feasibility.

RESULTS

Multitasking T₁ /T₂ /ADC maps were perfectly co-registered and free from image distortion. Phantom studies showed substantial quantitative agreement ( R 2 = 0.999 ) with reference protocols for T₁ /T₂ /ADC. In vivo studies showed nonsignificant T₁ (P = .248), T₂ (P = .97), ADC (P = .328) differences among the frontal, parietal, and occipital regions. Although Multitasking showed significant differences of T₁ (P = .03), T₂ (P < .001), and ADC (P = .001) biases against the references, the mean bias estimates were small (ΔT₁ % < 5%, ΔT₂ % < 7%, ΔADC% < 5%), with all intraclass correlation coefficients greater than 0.82 indicating "excellent" agreement. Patient studies showed that Multitasking T₁ /T₂ /ADC maps were consistent with the clinical qualitative images.

CONCLUSION

The Multitasking approach simultaneously quantifies T₁ /T₂ /ADC with substantial agreement with the references and is promising for clinical applications.

Isotropic resolution diffusion tensor imaging of lumbosacral and sciatic nerves using a phase-corrected diffusion-prepared 3D turbo spin echo

PURPOSE

To perform in vivo isotropic-resolution diffusion tensor imaging (DTI) of lumbosacral and sciatic nerves with a phase-navigated diffusion-prepared (DP) 3D turbo spin echo (TSE) acquisition and modified reconstruction incorporating intershot phase-error correction and to investigate the improvement on image quality and diffusion quantification with the proposed phase correction.

METHODS

Phase-navigated DP 3D TSE included magnitude stabilizers to minimize motion and eddy-current effects on the signal magnitude. Phase navigation of motion-induced phase errors was introduced before readout in 3D TSE. DTI of lower back nerves was performed in vivo using 3D TSE and single-shot echo planar imaging (ss-EPI) in 13 subjects. Diffusion data were phase-corrected per kz plane with respect to T2 -weighted data. The effects of motion-induced phase errors on DTI quantification was assessed for 3D TSE and compared with ss-EPI.

RESULTS

Non-phase-corrected 3D TSE resulted in artifacts in diffusion-weighted images and overestimated DTI parameters in the sciatic nerve (mean diffusivity [MD] = 2.06 ± 0.45). Phase correction of 3D TSE DTI data resulted in reductions in all DTI parameters (MD = 1.73 ± 0.26) of statistical significance (P ≤ 0.001) and in closer agreement with ss-EPI DTI parameters (MD = 1.62 ± 0.21).

CONCLUSION

DP 3D TSE with phase correction allows distortion-free isotropic diffusion imaging of lower back nerves with robustness to motion-induced artifacts and DTI quantification errors. Magn Reson Med 80:609-618, 2018. © 2018 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 License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited and is not used for commercial purposes.