MRSCloud: A cloud-based MRS tool for basis set simulation

GABA-edited MEGA-PRESS at 3 T: Does a measured macromolecule background improve linear combination modeling?

PURPOSE

The J-difference edited γ-aminobutyric acid (GABA) signal is contaminated by other co-edited signals-the largest of which originates from co-edited macromolecules (MMs)-and is consequently often reported as "GABA+." MM signals are broader and less well-characterized than the metabolites, and are commonly approximated using a Gaussian model parameterization. Experimentally measured MM signals are a consensus-recommended alternative to parameterized modeling; however, they are relatively under-studied in the context of edited MRS.

METHODS

To address this limitation in the literature, we have acquired GABA-edited MEGA-PRESS data with pre-inversion to null metabolite signals in 13 healthy controls. An experimental MM basis function was derived from the mean across subjects. We further derived a new parameterization of the MM signals from the experimental data, using multiple Gaussians to accurately represent their observed asymmetry. The previous single-Gaussian parameterization, mean experimental MM spectrum and new multi-Gaussian parameterization were compared in a three-way analysis of a public MEGA-PRESS dataset of 61 healthy participants.

RESULTS

Both the experimental MMs and the multi-Gaussian parameterization exhibited reduced fit residuals compared to the single-Gaussian approach (p = 0.034 and p = 0.031, respectively), suggesting they better represent the underlying data than the single-Gaussian parameterization. Furthermore, both experimentally derived models estimated larger MM fractional contribution to the GABA+ signal for the experimental MMs (58%) and multi-Gaussian parameterization (58%), compared to the single-Gaussian approach (50%).

CONCLUSIONS

Our results indicate that single-Gaussian parameterization of edited MM signals is insufficient and that both experimentally derived GABA+ spectra and their parameterized replicas improve the modeling of GABA+ spectra.

From Molecules to Behavior: Linking Brain Biochemistry in the Dorsal Anterior Cingulate Cortex to 'Anxiety to Eat' Responses in Anorexia Nervosa

Background

Anorexia nervosa (AN) is a mental and behavioral health condition characterized by an intense fear of weight or fat gain, severe restriction of food intake resulting in low body weight, and distorted self-perception of body shape or weight. While substantial research has focused on general anxiety in AN, less is known about eating-related anxiety and its underlying neural mechanisms. Therefore, we sought to characterize 'anxiety to eat' in AN and examine the neurometabolic profile within the dorsal anterior cingulate cortex (dACC), a brain region putatively involved in magnifying the threat response.

Methods

Women seeking inpatient treatment for AN and women of healthy weight without a lifetime history of an eating disorder (healthy controls; HC) completed a computer-based behavioral task assessing 'anxiety to eat' in response to images of higher (HED) and lower (LED) energy density foods. Participants also underwent magnetic resonance spectroscopy of the dACC in a 3T scanner.

Results

The AN group reported greater anxiety to eat both HED and LED foods relative to the HC group. Both groups reported greater anxiety to eat HED foods relative to LED foods. The neurometabolite myo-inositol (mI) was lower in the dACC in AN relative to HC, and mI levels negatively predicted anxiety to eat HED but not LED foods in the AN group only. mI levels in the dACC were independent of body weight, body mass, and general anxiety.

Conclusions

These findings provide critical new insight into the clinically challenging feature and underlying neural mechanisms of eating-related anxiety, and indicate mI levels in the dACC could serve as a novel biomarker of illness severity that is independent of body weight to identify individuals vulnerable to disordered eating or eating pathology as well as a potential therapeutic target.

Simultaneous multi-transient linear-combination modeling of MRS data improves uncertainty estimation

PURPOSE

The interest in applying and modeling dynamic MRS has recently grown. Two-dimensional modeling yields advantages for the precision of metabolite estimation in interrelated MRS data. However, it is unknown whether including all transients simultaneously in a 2D model without averaging (presuming a stable signal) performs similarly to one-dimensional (1D) modeling of the averaged spectrum. Therefore, we systematically investigated the accuracy, precision, and uncertainty estimation of both described model approaches.

