Frequency and phase correction for multiplexed edited MRS of GABA and glutathione

Simultaneous frequency and phase corrections of single-shot MRS data using cross-correlation

PURPOSE

The objective of this study was to propose a novel preprocessing approach to simultaneously correct for the frequency and phase drifts in MRS data using cross-correlation technique.

METHODS

The performance of the proposed method was first investigated at different SNR levels using simulation. Random frequency and phase offsets were added to a previously acquired STEAM human data at 7 T, simulating two different noise levels with and without baseline artifacts. Alongside the proposed spectral cross-correlation (SC) method, three other simultaneous alignment approaches were evaluated. Validation was performed on human brain data at 3 T and mouse brain data at 16.4 T.

RESULTS

The results showed that the SC technique effectively corrects for both small and large frequency and phase drifts, even at low SNR levels. Furthermore, the mean square measurement error of the SC algorithm was comparable to the other three methods used, with much faster processing time. The efficacy of the proposed technique was successfully demonstrated in both human brain MRS data and in a noisy MRS dataset acquired from a small volume-of-interest in the mouse brain.

CONCLUSION

The study demonstrated the availability of a fast and robust technique that accurately corrects for both small and large frequency and phase shifts in MRS.

Simultaneous Measurement of GABA, Glutathione, and Glutamate-Glutamine in the Thalamus using Edited MR Spectroscopy: Feasibility and Applications in Traumatic Brain Injury

BACKGROUND

MR spectroscopy (MRS) is a noninvasive tool for evaluating biochemical alterations, such as glutamate (Glu)/gamma-aminobutyric acid (GABA) imbalance and depletion of antioxidative glutathione (GSH) after traumatic brain injury (TBI). Thalamus, a critical and vulnerable region post-TBI, is challenging for MRS acquisitions, necessitating optimization to simultaneously measure GABA/Glu and GSH.

PURPOSE

To assess the feasibility and optimize acquisition and processing approaches for simultaneously measuring GABA, Glx (Glu + glutamine (Gln)), and GSH in the thalamus, employing Hadamard encoding and reconstruction of MEscher-GArwood (MEGA)-edited spectroscopy (HERMES).

STUDY TYPE

Prospective.

SUBJECTS

28 control subjects (age: 35.9 ± 15.1 years), and 17 mild TBI (mTBI) patients (age: 32.4 ± 11.3 years).

FIELD STRENGTH/SEQUENCE

3T/T1-weighted magnetization-prepared rapid gradient-echo (MP-RAGE), HERMES.

ASSESSMENT

We evaluated the impact of acquisition with spatial saturation bands and post-processing with spectral alignment on HERMES performance in the thalamus among controls. Within-subject variability was examined in five controls through repeated scans within a week. The HERMES spectra in the posterior cingulate cortex (PCC) of controls were used as a reference for assessing HERMES performance in a reliable target. Furthermore, we compared metabolite levels and fitting quality in the thalamus between mTBI patients and controls.

STATISTICAL TESTS

Unpaired t-tests and within-subject coefficient-of-variation (CV). A P-value <0.05 was deemed significant.

RESULTS

HERMES spectra, acquired with saturation bands and processed with spectral alignment, yielded reliable metabolite measurements in the thalamus. The mean within-subject CV for GABA, Glx, and GSH levels were 18%, 10%, and 16% in the thalamus (7%, 9%, and 16% in the PCC). GABA (3.20 ± 0.60 vs 2.51 ± 0.55, P < 0.01) and Glx (8.69 ± 1.23 vs 7.72 ± 1.19, P = 0.03) levels in the thalamus were significantly higher in mTBI patients than in controls, with GSH (1.27 ± 0.35 vs 1.22 ± 0.28, P = 0.65) levels showing no significant difference.

DATA CONCLUSION

Simultaneous measuring GABA/Glx and GSH using HERMES is feasible in the thalamus, providing valuable insight into TBI.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 2.

GABA and glutamate measurements in temporal cortex of autistic children

Autism spectrum disorder (ASD) is a neurodevelopmental disorder and presents with challenges in social communication. A hypothesized underlying contributing mechanism is the imbalance in excitation and inhibition (E/I), partly influenced by the levels of excitatory neurotransmitter glutamate (Glu) and inhibitory neurotransmitter γ-aminobutyric acid (GABA) in the brain. Although many have reported the levels of GABA and Glu in the brain, only a few reports address the temporal cortex and then only with a small sample of autistic children, and often only in one hemisphere. We used a macromolecular suppressed edited-magnetic resonance spectroscopy (MRS) sequence to study GABA and Glu (as potential key players influencing E/I) in a large sample of children with ASD in the right and left temporal cortices of children with (N = 56) and without (N = 30) ASD (7-18 years). As a group, children with ASD exhibited no differences in the left hemisphere (GABA and Glu Cohen's |d|: 0.24 and 0.03), but the right hemisphere showed higher GABA and lower Glu concentrations (GABA and Glu Cohen's |d|: 0.53 and 0.65) compared to neurotypicals. Furthermore, a negative association was found between the right hemisphere Glu levels of the ASD group and a clinical assessment tool (r = -0.361, p = 0.022), reflecting autism trait severity (social responsiveness scale). In conclusion, we highlight the chemical abnormalities in children with ASD through a cross-sectional measurement. Longitudinal studies are warranted to determine whether these chemical levels persist or resolve over development.

Frontal and occipital brain glutathione levels are unchanged in autistic adults

BACKGROUND

The neurobiological underpinnings of Autism Spectrum Disorder (ASD) are diverse and likely multifactorial. One possible mechanism is increased oxidative stress leading to altered neurodevelopment and brain function. However, this hypothesis has mostly been tested in post-mortem studies. So far, available in vivo studies in autistic individuals have reported no differences in glutathione (GSH) levels in frontal, occipital, and subcortical regions. However, these studies were limited by the technically challenging quantification of GSH, the main brain antioxidant molecule. This study aimed to overcome previous studies' limitations by using a GSH-tailored spectroscopy sequence and optimised quantification methodology to provide clarity on GSH levels in autistic adults.

