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-values
1,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
artifacts
3. 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
applications
4, 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 FoundationReferences
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).