METHODS

Monte Carlo simulations of synthetic MRS data were used to compare the accuracy and uncertainty estimation of simultaneous 2D multitransient linear-combination modeling (LCM) with 1D-LCM of the average. A total of 2,500 data sets per condition with different noise representations of a 64-transient MRS experiment at six signal-to-noise levels for two separate spin systems (scyllo-inositol and gamma-aminobutyric acid) were analyzed. Additional data sets with different levels of noise correlation were also analyzed. Modeling accuracy was assessed by determining the relative bias of the estimated amplitudes against the ground truth, and modeling precision was determined by SDs and Cramér-Rao lower bounds (CRLBs).

RESULTS

Amplitude estimates for 1D- and 2D-LCM agreed well and showed a similar level of bias compared with the ground truth. Estimated CRLBs agreed well between both models and with ground-truth CRLBs. For correlated noise, the estimated CRLBs increased with the correlation strength for the 1D-LCM but remained stable for the 2D-LCM.

CONCLUSION

Our results indicate that the model performance of 2D multitransient LCM is similar to averaged 1D-LCM. This validation on a simplified scenario serves as a necessary basis for further applications of 2D modeling.

Integrated Short-TE and Hadamard-edited Multi-Sequence (ISTHMUS) for Advanced MRS

Background

To examine data quality and reproducibility using ISTHMUS, which has been implemented as the standardized MR spectroscopy sequence for the multi-site Healthy Brain and Child Development (HBCD) study.

Methods

ISTHMUS is the consecutive acquisition of short-TE PRESS (32 transients) and long-TE HERCULES (224 transients) data with dual-TE water reference scans. Voxels were positioned in the centrum semiovale, dorsal anterior cingulate cortex, posterior cingulate cortex and bilateral thalamus regions. After acquisition, ISTHMUS data were separated into the PRESS and HERCULES portions for analysis and modeled separately using Osprey. In vivo experiments were performed in 10 healthy volunteers (6 female; 29.5±6.6 years). Each volunteer underwent two scans on the same day. Differences in metabolite measurements were examined. T2 correction based on the dual-TE water integrals were compared with: 1) T2 correction based the default white matter and gray matter T2 reference values in Osprey; 2) shorter WM and GM T2 values from recent literature; and 3) reduced CSF fractions.

Results

No significant difference in linewidth was observed between PRESS and HERCULES. Bilateral thalamus spectra had produced significantly higher (p<0.001) linewidth compared to the other three regions. Linewidth measurements were similar between scans, with scan-to-scan differences under 1 Hz for most subjects. Paired t-tests indicated a significant difference only in PRESS NAAG between the two thalamus scans (p=0.002). T2 correction based on shorter T2 values showed better agreement to the dual-TE water integral ratio.

Conclusions

ISTHMUS facilitated and standardized acquisition and post-processing and reduced operator workload to eliminate potential human error.

Investigating the brain's neurochemical profile at midlife in relation to dementia risk factors

