2435

Comparison of time-division multiplexing and multi-echo sequences for accelerated relaxation diffusion MRI
Qiang Liu1,2, Ante Zhu3, Imam Ahmed Shaik1, Maxim Zaitsev4, Thomas Foo3, Jon-Fredrik Nielsen5, Carl-Fredrik Westin1, Yogesh Rathi1, and Lipeng Ning1
1Brigham and Women's Hospital, Harvard Medical School, Boston, MA, United States, 2School of Biomedical Engineering, Southern Medical University, Guangzhou, China, 3Technology and Innovation Center, GE Healthcare, Niskayuna, NY, United States, 4Division of Medical Physics, Department of Radiology, University Medical Center Freiburg, Freiburg, Germany, 5fMRI Laboratory and Biomedical Engineering, University of Michigan, Ann Arbor, MI, United States

Synopsis

Keywords: Diffusion Acquisition, Quantitative Imaging

Motivation: Compare multi-echo and time-division multiplexing (TDM) sequences that accelerate multi-TE dMRI scans.

Goal(s): To investigate the accuracy of multi-TE diffusion measures using multi-echo and TDM sequences.

Approach: Standard single-TE spin-echo echo-planar imaging (EPI), dual-echo EPI, and TDM-EPI sequences were implemented using Pulseq with matched diffusion and readout gradients. A calibration phantom and a human subject were scanned on two scanners to examine the accuracy of R2 (1/T2) and apparent diffusivity coefficient (ADC) maps.

Results: TDM demonstrates up to 160% reduction in estimated R2 on phantom and 30% in in-vivo data compared to dual-echo sequence.

Impact: The proposed TDM sequence allows for fast and accurate T2 and ADC estimation compared to the conventional multi-echo sequence, making it a potential tool for relaxation-diffusion imaging.

Introduction

Joint modeling of diffusion magnetic resonance imaging (MRI) with multiple echo times (TEs) has been investigated to derive novel microstructural measures that improve standard quantitative relaxometry and diffusion imaging.1–3 However, the longer scan time is a limitation of this technique in clinical applications. A standard method to accelerate multi-TE scans is to use multi-echo sequences that apply multiple 180° refocusing pulses to generate a train of echoes. However, artifacts often arise in multi-echo images due to transmit field (B1+) inhomogeneity and imperfect refocusing pulses, leading to biases in T2 mapping.4,5 Time-division multiplexing (TDM) sequence6,7 is an alternative strategy to accelerate multi-TE scans by acquiring multiple slices at different TEs to avoid multiple 180° pulses per slice. It uses echo-time shifting gradients to remove interference between slices and coupling with gradient gradients. In this study, we aim to evaluate and compare the accuracy of multi-echo and TDM sequences for multi-TE diffusion scans.

Methods

This study follows approval from the local IRBs, and each human subject was scanned in MRI after signing consent forms. We implemented the standard single-TE EPI (reference), dual-echo, and TDM sequences using the vendor-neutral platform Pulseq with matched diffusion and readout gradients. A schematic illustration of the sequences is shown in Figure 1. The sequences were used to scan a phantom and a human subject on a Siemens Prisma whole-body scanner with the following parameters: voxel size=2x2, TE1/TE2/TR=60/90/4000ms, b-value=1000. The subject was scanned again on an advanced GE MAGNUS head-only scanner with voxel size=1.5x1.5, TE1/TE2/TR=45/90/4500ms, b-value=1000.
We computed the normalized root mean squared error (NRMSE) in signal intensity between TDM, dual-echo, and standard single echo sequences. We used exponential fitting to derive the R2 (1/T2) maps and calculated the relative error. ADC was fitted for all the sequences at the two echo times of TE1 and TE2 separately. For the phantom scan, we extracted the ADC measures from two groups of regions of interest (ROIs) categorized by low and high diffusivity values, with 8 ROIs in each group. Then, a Wilcoxon-signed rank test was performed to examine the difference in ADC measurements between the data from the three sequences. For the data from the MAGNUS scanner, we compared b0 image intensities at TE2 using relative error and derived R2 maps from the two echoes.

Results

