Comparison of spectral fitting methods for overlapping J-coupled metabolite resonances

Universal dynamic fitting of magnetic resonance spectroscopy

PURPOSE

Dynamic (2D) MRS is a collection of techniques where acquisitions of spectra are repeated under varying experimental or physiological conditions. Dynamic MRS comprises a rich set of contrasts, including diffusion-weighted, relaxation-weighted, functional, edited, or hyperpolarized spectroscopy, leading to quantitative insights into multiple physiological or microstructural processes. Conventional approaches to dynamic MRS analysis ignore the shared information between spectra, and instead proceed by independently fitting noisy individual spectra before modeling temporal changes in the parameters. Here, we propose a universal dynamic MRS toolbox which allows simultaneous fitting of dynamic spectra of arbitrary type.

METHODS

A simple user-interface allows information to be shared and precisely modeled across spectra to make inferences on both spectral and dynamic processes. We demonstrate and thoroughly evaluate our approach in three types of dynamic MRS techniques. Simulations of functional and edited MRS are used to demonstrate the advantages of dynamic fitting.

RESULTS

Analysis of synthetic functional 1H-MRS data shows a marked decrease in parameter uncertainty as predicted by prior work. Analysis with our tool replicates the results of two previously published studies using the original in vivo functional and diffusion-weighted data. Finally, joint spectral fitting with diffusion orientation models is demonstrated in synthetic data.

CONCLUSION

A toolbox for generalized and universal fitting of dynamic, interrelated MR spectra has been released and validated. The toolbox is shared as a fully open-source software with comprehensive documentation, example data, and tutorials.

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.

Localized MRS reliability of in vivo glutamate at 3 T in shortened scan times: a feasibility study

Glutamate is the prime excitatory neurotransmitter in the mammalian brain and has been implicated in a wide range of psychiatric conditions. To improve the applicability and clinical reach of magnetic resonance spectroscopy (MRS), research is needed to develop shortened, yet reliable, MRS scanning procedures for standard 1.5-3-T clinical magnetic resonance imaging (MRI) systems, particularly with young or vulnerable populations unable to tolerate longer protocols. To this end, we evaluated the test-retest reliability of a shortened J-resolved MRS sequence in healthy adolescents (n = 22) aged 12-14 years. Participants underwent a series of sequential 6-min MRS scans, with the participants remaining in situ between successive scans. Glutamate and other metabolites were acquired from the rostral anterior cingulate cortex, as glutamatergic function in this region has been implicated in a number of psychiatric illnesses. Thirteen neurochemicals were quantified as ratios to total creatine, and reliability scores were expressed as the percentage difference between the two scans for each metabolite. Test-retest assessment of glutamate was reliable, as scores were less than 10% different (7.1 ± 4.2%), and glutamate values across scans were significantly correlated (Pearson r = 0.680, p < 10^-4 ). Several other neurochemicals demonstrated satisfactory reliability, including choline (Cho) (7.4 ± 5.6%), glutathione (GSH) (8.6 ± 4.1%), myo-inositol (mI) (6.5 ± 7.1%) and N-acetylaspartate (NAA) (3.5 ± 3.6%), with test-retest correlations ranging from 0.747 to 0.953. A number of metabolites, however, did not demonstrate acceptable test-retest reliability using the current J-resolved MRS sequence, ranging from 13.8 ± 13.7% (aspartate, Asp) to 45.9 ± 38.3% (glycine, Gly). Collectively, test-retest analyses suggest that clinically viable quantitative data can be obtained on standard MRI systems for glutamate, as well as the other metabolites, during short scan times in a traditionally challenging brain region.

Edited 1 H magnetic resonance spectroscopy in vivo: Methods and metabolites

The Proton magnetic resonance (1 H-MRS) spectrum contains information about the concentration of tissue metabolites within a predefined region of interest (a voxel). The conventional spectrum in some cases obscures information about less abundant metabolites due to limited separation and complex splitting of the metabolite peaks. One method to detect these metabolites is to reduce the complexity of the spectrum using editing. This review provides an overview of the one-dimensional editing methods available to interrogate these obscured metabolite peaks. These methods include sequence optimizations, echo-time averaging, J-difference editing methods (single BASING, dual BASING, and MEGA-PRESS), constant-time PRESS, and multiple quantum filtering. It then provides an overview of the brain metabolites whose detection can benefit from one or more of these editing approaches, including ascorbic acid, γ-aminobutyric acid, lactate, aspartate, N-acetyl aspartyl glutamate, 2-hydroxyglutarate, glutathione, glutamate, glycine, and serine. Magn Reson Med 77:1377-1389, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Regional distributions of brain glutamate and glutamine in normal subjects

