T2* quantitation with chemical shift and multi-echo spiral imaging
Atsushi M Takahashi1

1McGovern Institute for Brain Research, Massachusetts Institute of Technology, Cambridge, MA, United States

Synopsis

Quantitation of T2* with spiral imaging sequences can be made in two distinct ways: Collecting data at various echo-times results in a measurement of the chemical shift after Fourier transformation along the echo-time dimension. Off resonance is intrinsically corrected by this processing. Alternately, multiple-echo spiral imaging can be used to quantitatively measure T2* as long as dephasing from B0 distortion is small over the duration of the spiral readout. Multiple spiral interleaves are used to reduce the readout time of the spiral. Both methods are demonstrated.

Introduction

Quantitation of various NMR properties of tissues elucidates information about the state of the underlying tissues. We would like to quantitatively measure T2* to attempt to understand the mechanisms underlying the BOLD effect and to tune and optimize T2* weighted fMRI sequences. Certain regions of the brain suffer from susceptibility artifacts from the nasal sinuses and ear canals. Certain neuro-biologically interesting regions of the brain are difficult, if not impossible to image with BOLD fMRI because of the "susceptibility holes" in pre-frontal cortex and along the inferior side of the temporal lobes.

Although echo planar imaging remains the work horse in the neuroscience community, other fast imaging techniques, like spiral, are frequently employed. Spiral k-space trajectories have been used to collect data in a variety of applications for a number of years.1,2 It was quickly noticed that spiral k-space trajectories resulted in localized blurring of the images in the presence of B0 inhomogeneities. Noll3,4 et al. found ways to correct this blurring with a number of method including multi-frequency reconstructing. Adalsteinsson5 et al. collected data at varying echo times to sample the FIDs and used an FFT in the time dimension to sample the chemical shift spectrum. Sarkar6 et al. refined the reconstruction to reduce side lobes from truncations of the data in echo-time and spatial k-space domains.

Methods

We implemented a multi-echo, multi-interleaf spiral pulse sequence in the IDEA pulse programming environment on our Siemens 3T Tim Trio MRI scanner (Siemens Medical Solutions). Two methods for measuring T2* will be discussed. A slow but reliable way of measuring T2* is to make repeated measurement of the FID but instead of identical measurements, the echo-time is stepped by a fixed step size, typically 500us. After Fourier transformation, this results in a spectral window of 1000Hz. Our pulse sequence allows us to collect these data with single shot or multi-interleaf spiral k-space trajectories. Hardware limits, fields-of-view, resolution and maximum data acquisition time, dictate the choice of parameters used to generate the trajectories at runtime. This type of spiral scan can be run with a single or multiple echoes but we typically collect these with only a single readout window. By making the step size much shorter than the readout duration, the blurring effects of off-resonance can be significantly reduced. It should be noted that intra-voxel dephasing due to the susceptibility gradients will still result in local signal loss.

The alternate method for measuring T2* is to collect multiple echoes. The number of echoes and times between echoes can be adjusted by using multiple interleaves. Increasing the number of interleaves shortens each readout period and the echo spacing can be reduced. The disadvantage of doing this is we violate Nyquist sampling in the echo-time dimension but a field map and multi-frequency reconstruction can be used to reduce blurring from off-resonance. Again, intra-voxel dephasing cannot be corrected and results in signal loss. However, shorter spiral readout times do reduce this signal loss.

Results

Fig. 1 shows the results of a 25 minute single shot spiral scan that was acquired by incrementing the echo time by 500usec for each of the 256 frames. The spiral was a 2D single shot spiral with a 192mm FOV, a resolution of 64 by 64 with a slice thickness of 3mm. The slew-rate was limited to 120 mT/m/msec to avoid peripheral nerve stimulation and resulted in a slew-rate limited spiral with a maximum gradient of only 25 mT/m. Fig. 2 shows one of the 256 frames. Figure 3 shows the data after Fourier transformation and the off-resonance blurring is inherently corrected.

Figures 4 and 5 show the fitted M0 and T2* from 1mm isotropic six-echo 32-interleaf spiral data acquired of the whole head.

Discussion and Conclusions

Various trajectories are possible for rapidly collecting data in k-space. Spiral has long been a reliable and robust method. We have explored a few variations of spiral imaging to make quantitative measurements of the brain in vivo in the presence of susceptibility gradients. Slow and tedious measurements are possible and future work will reduce the number of acquired echo-time steps.

Multi-echo, multi-shot imaging is also possible and interleaving the spirals reduces the time required for each readout. Readout time should be kept short to reduce the dephasing from off resonance when the echo time is not stepped. M0 and T2* were fitted to the multi-echo data but are not immediately comparable to the complete spectral measurements because T2* depends on the local susceptibility gradients across the voxel. Voxel size affects T2*.

Acknowledgements

The author wishes to acknowledge the the support of the McGovern Institute for Brain Research and the Athinoula A. Martinos Imaging Center at MIT for supporting this work. Technical support was provided by Siemens Medical Solutions.

References

1. CB Ahn, JH Kim, ZH Cho. High Speed Spiral-Scan Echo Planar NMR imaging. IEEE Trans Med Imaging 5(1):2-7, 1986

2. CH Meyer, BS Hu, DG Nishimura, A Macovski. Fast spiral coronary artery imaging. Magn Reson Med. 28(2):202-13, 1992

3. DC Noll, CH Meyer, JM Pauly, DG Nishimura, A Macovski. A homogeneity correction method for magnetic resonance imaging with time-varying gradients. IEEE Trans Med Imaging 10(4):629-37, 1991

4. DC Noll, JM Pauly, CH Meyer, DG Nishimura,A Macovski. Deblurring for non-2D Fourier transform magnetic resonance imaging. Magn Reson Med. 25(2):319-33, 1992

5. E Adalsteinsson, P Irarrazabal, DM Spielman, A Macovski. Three-dimensional spectroscopic imaging with time-varying gradients. Magn Reson Med. 33(4):461-6, 1995

6. S Sarkar, K Heberlein, X Hu. Truncation artifact reduction in spectroscopic imaging using a dual-density spiral k-space trajectory. Magn Reson Imaging. 20(10):743-57, 2002

Figures

Fig. 1. These are signal intensity values from a single voxel of a 3mm isotropic 2D multi-slice acquisition with a 256 point chemical shift image. Echo-times range from 2.3 to 129.3msec. A repetition time of 6 seconds was required to fit the 33 slices and the longest echo time delay. A single exponential function was fit to the data. These data are displayed without Fourier transformation in time but with multi-frequency reconstruction to correct blurring from B0 inhomogeneity.

Fig. 2. Images of frame 10 (TE=6.8msec) of the 256 TE data set.

Fig. 3. After Fourier transformation, we obtain a spectrum. Note that widths of the peaks vary inversely with T2*. There is also variation of B0 at various locations across the field of view as noted by the centroid of the peaks in the spectra. The underlying image is the height of the peaks of the spectra.

Fig. 4. Six echo, 1mm isotropic resolution spiral data was use to fit T2*. The repetition time was only 150msec therefore the longer T2* of CSF and other fluids was not reliably measured. Echo times of 5, 25, 45, 65, 85 and 105msec were acquired. The spiral trajectory had 32 interleaves, a readout time of 8.412msec and was restricted to a slew rate of 120 mT/m/msec and a maximum gradient of 33 mT/m.

Fig. 5. M0 map calculated from the same high resolution six echo data.



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