0491

Simultaneous ultra-short TE-MRS in two voxels using a SPECIAL sequence with Hadamard encoding
Masoumeh Dehghani1,2, Richard Edden 3, and Jamie Near1,2
1McGill University, Montreal, QC, Canada, 2Centre d'Imagerie Cérébrale, Montreal, QC, Canada, 3Johns Hopkins University, Baltimore, MD, United States

Synopsis

The spin echo, full Intensity acquired localized (SPECIAL) sequence consists of a localized spin-echo, preceded by an alternating ISIS pre-inversion for voxel localization. In this study we modified the SPECIAL sequence to simultaneously localize the signal at two different positions. The technique relies on a four-step inversion scheme involving two different inversion positions, followed by a Hadamard encoding reconstruction. Comparing the in vivo performance of dual-SPECIAL sequence to the conventional SPECIAL sequence demonstrated that the dual-Special sequence provides simultaneous metabolite profile from two different regions, reducing the acquisition time by a factor of two, and without any penalty in SNR.

Introduction

The 1H spin echo, full Intensity acquired localized (SPECIAL) sequence1, is one commonly used acquisition technique for localized MR spectroscopy in the brain. Like the STEAM2 sequence, SPECIAL allows acquisition of spectra at ultra-short echo times, but with the additional benefit of a two-fold increase in the signal to noise ratio (SNR). The SPECIAL sequence consists of a localized spin-echo (SE), preceded by an alternating ISIS pre-inversion module. The ISIS subtraction scheme, in which inversion-on scans are subtracted from inversion-off scans, is used for processing the raw SPECIAL data to yields fully localized spectra3. In this study we modified the SPECIAL sequence to simultaneously localize the signal at two different positions. The modified sequence, denoted as dual-SPECIAL, involves a four-step multiplexed inversion scheme with two different inversion positions, followed by a Hadamard reconstruction4 to yield localized short-TE spectra from two voxels in a single acquisition.

Methods

As with the conventional SPECIAL sequence, each acquisition of the dual-SPECIAL sequence consists of a slice-selective SE that is localized to a column of tissue. Localization of two regions was achieved by introducing two different inversion pulses, both perpendicular to the SE column. The two inversion locations correspond to the two voxel positions (voxel 1 and voxel 2) and were arranged in a sequence of four experiments in the dual-SPECIAL sequence (Fig.1). The first scan (OFF/OFF) is acquired from the entire SE column within the field of view of the RF receiver/transmitter coil. The other three experiments incorporated a slice‐selective full‐passage adiabatic inversion pulse followed by a spoiling gradient prior to the SE sequence. Scan ON/OFF is acquired using a single selective inversion pulse centered at the position of first voxel, scan OFF/ON is acquired using a single selective inversion pulse centered at the position of second voxel, scan ON/ON is acquired using a dual selective inversion pulse and two bands centered on the positions of first and second voxels, respectively. The performance of the dual-SPECIAL sequence was investigated in vivo and compared with the conventional SPECIAL sequence. One healthy female volunteer provided informed consent to participate in this study. Experiments were performed on a 3T Siemens Prisma MR scanner with a commercial body transmit-volume coil and 32-channel receive array. The SE column was oriented in the scanner y-direction (anterior-posterior). Bo field inhomogeneities were minimized within a section of the column measuring 2.0 × 10.0 × 2.0 cm3 (encompassing both regions of interest) using the GRE-shim procedure, resulting in a water linewidth of 6.3 Hz and 5.3 in posterior and anterior regions (Fig.2.a). Localized water suppressed 1H spectra were acquired using the dual-SPECIAL sequence and SPECIAL sequence from two equal-size voxels (2.0 × 2.0 × 2.0 cm3), one placed in the medial prefrontal cortex and one placed in the medial posterior cingulate cortex/precuneus posterior regions of VOI and separated by 8.0 cm. All acquisition parameters including TR/TE = 4000/8.5 ms and 128 averages were kept same between the two sequences. Spectral pre-processing steps, including coil combination, phase and frequency correction and averaging were performed in MATLAB using the FID-A toolkit5. To achieve the desired dual localization, the four subspectra from dual-Special sequence were aligned and combined using a Hadamard encoding approach (Fig. 1.b). The signal from posterior and anterior regions is acquired by linear combination of (OFF/OFF)+(ON/OFF)-(ON/ON)-(OFF/ON) and (OFF/OFF)+(OFF/ON)-(ON/ON)-(ON/OFF), respectively. The two subspectra in Special sequence were combined based on the 1D ISIS add‐subtract scheme3.

