Simultaneous mapping of metabolites and individual macromolecular components via ultra-short acquisition delay 1 H MRSI in the brain at 7T

Evidence for biexponential glutamate T2 relaxation in human visual cortex at 3T: A functional MRS study

Functional magnetic resonance spectroscopy (fMRS) measures dynamic changes in metabolite concentration in response to neural stimulation. The biophysical basis of these changes remains unclear. One hypothesis suggests that an increase or decrease in the glutamate signal detected by fMRS could be due to neurotransmitter movements between cellular compartments with different T2 relaxation times. Previous studies reporting glutamate (Glu) T2 values have generally sampled at echo times (TEs) within the range of 30-450 ms, which is not adequate to observe a component with short T2 (<20 ms). Here, we acquire MRS measurements for Glu, (t) total creatine (tCr) and total N-acetylaspartate (tNAA) from the visual cortex in 14 healthy participants at a range of TE values between 9.3-280 ms during short blocks (64 s) of flickering checkerboards and rest to examine both the short- and long-T2 components of the curve. We fit monoexponential and biexponential Glu, tCr and tNAA T2 relaxation curves for rest and stimulation and use Akaike information criterion to assess best model fit. We also include power calculations for detection of a 2% shift of Glu between compartments for each TE. Using pooled data over all participants at rest, we observed a short Glu T2-component with T2 = 10 ms and volume fraction of 0.35, a short tCr T2-component with T2 = 26 ms and volume fraction of 0.25 and a short tNAA T2-component around 15 ms with volume fraction of 0.34. No statistically significant change in Glu, tCr and tNAA signal during stimulation was detected at any TE. The volume fractions of short-T2 component between rest and active conditions were not statistically different. This study provides evidence for a short T2-component for Glu, tCr and tNAA but no evidence to support the hypothesis of task-related changes in glutamate distribution between short and long T2 compartments.

A Comparison of 7 Tesla MR Spectroscopic Imaging and 3 Tesla MR Fingerprinting for Tumor Localization in Glioma Patients

This paper investigated the correlation between magnetic resonance spectroscopic imaging (MRSI) and magnetic resonance fingerprinting (MRF) in glioma patients by comparing neuro-oncological markers obtained from MRSI to T1/T2 maps from MRF. Data from 12 consenting patients with gliomas were analyzed by defining hotspots for T1, T2, and various metabolic ratios, and comparing them using Sørensen-Dice similarity coefficients (DSCs) and the distances between their centers of intensity (COIDs). The median DSCs between MRF and the tumor segmentation were 0.73 (T1) and 0.79 (T2). The DSCs between MRSI and MRF were the highest for Gln/tNAA (T1: 0.75, T2: 0.80, tumor: 0.78), followed by Gly/tNAA (T1: 0.57, T2: 0.62, tumor: 0.54) and tCho/tNAA (T1: 0.61, T2: 0.58, tumor: 0.45). The median values in the tumor hotspot were T1 = 1724 ms, T2 = 86 ms, Gln/tNAA = 0.61, Gly/tNAA = 0.28, Ins/tNAA = 1.15, and tCho/tNAA = 0.48, and, in the peritumoral region, were T1 = 1756 ms, T2 = 102 ms, Gln/tNAA = 0.38, Gly/tNAA = 0.20, Ins/tNAA = 1.06, and tCho/tNAA = 0.38, and, in the NAWM, were T1 = 950 ms, T2 = 43 ms, Gln/tNAA = 0.16, Gly/tNAA = 0.07, Ins/tNAA = 0.54, and tCho/tNAA = 0.20. The results of this study constitute the first comparison of 7T MRSI and 3T MRF, showing a good correspondence between these methods.

Macromolecule modelling for improved metabolite quantification using short echo time brain 1 H MRS at 3 T and 7 T: The PRaMM Model

Purpose

To improve reliability of metabolite quantification at both, 3 T and 7 T, we propose a novel parametrized macromolecules quantification model (PRaMM) for brain 1 H MRS, in which the ratios of macromolecule peak intensities are used as soft constraints.

Methods

Full- and metabolite-nulled spectra were acquired in three different brain regions with different ratios of grey and white matter from six healthy volunteers, at both 3 T and 7 T. Metabolite-nulled spectra were used to identify highly correlated macromolecular signal contributions and estimate the ratios of their intensities. These ratios were then used as soft constraints in the proposed PRaMM model for quantification of full spectra. The PRaMM model was validated by comparison with a single component macromolecule model and a macromolecule subtraction technique. Moreover, the influence of the PRaMM model on the repeatability and reproducibility compared to those other methods was investigated.

Results

The developed PRaMM model performed better than the two other approaches in all three investigated brain regions. Several estimates of metabolite concentration and their Cramér-Rao lower bounds were affected by the PRaMM model reproducibility, and repeatability of the achieved concentrations were tested by evaluating the method on a second repeated acquisitions dataset. While the observed effects on both metrics were not significant, the fit quality metrics were improved for the PRaMM method (p≤0.0001). Minimally detectable changes are in the range 0.5 - 1.9 mM and percent coefficients of variations are lower than 10% for almost all the clinically relevant metabolites. Furthermore, potential overparameterization was ruled out.

Conclusion

Here, the PRaMM model, a method for an improved quantification of metabolites was developed, and a method to investigate the role of the MM background and its individual components from a clinical perspective is proposed.

Quantification of short echo time MRS signals with improved version of QUantitation based on quantum ESTimation algorithm

Magnetic resonance spectroscopy offers information about metabolite changes in the organism, which can be used in diagnosis. While short echo time proton spectra exhibit more distinguishable metabolites compared with proton spectra acquired with long echo times, their quantification (and providing estimates of metabolite concentrations) is more challenging. They are hampered by a background signal, which originates mainly from macromolecules (MM) and mobile lipids. An improved version of the quantification algorithm QUantitation based on quantum ESTimation (QUEST), with MM prior knowledge (QUEST-MM), dedicated to proton signals and invoking appropriate prior knowledge on MM, is proposed and tested. From a single acquisition, it enables better metabolite quantification, automatic estimation of the background, and additional automatic quantification of MM components, thus improving its applicability in the clinic. The proposed algorithm may facilitate studies that involve patients with pathological MM in the brain. QUEST-MM and three QUEST-based strategies for quantifying short echo time signals are compared in terms of bias-variance trade-off and Cramér-Rao lower bound estimates. The performances of the methods are evaluated through extensive Monte Carlo studies. In particular, the histograms of the metabolite and MM amplitude distributions demonstrate the performances of the estimators. They showed that QUEST-MM works better than QUEST (Subtract approach) and is a good alternative to QUEST when measured MM signal is unavailable or unsuitable. Quantification with QUEST-MM is shown for 1 H in vivo rat brain signals obtained with the SPECIAL pulse sequence at 9.4 T, and human brain signals obtained, respectively, with STEAM at 4 T and PRESS at 3 T. QUEST-MM is implemented in jMRUI and will be available for public use from version 7.1.