Changes in the brain's physiology in Alzheimer's disease are thought to occur early in the disease's trajectory. In this study our aim was to investigate the brain's neurochemical profile in a midlife cohort in relation to risk factors for future dementia using single voxel proton magnetic resonance spectroscopy. Participants in the multi-site PREVENT-Dementia study (age range 40-59 year old) underwent 3T magnetic resonance spectroscopy with the spectroscopy voxel placed in the posterior cingulate/precuneus region. Using LCModel, we quantified the absolute concentrations of myo-inositol, total N-acetylaspartate, total creatine, choline, glutathione and glutamate-glutamine for 406 participants (mean age 51.1; 65.3% female). Underlying partial volume effects were accounted for by applying a correction for the presence of cerebrospinal fluid in the magnetic resonance spectroscopy voxel. We investigated how metabolite concentrations related to apolipoprotein ɛ4 genotype, dementia family history, a risk score (Cardiovascular Risk Factors, Aging and Incidence of Dementia -CAIDE) for future dementia including non-modifiable and potentially-modifiable factors and dietary patterns (adherence to Mediterranean diet). Dementia family history was associated with decreased total N-acetylaspartate and no differences were found between apolipoprotein ɛ4 carriers and non-carriers. A higher Cardiovascular Risk Factors, Aging, and Incidence of Dementia score related to higher myo-inositol, choline, total creatine and glutamate-glutamine, an effect which was mainly driven by older age and a higher body mass index. Greater adherence to the Mediterranean diet was associated with lower choline, myo-inositol and total creatine; these effects did not survive correction for multiple comparisons. The observed associations suggest that at midlife the brain demonstrates subtle neurochemical changes in relation to both inherited and potentially modifiable risk factors for future dementia.

Metabolite T1 relaxation times decrease across the adult lifespan

Relaxation correction is an integral step in quantifying brain metabolite concentrations measured by in vivo magnetic resonance spectroscopy (MRS). While most quantification routines assume constant T1 relaxation across age, it is possible that aging alters T1 relaxation rates, as is seen for T2 relaxation. Here, we investigate the age dependence of metabolite T1 relaxation times at 3 T in both gray- and white-matter-rich voxels using publicly available metabolite and metabolite-nulled (single inversion recovery TI = 600 ms) spectra acquired at 3 T using Point RESolved Spectroscopy (PRESS) localization. Data were acquired from voxels in the posterior cingulate cortex (PCC) and centrum semiovale (CSO) in 102 healthy volunteers across 5 decades of life (aged 20-69 years). All spectra were analyzed in Osprey v.2.4.0. To estimate T1 relaxation times for total N-acetyl aspartate at 2.0 ppm (tNAA2.0) and total creatine at 3.0 ppm (tCr3.0), the ratio of modeled metabolite residual amplitudes in the metabolite-nulled spectrum to the full metabolite signal was calculated using the single-inversion-recovery signal equation. Correlations between T1 and subject age were evaluated. Spearman correlations revealed that estimated T1 relaxation times of tNAA2.0 (rs = -0.27; p < 0.006) and tCr3.0 (rs = -0.40; p < 0.001) decreased significantly with age in white-matter-rich CSO, and less steeply for tNAA2.0 (rs = -0.228; p = 0.005) and (not significantly for) tCr3.0 (rs = -0.13; p = 0.196) in graymatter-rich PCC. The analysis harnessed a large publicly available cross-sectional dataset to test an important hypothesis, that metabolite T1 relaxation times change with age. This preliminary study stresses the importance of further work to measure age-normed metabolite T1 relaxation times for accurate quantification of metabolite levels in studies of aging.

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.

Model-based frequency-and-phase correction of 1H MRS data with 2D linear-combination modeling

Purpose

Retrospective frequency-and-phase correction (FPC) methods attempt to remove frequency-and-phase variations between transients to improve the quality of the averaged MR spectrum. However, traditional FPC methods like spectral registration struggle at low SNR. Here, we propose a method that directly integrates FPC into a two-dimensional linear-combination model (2D-LCM) of individual transients ('model-based FPC'). We investigated how model-based FPC performs compared to the traditional approach, i.e., spectral registration followed by 1D-LCM in estimating frequency-and-phase drifts and, consequentially, metabolite level estimates.

Methods

We created synthetic in-vivo-like 64-transient short-TE sLASER datasets with 100 noise realizations at 5 SNR levels and added randomly sampled frequency and phase variations. We then used this synthetic dataset to compare the performance of 2D-LCM with the traditional approach (spectral registration, averaging, then 1D-LCM). Outcome measures were the frequency/phase/amplitude errors, the standard deviation of those ground-truth errors, and amplitude Cramér Rao Lower Bounds (CRLBs). We further tested the proposed method on publicly available in-vivo short-TE PRESS data.