Glutamate (Glu) and glutamine (Gln) play an important role in neuronal regulation and are of value as MRS-observable diagnostic biomarkers. In this study the relative concentrations of these metabolites have been measured in multiple regions in the normal brain using a short-TE whole-brain MRSI measurement at 3 T combined with a modified data analysis approach that used spatial averaging to obtain high-SNR spectra from atlas-registered anatomic regions or interest. By spectral fitting of high-SNR spectra this approach yielded reliable measurements across a wide volume of the brain. Spectral averaging also demonstrated increased SNR and improved fitting accuracy for the sum of Glu and Gln (Glx) compared with individual voxel fitting. Results in 26 healthy controls showed relatively constant Glu/Cr and Gln/Cr throughout the cerebrum, although with increased values in the anterior cingulum and paracentral lobule, and increased Gln/Cr in the superior motor area. The deep gray-matter regions of thalamus, putamen, and pallidum show lower Glu/Cr compared with cortical white-matter regions. Lobar measurements demonstrated reduced Glu/Cr and Gln/Cr in the cerebellum as compared with the cerebrum, where white-matter regions show significantly lower Glu/Cr and Gln/Cr as compared with gray-matter regions across multiple brain lobes. Regression analysis showed no significant effect of gender on Glu/Cr or Gln/Cr measurement; however, Glx/Cr ratio was found to be significantly negatively correlated with age in some lobar brain regions. In summary, this methodology provides the spectral quality necessary for reliable separation of Glu and Gln at 3 T from a single MRSI acquisition enabling generation of regional distributions of metabolites over a large volume of the brain, including cortical regions. Copyright © 2016 John Wiley & Sons, Ltd.

3D spatially encoded and accelerated TE-averaged echo planar spectroscopic imaging in healthy human brain

Several different pathologies, including many neurodegenerative disorders, affect the energy metabolism of the brain. Glutamate, a neurotransmitter in the brain, can be used as a biomarker to monitor these metabolic processes. One method that is capable of quantifying glutamate concentration reliably in several regions of the brain is TE-averaged (1) H spectroscopic imaging. However, this type of method requires the acquisition of multiple TE lines, resulting in long scan durations. The goal of this experiment was to use non-uniform sampling, compressed sensing reconstruction and an echo planar readout gradient to reduce the scan time by a factor of eight to acquire TE-averaged spectra in three spatial dimensions. Simulation of glutamate and glutamine showed that the 2.2-2.4 ppm spectral region contained 95% glutamate signal using the TE-averaged method. Peak integration of this spectral range and home-developed, prior-knowledge-based fitting were used for quantitation. Gray matter brain phantom measurements were acquired on a Siemens 3 T Trio scanner. Non-uniform sampling was applied retrospectively to these phantom measurements and quantitative results of glutamate with respect to creatine 3.0 (Glu/Cr) ratios showed a coefficient of variance of 16% for peak integration and 9% for peak fitting using eight-fold acceleration. In vivo scans of the human brain were acquired as well and five different brain regions were quantified using the prior-knowledge-based algorithm. Glu/Cr ratios from these regions agreed with previously reported results in the literature. The method described here, called accelerated TE-averaged echo planar spectroscopic imaging (TEA-EPSI), is a significant methodological advancement and may be a useful tool for categorizing glutamate changes in pathologies where affected brain regions are not known a priori. Copyright © 2016 John Wiley & Sons, Ltd.

Glutamate and glutamine: a review of in vivo MRS in the human brain

Our understanding of the roles that the amino acids glutamate (Glu) and glutamine (Gln) play in the mammalian central nervous system has increased rapidly in recent times. Many conditions are known to exhibit a disturbance in Glu-Gln equilibrium, and the exact relationships between these changed conditions and these amino acids are not fully understood. This has led to increased interest in Glu/Gln quantitation in the human brain in an array of conditions (e.g. mental illness, tumor, neuro-degeneration) as well as in normal brain function. Accordingly, this review has been undertaken to describe the increasing number of in vivo techniques available to study Glu and Gln separately, or pooled as 'Glx'. The present MRS methods used to assess Glu and Gln vary in approach, complexity, and outcome, thus the focus of this review is on a description of MRS acquisition approaches, and an indication of relative utility of each technique rather than brain pathologies associated with Glu and/or Gln perturbation. Consequently, this review focuses particularly on (1) one-dimensional (1)H MRS, (2) two-dimensional (1)H MRS, and (3) one-dimensional (13)C MRS techniques.

Magnetic resonance spectroscopy study of the glutamatergic system in adolescent males with high-functioning autistic disorder: a pilot study at 4T

The pilot study aimed at examining the neural glutamatergic activity in autism. Seven adolescent males (mean age: 14 ± 1.8; age range: 12-17 years) with intact intellectual capacity (mean IQ: 108 ± 14.26; IQ range: 85-127) suffering from autistic disorder and an equal number of age- and sex-matched healthy controls underwent a two-dimensional magnetic resonance spectroscopy scan at 4T. Results indicated significantly high glutamate (Glu) levels in the anterior cingulate cortex of autistic disorder versus control subjects (paired t test p = 0.01) and a trend for lower Glu in the right medial temporal lobe, which was not statistically different between the groups (paired t test p = 0.06). These preliminary findings support the glutamatergic dysregulation hypothesis in autism and need to be replicated in a larger sample.

The case of the missing glutamine

A theoretical study was performed to determine the accuracy and repeatability of multiple one-dimensional pulse sequences in the quantification of glutamine concentration at 3 T. Variable repeatability (12% to > 50%) and significant absolute error (-50% to +70%) were noted for the eight pulse sequences considered. Data acquired in vivo using three of the pulse sequences used for simulation matched the predicted repeatability well; among the pulse sequences considered, point-resolved spectroscopy (TE = 80 ms) offered minimal error and acceptable repeatability (12%) for brain glutamine measurements. Following correction for the expected bias of each pulse sequence, consistent glutamine measurements, in the 1-mM range, were reported with the three sequences. An explanation for the mismatch between in vivo (1)H MRS and in vitro (13)C/(1)H MRS at high field was attempted.