Results

Figure 1 illustrates the four-step inversion scheme and Hadamard‐encoding approach used for combination of subspectra. Figure 2. shows the VOI for shimming along with posterior and anterior voxels selected for spectra acquisition, and the four subspectra acquired using the dual-SPECIAL sequence. Figure 3 shows the localized spectra obtained from posterior and anterior regions by dual-SPECIAL sequence and SPECIAL sequence. In both anterior and posterior regions, the spectral appearance is very similar between the SPECIAL and dual-SPECIAL sequences. Moreover, the SNR of spectra acquired using dual-SPECIAL sequence is similar to the SPECIAL sequence, ~190 for posterior and ~170 for anterior regions.

Discussion

Using the described dual-SPECIAL sequence and Hadamard encoding approach, the simultaneous signals from anterior and posterior regions are achieved by coaddition of signal from within the voxels and cancellation of signals from all external regions. Compared to the conventional SPECIAL sequence, the quality of the spectra, linewidth and SNR, in dual-SPECIAL sequence is similar for both anterior and posterior regions while the total measurement time for acquisition of both voxels is reduced by half.

Conclusion

In summary, simultaneous measurement of signals from two regions of interest was shown to be feasible using the developed dual-SPECIAL sequence combined with a Hadamard reconstruction. The dual-SPECIAL sequence reduces the total measurement time by half relative to sequential measurements of each voxel using the SPECIAL sequence without any penalty in SNR. These preliminary results suggest that the newly developed dual-SPECIAL sequence can be a potential candidate for simultaneous measuring of metabolites profiles from two different regions in clinical and pre-clinical applications such as functional MRS or 13C MRS.

Acknowledgements

This work is supported by the National Engineering and Sciences Research Council (NSERC, RGPIN-2014-06072, J.N.), the Canadian Institutes of Health Research (PJT-165896, J.N.) and the Fonds de recherche du Quebec – Santé (FRQS, J.N.).

References

1. Mlynárik V, Gambarota G, Frenkel H, Gruetter R. Localized short-echo-time proton MR spectroscopy with full signal-intensity acquisition. Magnetic Resonance in Medicine. 2006; 56(5):965-970.

2. Frahm J, Merboldt K-D, Hänicke W. Localized proton spectroscopy using stimulated echoes. Journal of Magnetic Resonance. 1987; 72(3):502-508.

3. Ordidge R.J, Connelly A, Lohman J.A.B. Image-selected in Vivo spectroscopy (ISIS). A new technique for spatially selective nmr spectroscopy. Journal of Magnetic Resonance.1986; 66(2):283-294.

4. Hadamard, J.S. Résolution d’une question relative aux determinants. Bulletin des Sciences Mathematiques. 1893; 17: 240-246.

5. Simpson R, Devenyi GA, Jezzard P, Hennessy TJ, Near J. Advanced processing and simulation of MRS data using the FID appliance (FID-A)-An open source, MATLAB-based toolkit. Magnetic Resonance in Medicine. 2017; 77(1):23-33.

Figures

Figure 1. (a) Scheme of the dual-SPECIAL pulse sequence in four steps: The 180° pulses before SE are adiabatic hyperbolic (HS4, R =10, 5ms) and have B1 of 10 µT for single inversion (experiment B&C) and 20 µT for dual inversion (Experiment D). WS represents an additional water suppression pulse. (b) Schematic illustration of Hadamard‐encoding applied for reconstruction of signal from two voxels in dual-SPECIAL sequence.

Figure 2. Top: Three-dimensional T1-weighted anatomical image was acquired using MPRAGE. Localized VOI (2.0 × 10.0 × 2.0 cm3) was extended from posterior toward interior region of brain (yellow box). The red and green boxes show posterior and interior regions selected (2.0 × 1.0 × 2.0 cm3) for 1H MRS acquisition and separated by 8.0 cm. Bottom: Four different 1H MR subspectra acquired in vivo from the human brain at 3T using dual-SPECIAL sequence. Acquisition parameters are TR/TE = 4000/8.5 ms and 32×4 averages.

Figure 3. Localized 1H spectrum acquired from posterior and anterior region of brain by (a) dual-SPECIAL sequence and (b) SPECIAL sequence.

Proc. Intl. Soc. Mag. Reson. Med. 28 (2020)
0491