Meta-analysis and open-source database for in vivo brain Magnetic Resonance spectroscopy in health and disease

Proton (1H) Magnetic Resonance Spectroscopy (MRS) is a non-invasive tool capable of quantifying brain metabolite concentrations in vivo. Prioritization of standardization and accessibility in the field has led to the development of universal pulse sequences, methodological consensus recommendations, and the development of open-source analysis software packages. One on-going challenge is methodological validation with ground-truth data. As ground-truths are rarely available for in vivo measurements, data simulations have become an important tool. The diverse literature of metabolite measurements has made it challenging to define ranges to be used within simulations. Especially for the development of deep learning and machine learning algorithms, simulations must be able to produce accurate spectra capturing all the nuances of in vivo data. Therefore, we sought to determine the physiological ranges and relaxation rates of brain metabolites which can be used both in data simulations and as reference estimates. Using the Preferred Reporting Items for Systematic reviews and Meta-Analyses (PRISMA) guidelines, we've identified relevant MRS research articles and created an open-source database containing methods, results, and other article information as a resource. Using this database, expectation values and ranges for metabolite concentrations and T2 relaxation times are established based upon a meta-analyses of healthy and diseased brains.

High field brain proton magnetic resonance spectroscopy and volumetry in children with chronic, compensated liver disease - A pilot study

BACKGROUND

and rationale: There is increasing evidence that children or young adults having acquired liver disease in childhood display neurocognitive impairment which may become more apparent as they grow older. The molecular, cellular and morphological underpinnings of this clinical problem are incompletely understood.

AIM

Therefore, we used the advantages of highly-resolved proton magnetic resonance spectroscopy at ultra-high magnetic field to analyze the neurometabolic profile and brain morphometry of children with chronic, compensated liver disease, hypothesizing that with high field spectroscopy we would identify early evidence of rising brain glutamine and decreased myoinositol, such as has been described both in animals and humans with more significant liver disease.

METHODS

Patients (n = 5) and age-matched controls (n = 19) underwent 7T MR scans and short echo time ¹H MR spectra were acquired using the semi-adiabatic SPECIAL sequence in two voxels located in gray and white matter dominated prefrontal cortex, respectively. A 3D MP2RAGE sequence was also acquired for brain volumetry and T₁ mapping. Liver disease had to have developed at least 6 months before entering the study. Subjects underwent routine blood analysis and neurocognitive testing using validated methods within 3 months of MRI and MRS.

RESULTS

Five children currently aged 8-16 years with liver disease acquired in childhood were included. Baseline biological characteristics were similar among patients. There were no statistically significant differences between subjects and controls in brain metabolite levels or brain volumetry. Finally, there were minor neurocognitive fluctuations including attention deficit in one child, but none fell in the statistically significant range.

CONCLUSION

Children with chronic, compensated liver disease did not display an abnormal neurometabolic profile, neurocognitive abnormalities, or signal intensity changes in the globus pallidus. Despite the absence of neurometabolic changes, it is an opportunity to emphasize that it is only by developing the use of ¹H MRS at high field in the clinical arena that we will understand the significance and generalizability of these findings in children with CLD. Attention deficit was observed in one child. Healthy children displayed neurometabolic regional differences as previously reported in adult subjects.

Meta-analysis and Open-source Database for In Vivo Brain Magnetic Resonance Spectroscopy in Health and Disease

Proton ( 1 H) Magnetic Resonance Spectroscopy (MRS) is a non-invasive tool capable of quantifying brain metabolite concentrations in vivo . Prioritization of standardization and accessibility in the field has led to the development of universal pulse sequences, methodological consensus recommendations, and the development of open-source analysis software packages. One on-going challenge is methodological validation with ground-truth data. As ground-truths are rarely available for in vivo measurements, data simulations have become an important tool. The diverse literature of metabolite measurements has made it challenging to define ranges to be used within simulations. Especially for the development of deep learning and machine learning algorithms, simulations must be able to produce accurate spectra capturing all the nuances of in vivo data. Therefore, we sought to determine the physiological ranges and relaxation rates of brain metabolites which can be used both in data simulations and as reference estimates. Using the Preferred Reporting Items for Systematic reviews and Meta-Analyses (PRISMA) guidelines, we've identified relevant MRS research articles and created an open-source database containing methods, results, and other article information as a resource. Using this database, expectation values and ranges for metabolite concentrations and T 2 relaxation times are established based upon a meta-analyses of healthy and diseased brains.

Feasibility and implications of using subject-specific macromolecular spectra to model short echo time magnetic resonance spectroscopy data

Expert consensus recommends linear-combination modeling (LCM) of 1 H MR spectra with sequence-specific simulated metabolite basis function and experimentally derived macromolecular (MM) basis functions. Measured MM basis functions are usually derived from metabolite-nulled spectra averaged across a small cohort. The use of subject-specific instead of cohort-averaged measured MM basis functions has not been studied widely. Furthermore, measured MM basis functions are not widely available to non-expert users, who commonly rely on parameterized MM signals internally simulated by LCM software. To investigate the impact of the choice of MM modeling, this study, therefore, compares metabolite level estimates between different MM modeling strategies (cohort-mean measured; subject-specific measured; parameterized) in a lifespan cohort and characterizes its impact on metabolite-age associations. 100 conventional (TE = 30 ms) and metabolite-nulled (TI = 650 ms) PRESS datasets, acquired from the medial parietal lobe in a lifespan cohort (20-70 years of age), were analyzed in Osprey. Short-TE spectra were modeled in Osprey using six different strategies to consider the MM baseline. Fully tissue- and relaxation-corrected metabolite levels were compared between MM strategies. Model performance was evaluated by model residuals, the Akaike information criterion (AIC), and the impact on metabolite-age associations. The choice of MM strategy had a significant impact on the mean metabolite level estimates and no major impact on variance. Correlation analysis revealed moderate-to-strong agreement between different MM strategies (r > 0.6). The lowest relative model residuals and AIC values were found for the cohort-mean measured MM. Metabolite-age associations were consistently found for two major singlet signals (total creatine (tCr])and total choline (tCho)) for all MM strategies; however, findings for metabolites that are less distinguishable from the background signals associations depended on the MM strategy. A variance partition analysis indicated that up to 44% of the total variance was related to the choice of MM strategy. Additionally, the variance partition analysis reproduced the metabolite-age association for tCr and tCho found in the simpler correlation analysis. In summary, the inclusion of a single high signal-to-noise ratio MM basis function (cohort-mean) in the short-TE LCM leads to more lower model residuals and AIC values compared with MM strategies with more degrees of freedom (Gaussian parametrization) or subject-specific MM information. Integration of multiple LCM analyses into a single statistical model potentially allows to identify the robustness in the detection of underlying effects (e.g., metabolite vs. age), reduces algorithm-based bias, and estimates algorithm-related variance.