METHODS

We used spectral editing proton-magnetic resonance spectroscopy (1H-MRS) combined with linear combination model fitting to quantify GSH in the dorsomedial prefrontal cortex (DMPFC) and medial occipital cortex (mOCC) of autistic and non-autistic adults (male and female). We compared GSH levels between groups. We also examined correlations between GSH and current autism symptoms, measured using the Autism Quotient (AQ).

RESULTS

Data were available from 31 adult autistic participants (24 males, 7 females) and 40 non-autistic participants (21 males, 16 females); the largest sample to date. The GSH levels did not differ between groups in either region. No correlations with AQ were observed.

CONCLUSION

GSH levels as measured using 1H-MRS are unaltered in the DMPFC and mOCC regions of autistic adults, suggesting that oxidative stress in these cortical regions is not a marked neurobiological signature of ASD.

Impact of acquisition and modeling parameters on the test-retest reproducibility of edited GABA

Literature values vary widely for within-subject test-retest reproducibility of gamma-aminobutyric acid (GABA) measured with edited magnetic resonance spectroscopy (MRS). Reasons for this variation remain unclear. Here, we tested whether three acquisition parameters-(1) sequence complexity (two-experiment MEscher-GArwood Point RESolved Spectroscopy [MEGA-PRESS] vs. four-experiment Hadamard Encoding and Reconstruction of MEGA-Edited Spectroscopy [HERMES]); (2) editing pulse duration (14 vs. 20 ms); and (3) scanner frequency drift (interleaved water referencing [IWR] turned ON vs. OFF)-and two linear combination modeling variations-(1) three different coedited macromolecule models (called "1to1GABA", "1to1GABAsoft", and "3to2MM" in the Osprey software package); and (2) 0.55- versus 0.4-ppm spline baseline knot spacing-affected the within-subject coefficient of variation of GABA + macromolecules (GABA+). We collected edited MRS data from the dorsal anterior cingulate cortex from 20 participants (mean age: 30.8 ± 9.5 years; 10 males). Test and retest scans were separated by removing the participant from the scanner for 5-10 min. Each acquisition consisted of two MEGA-PRESS and two HERMES sequences with editing pulse durations of 14 and 20 ms (referred to here as MEGA-14, MEGA-20, HERMES-14, and HERMES-20; all TE = 80 ms, 224 averages). We identified the best test-retest reproducibility following postprocessing with a composite model of the 0.9- and 3-ppm macromolecules ("3to2MM"); this model performed particularly well for the HERMES data. Furthermore, sparser (0.55- compared with 0.4-ppm) spline baseline knot spacing yielded generally better test-retest reproducibility for GABA+. Replicating our prior results, linear combination modeling in Osprey compared with simple peak fitting in Gannet resulted in substantially better test-retest reproducibility. However, reproducibility did not consistently differ for MEGA-PRESS compared with HERMES, for 14- compared with 20-ms editing pulses, or for IWR-ON versus IWR-OFF. These results highlight the importance of model selection for edited MRS studies of GABA+, particularly for clinical studies that focus on individual patient differences in GABA+ or changes following an intervention.

A preliminary study of dynamic neurochemical changes in the dorsolateral prefrontal cortex during working memory

Working memory (WM) is one of the fundamental cognitive functions associated with the dorsolateral prefrontal cortex (DLPFC). However, the neurochemical mechanisms of WM, including the dynamic changes in neurometabolites such as glutamate and GABA in the DLPFC, remain unclear. Here, we investigated WM-related glutamate and GABA changes, alongside hemodynamic responses in the DLPFC, using a combination of functional magnetic resonance spectroscopy (fMRS) and functional magnetic resonance imaging (fMRI). During a WM task, we measured Glx (glutamate + glutamine) and GABA levels using GABA editing MEscher-GArwood Point REsolved Spectroscopy (MEGA-PRESS) sequence and blood-oxygen-level-dependent (BOLD) signal changes. In the DLPFC, we observed elevated Glx levels and increased BOLD signal changes during a 2-back task. Specifically, the Glx levels in the DLPFC were significantly higher during the 2-back task compared with fixation, although this difference was not significant when compared with a 0-back task. However, Glx levels during the 0-back task were higher than during fixation. Furthermore, there was a positive correlation between Glx levels in the DLPFC during the 2-back task and the corresponding BOLD signal changes. Notably, higher Glx increases were associated with increased DLPFC activation and lower WM task performance in individuals. No notable changes in DLPFC GABA levels were observed during WM processing. These findings suggest that the modulation of glutamatergic activity in the DLPFC may play a crucial role in both working memory processing and its associated performance outcomes.

Glutamate measurements using edited MRS

PURPOSE

To demonstrate J-difference coediting of glutamate using Hadamard encoding and reconstruction of Mescher-Garwood-edited spectroscopy (HERMES).

METHODS

Density-matrix simulations of HERMES (TE 80 ms) and 1D J-resolved (TE 31-229 ms) of glutamate (Glu), glutamine (Gln), γ-aminobutyric acid (GABA), and glutathione (GSH) were performed. HERMES comprised four sub-experiments with editing pulses applied as follows: (A) 1.9/4.56 ppm simultaneously (ONGABA /ONGSH ); (B) 1.9 ppm only (ONGABA /OFFGSH ); (C) 4.56 ppm only (OFFGABA /ONGSH ); and (D) 7.5 ppm (OFFGABA /OFFGSH ). Phantom HERMES and 1D J-resolved experiments of Glu were performed. Finally, in vivo HERMES (20-ms editing pulses) and 1D J-resolved (TE 31-229 ms) experiments were performed on 137 participants using 3 T MRI scanners. LCModel was used for quantification.

