Pushing the limits of speed and accuracy for 7T GABA MR spectroscopy to reveal GABA level fluctuations in resting brain
Arjan D. Hendriks1, Natalia Petridou1, Catalina S. Arteaga de Castro1, Mark W.J.M. Gosselink1, Alessio Fracasso1, and Dennis W.J. Klomp1

1Department of Radiology, University Medical Center Utrecht, Utrecht, Netherlands

Synopsis

By maximizing acquisition volume, using efficient semi LASER detection with editing and macro molecular nulling, GABA can be detected at high SNR within 3 minutes. GABA concentrations were measured in the visual cortex and in a phantom containing a known concentration of GABA. A bootstrapping method was used to determine the accuracy. The results indicate that the stability and fitting accuracy of the method is sufficient to detect concentration changes of GABA higher than 3% within a short scan time of less than 3 minutes. In resting state, GABA fluctuations up to 30% are found in in-vivo brain measurements.

Introduction

Knowledge about the main inhibitory neurotransmitter γ-aminobutyric acid (GABA) of the brain can give important insight in human brain function. This is valuable not only from a fundamental point of view, but also from a clinical perspective, since imbalances in excitatory and inhibitory processes are connected to several diseases like epilepsy, Parkinson, autoimmune inflammation and stroke [1].

Several studies have demonstrated significant changes in GABA levels between patient populations, however, with substantial variances of GABA levels within groups. In this work we analyzed the physiologic variations in GABA levels of the visual cortex in resting state using a fast and high SNR method for accurate detection of GABA.

Aim

The aim of this work is to assess the stability and fitting accuracy of the current fast GABA detection method to analyze variations in GABA concentration. This is done by bootstrapping the spectroscopy data and validation with a phantom.

Methods

Three healthy volunteers were scanned at a high field 7T system (Philips Healthcare, Cleveland, USA). For validation purposes, a phantom (ping-pong ball, with a diameter of 4 cm) containing water-dissolved GABA (6.6mM) and creatine (48mM) was constructed and scanned.

The set-up used was a half-cylinder open transmit coil with 8 transmit channels, combined with 2x16ch receive arrays covering the back of the head (Figure 1). Data acquisition was performed with a MEGA-sLASER [2] sequence with GABA editing, enhanced with macromolecule suppression and FOCI pulses for better B1 performance [3]. The sequence parameters were: single voxel, TE/TR= 74/5000 ms, in-vivo voxel size of 40x30x30 mm3, phantom voxel size of 17x20x17 mm3, spectral bandwidth: 4000 Hz, a total acquisition time of 2:46 min (16 averages) and 5:30 min (32 averages for bootstrapping). The spectroscopy volume inside the phantom was 6x smaller than the voxel volume in-vivo, due to phantom size limits. This is corrected, by selecting a GABA phantom concentration 6x larger than in-vivo concentration, maintaining a realistic SNR comparison between phantom and volunteers.

An in-house Matlab-based fitting model was used to fit the data. Bootstrapping was performed by selecting different sets of 16 odd-even pairs from a total of 32 pairs (32 averages). For each selected combination, the resulting spectrum was fitted. This was repeated 500 times.

Results

The spectral fitting of creatine and GABA is shown in Figure 2. The SNR and goodness of fit (assessed using the Cramer-Rao lower bounds) is calculated for increasing number of spectral averages. The SNR (Figure 2, middle column) increases with the number of averages. The goodness of fit (Figure 2, right column) stabilizes after 8 averages for creatine and after 16 averages for GABA. The time to acquire 16 averages is 2:46.

Bootstrapping shows the variation of multiple spectroscopy measurements in a phantom (Figure 3, left). The mean and standard deviation of the bootstrap for six different measurements is displayed. The bootstrapping shows a variation for creatine of less than 0.5%, which is similar the total system stability. GABA shows a variation of less than 3% both within bootstraps and between measurements. When measuring in volunteers (Figure 3, right), around 1% variation is found for creatine, which is a factor of 10 higher than the variation for the phantom data. For GABA, a variation within bootstraps is in the order of 8-12%. The difference in GABA concentration between individuals is up to 30%.

Discussion and conclusion

The SNR improves with the number of averages. The goodness of fit stabilizes for GABA at around 16 averages (2:46 min). For the phantom data, the variation of the GABA measurements is less than 3% within and between scans. However, for comparable measurements in volunteers a much larger variation is found. Also, the variation between consecutive measurements (10%-30%) is much larger than the variation within the scan (8-12%). This shows that fluctuations of GABA signals are substantially larger than those caused by intrinsic system noise and measurement instabilities. This implies that in the visual cortex at resting state, GABA physiology fluctuates up to 30%.

In conclusion, the results indicate that stability and fitting accuracy of the sensitivity optimized method employed here is sufficient to detect concentration changes of GABA in the visual cortex higher than 3%, within a fast scan time of less than 3 minutes. This brings us closer to more real-time measurements of GABA levels. The in-vivo measurements suggest that the large variation of GABA (up to 30% in this work) between individuals is mainly of physiologic origin. Moreover, individual GABA levels can vary within a temporal resolution of 3 minutes. This brings new insight in the interpretation of GABA spectroscopy measurements.

Acknowledgements

No acknowledgement found.

References

[1] M.J. Donahue and J. Near et al. 2010 NeuroImage; 53: 392–398

[2] A. Andreychenko et al. 2012 Magn Reson Med; 68:1018–1025

[3] C. S. Arteaga de Castro et al. 2013 NMR Biomed; 26: 1213–1219

[4] A.D. Hendriks et al. 2015 Proc. Intl. Soc. Mag. Reson. Med. 23. Abstract 3200

Figures

Figure 1: Overview of the half volume coil setup. The setup consists of an 8-channel transmit coil, 2x16ch receive coils and large screen for visual stimulation. Projection size is >60° visual angle, resulting in increased activation. Note how the large screen allows a big voxel of 36ml. Adapted from [4].

Figure 2: Spectral fitting is displayed (left) as well as the SNR (middle) and goodness of fit for different averages (right). Results are shown for GABA (top row) and creatine (bottom row). The goodness of fit (CRLB), stabilizes for the GABA fit after 16 averages, which are acquired in 2:46min.

Figure 3: Bootstrap of the fit for creatine (top) and GABA (bottom) for the phantom (left) and in-vivo (right). The mean, standard deviation and ratio (percentage) is shown. Despite the same intrinsic SNR of GABA, the in-vivo levels show up to 10-fold higher fluctuations (30%) than the phantom results (3%)



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