Variability by region and method in human brain sodium concentrations estimated by 23Na magnetic resonance imaging: a meta-analysis

Sodium imaging (23Na-MRI) is of interest in neurological conditions given potential sensitivity to the physiological and metabolic status of tissues. Benchmarks have so far been restricted to parenchyma or grey/white matter (GM/WM). We investigate (1) the availability of evidence, (2) regional pooled estimates and (3) variability attributable to region/methodology. MEDLINE literature search for tissue sodium concentration (TSC) measured in specified 'healthy' brain regions returned 127 reports, plus 278 retrieved from bibliographies. 28 studies met inclusion criteria, including 400 individuals. Reporting variability led to nested data structure, so we used multilevel meta-analysis and a random effects model to pool effect sizes. The pooled mean from 141 TSC estimates was 40.51 mM (95% CI 37.59-43.44; p < 0.001, I2Total=99.4%). Tissue as a moderator was significant (F214 = 65.34, p-val < .01). Six sub-regional pooled means with requisite statistical power were derived. We were unable to consider most methodological and demographic factors sought because of non-reporting, but each factor included beyond tissue improved model fit. Significant residual heterogeneity remained. The current estimates provide an empirical point of departure for better understanding in 23Na-MRI. Improving on current estimates supports: (1) larger, more representative data collection/sharing, including (2) regional data, and (3) agreement on full reporting standards.

Test-retest reproducibility of human brain multi-slice 1 H FID-MRSI data at 9.4T after optimization of lipid regularization, macromolecular model, and spline baseline stiffness

PURPOSE

This study analyzes the effects of retrospective lipid suppression, a simulated macromolecular prior knowledge and different spline baseline stiffness values on 9.4T multi-slice proton FID-MRSI data spanning the whole cerebrum of human brain and the reproducibility of respective metabolite ratio to total creatine (/tCr) maps for 10 brain metabolites.

METHODS

Measurements were performed twice on 5 volunteers using a short TR and TE FID MRSI 2D sequence at 9.4T. The effects of retrospective lipid L2-regularization, macromolecular spectrum and different LCModel baseline flexibilities on SNR, FWHM, fitting residual, Cramér-Rao lower bound, and metabolite ratio maps were investigated. Intra-subject, inter-session coefficient of variation and the test-retest reproducibility of the mean metabolite ratios (/tCr) of each slice was calculated.

RESULTS

Transversal, sagittal, and coronal slices of many metabolite ratio maps correspond to the anatomically expected concentration relations in gray and white matter for the majority of the cerebrum when using a flexible baseline in LCModel fit. Results from the second measurements of the same subjects show that slice positioning and data quality correlate significantly to the first measurement. L2-regularization provided effective suppression of lipid-artifacts, but should be avoided if no artifacts are detected.

CONCLUSION

Reproducible concentration ratio maps (/tCr) for 4 metabolites (total choline, N-acetylaspartate, glutamate, and myoinositol) spanning the majority of the cerebrum and 6 metabolites (N-acetylaspartylglutamate, γ-aminobutyric acid, glutathione, taurine, glutamine, and aspartate) covering 32 mm in the upper part of the brain were acquired at 9.4T using multi-slice FID MRSI with retrospective lipid suppression, a macromolecular spectrum and a flexible LCModel baseline.

Neuroimaging at 7 Tesla: a pictorial narrative review

Neuroimaging using the 7-Tesla (7T) human magnetic resonance (MR) system is rapidly gaining popularity after being approved for clinical use in the European Union and the USA. This trend is the same for functional MR imaging (MRI). The primary advantages of 7T over lower magnetic fields are its higher signal-to-noise and contrast-to-noise ratios, which provide high-resolution acquisitions and better contrast, making it easier to detect lesions and structural changes in brain disorders. Another advantage is the capability to measure a greater number of neurochemicals by virtue of the increased spectral resolution. Many structural and functional studies using 7T have been conducted to visualize details in the white matter and layers of the cortex and hippocampus, the subnucleus or regions of the putamen, the globus pallidus, thalamus and substantia nigra, and in small structures, such as the subthalamic nucleus, habenula, perforating arteries, and the perivascular space, that are difficult to observe at lower magnetic field strengths. The target disorders for 7T neuroimaging range from tumoral diseases to vascular, neurodegenerative, and psychiatric disorders, including Alzheimer's disease, Parkinson's disease, multiple sclerosis, epilepsy, major depressive disorder, and schizophrenia. MR spectroscopy has also been used for research because of its increased chemical shift that separates overlapping peaks and resolves neurochemicals more effectively at 7T than a lower magnetic field. This paper presents a narrative review of these topics and an illustrative presentation of images obtained at 7T. We expect 7T neuroimaging to provide a new imaging biomarker of various brain disorders.

In vivo Human MR Spectroscopy Using a Clinical Scanner: Development, Applications, and Future Prospects

MR spectroscopy (MRS) is a unique and useful method for noninvasively evaluating biochemical metabolism in human organs and tissues, but its clinical dissemination has been slow and often limited to specialized institutions or hospitals with experts in MRS technology. The number of 3-T clinical MR scanners is now increasing, representing a major opportunity to promote the use of clinical MRS. In this review, we summarize the theoretical background and basic knowledge required to understand the results obtained with MRS and introduce the general consensus on the clinical utility of proton MRS in routine clinical practice. In addition, we present updates to the consensus guidelines on proton MRS published by the members of a working committee of the Japan Society of Magnetic Resonance in Medicine in 2013. Recent research into multinuclear MRS equipped in clinical MR scanners is explained with an eye toward future development. This article seeks to provide an overview of the current status of clinical MRS and to promote the understanding of when it can be useful. In the coming years, MRS-mediated biochemical evaluation is expected to become available for even routine clinical practice.