RESULTS

HERMES simulation and phantom experiments show a Glu-edited signal at 2.34 ppm in the Hadamard sum combination A+B+C+D with no overlapping Gln signal. The J-resolved simulations and phantom experiments show substantial TE modulation of the Glu and Gln signals across the TEs, whose average yields a well-resolved Glu signal closely matching the Glu-edited signal from the HERMES sum spectrum. In vivo quantification of Glu show that the two methods are highly correlated (p < 0.001) with a bias of ∼10%, along with similar between-subject coefficients of variation (HERMES/TE-averaged: ∼7.3%/∼6.9%). Other Hadamard combinations produce the expected GABA-edited (A+B-C-D) or GSH-edited (A-B+C-D) signal.

CONCLUSION

HERMES simulation and phantom experiments show the separation of Glu from Gln. In vivo HERMES experiments yield Glu (without Gln), GABA, and GSH in a single MRS scan.

Assessing the Dose-Dependent Effects of tDCS on Neurometabolites using Magnetic Resonance Spectroscopy

Concurrent transcranial direct current stimulation (tDCS) and proton Magnetic Resonance Spectroscopy ( ¹ H MRS) experiments have shown up- or downregulation of neurotransmitter concentration. However, effects have been modest applying mostly lower current doses and not all studies found significant effects. Dose of stimulation might be an important variable in eliciting a consistent response. To investigate dose effects of tDCS on neurometabolites, we placed an electrode over the left supraorbital region (with a return electrode over the right mastoid bone) and utilized an MRS voxel (3x3x3cm) that was centered over the anterior cingulate/inferior mesial prefrontal region which is in the path of the current distribution. We conducted 5 epochs of acquisition, each one with a 9:18min acquisition time, and applied tDCS in the third epoch. We observed significant dose and polarity dependent modulation of GABA and to a lesser degree of Glutamine/Glutamate (GLX) with the highest and reliable changes seen with the highest current dose, 5mA (current density 0.39 mA/cm ² ), during and after the stimulation epoch compared with pre-stimulation baselines. The strong effect on GABA concentration (achieving a mean change of 63% from baseline, more than twice as much as reported with lower doses of stimulation) establishes tDCS-dose as an important parameter in eliciting a regional brain engagement and response. Furthermore, our experimental design in examining tDCS parameters and effects using shorter epochs of acquisitions might constitute a framework to explore the tDCS parameter space further and establish measures of regional engagement by non-invasive brain-stimulation.

Ongoing oxidative stress in individuals with post-acute sequelae of COVID-19

Objectives

Coronavirus disease 2019 (COVID-19) caused by SARS-CoV-2 infection is associated with lower plasma glutathione (GSH) levels due to oxidative stress. However, plasma levels may not reflect brain GSH levels. Individuals with post-acute sequelae of COVID-19 (PASC) have a higher prevalence of cognitive fatigue, which might be related to altered brain γ-aminobutyric-acid (GABA) levels. Hence, our study aims to measure the brain GSH and GABA levels in PASC.

Methods

29 PASC participants and 24 uninfected controls were recruited for this study. Each was evaluated with detailed neuropsychiatric assessments and an edited proton MRS (Hadamard Encoding and Reconstruction of Mega-Edited Spectroscopy, HERMES) method to measure GABA and GSH concentrations in predominantly grey matter (GM) and predominantly white matter (WM) brain frontal voxels.

Results

PASC participants were 219 ± 137 days since their COVID-19 diagnosis. Nine individuals with PASC were hospitalized. Compared to controls, individuals with PASC had similar levels of GABA in both brain regions, but lower GSH and greater age-related GSH decline in the frontal GM region.

Conclusions

The lower-than-normal frontal GM GSH level in participants with PASC suggest that they have ongoing oxidative stress in the brain, and that older individuals may be even more vulnerable to oxidative stress.

Impact of acquisition and modeling parameters on test-retest reproducibility of edited GABA

Literature values for within-subject test-retest reproducibility of gamma-aminobutyric acid (GABA), measured with edited magnetic resonance spectroscopy (MRS), vary widely. Reasons for this variation remain unclear. Here we tested whether sequence complexity (two-experiment MEGA-PRESS versus four-experiment HERMES), editing pulse duration (14 versus 20 ms), scanner frequency drift (interleaved water referencing (IWR) turned ON versus OFF), and linear combination modeling variations (three different co-edited macromolecule models and 0.55 versus 0.4 ppm spline baseline knot spacing) affected the within-subject coefficient of variation of GABA + macromolecules (GABA+). We collected edited MRS data from the dorsal anterior cingulate cortex from 20 participants (30.8 ± 9.5 years; 10 males). Test and retest scans were separated by removing the participant from the scanner for 5-10 minutes. Each acquisition consisted of two MEGA-PRESS and two HERMES sequences with editing pulse durations of 14 and 20 ms (referred to here as: MEGA-14, MEGA-20, HERMES-14, and HERMES-20; all TE = 80 ms, 224 averages). Reproducibility did not consistently differ for MEGA-PRESS compared with HERMES or for 14 compared with 20 ms editing pulses. A composite model of the 0.9 and 3 ppm macromolecules (particularly for HERMES) and sparser (0.55 compared with 0.4 ppm) spline baseline knot spacing yielded generally better test-retest reproducibility for GABA+. Replicating our prior results, linear combination modeling in Osprey compared with simple peak fitting in Gannet resulted in substantially better test-retest reproducibility. These results highlight the importance of model selection for edited MRS studies of GABA+, particularly for clinical studies which focus on individual patient differences in GABA+ or changes following an intervention.