Results

2D-LCM estimates (and accounts for) frequency-and-phase variations directly from uncorrected data with equivalent or better fidelity than the conventional approach. Furthermore, 2D-LCM metabolite amplitude estimates were at least as accurate, precise, and certain as the conventionally derived estimates. 2D-LCM estimation of frequency and phase correction and amplitudes performed substantially better at low-to-very-low SNR.

Conclusion

Model-based FPC with 2D linear-combination modeling is feasible and has great potential to improve metabolite level estimation for conventional and dynamic MRS data, especially for low-SNR conditions, e.g., long TEs or strong diffusion weighting.

Brain glutathione and GABA+ levels in autistic children

Autism spectrum disorder (ASD) is a neurodevelopmental condition characterized by social communication challenges and repetitive behaviors. Altered neurometabolite levels, including glutathione (GSH) and gamma-aminobutyric acid (GABA), have been proposed as potential contributors to the biology underlying ASD. This study investigated whether cerebral GSH or GABA levels differ between a cohort of children aged 8-12 years with ASD (n = 52) and typically developing children (TDC, n = 49). A comprehensive analysis of GSH and GABA levels in multiple brain regions, including the primary motor cortex (SM1), thalamus (Thal), medial prefrontal cortex (mPFC), and supplementary motor area (SMA), was conducted using single-voxel HERMES MR spectroscopy at 3T. The results revealed no significant differences in cerebral GSH or GABA levels between the ASD and TDC groups across all examined regions. These findings suggest that the concentrations of GSH (an important antioxidant and neuromodulator) and GABA (a major inhibitory neurotransmitter) do not exhibit marked alterations in children with ASD compared to TDC. A statistically significant positive correlation was observed between GABA levels in the SM1 and Thal regions with ADHD inattention scores. No significant correlation was found between metabolite levels and hyper/impulsive scores of ADHD, measures of core ASD symptoms (ADOS-2, SRS-P) or adaptive behavior (ABAS-2). While both GSH and GABA have been implicated in various neurological disorders, the current study provides valuable insights into the specific context of ASD and highlights the need for further research to explore other neurochemical alterations that may contribute to the pathophysiology of this complex disorder.

sLASER and PRESS perform similarly at revealing metabolite-age correlations at 3 T

PURPOSE

To compare the respective ability of PRESS and sLASER to reveal biological relationships, using age as a validation covariate at 3 T.

METHODS

MRS data were acquired from 102 healthy volunteers using PRESS and sLASER in centrum semiovale and posterior cingulate cortex (PCC). Acquisition parameters included TR/TE = 2000/30 ms, 96 transients, and 2048 datapoints sampled at 2 kHz. Spectra were analyzed using Osprey. SNR, FWHM linewidth of total creatine, and metabolite concentrations were extracted. A linear model was used to compare SNR and linewidth. Paired t-tests were used to assess differences in metabolite measurements between PRESS and sLASER. Correlations were used to evaluate the relationship between PRESS and sLASER metabolite estimates, as well as the strength of each metabolite-age relationship. Coefficients of variation were calculated to assess inter-subject variability in each metabolite measurement.

RESULTS

SNR and linewidth were significantly higher (p < 0.01) for sLASER than PRESS in PCC. Paired t-tests showed significant differences between PRESS and sLASER in most metabolite measurements. PRESS-sLASER measurements were significantly correlated (p < 0.05) for most metabolites. Metabolite-age relationships were consistently identified using both methods. Similar coefficients of variation were observed for most metabolites.

CONCLUSION

The study results suggest strong agreement between PRESS and sLASER in identifying relationships between brain metabolites and age in centrum semiovale and PCC data acquired at 3 T. sLASER is technically desirable due to the reduced chemical shift displacement artifact; however, PRESS performed similarly in homogeneous brain regions at clinical field strength.

