Diffusion weighted MR spectroscopy without water suppression allows to use water as inherent reference signal to correct for motion-related signal drop
André Döring1, Victor Adalid Lopez1, Vaclav Brandejsky1, Roland Kreis1, and Chris Boesch1

1Depts. Radiology and Clinical Research, University Bern, Bern, Switzerland

Synopsis

A non-water suppressed diffusion-weighting MR spectroscopy sequence based on metabolite-cycling and STEAM is presented and tested in-vitro and in-vivo. The water peak as an inherent reference facilitates a post processing correction of the signal drop induced in individual acquisitions by cardiac and other motion. The correction leads to improved spectral resolution on one hand, but more importantly also to more accurate fitting of ADC values that are found to be smaller than without correction and most likely closer to the true values - and hence better suited for physiological interpretation.

Purpose

Diffusion weighted spectroscopy (DWS) is limited by a low SNR, especially for high b-values1,2. Hence, averaging over many measurements is inevitable, which prolongs the measurement time and, thus, leads to a high susceptibly for motion and line broadening artifacts3. The usage of a non-water-suppressed (nWS) sequence could alleviate some of the problems allowing the use of the water signal as inherent reference, not only for optimal frequency-, eddy current-, and phase-correction as in previous applications4, but also for compensation of motion-related signal drops. This work aims to demonstrate that the measurement of diffusion coefficients (ADC) of brain metabolites substantially profits from using water as internal reference for signal correction in MRS with metabolite cycling (MC) instead of water presaturation.

Methods

A DWS sequence based on STEAM5 was extended to include: (i) adiabatic MC pulses4,6 implemented in the TM period and (ii) adiabatic inversion recovery for FLAIR CSF suppression. The sequence is illustrated schematically in Fig. 1. The sequence was tested in a “braino” phantom with aqueous solutions of typical brain metabolites and in gray matter (GM) of three healthy volunteers (VOL). Acquisition parameters included: TE/TM/TR=37ms/150ms/2500ms; diffusion gradient length δ=11ms; diffusion time Δ=168.2ms and maximum amplitude Gmax=38mT/m. The crusher gradient strength and duration was optimized to avoid interfering spurious echoes and to minimize eddy currents. Maximum b-value in a single direction was 5236s/mm² with 32-128 acquisition per spectrum depending on the b-value. Spectra were acquired on a 3T Siemens TRIO scanner using a multi-channel headcoil. Single shots were stored, though already combined into a single FID (phase-corrected based on 1st data point). Spectral fitting was done in a 2D-fashion in FiTAID7 (simulated metabolite basis sets, experimental macromolecule base spectrum), including ADC determination with a mono-exponential model.

Motion-compensation: Since spectra were recorded without triggering, single acquisitions are affected by various degrees of brain motion, evidenced by a signal drop. This drop was quantified in the non-suppressed water signal (using the median of the top quartile of all shots as reference level) and, thus, all single acquisitions were scaled up to this level before frequency alignment and signal averaging. Only acquisitions that increased the overall SNR were added in. The resulting data was used to determine the water and metabolite signals by adding or subtracting the up- and downfield inverted acquisitions, followed by eddy current correction4.

Results

Fig. 2 summarizes the results obtain in the phantom, in particular the fact that ADCs for metabolites and water are obtained simultaneously – with all components showing mono-exponential decay as expected and the ADCs agreeing well with the literature (Fig. 2b,c and 4). In the in-vivo study, the water signal could be used as reference line for spectral correction and motion-related intensity compensation up to the highest b-value. The effect of eliminating the scans with largest motion-distortion and of intensity scaling is illustrated in Fig. 3. Besides the increase in signal intensity at high b, it can be appreciated from Fig. 3a,b that the lineshape becomes more symmetric and that ghosting artifacts are reduced. Fig. 3c illustrates that the motion-corrected water signal exhibits a non-mono-exponential decay. However, the metabolic attenuation is almost mono-exponential up to b=5236s/mm² (cf. Fig. 3c inset). Estimated ADCs for our initial 3 subjects are smaller after the correction for motion-related signal drops (Fig. 4).

Discussion