Connectivity-Guided Theta Burst Transcranial Magnetic Stimulation Versus Repetitive Transcranial Magnetic Stimulation for Treatment-Resistant Moderate to Severe Depression: Magnetic Resonance Imaging Protocol and SARS-CoV-2-Induced Changes for a Randomized Double-blind Controlled Trial

Background

Depression is a substantial health and economic burden. In approximately one-third of patients, depression is resistant to first-line treatment; therefore, it is essential to find alternative treatments. Transcranial magnetic stimulation (TMS) is a neuromodulatory treatment involving the application of magnetic pulses to the brain that is approved in the United Kingdom and the United States in treatment-resistant depression. This trial aims to compare the clinical effectiveness, cost-effectiveness, and mechanism of action of standard treatment repetitive TMS (rTMS) targeted at the F3 electroencephalogram site with a newer treatment—a type of TMS called theta burst stimulation (TBS) targeted based on measures of functional brain connectivity. This protocol outlines brain imaging acquisition and analysis for the Brain Imaging Guided Transcranial Magnetic Stimulation in Depression (BRIGhTMIND) study trial that is used to create personalized TMS targets and answer the proposed mechanistic hypotheses.

Objective

The aims of the imaging arm of the BRIGhTMIND study are to identify functional and neurochemical brain signatures indexing the treatment mechanisms of rTMS and connectivity-guided intermittent theta burst TMS and to identify imaging-based markers predicting response to treatment.

Methods

The study is a randomized double-blind controlled trial with 1:1 allocation to either 20 sessions of TBS or standard rTMS. Multimodal magnetic resonance imaging (MRI) is acquired for each participant at baseline (before TMS treatment) with T1-weighted and task-free functional MRI during rest used to estimate TMS targets. For participants enrolled in the mechanistic substudy, additional diffusion-weighted sequences are acquired at baseline and at posttreatment follow-up 16 weeks after treatment randomization. Core data sets of T1-weighted and task-free functional MRI during rest are acquired for all participants and are used to estimate TMS targets. Additional sequences of arterial spin labeling, magnetic resonance spectroscopy, and diffusion-weighted images are acquired depending on the recruitment site for mechanistic evaluation. Standard rTMS treatment is targeted at the F3 electrode site over the left dorsolateral prefrontal cortex, whereas TBS treatment is guided using the coordinate of peak effective connectivity from the right anterior insula to the left dorsolateral prefrontal cortex. Both treatment targets benefit from the level of MRI guidance, but only TBS is provided with precision targeting based on functional brain connectivity.

Results

Recruitment began in January 2019 and is ongoing. Data collection is expected to continue until January 2023.

Conclusions

This trial will determine the impact of precision MRI guidance on rTMS treatment and assess the neural mechanisms underlying this treatment in treatment-resistant depressed patients.

Trial Registration

ISRCTN Registry ISRCTN19674644; https://www.isrctn.com/ISRCTN19674644

International Registered Report Identifier (IRRID)

DERR1-10.2196/31925

Prediction learning in adults with autism and its molecular correlates

Background

According to Bayesian hypotheses, individuals with Autism Spectrum Disorder (ASD) have difficulties making accurate predictions about their environment. In particular, the mechanisms by which they assign precision to predictions or sensory inputs would be suboptimal in ASD. These mechanisms are thought to be mostly mediated by glutamate and GABA. Here, we aimed to shed light on prediction learning in ASD and on its neurobiological correlates.

Methods

Twenty-six neurotypical and 26 autistic adults participated in an associative learning task where they had to learn a probabilistic association between a tone and the rotation direction of two dots, in a volatile context. They also took part in magnetic resonance spectroscopy (MRS) measurements to quantify Glx (glutamate and glutamine), GABA + and glutathione in a low-level perceptual region (occipital cortex) and in a higher-level region involved in prediction learning (inferior frontal gyrus).

Results

Neurotypical and autistic adults had their percepts biased by their expectations, and this bias was smaller for individuals with a more atypical sensory sensitivity. Both groups were able to learn the association and to update their beliefs after a change in contingency. Interestingly, the percentage of correct predictions was correlated with the Glx/GABA + ratio in the occipital cortex (positive correlation) and in the right inferior frontal gyrus (negative correlation). In this region, MRS results also showed an increased concentration of Glx in the ASD group compared to the neurotypical group.

Limitations

We used a quite restrictive approach to select the MR spectra showing a good fit, which led to the exclusion of some MRS datasets and therefore to the reduction of the sample size for certain metabolites/regions.

Conclusions

Autistic adults appeared to have intact abilities to make predictions in this task, in contrast with the Bayesian hypotheses of ASD. Yet, higher ratios of Glx/GABA + in a frontal region were associated with decreased predictive abilities, and ASD individuals tended to have more Glx in this region. This neurobiological difference might contribute to suboptimal predictive mechanisms in ASD in certain contexts.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13229-021-00470-6.

In Vivo Brain GSH: MRS Methods and Clinical Applications

Glutathione (GSH) is an important antioxidant implicated in several physiological functions, including the oxidation−reduction reaction balance and brain antioxidant defense against endogenous and exogenous toxic agents. Altered brain GSH levels may reflect inflammatory processes associated with several neurologic disorders. An accurate and reliable estimation of cerebral GSH concentrations could give a clear and thorough understanding of its metabolism within the brain, thus providing a valuable benchmark for clinical applications. In this context, we aimed to provide an overview of the different magnetic resonance spectroscopy (MRS) technologies introduced for in vivo human brain GSH quantification both in healthy control (HC) volunteers and in subjects affected by different neurological disorders (e.g., brain tumors, and psychiatric and degenerative disorders). Additionally, we aimed to provide an exhaustive list of normal GSH concentrations within different brain areas. The definition of standard reference values for different brain areas could lead to a better interpretation of the altered GSH levels recorded in subjects with neurological disorders, with insights into the possible role of GSH as a biomarker and therapeutic target.