Simultaneous multi-transient linear-combination modeling of MRS data improves uncertainty estimation

Purpose

The interest in applying and modeling dynamic MRS has recently grown. 2D modeling yields advantages for the precision of metabolite estimation in interrelated MRS data. However, it is unknown whether including all transients simultaneously in a 2D model without averaging (presuming a stable signal) performs similarly to 1D modeling of the averaged spectrum. Therefore, we systematically investigated the accuracy, precision, and uncertainty estimation of both described model approaches.

Methods

Monte Carlo simulations of synthetic MRS data were used to compare the accuracy and uncertainty estimation of simultaneous 2D multi-transient LCM with 1D-LCM of the average. 2,500 datasets per condition with different noise representations of a 64-transient MRS experiment at 6 signal-to-noise levels for two separate spin systems (scyllo-inositol and GABA) were analyzed. Additional datasets with different levels of noise correlation were also analyzed. Modeling accuracy was assessed by determining the relative bias of the estimated amplitudes against the ground truth, and modeling precision was determined by standard deviations and Cramér-Rao Lower Bounds (CRLB).

Results

Amplitude estimates for 1D- and 2D-LCM agreed well and showed similar level of bias compared to the ground truth. Estimated CRLBs agreed well between both models and with ground truth CRLBs. For correlated noise the estimated CRLBs increased with the correlation strength for the 1D-LCM but remained stable for the 2D-LCM.

Conclusion

Our results indicate that the model performance of 2D multi-transient LCM is similar to averaged 1D-LCM. This validation on a simplified scenario serves as necessary basis for further applications of 2D modeling.

CloudBrain-MRS: An intelligent cloud computing platform for in vivo magnetic resonance spectroscopy preprocessing, quantification, and analysis

Magnetic resonance spectroscopy (MRS) is an important clinical imaging method for diagnosis of diseases. MRS spectrum is used to observe the signal intensity of metabolites or further infer their concentrations. Although the magnetic resonance vendors commonly provide basic functions of spectrum plots and metabolite quantification, the spread of clinical research of MRS is still limited due to the lack of easy-to-use processing software or platform. To address this issue, we have developed CloudBrain-MRS, a cloud-based online platform that provides powerful hardware and advanced algorithms. The platform can be accessed simply through a web browser, without the need of any program installation on the user side. CloudBrain-MRS also integrates the classic LCModel and advanced artificial intelligence algorithms and supports batch preprocessing, quantification, and analysis of MRS data from different vendors. Additionally, the platform offers useful functions: (1) Automatically statistical analysis to find biomarkers for diseases; (2) Consistency verification between the classic and artificial intelligence quantification algorithms; (3) Colorful three-dimensional visualization for easy observation of individual metabolite spectrum. Last, data of both healthy subjects and patients with mild cognitive impairment are used to demonstrate the functions of the platform. To the best of our knowledge, this is the first cloud computing platform for in vivo MRS with artificial intelligence processing. We have shared our cloud platform at MRSHub, providing at least two years of free access and service. If you are interested, please visit https://mrshub.org/software_all/#CloudBrain-MRS or https://csrc.xmu.edu.cn/CloudBrain.html.

Neurochemical Ameliorating of the Hippocampus in Dyslipidemic Alzheimer Patients Following Silymarin; a Double-Blind Placebo-Controlled Randomized Clinical Trial

Background

Amyloid-beta (Aβ) production is a normal physiological process, and an imbalance in Aβ production/excretion rate is the basis of the plaque load increase in AD. LRP1 is involved in both central clearance of Aβ from the CNS and transport of Aβ toward peripheral organs. In this study, the effect of silymarin combination compared to rosuvastatin and placebo on neuro-metabolites and serum levels of LRP1 and Aβ1-42 proteins and oxidative stress enzymes and lipid and cognitive tests of Iranian AD patients.

