Fully automated macromolecule suppressed single voxel glutamate spectroscopy (FAMOUS SVGS)

Purpose The aim of the study was to develop and validate a new localized 1H MRS pulse sequence and automated post-processing software for the quantification of brain Glutamate (Glu) in clinical conditions at 7.0T in order to get reliable and reproducible results for acute intervention studies. Methods Here we describe a new localized proton MRS method “Fully Automated MacrOmolecUle Suppressed Single Voxel Glutamate Spectroscopy (FAMOUS SVGS)” for measuring Glu. FAMOUS SVGS method consists of a new pulse sequence with optimized switchable water, metabolites and outer volume suppression modules, as well as a frequency selective inversion pulse and automated post-processing of the five spectra obtained. FAMOUS SVGS method was first validated with glutamate phantoms and then validated with test–retest repeatability studies in the occipital cortex of five normal volunteers at 7.0T. Results Glutamate concentrations estimated from phantoms with FAMOUS SVGS method correlated well with actual concentrations. Test–retest repeatability studies in human brain in vivo yielded less than 0.3 mM intra-subject variations in Glu concentrations. Conclusions FAMOUS SVGS method enables Glu quantification in vivo at 7.0T with test–retest variability of less than 0.3 mM. We expect that we can reliably measure ≥0.5 mM change in glutamate due to any acute intervention.


Background
Glutamate (Glu) is a major excitatory neurotransmitter that is present in abundance in the synaptic vesicles of the human brain [1]. It plays an important role in the day-to-day functions such as cognitive, learning and memory, a key metabolite in cellular metabolism, and helps in removal of ammonia from the body. Its quantification has gained significant importance due to its implications in many of the neuropsychological as well as neurodegenerative disorders. In these brain disorders, as the changes in the Glu levels are expected to be relatively small, it is important to have a non-invasive method that measures these small changes in Glu concentration with high degree of accuracy and repeatability.
Single voxel proton magnetic resonance spectroscopy ( 1 H MRS) which is a safe, non-invasive method for the detection of human brain metabolites in vivo, has been widely applied for the study of Glu levels in many of the neurological disorders as well as in cases of insomnia and addiction disorders [2][3][4][5][6][7][8][9][10][11][12][13][14]. With the availability of 7.0T ultra-high field magnets, and the phased-array head coils, there is an improved resolution in the 1 H MRS and appreciable gain in signal-to-noise ratio (SNR) over the lower field strength magnets [15]. In a previous acute gabapentin drug study done at 7.0T, the effect of drug produced detectable change in GABA but not for glutamate [16]. Based on literature reports and from our own experience, we feel the following issues might contribute to the glutamate quantification errors: (i) Eddy currents, Nanga et al. J Transl Med (2016) 14:220 (ii) Phased-array combination from the multi-channel receive head coils, (iii) Chemical shift artifacts, (iv) Phase correction, (v) Baseline correction, and (vi) Macromolecular contaminations. There are few test-retest studies done at 7.0T, where mean inter-subject variations in glutamate estimation are reported to be around 5 % [17,18], but intra-subject test-retest variation for each individual in glutamate estimation is not reported. Our interest is in detecting the absolute quantification of glutamate for pre-and post acute interventions of an individual with an aim to detect ≥0.5 mM change in glutamate concentrations in clinical conditions such as temporal lobe epilepsy, cortical dysplasia, Alzheimer's, Psychiatric disorders, etc.
This has motivated us to develop a new strategy to design a new pulse sequence as well as develop an automated processing pipeline with no user intervention to mitigate the above-mentioned issues with the single voxel 1 H MR spectra, for improved repeatability of the Glu measurement.
In this study, we demonstrate a new method for quantification of the Glu at 2.35 ppm from the Fully Automated MacrOmolecUle Suppressed Single Voxel Glutamate Spectroscopy (FAMOUS SVGS). FAMOUS SVGS 1 H MRS with an automated post-processing provides a simple spectrum with a flat baseline. FAMOUS SVGS method was validated in vitro using glutamate monosodium phantoms. We investigated in vivo efficacy of this method by acquiring the data from the voxel in occipital cortex from multiple subjects, at two-time points.