The trajectory of cortical GABA across the lifespan, an individual participant data meta-analysis of edited MRS studies

γ-Aminobutyric acid (GABA) is the principal inhibitory neurotransmitter in the human brain and can be measured with magnetic resonance spectroscopy (MRS). Conflicting accounts report decreases and increases in cortical GABA levels across the lifespan. This incompatibility may be an artifact of the size and age range of the samples utilized in these studies. No single study to date has included the entire lifespan. In this study, eight suitable datasets were integrated to generate a model of the trajectory of frontal GABA estimates (as reported through edited MRS; both expressed as ratios and in institutional units) across the lifespan. Data were fit using both a log-normal curve and a nonparametric spline as regression models using a multi-level Bayesian model utilizing the Stan language. Integrated data show that an asymmetric lifespan trajectory of frontal GABA measures involves an early period of increase, followed by a period of stability during early adulthood, with a gradual decrease during adulthood and aging that is described well by both spline and log-normal models. The information gained will provide a general framework to inform expectations of future studies based on the age of the population being studied.

Comparison of methods for spectral alignment and signal modelling of GABA-edited MR spectroscopy data

Many methods exist for aligning and quantifying magnetic resonance spectroscopy (MRS) data to measure in vivo γ-aminobutyric acid (GABA). Research comparing the performance of these methods is scarce partly due to the lack of ground-truth measurements. The concentration of GABA is approximately two times higher in grey matter than in white matter. Here we use the proportion of grey matter within the MRS voxel as a proxy for ground-truth GABA concentration to compare the performance of four spectral alignment methods (i.e., retrospective frequency and phase drift correction) and six GABA signal modelling methods. We analyse a diverse dataset of 432 MEGA-PRESS scans targeting multiple brain regions and find that alignment to the creatine (Cr) signal produces GABA+ estimates that account for approximately twice as much of the variance in grey matter as the next best performing alignment method. Further, Cr alignment was the most robust, producing the fewest outliers. By contrast, all signal modelling methods, except for the single-Lorentzian model, performed similarly well. Our results suggest that variability in performance is primarily caused by differences in the zero-order phase estimated by each alignment method, rather than frequency, resulting from first-order phase offsets within subspectra. These results provide support for Cr alignment as the optimal method of processing MEGA-PRESS to quantify GABA. However, more broadly, they demonstrate a method of benchmarking quantification of in vivo metabolite concentration from other MRS sequences.

Preprocessing, analysis and quantification in single-voxel magnetic resonance spectroscopy: experts' consensus recommendations

Once an MRS dataset has been acquired, several important steps must be taken to obtain the desired metabolite concentration measures. First, the data must be preprocessed to prepare them for analysis. Next, the intensity of the metabolite signal(s) of interest must be estimated. Finally, the measured metabolite signal intensities must be converted into scaled concentration units employing a quantitative reference signal to allow meaningful interpretation. In this paper, we review these three main steps in the post-acquisition workflow of a single-voxel MRS experiment (preprocessing, analysis and quantification) and provide recommendations for best practices at each step.

Spectral editing in 1 H magnetic resonance spectroscopy: Experts' consensus recommendations

Spectral editing in in vivo ¹ H-MRS provides an effective means to measure low-concentration metabolite signals that cannot be reliably measured by conventional MRS techniques due to signal overlap, for example, γ-aminobutyric acid, glutathione and D-2-hydroxyglutarate. Spectral editing strategies utilize known J-coupling relationships within the metabolite of interest to discriminate their resonances from overlying signals. This consensus recommendation paper provides a brief overview of commonly used homonuclear editing techniques and considerations for data acquisition, processing and quantification. Also, we have listed the experts' recommendations for minimum requirements to achieve adequate spectral editing and reliable quantification. These include selecting the right editing sequence, dealing with frequency drift, handling unwanted coedited resonances, spectral fitting of edited spectra, setting up multicenter clinical trials and recommending sequence parameters to be reported in publications.

GSH and GABA decreases in IDH1-mutated low-grade gliomas detected by HERMES spectral editing at 3 T in vivo

Isocitrate dehydrogenase 1 (IDH1) mutational status is an important prognostic biomarker in gliomas. γ-aminobutyric acid (GABA) and reduced glutathione (GSH) play an important role in energy production, which is related to tumor progression. Hadamard Encoding and Reconstruction of Mega-Edited Spectroscopy (HERMES) is able to detect GABA and GSH in healthy controls. This study aims to examine GABA and GSH alterations in IDH1-mutated low-grade gliomas using HERMES. We prospectively enrolled 14 suspected low-grade gliomas and 6 healthy control patients in this study, all cases underwent a 3 T MRI scan, including T1-weighted imaging and HERMES acquisition with a volume of interest 3 × 3 × 3 cm3. HERMES detects a "GABA+" signal that includes contributions from macromolecules and homocarnosine. GABA+ and GSH in tumor foci (group 1), contralateral cerebral regions (group 2) and healthy controls (group 3) were quantified using Gannet. The fitting errors and SNR of HERMES for GABA+ and GSH were analyzed; FWHM of the unsuppressed water signal was also recorded. The Wilcoxon signed-rank test was performed to test for differences between contralateral GABA+ and GSH levels, and differences in GABA+, GSH and fitting errors/SNR between the three groups were analyzed using analysis of variance (ANOVA). Eleven IDH1-mutant low-grade gliomas (5 Female and 6 Male, age 33-69) and 6 healthy subjects (2 Female and 4 Male, age 35-60) were finally enrolled this study. The mean water linewidth across all subjects was 9.67 ± 2.28 Hz. The Wilcoxon signed-rank test revealed that GABA+ and GSH were decreased significantly in glioma foci compared with contralateral regions, whereas no differences were seen between the left and right regions in healthy controls. ANOVA showed that GABA+ and GSH levels in tumor were lower than contralaterally and in healthy controls, while no differences were observed between the contralateral healthy tissue and healthy controls. No differences of fitting errors or SNR were found between tumors, contralateral regions or healthy controls. Our results suggest that HERMES is a reliable tool to simultaneously measure GABA and GSH alterations in low-grade gliomas with IDH1 mutations.

