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

DTI-MR fingerprinting for rapid high-resolution whole-brain T1 , T2 , proton density, ADC, and fractional anisotropy mapping

PURPOSE

This study aims to develop a high-efficiency and high-resolution 3D imaging approach for simultaneous mapping of multiple key tissue parameters for routine brain imaging, including T1 , T2 , proton density (PD), ADC, and fractional anisotropy (FA). The proposed method is intended for pushing routine clinical brain imaging from weighted imaging to quantitative imaging and can also be particularly useful for diffusion-relaxometry studies, which typically suffer from lengthy acquisition time.

METHODS

To address challenges associated with diffusion weighting, such as shot-to-shot phase variation and low SNR, we integrated several innovative data acquisition and reconstruction techniques. Specifically, we used M1-compensated diffusion gradients, cardiac gating, and navigators to mitigate phase variations caused by cardiac motion. We also introduced a data-driven pre-pulse gradient to cancel out eddy currents induced by diffusion gradients. Additionally, to enhance image quality within a limited acquisition time, we proposed a data-sharing joint reconstruction approach coupled with a corresponding sequence design.

RESULTS

The phantom and in vivo studies indicated that the T1 and T2 values measured by the proposed method are consistent with a conventional MR fingerprinting sequence and the diffusion results (including diffusivity, ADC, and FA) are consistent with the spin-echo EPI DWI sequence.

CONCLUSION

The proposed method can achieve whole-brain T1 , T2 , diffusivity, ADC, and FA maps at 1-mm isotropic resolution within 10 min, providing a powerful tool for investigating the microstructural properties of brain tissue, with potential applications in clinical and research settings.

3D diffusion MRI with twin navigator-based GRASE and comparison with 2D EPI for tractography in the human brain

PURPOSE

3D pulse sequences enable high-resolution acquisition with a high SNR and ideal slice profiles, which, however, is particularly difficult for diffusion MRI (dMRI) due to the additional phase errors from diffusion encoding.

METHODS

We proposed a twin navigator-based 3D diffusion-weighted gradient spin-echo (GRASE) sequence to correct the phase errors between shots and between odd and even spin echoes for human whole-brain acquisition. We then compared the SNR of 3D GRASE and 2D simultaneous multi-slice EPI within the same acquisition time. We further tested the performance of 2D versus 3D acquisition at equivalent SNR on fiber tracking and microstructural mapping, using the diffusion tensor and high-order fiber orientation density-based metrics.

RESULTS

The proposed twin navigator approach removed multi-shot phase errors to some extent in the whole brain dMRI, and the 2D navigator performed better than the 1D navigator. Comparisons of SNR between the 2D simultaneous multi-slice EPI and 3D GRASE sequences demonstrated that the SNR of the GRASE sequence was 1.4-1.5-fold higher than the EPI sequence at an equivalent scan time. More importantly, we found a significantly higher fiber cross-section in the cerebrospinal tract, as well as richer subcortical fibers (U-fibers) using the 3D GRASE sequence compared to 2D EPI.

CONCLUSION

The twin navigator-based 3D diffusion-weighted-GRASE sequence minimized the multishot phase error and effectively improved the SNR for whole-brain dMRI acquisition. We found differences in fiber tracking and microstructural mapping between 2D and 3D acquisitions, possibly due to the different slice profiles.

Inversion-Recovery-Prepared Oscillating Gradient Sequence Improves Diffusion-Time Dependency Measurements in the Human Brain

BACKGROUND

Oscillating gradient diffusion MRI (dMRI) enables measurements at a short diffusion-time (t_d ), but it is challenging for clinical systems. Particularly, the low b-value and low resolution may give rise to cerebrospinal fluid (CSF) contamination.

PURPOSE

To assess the effect of CSF partial volume on t_d -dMRI measurements and efficacy of inversion-recovery (IR) prepared oscillating and pulsed gradient dMRI sequence to improve t_d -dMRI measurements in the human brain.

STUDY TYPE

Prospective.

SUBJECTS

Ten normal volunteers and six glioma patients.

FIELD STRENGTH/SEQUENCE

A 3 T; three-dimensional (3D) IR-prepared oscillating gradient-prepared gradient spin-echo (GRASE) and two-dimensional (2D) IR-prepared oscillating gradient echo-planar imaging (EPI) sequences.

ASSESSMENT