Comparison of linear combination modeling strategies for edited magnetic resonance spectroscopy at 3 T

J-difference-edited spectroscopy is a valuable approach for the in vivo detection of γ-aminobutyric-acid (GABA) with magnetic resonance spectroscopy (MRS). A recent expert consensus article recommends linear combination modeling (LCM) of edited MRS but does not give specific details regarding implementation. This study explores different modeling strategies to adapt LCM for GABA-edited MRS. Sixty-one medial parietal lobe GABA-edited MEGA-PRESS spectra from a recent 3-T multisite study were modeled using 102 different strategies combining six different approaches to account for co-edited macromolecules (MMs), three modeling ranges, three baseline knot spacings, and the use of basis sets with or without homocarnosine. The resulting GABA and GABA+ estimates (quantified relative to total creatine), the residuals at different ranges, standard deviations and coefficients of variation (CVs), and Akaike information criteria, were used to evaluate the models' performance. Significantly different GABA+ and GABA estimates were found when a well-parameterized MM_3co basis function was included in the model. The mean GABA estimates were significantly lower when modeling MM_3co , while the CVs were similar. A sparser spline knot spacing led to lower variation in the GABA and GABA+ estimates, and a narrower modeling range-only including the signals of interest-did not substantially improve or degrade modeling performance. Additionally, the results suggest that LCM can separate GABA and the underlying co-edited MM_3co . Incorporating homocarnosine into the modeling did not significantly improve variance in GABA+ estimates. In conclusion, GABA-edited MRS is most appropriately quantified by LCM with a well-parameterized co-edited MM_3co basis function with a constraint to the nonoverlapped MM_0.93 , in combination with a sparse spline knot spacing (0.55 ppm) and a modeling range of 0.5-4 ppm.

Inter-subject stability and regional concentration estimates of 3D-FID-MRSI in the human brain at 7 T

PURPOSE

Recently, a 3D-concentric ring trajectory (CRT)-based free induction decay (FID)-MRSI sequence was introduced for fast high-resolution metabolic imaging at 7 T. This technique provides metabolic ratio maps of almost the entire brain within clinically feasible scan times, but its robustness has not yet been thoroughly investigated. Therefore, we have assessed quantitative concentration estimates and their variability in healthy volunteers using this approach.

METHODS

We acquired whole-brain 3D-CRT-FID-MRSI at 7 T in 15 min with 3.4 mm nominal isometric resolution in 24 volunteers (12 male, 12 female, mean age 27 ± 6 years). Concentration estimate maps were calculated for 15 metabolites using internal water referencing and evaluated in 55 different regions of interest (ROIs) in the brain. Data quality, mean metabolite concentrations, and their inter-subject coefficients of variation (CVs) were compared for all ROIs.

RESULTS

Of 24 datasets, one was excluded due to motion artifacts. The concentrations of total choline, total creatine, glutamate, myo-inositol, and N-acetylaspartate in 44 regions were estimated within quality thresholds. Inter-subject CVs (mean over 44 ROIs/minimum/maximum) were 9%/5%/19% for total choline, 10%/6%/20% for total creatine, 11%/7%/24% for glutamate, 10%/6%/19% for myo-inositol, and 9%/6%/19% for N-acetylaspartate.

DISCUSSION

We defined the performance of 3D-CRT-based FID-MRSI for metabolite concentration estimate mapping, showing which metabolites could be robustly quantified in which ROIs with which inter-subject CVs expected. However, the basal brain regions and lesser-signal metabolites in particular remain as a challenge due susceptibility effects from the proximity to nasal and auditory cavities. Further improvement in quantification and the mitigation of B₀ /B₁ -field inhomogeneities will be necessary to achieve reliable whole-brain coverage.

ProFit-1D-A 1D fitting software and open-source validation data sets

PURPOSE

Accurate and precise MRS fitting is crucial for metabolite concentration quantification of ¹ H-MRS spectra. LCModel, a spectral fitting software, has shown to have certain limitations to perform advanced spectral fitting by previous literature. Herein, we propose an open-source spectral fitting algorithm with adaptive spectral baseline determination and more complex cost functions.

THEORY

The MRS spectra are characterized by several parameters, which reflect the environment of the contributing metabolites, properties of the acquisition sequence, or additional disturbances. Fitting parameters should accurately describe these parameters. Baselines are also a major contributor to MRS spectra, in which smoothness of the spline baselines used for fitting can be adjusted based on the properties of the spectra. Three different cost functions used for the minimization problem were also investigated.

METHODS

The newly developed ProFit-1D fitting algorithm is systematically evaluated for simulations of several types of possible in vivo parameter variations. Although accuracy and precision are tested with simulated spectra, spectra measured in vivo at 9.4 T are used for testing precision using subsets of averages. ProFit-1D fitting results are also compared with LCModel.

RESULTS

Both ProFit-1D and LCModel fitted the spectra well with induced parameter and baseline variations. ProFit-1D proved to be more accurate than LCModel for simulated spectra. However, LCModel showed a somewhat increased precision for some spectral simulations and for in vivo data.

CONCLUSION

The open-source ProFit-1D fitting algorithm demonstrated high accuracy while maintaining precise metabolite concentration quantification. Finally, through the newly proposed cost functions, new ways to improve fitting were shown.

In vivo macromolecule signals in rat brain 1 H-MR spectra at 9.4T: Parametrization, spline baseline estimation, and T2 relaxation times

PURPOSE

Reliable detection and fitting of macromolecules (MM) are crucial for accurate quantification of brain short-echo time (TE) ¹ H-MR spectra. An experimentally acquired single MM spectrum is commonly used. Higher spectral resolution at ultra-high field (UHF) led to increased interest in using a parametrized MM spectrum together with flexible spline baselines to address unpredicted spectroscopic components. Herein, we aimed to: (1) implement an advanced methodological approach for post-processing, fitting, and parametrization of 9.4T rat brain MM spectra; (2) assess the concomitant impact of the LCModel baseline and MM model (ie, single vs parametrized); and (3) estimate the apparent T₂ relaxation times for seven MM components.

METHODS