Motion correction in magnetic resonance spectroscopy

In vivo proton magnetic resonance spectroscopy and spectroscopic imaging (MRS/MRSI) are valuable tools to study normal and abnormal human brain physiology. However, they are sensitive to motion, due to strong crusher gradients, long acquisition times, reliance on high magnetic field homogeneity, and particular acquisition methods such as spectral editing. The effects of motion include incorrect spatial localization, phase fluctuations, incoherent averaging, line broadening, and ultimately quantitation errors. Several retrospective methods have been proposed to correct motion-related artifacts. Recent advances in hardware also allow prospective (real-time) correction of the effects of motion, including adjusting voxel location, center frequency, and magnetic field homogeneity. This article reviews prospective and retrospective methods available in the literature and their implications for clinical MRS/MRSI. In combination, these methods can attenuate or eliminate most motion-related artifacts and facilitate the acquisition of high-quality data in the clinical research setting.

Correcting frequency and phase offsets in MRS data using robust spectral registration

An algorithm for retrospective correction of frequency and phase offsets in MRS data is presented. The algorithm, termed robust spectral registration (rSR), contains a set of subroutines designed to robustly align individual transients in a given dataset even in cases of significant frequency and phase offsets or unstable lipid contamination and residual water signals. Data acquired by complex multiplexed editing approaches with distinct subspectral profiles are also accurately aligned. Automated removal of unstable lipid contamination and residual water signals is applied first, when needed. Frequency and phase offsets are corrected in the time domain by aligning each transient to a weighted average reference in a statistically optimal order using nonlinear least-squares optimization. The alignment of subspectra in edited datasets is performed using an approach that specifically targets subtraction artifacts in the frequency domain. Weighted averaging is then used for signal averaging to down-weight poorer-quality transients. Algorithm performance was assessed on one simulated and 67 in vivo pediatric GABA-/GSH-edited HERMES datasets and compared with the performance of a multistep correction method previously developed for aligning HERMES data. The performance of the novel approach was quantitatively assessed by comparing the estimated frequency/phase offsets against the known values for the simulated dataset or by examining the presence of subtraction artifacts in the in vivo data. Spectral quality was improved following robust alignment, especially in cases of significant spectral distortion. rSR reduced more subtraction artifacts than the multistep method in 64% of the GABA difference spectra and 75% of the GSH difference spectra. rSR overcomes the major challenges of frequency and phase correction.

Separate N-acetyl aspartyl glutamate, N-acetyl aspartate, aspartate, and glutamate quantification after pediatric mild traumatic brain injury in the acute phase

PURPOSE

To separately measure N-acetyl aspartul glutamate (NAAG), N-acetyl aspartate (NAA), aspartate (Asp), and glutamate (Glu) concentrations in white matter (WM) using J-editing techniques in patients with mild traumatic brain injury (mTBI) in the acute phase.

METHODS

Twenty-four patients with closed concussive head injury and 29 healthy volunteers were enrolled in the current study. For extended ¹ H MRS examination, patients and controls were equally divided into two subgroups. In subgroup 1 (12 patients/15 controls), NAAG and NAA concentrations were measured in WM separately with MEGA-PRESS (echo time/repetition time [TE/TR] = 140/2000 ms; δ ON NAA / δ OFF NAA = 4.84/4.38 ppm, δ ON NAAG / δ OFF NAAG = 4.61/4.15 ppm). In subgroup 2 (12 patients/14 controls), Asp and Glu concentrations were acquired with MEGA-PRESS (TE/TR = 90/2000 ms; δ ON Asp / δ OFF Asp = 3.89/5.21 ppm) and TE-averaged PRESS (TE from 35 ms to 185 ms with 2.5-ms increments; TR = 2000 ms) pulse sequences, respectively.

RESULTS

tNAA and NAAG concentrations were found to be reduced, while NAA concentrations were unchanged, after mild mTBI. Reduced Asp and elevated myo-inositol (mI) concentrations were also found.

CONCLUSION

The main finding of the study is that the tNAA signal reduction in WM after mTBI is associated with a decrease in the NAAG concentration rather than a decrease in the NAA concentration, as was thought previously. This finding highlights the importance of separating these signals, at least for WM studies, to avoid misinterpretation of the results. NAAG plays an important role in selectively activating mGluR3 receptors, thus providing neuroprotective and neuroreparative functions immediately after mTBI. NAAG shows potential for the development of new therapeutic strategies for patients with injuries of varying severity.

Effect of Age on GABA+ and Glutathione in a Pediatric Sample

BACKGROUND AND PURPOSE

Gamma-aminobutyric acid (GABA) is the primary inhibitory neurotransmitter in the human brain and is implicated in several neuropathologies. Glutathione is a major antioxidant in the brain and is considered a marker of oxidative stress. Several studies have reported age-related declines in GABA levels in adulthood, but the trajectory of both GABA and glutathione during childhood has not been well explored. The aim of this study is to establish how GABA and glutathione vary with age during early development.

MATERIALS AND METHODS