In DWS, even weak motion causes enhanced signal attenuation at high b-values, which can easily be misinterpreted as faster diffusion. Using the co-recorded high SNR peak of water, this attenuation can be quantified in single shots and a correction applied to the spectra. As seen from Fig. 3a the signal attenuation for the motion-corrected spectra is thus much lower at high b-values than for the non-motion-corrected case. Furthermore, the fact that the lineshape becomes more symmetric with elimination of outliers (and also spectral reference-based corrections like frequency-alignment) allows for better resolution of metabolites with multiplet patterns like mI or Glu and thus improved model fitting. The current fitting model is preliminary, possibly accounting for apparent intersubject differences.

Conclusion

It is demonstrated, that DWS in conjunction with MC can provide spectra with equivalent quality to those with common water suppression techniques. In addition to eliminating any influences of water-exchange, it is shown, that the information provided by the water signal can be used for motion-correction, not only in terms of easy zero-order phasing, spectral re-alignment and eddy-current correction, but also compensation for motion-related signal drop. Overall, this leads to improved spectral resolution, more stable fitting, and determination of ADC values that are closer to the true values, and hence better physiological interpretation.

Acknowledgements

Supported by the Swiss National Science Foundation

References

1. Valette, J., Guillermier, M., Besret, L., Boumezbeur, F., Hantraye, P., and Lebon, V. Optimized diffusion-weighted spectroscopy for measuring brain glutamate apparent diffusion coefficient on a whole-body MR system. NMR in Biomedicine, 18(8), 527–533. (2005).

2. Kan, H. E., Techawiboonwong, A., Van Osch, M. J. P., Versluis, M. J., Deelchand, D. K., Henry, P. G., Marjanska, M., Van Buchem, M., Webb, A. G., and Ronen, I. Differences in apparent diffusion coefficients of brain metabolites between grey and white matter in the human brain measured at 7 T. Magn. Reson. Med. 67:1203–1209 (2012).

3. Keating, B., and Ernst, T. Real-time dynamic frequency and shim correction for single-voxel magnetic resonance spectroscopy. Magn. Reson. Med. 68:1339–1345 (2012).

4. E. L. MacMillan, D. G. Q. Chong, W. Dreher, A. Henning, C. Boesch, and R. Kreis. Magnetization exchange with water and T1 relaxation of the downfield resonances in human brain spectra at 3.0 T. Magn Reson Med 65:1239-1246 (2011).

5. Brandejsky, V., Boesch, C., and Kreis, R. Proton diffusion tensor spectroscopy of metabolites in human muscle in vivo. Magn. Reson. Med. 73:481–487 (2014).

6. T. L. Hwang, P. C. van Zijl, and M. Garwood. Asymmetric adiabatic pulses for NH selection. J Magn Reson 138:173-177, (1999).

7. D. G. Chong, R. Kreis, C. S. Bolliger, C. Boesch, and J. Slotboom. Two-dimensional linear-combination model fitting of magnetic resonance spectra to define the macromolecule baseline using FiTAID, a Fitting Tool for Arrays of Interrelated Datasets. MAGMA 24:147-164 (2011).

Figures

Fig.1: The DWS STEAM sequence is optimized for minimization of TE to resolve metabolites even with small T2 and comprises an adiabatic FLAIR pulse with inversion time TI adjusted for CSF suppression and an adiabatic MC pulse placed into the mixing period TM.

Fig.2: (a) Typical results from a “braino” phantom showing simultaneously acquired metabolite and water spectra for 11 b-values. (b) and (c) present the fitted water/metabolite areas as function of b-value (single direction along room diagonal). Water and metabolites exhibit nearly mono-exponential decays over the whole b-value range indicating non-restricted diffusion.

Fig.3: (a) Non-water-suppressed diffusion-weighted spectra from human brain with/without correction for motion-related signaldrop based on the co-acquired water. (b) Effect of signal correction for the water peak intensity and lineshape (logarithmic scale). (c) In vivo diffusion-related signal decay for water (non-mono-exponential) and metabolites (mono-exponential, i.e. linear on the logarithmic scale).

Fig.4: ADC values obtained in vitro and in vivo in preliminary applications to 3 subjects. ADCs are in agreement with literature. Without correction, they are overestimated due to motion. The co-acquired water signal was used as internal reference in post-processing, improving. fit quality and accuracy of the determined ADC values.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
2395