A single inversion recovery sequence combined with advanced AMARES prior knowledge was used to eliminate the metabolite residuals, fit, and parametrize 10 MM components directly from 9.4T rat brain in vivo ¹ H-MR spectra at different TEs. Monte Carlo simulations were also used to assess the concomitant influence of parametrized MM and DKNTMN parameter in LCModel.

RESULTS

A very stiff baseline (DKNTMN ≥ 1 ppm) in combination with a single MM spectrum led to deviations in metabolite concentrations. For some metabolites the parametrized MM showed deviations from the ground truth for all DKNTMN values. Adding prior knowledge on parametrized MM improved MM and metabolite quantification. The apparent T₂ ranged between 12 and 24 ms for seven MM peaks.

CONCLUSION

Moderate flexibility in the spline baseline was required for reliable quantification of real/experimental spectra based on in vivo and Monte Carlo data. Prior knowledge on parametrized MM improved MM and metabolite quantification.

Contribution of macromolecules to brain 1 H MR spectra: Experts' consensus recommendations

Proton MR spectra of the brain, especially those measured at short and intermediate echo times, contain signals from mobile macromolecules (MM). A description of the main MM is provided in this consensus paper. These broad peaks of MM underlie the narrower peaks of metabolites and often complicate their quantification but they also may have potential importance as biomarkers in specific diseases. Thus, separation of broad MM signals from low molecular weight metabolites enables accurate determination of metabolite concentrations and is of primary interest in many studies. Other studies attempt to understand the origin of the MM spectrum, to decompose it into individual spectral regions or peaks and to use the components of the MM spectrum as markers of various physiological or pathological conditions in biomedical research or clinical practice. The aim of this consensus paper is to provide an overview and some recommendations on how to handle the MM signals in different types of studies together with a list of open issues in the field, which are all summarized at the end of the paper.

Accelerated MR spectroscopic imaging-a review of current and emerging techniques

Over more than 30 years in vivo MR spectroscopic imaging (MRSI) has undergone an enormous evolution from theoretical concepts in the early 1980s to the robust imaging technique that it is today. The development of both fast and efficient sampling and reconstruction techniques has played a fundamental role in this process. State-of-the-art MRSI has grown from a slow purely phase-encoded acquisition technique to a method that today combines the benefits of different acceleration techniques. These include shortening of repetition times, spatial-spectral encoding, undersampling of k-space and time domain, and use of spatial-spectral prior knowledge in the reconstruction. In this way in vivo MRSI has considerably advanced in terms of spatial coverage, spatial resolution, acquisition speed, artifact suppression, number of detectable metabolites and quantification precision. Acceleration not only has been the enabling factor in high-resolution whole-brain 1 H-MRSI, but today is also common in non-proton MRSI (31 P, 2 H and 13 C) and applied in many different organs. In this process, MRSI techniques had to constantly adapt, but have also benefitted from the significant increase of magnetic field strength boosting the signal-to-noise ratio along with high gradient fidelity and high-density receive arrays. In combination with recent trends in image reconstruction and much improved computation power, these advances led to a number of novel developments with respect to MRSI acceleration. Today MRSI allows for non-invasive and non-ionizing mapping of the spatial distribution of various metabolites' tissue concentrations in animals or humans, is applied for clinical diagnostics and has been established as an important tool for neuro-scientific and metabolism research. This review highlights the developments of the last five years and puts them into the context of earlier MRSI acceleration techniques. In addition to 1 H-MRSI it also includes other relevant nuclei and is not limited to certain body regions or specific applications.

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.

Longitudinal quantification of metabolites and macromolecules reveals age- and sex-related changes in the healthy Fischer 344 rat brain

Normal aging is associated with numerous biological changes, including altered brain metabolism and tissue chemistry. In vivo characterization of the neurochemical profile during aging is possible using magnetic resonance spectroscopy, a powerful noninvasive technique capable of quantifying brain metabolites involved in physiological processes that become impaired with age. A prominent macromolecular signal underlies those of brain metabolites and is particularly visible at high fields; parameterization of this signal into components improves quantification and expands the number of biomarkers comprising the neurochemical profile. The present study reports, for the first time, the simultaneous absolute quantification of brain metabolites and individual macromolecules in aging male and female Fischer 344 rats, measured longitudinally using proton magnetic resonance spectroscopy at 7 T. We identified age- and sex-related changes in neurochemistry, with prominent differences in metabolites implicated in anaerobic energy metabolism, antioxidant defenses, and neuroprotection, as well as numerous macromolecule changes. These findings contribute to our understanding of the neurobiological processes associated with healthy aging, critical for the proper identification and management of pathologic aging trajectories. This article is part of the Virtual Special Issue titled COGNITIVE NEUROSCIENCE OF HEALTHY AND PATHOLOGICAL AGING. The full issue can be found on ScienceDirect athttps://www.sciencedirect.com/journal/neurobiology-of-aging/special-issue/105379XPWJP.

Water and lipid suppression techniques for advanced 1 H MRS and MRSI of the human brain: Experts' consensus recommendations

The neurochemical information provided by proton magnetic resonance spectroscopy (MRS) or MR spectroscopic imaging (MRSI) can be severely compromised if strong signals originating from brain water and extracranial lipids are not properly suppressed. The authors of this paper present an overview of advanced water/lipid-suppression techniques and describe their advantages and disadvantages. Moreover, they provide recommendations for choosing the most appropriate techniques for proper use. Methods of water signal handling are primarily focused on the VAPOR technique and on MRS without water suppression (metabolite cycling). The section on lipid-suppression methods in MRSI is divided into three parts. First, lipid-suppression techniques that can be implemented on most clinical MR scanners (volume preselection, outer-volume suppression, selective lipid suppression) are described. Second, lipid-suppression techniques utilizing the combination of k-space filtering, high spatial resolutions and lipid regularization are presented. Finally, three promising new lipid-suppression techniques, which require special hardware (a multi-channel transmit system for dynamic B₁⁺ shimming, a dedicated second-order gradient system or an outer volume crusher coil) are introduced.

Influence of fitting approaches in LCModel on MRS quantification focusing on age-specific macromolecules and the spline baseline