Twenty-three healthy children (5.6-13.9 years of age) were recruited for this study. MR imaging/MR spectroscopy experiments were conducted on a 3T MR scanner. A 27-mL MR spectroscopy voxel was positioned in the frontal lobe. J-difference edited MR spectroscopy was used to spectrally edit GABA and glutathione. Data were analyzed using the Gannet software, and GABA+ (GABA + macromolecules/homocarnosine) and glutathione were quantified using water (GABA+_H2O and Glutathione_H2O) and Cr (GABA+/Cr and glutathione/Cr) as concentration references. Also, the relative gray matter contribution to the voxel volume (GM_ratio) was estimated from structural images. Pearson correlation coefficients were used to examine the association between age and GABA+_H2O (and glutathione_H2O), between age and GABA+/Cr (and glutathione/Cr), and between age and GM_ratio.

RESULTS

Both GABA+_H2O (r = 0.63, P = .002) and GABA+/Cr (r = 0.48, P = .026) significantly correlated with age, whereas glutathione measurements and GM_ratio did not.

CONCLUSIONS

We demonstrate increases in GABA and no differences in glutathione with age in a healthy pediatric sample. This study provides insight into neuronal maturation in children and may facilitate better understanding of normative behavioral development and the pathophysiology of developmental disorders.

Simultaneous edited MRS of GABA, glutathione, and ethanol

The aim of this work was to develop simultaneous edited MRS of γ-aminobutyric acid (GABA), glutathione (GSH), and ethanol (EtOH) using Hadamard encoding and reconstruction of MEGA-edited spectroscopy (HERMES) at 3T. Density-matrix simulations of HERMES were carried out and compared with phantom experiments. In vivo experiments were performed in six healthy volunteers about 30 min after alcohol consumption. Simulations of HERMES showed GABA-, GSH-, and EtOH-edited spectra with low levels of crosstalk and excellent agreement with phantom spectra. In vivo experiments showed well edited GABA signals at 3.0 ppm, GSH at 2.95 ppm, and EtOH at 1.18 ppm in the respective Hadamard combination spectra. Measured integral ratios were 0.082 ± 0.012 for GABA/Cr, 0.037 ± 0.006 for GSH/Cr, and 0.305 ± 0.129 for EtOH/Cr. Simulated, phantom, and in vivo measurements of HERMES show excellent separation of GABA-, GSH-, and EtOH-edited signals with negligible levels of crosstalk. HERMES allows a threefold acceleration of editing while maintaining spectral quality compared with sequentially acquired MEGA-PRESS measurements.

An evaluation of the reproducibility of 1H-MRS GABA and GSH levels acquired in healthy volunteers with J-difference editing sequences at varying echo times

Recent advances in J-difference-edited proton magnetic resonance spectroscopy (¹H MRS) data acquisition and processing have led to the development of Hadamard Encoding and Reconstruction of MEGA-Edited Spectroscopy (HERMES) techniques, which enable the simultaneous measurement of ɣ-aminobutyric acid (GABA), the primary inhibitory amino acid neurotransmitter in the central nervous system, and of glutathione (GSH), the most abundant antioxidant in living tissue, at the commonly available magnetic field strength of 3 T. However, the reproducibility of brain levels of GABA and GSH measured across multiple scans in human subjects using HERMES remains to be established. In the present study, twelve healthy volunteers completed two consecutive HERMES scans of the dorsal anterior cingulate cortex (dACC) to assess the test-retest reproducibility of the technique for GABA and GSH measurements at TE = 80 ms. Eleven of the twelve participants additionally completed two consecutive MEGA-PRESS scans at TE = 120 ms, with editing pulses configured for GSH acquisition, to compare the reliability of GSH in the same voxel measured using the standard MEGA-PRESS at TE = 120 ms. The primary findings of study were that, 1) the coefficient of variation (CV) of measuring GABA with HERMES was 16.7%, which is in agreement with the reliability we previously reported for measuring GABA using MEGA-PRESS; and 2) the reliability of measuring GSH with MEGA-PRESS at TE = 120 ms was more than twice as high as that for measuring the antioxidant with HERMES at TE = 80 ms (CV = 7.3% vs. 19.0% respectively). These findings suggest that HERMES and MEGA-PRESS offer similar reliabilities for measuring GABA, while MEGA-PRESS at TE = 120 ms is more reliable for measuring GSH relative to HERMES at TE = 80 ms.

Simultaneous editing of GABA and GSH with Hadamard-encoded MR spectroscopic imaging

PURPOSE

To evaluate the feasibility of simultaneous MR spectroscopic imaging (MRSI) of gamma-aminobutyric acid (GABA) and glutathione (GSH) in the human brain using Hadamard Encoding and Reconstruction of MEGA-Edited Spectroscopy (HERMES).

METHODS

Point RESolved Spectroscopy (PRESS)-localized MRSI was performed in GABA and GSH phantoms and in the human brain (n = 3) using HERMES editing and compared to conventional MEGA editing of each metabolite. Multiplet patterns, signal intensities, and metabolite crosstalk were compared between methods. GABA+ and GSH levels were compared between methods for bias and variability. Linear regression of HERMES-MRSI GABA+/H₂ O and GSH/H₂ O versus gray matter (GM) fraction were performed to assess differences between GM and white matter (WM).

RESULTS

Phantom HERMES-MRSI scans gave comparable GABA+ and GSH signals to MEGA-MRSI across the PRESS-localized volume. In vivo, HERMES-reconstructed GABA+ and GSH values had minimal measurement bias and variability relative to MEGA-MRSI. Intersubject coefficients of variation (CV) from two regions within the PRESS-localized volume for HERMES and MEGA were 6-12% for GABA+ and 6-19% for GSH. Interregion CVs were 5-15% for GABA+ and 3-17% for GSH. The GABA+/H₂ O and GSH/H₂ O ratios were ~1.8 times higher and ~1.9 times higher, respectively, in GM than in WM.

CONCLUSION

HERMES-MRSI of GABA+ and GSH was found to be practical in the human brain with minimal measurement bias and comparable variability to separate MEGA-edited acquisitions of each metabolite performed in double the scan time. The HERMES-MRSI is a promising method for simultaneously mapping the distribution of multiple low-concentration metabolites.