Methods

In this double-blind placebo-controlled study, thirty-six mild AD patients were divided into groups (n=12) of silymarin 140mg, placebo, and rosuvastatin 10mg. Medications were administered 3 times a day for 6 months. Clinical tests, lipid profile (TG, HDL, TC, and LDL), Aβ1-42, and LRP1 markers were measured at the beginning and end of the intervention. Magnetic resonance spectroscopy (MRS) was used to measure metabolites. Using SPSS software a one-way ANOVA test was used to compare the means of the quantitative variables and Pearson and Spearman's correlations to measure the correlation. GraphPad Prism software was used for drawing graphs. P < 0.05 was considered a significant.

Results

The levels of LRP1 and Aβ1-42 in the silymarin group were significantly increased compared to the other groups (P < 0.05). NAA/mI in the silymarin group had a significant increase compared to both placebo and rosuvastatin groups (P < 0.05). Right and left hippocampal mI/Cr directly correlated with TG (r = 0.603, P = 0.003 and r = 0.595, P = 0.004, respectively). NAA/Cr of the right and left hippocampus was inversely related to TG (r = -0.511, P = 0.0033, and r = -0.532, P = 0.0021, respectively). NAA/Cr and NAA/mI of bilateral hippocampi directly correlated with HDL (P < 0.05). An inverse correlation was observed between the Aβ1-42 and mI/Cr of the right and left hippocampus (r = -0.661, P = 0.000 and r = -0.638, P = 0.000, respectively).

Conclusion

Donepezil and silymarin improved lipid profile associated with increased NAA/Cr, and decreased mI/Cr, in AD patients. Biomarker NAA/mI can be clinically significant in examining AD pathology. Measurement of the lipid factors and neurometabolites can be a suitable method for monitoring this disease.

Brain Glutathione and GABA+ levels in autistic children

Autism spectrum disorder (ASD) is a neurodevelopmental condition characterized by social communication challenges and repetitive behaviors. Altered neurometabolite levels, including glutathione (GSH) and gamma-aminobutyric acid (GABA), have been proposed as potential contributors to the biology underlying ASD. This study investigated whether cerebral GSH or GABA levels differ between a large cohort of children aged 8-12 years with ASD (n=52) and typically developing children (TDC, n=49). A comprehensive analysis of GSH and GABA levels in multiple brain regions, including the primary motor cortex (SM1), thalamus (Thal), medial prefrontal cortex (mPFC), and supplementary motor area (SMA), was conducted using single-voxel HERMES MR spectroscopy at 3T. The results revealed no significant differences in cerebral GSH or GABA levels between the ASD and TDC groups across all examined regions. These findings suggest that the concentrations of GSH (an important antioxidant and neuromodulator) and GABA (a major inhibitory neurotransmitter) do not exhibit marked alterations in children with ASD compared to TDC. A statistically significant positive correlation was observed between GABA levels in the SM1 and Thal regions with ADHD inattention scores. No significant correlation was found between metabolite levels and hyper/impulsive scores of ADHD, measures of core ASD symptoms (ADOS-2, SRS-P) or adaptive behavior (ABAS-2). While both GSH and GABA have been implicated in various neurological disorders, the current study provides valuable insights into the specific context of ASD and highlights the need for further research to explore other neurochemical alterations that may contribute to the pathophysiology of this complex disorder.

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.

Continuous Automated Analysis Workflow for MRS Studies