Quantification of neurochemical concentrations from ¹ H MR spectra is challenged by incomplete knowledge of contributing signals. Some experimental conditions hinder the acquisition of artifact-free spectra and impede the acquisition of condition-specific macromolecule (MM) spectra. This work studies differences caused by fitting solutions routinely employed to manage resonances from MM and lipids. High quality spectra (free of residual water and lipid artifacts and for which condition-specific MM spectra are available) are used to understand the influences of spline baseline flexibility and noncondition-specific MM on neurochemical quantification. Fitting with moderate spline flexibility or using noncondition-specific MM led to quantification that differed from when an appropriate, fully specified model was used. This occurred for all neurochemicals to an extent that varied in magnitude among and within approaches. The spline baseline was more tortuous when less constrained and when used in combination with noncondition-specific MM. Increasing baseline flexibility did not reproduce concentrations quantified under appropriate conditions when spectra were fitted using a MM spectrum measured from a mismatched cohort. Using the noncondition-specific MM spectrum led to quantification differences that were comparable in size with using a fitting model that had moderate freedom, and these influences were additive. Although goodness of fit was better with greater fitting flexibility, quantification differed from when fitting with a fully specified model that is appropriate for low noise data. Notable GABA and PE concentration differences occurred with lower estimates of measurement error when fitting with greater spline flexibility or noncondition-specific MM. These data support the need for improved metrics of goodness of fit. Attempting to correct for artifacts or absence of a condition-specific MM spectrum via increased spline flexibility and usage of noncondition-specific MM spectra cannot replace artifact-free data quantified with a condition-specific MM spectrum.

Comparison of different linear-combination modeling algorithms for short-TE proton spectra

Short-TE proton MRS is used to study metabolism in the human brain. Common analysis methods model the data as a linear combination of metabolite basis spectra. This large-scale multi-site study compares the levels of the four major metabolite complexes in short-TE spectra estimated by three linear-combination modeling (LCM) algorithms. 277 medial parietal lobe short-TE PRESS spectra (TE = 35 ms) from a recent 3 T multi-site study were preprocessed with the Osprey software. The resulting spectra were modeled with Osprey, Tarquin and LCModel, using the same three vendor-specific basis sets (GE, Philips and Siemens) for each algorithm. Levels of total N-acetylaspartate (tNAA), total choline (tCho), myo-inositol (mI) and glutamate + glutamine (Glx) were quantified with respect to total creatine (tCr). Group means and coefficient of variations of metabolite estimates agreed well for tNAA and tCho across vendors and algorithms, but substantially less so for Glx and mI, with mI systematically estimated as lower by Tarquin. The cohort mean coefficient of determination for all pairs of LCM algorithms across all datasets and metabolites was R 2 ¯ = 0.39, indicating generally only moderate agreement of individual metabolite estimates between algorithms. There was a significant correlation between local baseline amplitude and metabolite estimates (cohort mean R 2 ¯ = 0.10). While mean estimates of major metabolite complexes broadly agree between linear-combination modeling algorithms at group level, correlations between algorithms are only weak-to-moderate, despite standardized preprocessing, a large sample of young, healthy and cooperative subjects, and high spectral quality. These findings raise concerns about the comparability of MRS studies, which typically use one LCM software and much smaller sample sizes.

Macromolecule suppressed GABA levels show no relationship with age in a pediatric sample

The inhibitory neurotransmitter γ-Aminobutyric acid (GABA) plays a crucial role in cortical development. Therefore, characterizing changes in GABA levels during development has important implications for the study of healthy development and developmental disorders. Brain GABA levels can be measured non-invasively using GABA-edited magnetic resonance spectroscopy (MRS). However, the most commonly used editing technique to measure GABA results in contamination of the GABA signal with macromolecules (MM). Therefore, GABA measured using this technique is often referred to as GABA+ . While few in number, previous studies have shown GABA+ levels increase with age during development. However, these studies are unable to specify whether it is specifically GABA that is increasing or, instead, if levels of MM increase. In this study, we use a GABA-editing technique specifically designed to suppress the MM signal (MM-supp GABA). We find no relationship between MM-supp GABA and age in healthy children aged 7-14 years. These findings suggest that the relationship between GABA+ and age is driven by changes in MM levels, not by changes in GABA levels. Moreover, these findings highlight the importance of accounting for MM levels in MRS quantification.

Realistic head-shaped phantom with brain-mimicking metabolites for 7 T spectroscopy and spectroscopic imaging

PURPOSE

Moving to ultra-high fields (≥7 T), the inhomogeneity of both RF (B1 ) and static (B0 ) magnetic fields increases, which further motivates us to design a realistic head-shaped phantom, especially for spectroscopic imaging. Such phantoms provide images similar to the human brain and serve as a reliable tool for developing and examining methods in MRI. This study aims to develop and characterize a realistic head-shaped phantom filled with brain-mimicking metabolites for MRS and magnetic resonance spectroscopic imaging in a 7 T MRI scanner.

METHODS

A 3D head-shaped container with three sections-mimicking brain, muscle and precranial lipid-was constructed. The phantom was designed to provide robustness to heating, mechanical damage and leakage, with easy refilling. The head's shape and the agarose mixture were optimized to provide B0 and B1 distributions and T1 /T2 relaxation values similar to those of human brain. Eight brain-tissue-mimicking metabolites were included for spectroscopy. The phantom was evaluated for localized spectroscopy, fast spectroscopic imaging and fat suppression.

RESULTS

The B0 and B1 maps showed distribution similar to that of human brain, with increased B0 inhomogeneity near the nasal and ear areas and reduced B1 in the temporal lobe and brain stem regions, as expected in vivo. The metabolites' concentrations were verified by single-voxel spectroscopy, showing an average deviation of 11%. Fast spectroscopic imaging and imaging with fat suppression were demonstrated.

CONCLUSION

A 3D head-shaped phantom for human brain imaging and spectroscopic imaging in 7 T MRI was demonstrated, making it a realistic phantom for methodology development at 7 T.

High-resolution metabolic imaging of high-grade gliomas using 7T-CRT-FID-MRSI

OBJECTIVES

Successful neurosurgical intervention in gliomas depends on the precision of the preoperative definition of the tumor and its margins since a safe maximum resection translates into a better patient outcome. Metabolic high-resolution imaging might result in improved presurgical tumor characterization, and thus optimized glioma resection. To this end, we validated the performance of a fast high-resolution whole-brain 3D-magnetic resonance spectroscopic imaging (MRSI) method at 7T in a patient cohort of 23 high-grade gliomas (HGG).

MATERIALS AND METHODS

