We investigate the relationships between scan duration, signal-to-noise ratio (SNR) and group-level variance in GABA-edited MRS. Typically, GABA editing takes ~10 min for a 27-mL voxel. GABA+/Cr measurements from five voxels from 18 participants were analyzed by cumulatively binning the averages within each dataset to determine the effects on SNR and group-level variance. Sample size calculations estimated the required sample sizes needed for different predicted effect sizes in GABA-edited MRS studies. We show that the duration of GABA-edited acquisitions can be reduced if taking into account a statistically acceptable amount of group-level variance and the magnitudes of predicted effects.
Eighteen GABA-edited datasets were collected in vivo using MEGA-PRESS (3) on a 3T Philips Achieva scanner using a 32-channel head coil in voxels placed in the auditory cortex (AUD), dorsolateral prefrontal cortex (DLPFC), frontal eye field (FEF), occipital cortex (OCC) and sensorimotor cortex (SM). Voxels were 3×3×3 cm3, except the AUD voxel, which was 4×3×2 cm3. Acquisition parameters were: TE/TR = 68/2000 ms, 14-ms editing pulses (ON/OFF = 1.9/7.46 ppm), 320 averages (~10 min scan time), interleaving ON/OFF every 2 TRs, 2048 data points, 2 kHz spectral width.
Data were analyzed in Gannet (4). Due to excessive B0 field drift or participant head motion that could not be rectified with retrospective frequency/phase correction, 10 datasets were excluded from further analysis (2 AUD, 2 DLPFC, 1 FEF, 3 OCC, 2 SM). In the remaining datasets, up to 20 outlier averages were removed prior to averaging.
To test the SNR dependence of edited GABA measurements, for each dataset, difference spectra were generated by binning the first 4x averages successively, where 1 ≤ x ≤ 75. For consistency, only the first 300 (non-excluded) averages were considered for all datasets. For every difference spectrum, the GABA and Glx signals were fit with a three-Gaussian model and GABA+/Cr was quantified. SNR was measured as the ratio between the amplitude of the modeled GABA signal and twice the standard deviation of noise based on a recently reported algorithm (5). Fit error was calculated as the standard deviation of fit residuals normalized to modeled GABA signal amplitude.
Group-level variance in GABA+/Cr was measured as
the coefficient of variation (CV) at each signal-average bin. The relationship between number of averages and SNR was modeled using a
square-root function $$$a\sqrt{n}$$$. The relationship
between number of averages and fit error was modeled using an
inverse-square-root function $$$a\sqrt{n^{-1}}+b$$$,
where b accounts for the
error arising from imperfect modeling of the GABA signal. Based on the modeled
relationship between SNR and CV, sample sizes required for detecting group differences
of 10–30% in GABA+/Cr (power = 80%; alpha = 0.05) were calculated as a function
of SNR.
1. Mullins PG, McGonigle DJ, O’Gorman RL, Puts NAJ, Vidyasagar R, Evans CJ, Cardiff Symposium on MRS of GABA, Edden RAE. Current practice in the use of MEGA-PRESS spectroscopy for the detection of GABA. Neuroimage 2014;86:43–52. doi: 10.1016/j.neuroimage.2012.12.004.
2. Puts NAJ, Edden RAE. In vivo magnetic resonance spectroscopy of GABA: A methodological review. Prog. Nucl. Magn. Reson. Spectrosc. 2012;60:29–41. doi: 10.1016/j.pnmrs.2011.06.001.
3. Mescher M, Merkle H, Kirsch J, Garwood M, Gruetter R. Simultaneous in vivo spectral editing and water suppression. NMR Biomed. 1998;11:266–272. doi: 10.1002/(SICI)1099-1492(199810)11:6<266::AID-NBM530>3.0.CO;2-J.
4. Edden RAE, Puts NAJ, Harris AD, Barker PB, Evans CJ. Gannet: A batch-processing tool for the quantitative analysis of gamma-aminobutyric acid-edited MR spectroscopy spectra. J. Magn. Reson. Imaging 2014;40:1445–1452. doi: 10.1002/jmri.24478.
5. Mikkelsen M, Barker PB, Bhattacharyya PK, et al. Big GABA: Edited MR spectroscopy at 24 research sites. Neuroimage 2017;159:32–45. doi: 10.1016/j.neuroimage.2017.07.021.