Multi-vendor standardized sequence for edited magnetic resonance spectroscopy

Spectral editing allows direct measurement of low-concentration metabolites, such as GABA, glutathione (GSH) and lactate (Lac), relevant for understanding brain (patho)physiology. The most widely used spectral editing technique is MEGA-PRESS, which has been diversely implemented across research sites and vendors, resulting in variations in the final resolved edited signal. In this paper, we describe an effort to develop a new universal MEGA-PRESS sequence with HERMES functionality for the major MR vendor platforms with standardized RF pulse shapes, durations, amplitudes and timings. New RF pulses were generated for the universal sequence. Phantom experiments were conducted on Philips, Siemens, GE and Canon 3 T MRI scanners using 32-channel head coils. In vivo experiments were performed on the same six subjects on Philips and Siemens scanners, and on two additional subjects, one on GE and one on Canon scanners. On each platform, edited MRS experiments were conducted with the vendor-native and universal MEGA-PRESS sequences for GABA (TE = 68 ms) and Lac editing (TE = 140 ms). Additionally, HERMES for GABA and GSH was performed using the universal sequence at TE = 80 ms. The universal sequence improves inter-vendor similarity of GABA-edited and Lac-edited MEGA-PRESS spectra. The universal HERMES sequence yields both GABA- and GSH-edited spectra with negligible levels of crosstalk on all four platforms, and with strong agreement among vendors for both edited spectra. In vivo GABA+/Cr, Lac/Cr and GSH/Cr ratios showed relatively low variation between scanners using the universal sequence. In conclusion, phantom and in vivo experiments demonstrate successful implementation of the universal sequence across all four major vendors, allowing editing of several metabolites across a range of TEs.

Advanced Hadamard-encoded editing of seven low-concentration brain metabolites: Principles of HERCULES

PURPOSE

To demonstrate the framework of a novel Hadamard-encoded spectral editing approach for simultaneously detecting multiple low-concentration brain metabolites in vivo at 3T.

METHODS

HERCULES (Hadamard Editing Resolves Chemicals Using Linear-combination Estimation of Spectra) is a four-step Hadamard-encoded editing scheme. 20-ms editing pulses are applied at: (A) 4.58 and 1.9 ppm; (B) 4.18 and 1.9 ppm; (C) 4.58 ppm; and (D) 4.18 ppm. Edited signals from γ-aminobutyric acid (GABA), glutathione (GSH), ascorbate (Asc), N-acetylaspartate (NAA), N-acetylaspartylglutamate (NAAG), aspartate (Asp), lactate (Lac), and likely 2-hydroxyglutarate (2-HG) are separated with reduced signal overlap into distinct Hadamard combinations: (A+B+C+D); (A+B-C-D); and (A-B+C-D). HERCULES uses a novel multiplexed linear-combination modeling approach, fitting all three Hadamard combinations at the same time, maximizing the amount of information used for model parameter estimation, in order to quantify the levels of these compounds. Fitting also allows estimation of the levels of total choline (tCho), myo-inositol (Ins), glutamate (Glu), and glutamine (Gln). Quantitative HERCULES results were compared between two grey- and white-matter-rich brain regions (11 min acquisition time each) in 10 healthy volunteers. Coefficients of variation (CV) of quantified measurements from the HERCULES fitting approach were compared against those from a single-spectrum fitting approach, and against estimates from short-TE PRESS data.

RESULTS

HERCULES successfully segregates overlapping resonances into separate Hadamard combinations, allowing for the estimation of levels of seven coupled metabolites that would usually require a single 11-min editing experiment each. Metabolite levels and CVs agree well with published values. CVs of quantified measurements from the multiplexed HERCULES fitting approach outperform single-spectrum fitting and short-TE PRESS for most of the edited metabolites, performing only slightly to moderately worse than the fitting method that gives the lowest CVs for tCho, NAA, NAAG, and Asp.

CONCLUSION

HERCULES is a new experimental approach with the potential for simultaneous editing and multiplexed fitting of up to seven coupled low-concentration and six high-concentration metabolites within a single 11-min acquisition at 3T.

Testosterone is related to GABA+ levels in the posterior-cingulate in unmedicated depressed women during reproductive life

BACKGROUND

The role of testosterone (T) in the pathophysiology of affective disorders and anxiety is broadly supported. Evidence suggests that T has anxiolytic and antidepressant properties. One proposed route for the central effects of T is its interaction with the gamma-aminobutyric acid (GABA) system. We explored the relationship between T levels and GABA+ levels in anterior-cingulate (ACC) and the posterior-cingulate (PCC) regions in depressed women, using magnetic resonance spectroscopy (1H-MRS).

METHODS

Twenty-one depressed patients with regularly cycling who were not taking hormonal or psychotropic drugs were recruited. We assessed severity of depression using the Hamilton Depression Rating Scale (HDRS). Blood samples were taken for quantification of free (FT) and total testosterone (TT) on the day of the magnetic resonance (MR) scan. We evaluated GABA+ levels in the PCC and ACC, using the Hadamard Encoding and Reconstruction of MEGA-Edited Spectroscopy (HERMES) sequence. Pearson correlations were used to evaluate the association between FT, TT, GABA+ concentrations, and HDRS scores.

RESULTS

TT and FT levels were positively correlated with GABA+ levels in the PCC. No correlation was observed between T levels and GABA+ levels in the ACC. The HDRS total scores correlated negatively with FT levels.

LIMITATIONS

Limitations include the cross-sectional evaluation and the lack of a comparative healthy group.

CONCLUSIONS

Our findings suggest that the potential anxiolytic and antidepressant properties of T are related to increased GABA+ levels in the PCC. This observation may contribute to increased understanding of the role of T in depressive and anxiety symptoms in women.