We preoperatively measured 23 patients with histologically verified HGGs (17 male, 8 female, age 53 ± 15) with an MRSI sequence based on concentric ring trajectories with a 64 × 64 × 39 measurement matrix, and a 3.4 × 3.4 × 3.4 mm3 nominal voxel volume in 15 min. Quantification used a basis-set of 17 components including N-acetyl-aspartate (NAA), total choline (tCho), total creatine (tCr), glutamate (Glu), glutamine (Gln), glycine (Gly) and 2-hydroxyglutarate (2HG). The resultant metabolic images were evaluated for their reliability as well as their quality and compared to spatially segmented tumor regions-of-interest (necrosis, contrast-enhanced, non-contrast enhanced + edema, peritumoral) based on clinical data and also compared to histopathology (e.g., grade, IDH-status).

RESULTS

Eighteen of the patient measurements were considered usable. In these patients, ten metabolites were quantified with acceptable quality. Gln, Gly, and tCho were increased and NAA and tCr decreased in nearly all tumor regions, with other metabolites such as serine, showing mixed trends. Overall, there was a reliable characterization of metabolic tumor areas. We also found heterogeneity in the metabolic images often continued into the peritumoral region. While 2HG could not be satisfyingly quantified, we found an increase of Glu in the contrast-enhancing region of IDH-wildtype HGGs and a decrease of Glu in IDH1-mutant HGGs.

CONCLUSIONS

We successfully demonstrated high-resolution 7T 3D-MRSI in HGG patients, showing metabolic differences between tumor regions and peritumoral tissue for multiple metabolites. Increases of tCho, Gln (related to tumor metabolism), Gly (related to tumor proliferation), as well as decreases in NAA, tCr, and others, corresponded very well to clinical tumor segmentation, but were more heterogeneous and often extended into the peritumoral region.

A novel method to measure T1 -relaxation times of macromolecules and quantification of the macromolecular resonances

PURPOSE

Macromolecular peaks underlying metabolite spectra influence the quantification of metabolites. Therefore, it is important to understand the extent of contribution from macromolecules (MMs) in metabolite quantification. However, to model MMs more accurately in spectral fitting, differences in T₁ relaxation times among individual MM peaks must be considered. Characterization of T₁ -relaxation times for all individual MM peaks using a single inversion recovery technique is difficult due to eventual contributions from metabolites. On the contrary, a double inversion recovery (DIR) technique provided flexibility to acquire MM spectra spanning a range of longitudinal magnetizations with minimal metabolite influence. Thus, a novel method to determine T₁ -relaxation times of individual MM peaks is reported in this work.

METHODS

Extensive Bloch simulations were performed to determine inversion time combinations for a DIR technique that yielded adequate MM signal with varying longitudinal magnetizations while minimizing metabolite contributions. MM spectra were acquired using DIR-metabolite-cycled semi-LASER sequence. LCModel concentrations were fitted to the DIR signal equation to calculate T₁ -relaxation times.

RESULTS

T₁ -relaxation times of MMs range from 204 to 510 ms and 253 to 564 ms in gray- and white-matter rich voxels respectively at 9.4T. Additionally, concentrations of 13 MM peaks are reported.

CONCLUSION

A novel DIR method is reported in this work to calculate T₁ -relaxation times of MMs in the human brain. T₁ -relaxation times and relaxation time corrected concentrations of individual MMs are reported in gray- and white-matter rich voxels for the first time at 9.4T.

Effects of different macromolecular models on reproducibility of FID-MRSI at 7T

PURPOSE

A properly characterized macromolecular (MM) contribution is essential for accurate metabolite quantification in FID-MRSI. MM information can be included into the fitting model as a single component or parameterized and included over several individual MM resonances, which adds flexibility when pathologic changes are present but is prone to potential overfitting. This study investigates the effects of different MM models on MRSI reproducibility.

METHODS

Clinically feasible, high-resolution FID-MRSI data were collected in ~5 min at 7 Tesla from 10 healthy volunteers and quantified via LCModel (version 6.3) with 3 basis sets, each with a different approach for how the MM signal was handled: averaged measured whole spectrum (full MM), 9 parameterized components (param MM) with soft constraints to avoid overparameterization, or without any MM information included in the fitting prior knowledge. The test-retest reproducibility of MRSI scans was assessed voxel-wise using metabolite coefficients of variation and intraclass correlation coefficients and compared between the basis sets. Correlations of concentration estimates were investigated for the param MM fitting model.

RESULTS

The full MM model provided the most reproducible quantification of total NAA, total Cho, myo-inositol, and glutamate + glutamine ratios to total Cr (coefficients of variations ≤ 8%, intraclass correlation coefficients ≥ 0.76). Using the param MM model resulted in slightly lower reproducibility (up to +3% higher coefficients of variations, up to -0.1 decreased intraclass correlation coefficients). The quantification of the parameterized macromolecules did not affect quantification of the overlapping metabolites.

CONCLUSION

Clinically feasible FID-MRSI with an experimentally acquired MM spectrum included in prior knowledge provides highly reproducible quantification for the most common neurometabolites in healthy volunteers. Parameterization of the MM spectrum may be preferred as a compromise between quantification accuracy and reproducibility when the MM content is expected to be pathologically altered.

Overdiscrete echo-planar spectroscopic imaging with correlated higher-order phase correction

PURPOSE

To introduce a robust methodology for fast ¹ H MRSI of the brain at 3T with improved SNR and reduced phase-related artifacts.

METHOD

An accelerated acquisition scheme using echo-planar spectroscopic imaging (EPSI) was combined with the overdiscrete reconstruction framework. This approach enables the interleaved acquisition of a water reference scan at each phase encoding step, maximizing its correlation with the water-suppressed measurement. Moreover, a generalized high-order phase correction was incorporated into the reconstruction pipeline. The spatial-temporal phase correction term was estimated from the reference scan and interpolated to high resolution using a polynomial basis. The method was implemented at 3T and validated with phantom and in vivo experiments.

RESULTS

The methodology showed the elimination of spectral artifacts generated by phase disturbances and achieved mean SNR gains in vivo of 3.18 and 1.19 compared to standard reconstructions with corrections performed at nominal and high resolution, respectively. EPSI scans with interleaved water acquisition showed to be robust to system instabilities and potentially to patient motion. Moreover, phase distortions were effectively corrected in a single step, avoiding additional reference measurements and post-processing steps.

CONCLUSION

The overdiscrete reconstruction framework with high-order phase correction allowed to effectively correct for distortions, related to B₀ inhomogeneities, B₀ drift, eddy currents, and system vibrations. Furthermore, the presented reconstruction method, combined with EPSI acquisitions, demonstrated improved measurement stability, substantial SNR enhancement, better spectral linewidth, and effective artifact removal.