The top two rows of Figure 2 compare the signal intensities and ADC maps from three sequences at TE1 (45ms). As expected, all three sequences show similar image intensities and both dual-echo and TDM have low NRMSE (~0.16) and absolute difference in ADC (~0.23) compared to the standard sequence. The 3rd-4th row of Figure 2 compares signal intensities and ADC maps at TE2 (90ms). The dual-echo sequence showed a much higher NRMSE (0.3 v.s. 0.16; 87.5% lower error with TDM) and higher absolute difference (0.32 v.s. 0.24, 35% lower difference with TDM) in ADC maps compared to the TDM sequence. The last row shows the R2 maps. The estimated R2 map from the TDM sequence closely resembles that from the single echo gold standard. However, the R2 map obtained from the dual-echo sequence exhibits a noticeable bias (~40% difference in longer T2 and ~160% in shorter T2 ROIs).
The blue and red boxes in Figure 3A show the two groups of ROIs on the phantom of which ADC values are displayed in Figure 3B. Notably, a significant difference is observed between dual-echo and standard single-echo sequences at TE2 for the group with lower ADC (p=0.0391). No statistical difference is seen using the TDM sequence.
Figure 4 illustrates the evident bias in the TE2 intensity for the dual-echo sequence in in-vivo data. Comparatively, the R2 map obtained from TDM closely matches the single-echo sequence, whereas the R2 map from the dual-echo sequence shows a ~30% bias.
Figure 5 shows that the dual-echo sequence also induces ~30% higher relative error in signal intensity and biased R2 maps, on the MAGNUS head-only MRI system.

Discussion and Conclusion

In this study, with Pulseq, we implemented two sequences designed for fast multi-TE imaging and compared them with the standard single-echo sequence. For both phantom and in-vivo data, our proposed TDM sequence accurately produced T2 maps, while the traditionally used multi-echo sequence exhibited strong bias. Interestingly, the ADC map calculated from dual-echo showed a higher bias at low ADC in the phantom. Leveraging the high-performance MAGNUS gradient system, we achieved a higher resolution DWI with shorter TEs, highlighting the robust advantage of T2 estimation from our proposed TDM sequence.

Acknowledgements

This study is supported by Massachusetts Life Sciences Center, and NIH grants K01MH117346, R21MH126396, R01MH116173, R01MH125860, R01EB032378, and R01NS125781.

References

1 Ning, L., Gagoski, B., Szczepankiewicz, F., Westin, C. F. & Rathi, Y. Joint RElaxation-Diffusion Imaging Moments to Probe Neurite Microstructure. IEEE Trans Med Imaging 39, 668–677 (2020).

2. Kim, D., Doyle, E. K., Wisnowski, J. L., Kim, J. H. & Haldar, J. P. Diffusion-relaxation correlation spectroscopic imaging: A multidimensional approach for probing microstructure. Magn Reson Med 78, 2236–2249 (2017).

3. Afzali, M. et al. MR Fingerprinting with b-Tensor Encoding for Simultaneous Quantification of Relaxation and Diffusion in a Single Scan. Magn Reson Med 88, 2043–2057 (2022).

4. Petrovic, A., Aigner, C. S., Rund, A. & Stollberger, R. A time domain signal equation for multi-echo spin-echo sequences with arbitrary excitation and refocusing angle and phase. Journal of Magnetic Resonance 309, (2019).

5. Klupp, E. et al. B1-insensitive T2 mapping of healthy thigh muscles using a T2-prepared 3D TSE sequence. PLoS One 12, (2017).

6. Ji, Y., Gagoski, B., Hoge, W. S., Rathi, Y. & Ning, L. Accelerated diffusion and relaxation-diffusion MRI using time-division multiplexing EPI. Magn Reson Med 86, 2528–2541 (2021).

7. Ji, Y. et al. Accelerating joint relaxation-diffusion MRI by integrating time division multiplexing and simultaneous multi-slice (TDM-SMS) strategies. Magn Reson Med 87, 2697–2709 (2022).

Figures

Figure 1. Sequence diagrams of single-TE (Standard), Dual-echo, and Time-division multiplexing (TDM) sequences. The scan protocols for the diffusion phantom and in vivo scan are also shown below the images.

Figure 2. A representative case of the phantom scan. As expected, the TDM and dual-echo images are similar to the Standard for TE1, as shown in the first two rows. The 4th row shows that the TE2 of the dual-echo leads to a larger bias (NRMSE=0.3) compared to TDM. The last row shows that the R2 map derived from TDM closely resembles the Standard, but the Dual-echo sequence leads to ~ 40% bias in R2 and much higher error in the ADC map of TE2.

Figure 3. The blue and red boxes in (A) indicate two groups of regions of interest (ROIs). (B) shows the comparison results between the ADC measurements obtained with different sequences. There is a significant difference between the ADC derived from Dual-echo at TE2 compared to the ADC values from Standard (p=0.0391) in the group of lower ADC values.

Figure 4. A representative case of in-vivo scan on Prisma. The intensity of DWI of Dual-echo TE2 shows a larger bias than the TDM of TE2. As highlighted in the R2 map, TDM obtains an R2 map that matches the R2 map from Standard, while the R2 map derived from Dual-echo shows approximately 30% bias.

Figure 5. Comparison of signal intensity (Rows 1-3) and R2 maps (Rows 4-5) for the data acquired from the GE MAGNUS scanner. The Dual-echo sequence shows a much higher bias in the signal intensity and R2 maps than the TDM sequence.

Proc. Intl. Soc. Mag. Reson. Med. 32 (2024)
2435
DOI: https://doi.org/10.58530/2024/2435