We assessed the t_d -dependent patterns of apparent diffusion coefficient (ADC) in several gray and white matter structures, including the hippocampal subfields (head, body, and tail), cortical gray matter, thalamus, and posterior white matter in normal volunteers. Pulsed gradient (0 Hz) and oscillating gradients at frequencies of 20 Hz, 40 Hz, and 60 Hz dMRI were acquired with GRASE and EPI sequences with or without the IR module. We also tested the t_d -dependency patterns in glioma patients using the EPI sequence with or without the IR module.

STATISTICAL TESTS

The differences in ADC across the different t_d s were compared by one-way ANOVA followed by post hoc pairwise t-tests with Bonferroni correction.

RESULTS

In the healthy subjects, brain regions that were possibly contaminated by CSF signals, such as the hippocampus (head, body, and tail) and cortical gray matter, t_d -dependent ADC changes were only significant with the IR-prepared 2D and 3D sequences but not with the non-IR sequences. In brain glioblastomas patients, significantly higher t_d -dependence was observed in the tumor region with the IR module than that without IR (slope = 0.0196 μm² /msec² vs. 0.0034 μm² /msec² ).

CONCLUSION

The IR-prepared sequence effectively suppressed the CSF partial volume effect and significantly improved the t_d -dependent measurements in the human brain.

EVIDENCE LEVEL

1 TECHNICAL EFFICACY: Stage 1.

Tuned bipolar oscillating gradients for mapping frequency dispersion of diffusion kurtosis in the human brain

PURPOSE

Oscillating gradient spin-echo (OGSE) sequences have demonstrated an ability to probe time-dependent microstructural features, although they often suffer from low SNR due to increased TEs. In this work we introduce frequency-tuned bipolar (FTB) gradients as a variation of oscillating gradients with reduced TE and demonstrate their utility by mapping the frequency dispersion of kurtosis in human subjects.

METHODS

An FTB oscillating gradient waveform is presented that provides encoding of 1.5 net oscillation periods, thereby reducing the TE of the acquisition. Simulations were performed to determine an optimal protocol based on the SNR of kurtosis frequency dispersion-defined as the difference in kurtosis between pulsed and oscillating gradient acquisitions. Healthy human subjects were scanned at 7T using pulsed gradient and an optimized 23 Hz FTB protocol, which featured a maximum b-value of 2500 s/mm² . In addition, to directly compare existing methods, measurements using traditional cosine OGSE were also acquired.

RESULTS

FTB oscillating gradients demonstrated equivalent frequency-dependent diffusion measurements compared with cosine-modulated OGSE while enabling a significant reduction in TE. Optimization and in vivo results suggest that FTB gradients provide increased SNR of kurtosis dispersion maps compared with traditional cosine OGSE. The optimized FTB gradient protocol demonstrated consistent reductions in apparent kurtosis values and increased diffusivity in generated frequency dispersion maps.

CONCLUSIONS

This work presents an alternative to traditional cosine OGSE sequences, enabling more time-efficient acquisitions of frequency-dependent diffusion quantities as demonstrated through in vivo kurtosis frequency dispersion maps.

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.

Nonparametric 5D D-R2 distribution imaging with single-shot EPI at 21.1 T: Initial results for in vivo rat brain

In vivo human diffusion MRI is by default performed using single-shot EPI with greater than 50-ms echo times and associated signal loss from transverse relaxation. The individual benefits of the current trends of increasing B₀ to boost SNR and employing more advanced signal preparation schemes to improve the specificity for selected microstructural properties eventually may be cancelled by increased relaxation rates at high B₀ and echo times with advanced encoding. Here, initial attempts to translate state-of-the-art diffusion-relaxation correlation methods from 3 T to 21.1 T are made to identify hurdles that need to be overcome to fulfill the promises of both high SNR and readily interpretable microstructural information.

VRfRNet: Volumetric ROI fODF reconstruction network for estimation of multi-tissue constrained spherical deconvolution with only single shell dMRI

Diffusion MRI (dMRI) is one of the most popular techniques for studying the brain structure, mainly the white matter region. Among several sampling methods in dMRI, the high angular resolution diffusion imaging (HARDI) technique has attracted researchers due to its more accurate fiber orientation estimation. However, the current single-shell HARDI makes the intravoxel structure challenging to estimate accurately. While multi-shell acquisition can address this problem, it takes a longer scanning time, restricting its use in clinical applications. In addition, most existing dMRI scanners with low gradient-strengths often acquire single-shell up to b=1000s/mm² because of signal-to-noise ratio issues and severe image artefacts. Hence, we propose a novel generative adversarial network, VRfRNet, for the reconstruction of multi-shell multi-tissue fiber orientation distribution function from single-shell HARDI volumes. Such a transformation learning is performed in the spherical harmonics (SH) space, as raw input HARDI volume is transformed to SH coefficients to soften gradient directions. The proposed VRfRNet consists of several modules, such as multi-context feature enrichment module, feature level attention, and softmax level attention. In addition, three loss functions have been used to optimize network learning, including L₁, adversarial, and total variation. The network is trained and tested using standard qualitative and quantitative performance metrics on the publicly available HCP data-set.