Effect of linewidth on estimation of metabolic concentration when using water lineshape spectral model fitting for single voxel proton spectroscopy at 7 T

Good B₀ field homogeneity is considered an essential requirement to obtain high-quality MRS data. Many commonly used spectral fitting methods assume that all metabolite signals have Lorentzian or Gaussian shapes. However, B₀ inhomogeneity can both broaden the linewidth and modify the lineshape. In this study, it is hypothesized that a realistic metabolite fitting model, which accounts for B₀ homogeneity on the basis of the water lineshape, will improve the accuracy of estimation of metabolite concentrations. In-vivo water suppressed/unsuppressed single voxel spectroscopy signals were acquired under three different B₀ field homogeneity regimes. Individual realistic basis sets were created for each acquisition. Frequency-domain spectral fitting with LCModel was used to quantify the metabolite concentrations with fitting uncertainties given in terms of the Cramer-Rao lower bound. The quantification results obtained using the water lineshape basis set yielded similar concentrations independent of linewidth and showed a larger fitting error as the linewidth increased. The conventional approach, however quantifies metabolite concentrations with greater variations despite showing a supposedly improved fitting quality. The water lineshape basis set achieved single voxel spectroscopy accuracy that is less sensitive to the linewidth compared to the conventional spectral fitting method for the range of linewidths used in this study, but the precision deteriorated with worsening B₀ field inhomogeneity. The beneficial effect was ascribed to a reduction in the number of degrees of freedom when using the water lineshape to generate the basis set.

High-resolution metabolic mapping of gliomas via patch-based super-resolution magnetic resonance spectroscopic imaging at 7T

OBJECTIVES

To demonstrate the feasibility of 7 T magnetic resonance spectroscopic imaging (MRSI), combined with patch-based super-resolution (PBSR) reconstruction, for high-resolution multi-metabolite mapping of gliomas.

MATERIALS AND METHODS

Ten patients with WHO grade II, III and IV gliomas (6/4, male/female; 45 ± 9 years old) were prospectively measured between 2014 and 2018 on a 7 T whole-body MR imager after routine 3 T magnetic resonance imaging (MRI) and positron emission tomography (PET). Free induction decay MRSI with a 64 × 64-matrix and a nominal voxel size of 3.4 × 3.4 × 8 mm³ was acquired in six minutes, along with standard T1/T2-weighted MRI. Metabolic maps were obtained via spectral LCmodel processing and reconstructed to 0.9 × 0.9 × 8 mm³ resolutions via PBSR.

RESULTS

Metabolite maps obtained from combined 7 T MRSI and PBSR resolved the density of metabolic activity in the gliomas in unprecedented detail. Particularly in the more heterogeneous cases (e.g. post resection), metabolite maps enabled the identification of complex metabolic activities, which were in topographic agreement with PET enhancement.

CONCLUSIONS

PBSR-MRSI combines the benefits of ultra-high-field MR systems, cutting-edge MRSI, and advanced postprocessing to allow millimetric resolution molecular imaging of glioma tissue beyond standard methods. An ideal example is the accurate imaging of glutamine, which is a prime target of modern therapeutic approaches, made possible due to the higher spectral resolution of 7 T systems.

High-resolution echo-planar spectroscopic imaging at ultra-high field

MR spectroscopic imaging (MRSI) at ultra-high field (≥7 T) benefits from improved sensitivity that allows the detection of low-concentration metabolites in the brain. However, optimized acquisition techniques are required to overcome inherent limitations of MRSI at ultra-high field. This work describes an optimized method for fast high-resolution ¹ H-MRSI of the brain at 7 T. The proposed acquisition sequence combines precise volume localization using semi-localization by adiabatic selective refocusing, fast spatial encoding using high-bandwidth symmetric echo-planar spectroscopic imaging (EPSI), and robust water suppression with variable power and optimized relaxation delays. This showed improved robustness to B₀ and B₁⁺ inhomogeneities, eddy currents, nuisance signal contamination and system instabilities. Furthermore, a method for correction of phase inconsistencies in symmetric EPSI enabled high-bandwidth measurements at 7 T. The proposed correction effectively removed spectral ghosting using a single-shot water reference scan. This framework was tested in healthy volunteers at 7 T and spectral quality was compared with lower-spatial-resolution scans, measured at 3 T using the same methodology. A gain in the signal-to-noise ratio (SNR) per unit volume and unit time of 1.57 was achieved, keeping acquisition time short (5 min) and the specific absorption rate within the permitted limits. This SNR enhancement obtained at ultra-high field enabled high-resolution (0.25-0.375 mL) metabolite mapping of the brain within a clinically feasible scan time. The correlation of the reconstructed maps with anatomical structures was observed, showing the diagnostic potential of the technique.

Macromolecule-suppressed GABA measurements correlate more strongly with behavior than macromolecule-contaminated GABA+ measurements

The inhibitory neurotransmitter γ-aminobutyric acid (GABA) is known to be fundamental to the neuronal processes underlying visual orientation and vibrotactile frequency and amplitude discrimination. Previous studies have demonstrated that performance on visual and vibrotactile psychophysics tasks is associated with in vivo measurements of "GABA+" levels - a measure of GABA substantially contaminated by a macromolecular (MM) signal. Here, we establish that these prior findings are indeed driven by the GABA fraction of that signal. Edited magnetic resonance spectroscopy (MRS) was used to measure GABA with and without MM suppression in the sensorimotor (SM1) and occipital cortices in 14 healthy male adults. Volunteers also underwent psychophysical experiments to assess their performance on visual orientation discrimination and vibrotactile amplitude and frequency discrimination. We show that MM-suppressed GABA levels correlate more strongly with individual differences in vibrotactile (in the case of SM1 GABA; amplitude: r = -0.63, p = 0.03; frequency: r = -0.62, p = 0.02) and visual orientation (in the case of occipital GABA; r = -0.59, p = 0.05) discrimination thresholds than GABA levels contaminated by MM (vibrotactile amplitude: r = -0.36, p = 0.30; vibrotactile frequency: r = -0.53, p = 0.09; visual orientation: r = 0.21, p = 0.55). These findings further support the view that measurements of endogenous GABA acquired with edited MRS can usefully probe neurochemical-behavioral relationships in humans. Moreover, the more specific measurement of GABA used in this study provides increased statistical power to observe these regionally specific relationships.