Magnetic resonance spectroscopy (MRS) can non-invasively measure levels of endogenous metabolites in living tissue and is of great interest to neuroscience and clinical research. To this day, MRS data analysis workflows differ substantially between groups, frequently requiring many manual steps to be performed on individual datasets, e.g., data renaming/sorting, manual execution of analysis scripts, and manual assessment of success/failure. Manual analysis practices are a substantial barrier to wider uptake of MRS. They also increase the likelihood of human error and prevent deployment of MRS at large scale. Here, we demonstrate an end-to-end workflow for fully automated data uptake, processing, and quality review.The proposed continuous automated MRS analysis workflow integrates several recent innovations in MRS data and file storage conventions. They are efficiently deployed by a directory monitoring service that automatically triggers the following steps upon arrival of a new raw MRS dataset in a project folder: (1) conversion from proprietary manufacturer file formats into the universal format NIfTI-MRS; (2) consistent file system organization according to the data accumulation logic standard BIDS-MRS; (3) executing a command-line executable of our open-source end-to-end analysis software Osprey; (4) e-mail delivery of a quality control summary report for all analysis steps.The automated architecture successfully completed for a demonstration dataset. The only manual step required was to copy a raw data folder into a monitored directory.Continuous automated analysis of MRS data can reduce the burden of manual data analysis and quality control, particularly for non-expert users and multi-center or large-scale studies and offers considerable economic advantages.

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.

Vespa: Integrated applications for RF pulse design, spectral simulation and MRS data analysis

PURPOSE

The Vespa package (Versatile Simulation, Pulses, and Analysis) is described and demonstrated. It provides workflows for developing and optimizing linear combination modeling (LCM) fitting for ¹ H MRS data using intuitive graphical user interface interfaces for RF pulse design, spectral simulation, and MRS data analysis. Command line interfaces for embedding workflows in MR manufacturer platforms and utilities for synthetic dataset creation are included. Complete provenance is maintained for all steps in workflows.

THEORY AND METHODS

Vespa is written in Python for compatibility across operating systems. It embeds the PyGAMMA spectral simulation library for spectral simulation. Multiprocessing methods accelerate processing and visualization. Applications use the Vespa database for results storage and cross-application access. Three projects demonstrate pulse, sequence, simulation, and data analysis workflows: (1) short TE semi-LASER single-voxel spectroscopy (SVS) LCM fitting, (2) optimizing MEGA-PRESS (MEscher-GArwood Point RESolved Spectroscopy) flip angle and LCM fitting, and (3) creating a synthetic short TE dataset.

RESULTS

The LCM workflows for in vivo basis set creation and spectral analysis showed reasonable results for both the short TE semi-LASER and MEGA-PRESS. Examples of pulses, simulations, and data fitting are shown in Vespa application interfaces for various steps to demonstrate the interactive workflow.

CONCLUSION

Vespa provides an efficient and extensible platform for characterizing RF pulses, pulse design, spectral simulation optimization, and automated LCM fitting via an interactive platform. Modular design and command line interface make it easy to embed in other platforms. As open source, it is free to the MRS community for use and extension. Vespa source code and documentation are available through GitHub.

J-difference GABA-edited MRS reveals altered cerebello-thalamo-cortical metabolism in patients with hepatic encephalopathy

Hepatic encephalopathy (HE) is a common neurological manifestation of liver cirrhosis and is characterized by an increase of ammonia in the brain accompanied by a disrupted neurotransmitter balance, including the GABAergic and glutamatergic systems. The aim of this study is to investigate metabolic abnormalities in the cerebello-thalamo-cortical system of HE patients using GABA-edited MRS and links between metabolite levels, disease severity, critical flicker frequency (CFF), motor performance scores, and blood ammonia levels. GABA-edited MRS was performed in 35 participants (16 controls, 19 HE patients) on a clinical 3 T MRI system. MRS voxels were placed in the right cerebellum, left thalamus, and left motor cortex. Levels of GABA+ and of other metabolites of interest (glutamine, glutamate, myo-inositol, glutathione, total choline, total NAA, and total creatine) were assessed. Group differences in metabolite levels and associations with clinical metrics were tested. GABA+ levels were significantly increased in the cerebellum of patients with HE. GABA+ levels in the motor cortex were significantly decreased in HE patients, and correlated with the CFF (r = 0.73; p < .05) and motor performance scores (r = -0.65; p < .05). Well-established HE-typical metabolite patterns (increased glutamine, decreased myo-inositol and total choline) were confirmed in all three regions and were closely linked to clinical metrics. In summary, our findings provide further evidence for alterations in the GABAergic system in the cerebellum and motor cortex in HE. These changes were accompanied by characteristic patterns of osmolytes and oxidative stress markers in the cerebello-thalamo-cortical system. These metabolic disturbances are a likely contributor to HE motor symptoms in HE. In patients with hepatic encephalopathy, GABA+ levels in the cerebello-thalamo-cortical loop are significantly increased in the cerebellum and significantly decreased in the motor cortex. GABA+ levels in the motor cortex strongly correlate with critical flicker frequency (CFF) and motor performance score (pegboard test tPEG), but not blood ammonia levels (NH3).