Time-Dependent Diffusion MRI for Quantitative Microstructural Mapping of Prostate Cancer

Background Recently developed time-dependent diffusion MRI has potential in characterizing cellular tissue microstructures; however, its value in imaging prostate cancer (PCa) remains unknown. Purpose To investigate the feasibility of time-dependent diffusion MRI-based microstructural mapping for noninvasively characterizing cellular properties of PCa and for discriminating between clinically significant PCa and clinically insignificant disease. Materials and Methods Men with a clinical suspicion of PCa were enrolled prospectively between October 2019 and August 2020. Time-dependent diffusion MRI data were acquired with pulsed and oscillating gradient diffusion MRI sequences at an equivalent diffusion time of 7.5-30 msec on a 3.0-T scanner. Time-dependent diffusion MRI-based microstructural parameters, including cell diameter, intracellular volume fraction, cellularity, and diffusivities, were estimated with a two-compartment model. These were compared for different International Society of Urological Pathology grade groups (GGs), and their performance in discriminating clinically significant PCa (GG >1) from clinically insignificant disease (benign and GG 1) was determined with a linear discriminant analysis. The fitted microstructural parameters were validated by means of correlation with histopathologic measurements. Results In the 48 enrolled men, the time-dependent diffusion MRI measurements showed that higher GG was correlated with higher intracellular volume fraction and higher cellularity (intracellular volume fraction = 0.22, 0.36, 0.34, 0.37, and 0.40 in GGs 1-5, respectively; P < .001 at one-way analysis of variance), while lower cell diameter was found at higher GGs (diameter = 23.4, 18.3, 19.2, 17.9, and 18.5 μm in GGs 1-5, respectively; P = .002). Among all measurements derived from time-dependent diffusion MRI, cellularity achieved the highest diagnostic performance, with an accuracy of 92% (44 of 48 participants) and area under the receiver operating characteristic curve of 0.96 (95% CI: 0.87, 0.99) in discriminating clinically significant PCa from clinically insignificant disease. Microstructural mapping was supported by positive correlations between time-dependent diffusion MRI-based and pathologic examination-based intracellular volume fraction (r = 0.83; P < .001). Conclusion Time-dependent diffusion MRI-based microstructural mapping correlates with pathologic findings and demonstrates promise for characterizing prostate cancer. © RSNA, 2022 Online supplemental material is available for this article. See also the editorial by Chatterjee and Oto in this issue.

Evaluating diffusion dispersion across an extended range of b-values and frequencies: Exploiting gap-filled OGSE shapes, strong gradients, and spiral readouts

Purpose

To address the long echo times and relatively weak diffusion sensitization that typically limit oscillating gradient spin‐echo (OGSE) experiments, an OGSE implementation combining spiral readouts, gap‐filled oscillating gradient shapes providing stronger diffusion encoding, and a high‐performance gradient system is developed here and utilized to investigate the tradeoff between b‐value and maximum OGSE frequency in measurements of diffusion dispersion (i.e., the frequency dependence of diffusivity) in the in vivo human brain. In addition, to assess the effects of the marginal flow sensitivity introduced by these OGSE waveforms, flow‐compensated variants are devised for experimental comparison.

Methods

Using DTI sequences, OGSE acquisitions were performed on three volunteers at b‐values of 300, 500, and 1000 s/mm² and frequencies up to 125, 100, and 75 Hz, respectively; scans were performed for gap‐filled oscillating gradient shapes with and without flow sensitivity. Pulsed gradient spin‐echo DTI acquisitions were also performed at each b‐value. Upon reconstruction, mean diffusivity (MD) maps and maps of the diffusion dispersion rate were computed.

Results

The power law diffusion dispersion model was found to fit best to MD measurements acquired at b = 1000 s/mm² despite the associated reduction of the spectral range; this observation was consistent with Monte Carlo simulations. Furthermore, diffusion dispersion rates without flow sensitivity were slightly higher than flow‐sensitive measurements.

Conclusion

The presented OGSE implementation provided an improved depiction of diffusion dispersion and demonstrated the advantages of measuring dispersion at higher b‐values rather than higher frequencies within the regimes employed in this study.

Time-dependent diffusion MRI probes cerebellar microstructural alterations in a mouse model of Down syndrome