FAMOUS SVGS sequence
FAMOUS SVGS pulse sequence is shown in Fig. 1. The sequence consists of a water reference acquisition, two control acquisitions with selective inversion set at 7.05 ppm ["water suppressed" (Control1) and "water and metabolite suppressed" (Control2)], and two selective inversion acquisitions set at 2.35 ppm ["water suppressed" (Selinv1) and "water and metabolite suppressed" (Selinv2)]. The number of averages for the water reference acquisition and water-suppressed acquisitions can be independently set. The basic sequence has a magnetization preparation and single voxel acquisition blocks. The preparation block consists of double inversion for metabolite suppression, water suppression, outer volume suppression (OVS) and frequency selective inversion (FSI). The acquisition block consists of the standard PRESS sequence. The double inversion segment is designed with broad band adiabatic inversion pulses that do not excite water (pulse width = 18 ms, bandwidth = 3000 Hz (10 ppm centered at −1 ppm). The inversion delays are optimized to provide good inversion for metabolite T 1 values between 1 and 2 s. Water suppression is performed using seven water suppression pulses (α1-α7); Gaussian shape with user-defined bandwidth and seven inter-pulse delays (d1-d7). Delay d4 is designed to be long enough (>20 ms) to accommodate up to 4 OVS pulses (3.5 ms radio frequency (RF) pulse (bandwidth = 5600 Hz) and 1.3 ms gradient spoiler). Delay d7 is designed to be long enough (>45 ms) to accommodate both the OVS and FSI pulses. The FSI pulse is a Gaussian pulse with its width calculated to provide 50 % inversion over a user defined inversion bandwidth (typically 60 Hz to cover the distance between the center peak of the triplet to the either side (7 Hz + 7 Hz) and an additional 45 Hz for shim inhomogeneities). The water reference acquisition is performed with all spoiler gradients and OVS RF pulses without any water suppression, double inversion or FSI pulses. This is essential since we use the water reference spectrum to correct eddy current induced dynamic phase shifts for all the five spectra. "Water suppressed" spectra are acquired with water suppression, OVS and FSI pulses. "Water and metabolite suppressed" spectra are obtained with double inversion, water suppression, OVS and FSI pulses. The FSI chemical shift (invppm) is user selectable (for Glu, invppm is 2.35 ppm) and the sequence automatically calculates the offset frequencies for the control and Selinv acquisitions to be symmetrically placed around water resonance (at 4.7 + 2.35 ppm) and at 2.35 ppm respectively. Inversion at (4.7 + 2.35 ppm) is done to remove any residual baseline water difference. For the SVS PRESS block, custom designed Shinnar-Le Roux (SLR) RF pulses, were used for uniform excitation and refocusing profile [19]. Relevant SLR pulse design parameters are: (a) excitation pulse: pulse width = 2.7 ms, time bandwidth product = 9, pass band ripple = 1 %, stop band ripple = 0.2 %, flip angle = 90°, and (b) refocusing pulse: pulse width = 6 ms, time bandwidth product = 12, pass band ripple = 1 %, stop band ripple = 0.1 %, flip angle = 150°. In the PRESS part of the sequence, for the water acquisition the excitation and refocusing pulses were set on resonance with water peak at 4.7 ppm while for all "water-suppressed" and "water and metabolite" spectra, the excitation and refocusing pulses were set at the glutamate chemical shift of 2.35 ppm to avoid chemical shift artifact for glutamate. This chemical shift based offset was added to the gradient strength and positionbased frequency offset for all PRESS pulses. This ensures that the signal for water reference, Control1, Control2, Selinv1 and Selinv2 spectra all originate from the same location. The scan protocol uses a smaller number of averages for the water reference acquisition block and a larger number of averages for the water-suppressed and water-suppressed metabolite nulled acquisition blocks to optimize scan time. Nanga et al. J Transl Med (2016) 14:220

MR scan protocol
All the single voxel 1 H MRS experiments were performed on a 7.0T whole body MRI scanner (Siemens Medical Systems, Erlangen, Germany) with a vendor supplied volume coil transmit/32-Channel receive proton head phased array coil. The acquisition protocol consisted of a gradient echo localizer, a T 1 -weighted three-dimensional (3D) magnetization prepared rapid gradient echo (MPRAGE) of whole brain. Localized automated first order and second order shimming of the B 0 field was performed on the voxel of interest in order to obtain a water line width of <15 Hz (0.05 ppm) or less using FASTMAP shim method [20,21] provided by Siemens as work in progress (WIP) package. Optimal reference voltage for the region of interest was calculated by acquiring two water spectra with reference voltages of 100 and 200 V using stimulated echo acquisition mode (STEAM) sequence. Assuming the peak amplitudes of water signal to be S1 and S2 for the above two reference voltages, the optimal reference voltage "B1" for the region of interest was calculated using the following equation: Single voxel 1 H MRS was obtained with FAMOUS SVGS sequence.

Phantom studies
We verified the accuracy of the FAMOUS SVGS method using the glutamate monosodium phantoms. A stock solution of 100 mM Glutamate monosodium phantom was prepared using PBS buffer. From this stock solution, four phantoms were prepared at physiological pH to give final concentrations of 6, 9, 12, 15 mM. Even though with phantoms we could optimize shims to get line widths <5 Hz, to mimic the in vivo conditions, shim currents were mis-set to obtain a line width comparable to that of in vivo (15 Hz). Single voxel 1 H MR spectra were obtained on the phantom with FAMOUS SVGS at room temperature using following parameters: voxel-size = 20 × 30 × 20 mm 3 , repetition time

Human studies Ethics, consent and permission
All the human studies were conducted under an approved University of Pennsylvania Institutional Review Board protocol with written informed consent obtained from each volunteer after explaining the study protocol. Voxel was positioned in the occipital cortex (OCC, 20 × 30 × 20 mm 3 ) with saturation band covering the skull area and the spectra was obtained with the following parameters: repetition time (TR) = 3000 ms and time to echo (TE) = 21 ms, number of points = 4096, sweep width = 6000 Hz, dummy scans = 4, averages = 8 (water), 32 (Control 1), 32 (Selinv1), 32 (Control 2), 32 (Selinv2).
Total acquisition time to obtain the spectra was 7 min:48 s. All the five volunteers (all males, aged 19-62 years) were scanned twice within day for repeatability studies.

Voxel placement for studies involving multiple scans
We used our in-house developed software, "Imscribe", described elsewhere [22] to maintain the consistency in voxel placement for multiple acquisitions of 1 H MRS from the same voxel.
Imscribe (accessible at: cmroi.med.upenn.edu/ imscribe), uses high resolution T 1 -weighted images and the voxel information from the spectroscopy file acquired from first scan as target template for subsequent scans and performs affine co-registration giving the information for the new voxel placement as illustrated in Fig. 2.

MRS data processing
For spectroscopy post-processing, we have an in-house developed automated post-processing package using MATLAB (R2014a, MathWorks, Natick, MA). Postprocessing was done on the raw data from the system before any phased array channel combination and the flow chart is depicted in Fig. 3. For the phantom spectra, the FAMOUS Glu spectrum was obtained as (Control 1-Selinv1) and fitted as Lorentzian. Water content was assumed to be 100 % for phantoms.
In-house developed spectral processing method for FAMOUS SVGS works in the frequency domain. In this method, an exponential apodization of ~7 Hz was applied to the final time domain data after taking care of eddy current induced phase shifts, followed by automated Fourier transformation, phase correction and macromolecule removal. For fitting of water reference spectrum, we used non-linear least squares sense fitting (MATLAB "lsqcurvefit" routine). The integral of the Glu peak from 2.25 to 2.45 ppm (based on the inversion pulse bandwidth of 60 Hz) was normalized by water reference integral for absolute quantification of Glu. Water content was assumed to be 80 % for human brains.

Statistics
Repeatability intra-class correlation coefficient (ICC) with 95 % confidence intervals (CI) values of single measures for the test-retest data was calculated with one-way random model using the SPSS software [version 18]. ICC value of <0.0 is interpreted as poor agreement, 0.0-0.2 indicates slight agreement, 0.21-0.4 indicates fair agreement, 0.41-0.6 indicates moderate agreement, 0.61-0.8 indicates substantial agreement, whereas ICC >0.81 indicates an almost perfect agreement [23]. ICC with p value less than 0.05 is considered as statistically significant.

Results
The concentration of glutamate from the voxel positioned on the glutamate monosodium phantom, as shown in Fig. 4, from the FAMOUS SVGS shows a good linear correlation between the actual and measured glutamate concentrations (slope = 1.03, goodness of fit as measured by R 2 = 0.997). Figure 5 shows the typical water spectrum obtained from the volunteer in OCC voxel and the typical Con-trol1 and Control2 spectra. One can clearly see a 140fold reduction in water signal in the later two spectra and metabolite suppression has very little effect on the quality of water suppression. Figure 6 shows the OCC voxel position and the typical Control1, Control2, (Control1-Selinv1), (Control2-Selinv2) and FAMOUS Glu spectra from a volunteer. In the inset is shown the FAMOUS Glu spectrum along with the inversion profile of the FSI pulse. FAMOUS Glu spectrum is integrated from 2.25 to 2.45 ppm, which is based on the inversion bandwidth of the FSI pulse as shown in inset of Fig. 6.   Fig. 4 Water-suppressed control spectrum (a) and water-suppressed inversion spectrum (b) along with the FAMOUS Glu spectrum (c) obtained on region of interest shown on glutamate monosodium phantom image (d) along with the FAMOUS Glu spectra obtained from different concentrations of glutamate monosodium phantom (e) and a plot of actual vs measured Glu concentration (f) from FAMOUS Glu spectra FAMOUS Glu spectra of test-retest of all the five subjects for each scan are shown in Fig. 7. The mean and standard deviation of Glu concentration from FAMOUS Glu spectra for the ten scans obtained from the five volunteers in the OCC voxel is 9.57 ± 0.63 mM; [CoV = 6.6 % calculated as (standard deviation/mean) × 100].
Absolute changes in Glu variations for each volunteer on OCC voxel are summarized in Table 1. The maximal Fig. 5 Typical water reference spectrum from one of the volunteer on OCC voxel along with the Control1 (black) and Control2 (grey) spectra demonstrating a good water suppression in the later two spectra. This suppression of water peak was 140-fold

Discussion
We have successfully demonstrated the feasibility of FAMOUS SVGS with studies on glutamate phantoms at 7.0T. The major advantage of this method is the fully automated processing with no user intervention for either phasing or baseline correction. We have also shown that we can estimate Glu concentrations in vivo from the data obtained from healthy subjects by FAMOUS SVGS with an intra-subject variability of less than 0.3 mM. This is most likely due to the simplification of the spectrum with negligible baseline and reduction of lipid/macromolecular contamination. We do find that the Glu concentrations in the occipital cortex from the FAMOUS SVGS are comparable to the previously published results from human brain autopsy study [24].
One of our goals in the current study was to bring down the scan-to-scan variation of Glu measurements to  Although the FAMOUS SVGS method seems to be simplistic with the incorporation of a frequency selective inversion pulse, one might raise the following concern: (i) Effect of T 1 recovery of Glu due to delay between inversion pulse and PRESS sequence: The reported T 1 values of the metabolites at 7.0T were on the order of 1.21-1.83 s [25]. Therefore, the T 1 -weighting of the MR signals from the 12 ms delay (half the inversion pulse width = 10 ms and gradient crusher = 2 ms) between the inversion pulse and the PRESS sequence is negligible (<1 %) (ii) Macromolecule/lipid contamination at 2.35 ppm: This has been accounted in our method.
For the voxels of interest where local B 0 inhomogeneity is large, we expect the accuracy of quantification to deteriorate due to poor shimming. Additionally, this method is highly sensitive to subject motion and hence we needed to immobilize the head with proper padding.
An alternate strategy we have considered was the use of long echo time without any inversion to reduce the macromolecule/lipid content from the control in vivo spectra. But the complicated coupling pattern of the glutamate peak poses additional quantification issues due to which we preferred to go along with the short echo spectroscopy with FAMOUS SVGS.
Although we have demonstrated the performance of FAMOUS SVGS in combination with in-house developed automated post-processing method at 7.0T, it can be easily implemented on animal imaging scanners ≥7.0T.

Conclusions
In conclusion, FAMOUS SVGS can be one of the useful 1 H MRS tools in any studies that measure regional glutamate changes in brain. The significant improvement in repeatability of the measurement has implications in quantifying relatively smaller change in Glu, reducing the number of subjects and time points needed in a clinical trial and thereby contributing to significant savings in study cost and time. Authors' contributions RPRN contributed to conception and design of the study, performed the experiments and wrote the manuscript. HH contributed to coding the pulse sequence, manuscript writing and editing. RR provided conception and contributed to manuscript editing. All authors read and approved the final manuscript.