Brain total creatine differs between primary progressive aphasia (PPA) subtypes and correlates with disease severity

Primary progressive aphasia (PPA) is comprised of three subtypes: logopenic (lvPPA), non-fluent (nfvPPA), and semantic (svPPA). We used magnetic resonance spectroscopy (MRS) to measure tissue-corrected metabolite levels in the left inferior frontal gyrus (IFG) and right sensorimotor cortex (SMC) from 61 PPA patients. We aimed to: (1) characterize subtype differences in metabolites; and (2) test for metabolite associations with symptom severity. tCr differed by subtype across the left IFG and right SMC. tCr levels were lowest in lvPPA and highest in svPPA. tCr levels predicted lvPPA versus svPPA diagnosis. Higher IFG tCr and lower Glx correlated with greater disease severity. As tCr is involved in brain energy metabolism, svPPA pathology might involve changes in specific cellular energy processes. Perturbations to cellular energy homeostasis in language areas may contribute to symptoms. Reduced cortical excitatory capacity (i.e. lower Glx) in language regions may also contribute to symptoms. Thus, tCr may be useful for differentiating between PPA subtypes, and both tCr and Glx might have utility in understanding PPA mechanisms and tracking progression.

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.

sLASER and PRESS Perform Similarly at Revealing Metabolite-Age Correlations

Purpose

To compare the respective ability of PRESS and sLASER to reveal biological relationships, using age as a validation covariate.

Methods

MRS data were acquired from 102 healthy volunteers using PRESS and sLASER in centrum semiovale (CSO) and posterior cingulate cortex (PCC) regions. Acquisition parameters included TR/TE 2000/30 ms; 96 transients; 2048 datapoints sampled at 2 kHz.Spectra were analyzed using Osprey. Signal-to-noise ratio (SNR), full-width-half-maximum linewidth of tCr, and metabolite concentrations were extracted. A linear model was used to compare SNR and linewidth. Paired t-tests were used to assess differences in metabolite measurements between PRESS and sLASER. Correlations were used to evaluate the relationship between PRESS and sLASER metabolite estimates, as well as the strength of each metabolite-age relationship. Coefficients of variation were calculated to assess inter-subject variability in each metabolite measurement.

Results

SNR and linewidth were significantly higher (p<0.05) for sLASER than PRESS. Paired t-tests showed significant differences between PRESS and sLASER in most metabolite measurements. Metabolite measures were significantly correlated (p<0.05) for most metabolites between the two methods except GABA, Gln and Lac in CSO and GSH, Lac and NAAG in PCC. Metabolite-age relationships were consistently identified using both PRESS and sLASER. Similar CVs were observed for most metabolites.

Conclusion

The study results suggest strong agreement between PRESS and sLASER in identifying relationships between brain metabolites and age in CSO and PCC data acquired at 3T. sLASER is technically desirable due to the reduced chemical shift displacement artifact; however, PRESS performed similarly in 'good' brain regions at clinical field strength.