The cerebellum is a complex system with distinct cortical laminar organization. Alterations in cerebellar microstructure are common and associated with many factors such as genetics, cancer and ageing. Diffusion MRI (dMRI) provides a non-invasive tool to map the brain structural organization, and the recently proposed diffusion-time (td )-dependent dMRI further improves its capability to probe the cellular and axonal/dendritic microstructures by measuring water diffusion at multiple spatial scales. The td -dependent diffusion profile in the cerebellum and its utility in detecting cerebellar disorders, however, are not yet elucidated. Here, we first deciphered the spatial correspondence between dMRI contrast and cerebellar layers, based on which the cerebellar layer-specific td -dependent dMRI patterns were characterized in both euploid and Ts65Dn mice, a mouse model of Down syndrome. Using oscillating gradient dMRI, which accesses diffusion at short td 's by modulating the oscillating frequency, we detected subtle changes in the apparent diffusivity coefficient of the cerebellar internal granular layer and Purkinje cell layer of Ts65Dn mice that were not detectable by conventional pulsed gradient dMRI. The detection sensitivity of oscillating gradient dMRI increased with the oscillating frequency at both the neonatal and adult stages. The td -dependence, quantified by ΔADC map, was reduced in Ts65Dn mice, likely associated with the reduced granule cell density and abnormal dendritic arborization of Purkinje cells as revealed from histological evidence. Our study demonstrates superior sensitivity of short-td diffusion using oscillating gradient dMRI to detect cerebellar microstructural changes in Down syndrome, suggesting the potential application of this technique in cerebellar disorders.

MR cell size imaging with temporal diffusion spectroscopy

Cytological features such as cell size and intracellular morphology provide fundamental information on cell status and hence may provide specific information on changes that arise within biological tissues. Such information is usually obtained by invasive biopsy in current clinical practice, which suffers several well-known disadvantages. Recently, novel MRI methods such as IMPULSED (imaging microstructural parameters using limited spectrally edited diffusion) have been developed for direct measurements of mean cell size non-invasively. The IMPULSED protocol is based on using temporal diffusion spectroscopy (TDS) to combine measurements of water diffusion over a wide range of diffusion times to probe cellular microstructure over varying length scales. IMPULSED has been shown to provide rapid, robust, and reliable mapping of mean cell size and is suitable for clinical imaging. More recently, cell size distributions have also been derived by appropriate analyses of data acquired with IMPULSED or similar sequences, which thus provides MRI-cytometry. This review summarizes the basic principles, practical implementations, validations, and example applications of MR cell size imaging based on TDS and demonstrates how cytometric information can be used in various applications. In addition, the limitations and potential future directions of MR cytometry are identified including the diagnosis of nonalcoholic steatohepatitis of the liver and the assessment of treatment response of cancers.

A novel spectrally selective fat saturation pulse design with robustness to B0 and B1 inhomogeneities: A demonstration on 3D T1-weighted breast MRI at 3 T

PURPOSE

Spectrally selective fat saturation (FatSat) sequence is commonly used to suppress signal from adipose tissue. Conventional SINC-shaped pulses are sensitive to B₀ off-resonance and B₁⁺ offset. Uniform fat saturation with large spatial coverage is especially challenging for the body and breast MRI. The aim of this study is to develop spectrally selective FatSat pulses that offer more immunity to B₀/B₁⁺ field inhomogeneities than SINC pulses and evaluate them in bilateral breast imaging at 3 T.

MATERIALS AND METHODS

Optimized composite pulses (OCP) were designed based on the optimal control theory with robustness to a targeted B₀/ B₁⁺ conditions. OCP pulses also allows flexible flip angles to meet different requirements. Comparisons with the vendor-provided SINC pulses were conducted by numerical simulation and in vivo scans using a 3D T₁-weighted (T₁w) gradient-echo (GRE) sequence with coverage of the whole-breast.

RESULTS

Simulation revealed that OCP pulses yielded almost half of the transition band and much less sensitivity to B₁⁺ inhomogeneity compared to SINC pulses with B₀ off-resonance within ±200 Hz and B₁⁺ scale error within ±0.3 (P < 0.001). Across five normal subjects, OCP FatSat pulses produced 25-41% lower residual fat signals (P < 0.05) with 27-36% less spatial variation (P < 0.05) than SINC.

CONCLUSION

In contrast to conventional SINC-shaped pulses, the newly designed OCP FatSat pulses mitigated challenges of wide range of B₀/ B₁⁺ field inhomogeneities and achieved more uniform fat suppression in bilateral breast T₁w imaging at 3 T.