INSPECTOR: free software for magnetic resonance spectroscopy data inspection, processing, simulation and analysis

Retrospective analysis of Braak stage- and APOE4 allele-dependent associations between MR spectroscopy and markers of tau and neurodegeneration in cognitively unimpaired elderly

PURPOSE

The pathological hallmarks of Alzheimer's disease (AD), amyloid, tau, and associated neurodegeneration, are present in the cortical gray matter (GM) years before symptom onset, and at significantly greater levels in carriers of the apolipoprotein E4 (APOE4) allele. Their respective biomarkers, A/T/N, have been found to correlate with aspects of brain biochemistry, measured with magnetic resonance spectroscopy (MRS), indicating a potential for MRS to augment the A/T/N framework for staging and prediction of AD. Unfortunately, the relationships between MRS and A/T/N biomarkers are unclear, largely due to a lack of studies examining them in the context of the spatial and temporal model of T/N progression. Advanced MRS acquisition and post-processing approaches have enabled us to address this knowledge gap and test the hypotheses, that glutamate-plus-glutamine (Glx) and N-acetyl-aspartate (NAA), metabolites reflecting synaptic and neuronal health, respectively, measured from regions on the Braak stage continuum, correlate with: (i) cerebrospinal fluid (CSF) p-tau181 level (T), and (ii) hippocampal volume or cortical thickness of parietal lobe GM (N). We hypothesized that these correlations will be moderated by Braak stage and APOE4 genotype.

METHODS

We conducted a retrospective imaging study of 34 cognitively unimpaired elderly individuals who received APOE4 genotyping and lumbar puncture from pre-existing prospective studies at the NYU Grossman School of Medicine between October 2014 and January 2019. Subjects returned for their imaging exam between April 2018 and February 2020. Metabolites were measured from the left hippocampus (Braak II) using a single-voxel semi-adiabatic localization by adiabatic selective refocusing sequence; and from the bilateral posterior cingulate cortex (PCC; Braak IV), bilateral precuneus (Braak V), and bilateral precentral gyrus (Braak VI) using a multi-voxel echo-planar spectroscopic imaging sequence. Pearson and Spearman correlations were used to examine the relationships between absolute levels of choline, creatine, myo-inositol, Glx, and NAA and CSF p-tau181, and between these metabolites and hippocampal volume or parietal cortical thicknesses. Covariates included age, sex, years of education, Fazekas score, and months between CSF collection and MRI exam.

RESULTS

There was a direct correlation between hippocampal Glx and CSF p-tau181 in APOE4 carriers (Pearson's r = 0.76, p = 0.02), but not after adjusting for covariates. In the entire cohort, there was a direct correlation between hippocampal NAA and hippocampal volume (Spearman's r = 0.55, p = 0.001), even after adjusting for age and Fazekas score (Spearman's r = 0.48, p = 0.006). This relationship was observed only in APOE4 carriers (Pearson's r = 0.66, p = 0.017), and was also retained after adjustment (Pearson's r = 0.76, p = 0.008; metabolite-by-carrier interaction p = 0.03). There were no findings in the PCC, nor in the negative control (late Braak stage) regions of the precuneus and precentral gyrus.

CONCLUSIONS

Our findings are in line with the spatially- and temporally-resolved Braak staging model of pathological severity in which the hippocampus is affected earlier than the PCC. The correlations, between MRS markers of synaptic and neuronal health and, respectively, T and N pathology, were found exclusively within APOE4 carriers, suggesting a connection with AD pathological change, rather than with normal aging. We therefore conclude that MRS has the potential to augment early A/T/N staging, with the hippocampus serving as a more sensitive MRS target compared to the PCC.

The impact of spectral basis set composition on estimated levels of cingulate glutamate and its associations with different personality traits

BACKGROUND

1H-MRS is increasingly used in basic and clinical research to explain brain function and alterations respectively. In psychosis research it is now one of the main tools to investigate imbalances in the glutamatergic system. Interestingly, however, the findings are extremely variable even within patients of similar disease states. One reason may be the variability in analysis strategies, despite suggestions for standardization. Therefore, our study aimed to investigate the extent to which the basis set configuration- which metabolites are included in the basis set used for analysis- would affect the spectral fit and estimated glutamate (Glu) concentrations in the anterior cingulate cortex (ACC), and whether any changes in levels of glutamate would be associated with psychotic-like experiences and autistic traits.

METHODS

To ensure comparability, we utilized five different exemplar basis sets, used in research, and two different analysis tools, r-based spant applying the ABfit method and Osprey using the LCModel.

RESULTS

Our findings revealed that the types of metabolites included in the basis set significantly affected the glutamate concentration. We observed that three basis sets led to more consistent results across different concentration types (i.e., absolute Glu in mol/kg, Glx (glutamate + glutamine), Glu/tCr), spectral fit and quality measurements. Interestingly, all three basis sets included phosphocreatine. Importantly, our findings also revealed that glutamate levels were differently associated with both schizotypal and autistic traits depending on basis set configuration and analysis tool, with the same three basis sets showing more consistent results.

CONCLUSIONS

Our study highlights that scientific results may be significantly altered depending on the choices of metabolites included in the basis set, and with that emphasizes the importance of carefully selecting the configuration of the basis set to ensure accurate and consistent results, when using MR spectroscopy. Overall, our study points out the need for standardized analysis pipelines and reporting.

Uncertainty propagation in absolute metabolite quantification for in vivo MRS of the human brain

PURPOSE

Absolute spectral quantification is the standard method for deriving estimates of the concentration from metabolite signals measured using in vivo proton MRS (1 H-MRS). This method is often reported with minimum variance estimators, specifically the Cramér-Rao lower bound (CRLB) of the metabolite signal amplitude's scaling factor from linear combination modeling. This value serves as a proxy for SD and is commonly reported in MRS experiments. Characterizing the uncertainty of absolute quantification, however, depends on more than simply the CRLB. The uncertainties of metabolite-specific (T1m , T2m ), reference-specific (T1ref , T2ref ), and sequence-specific (TR , TE ) parameters are generally ignored, potentially leading to an overestimation of precision. In this study, the propagation of uncertainty is used to derive a comprehensive estimate of the overall precision of concentrations from an internal reference.

METHODS

The propagated uncertainty is calculated using analytical derivations and Monte Carlo simulations and subsequently analyzed across a set of commonly measured metabolites and macromolecules. The effect of measurement error from experimentally obtained quantification parameters is estimated using published uncertainties and CRLBs from in vivo 1 H-MRS literature.

RESULTS

The additive effect of propagated measurement uncertainty from applied quantification correction factors can result in up to a fourfold increase in the concentration estimate's coefficient of variation compared to the CRLB alone. A case study analysis reveals similar multifold increases across both metabolites and macromolecules.

CONCLUSION

The precision of absolute metabolite concentrations derived from 1 H-MRS experiments is systematically overestimated if the uncertainties of commonly applied corrections are neglected as sources of error.

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.

Linewidth-related bias in modelled concentration estimates from GABA-edited 1H-MRS

J-difference-edited MRS is widely used to study GABA in the human brain. Editing for low-concentration target molecules (such as GABA) typically exhibits lower signal-to-noise ratio (SNR) than conventional non-edited MRS, varying with acquisition region, volume and duration. Moreover, spectral lineshape may be influenced by age-, pathology-, or brain-region-specific effects of metabolite T2, or by task-related blood-oxygen level dependent (BOLD) changes in functional MRS contexts. Differences in both SNR and lineshape may have systematic effects on concentration estimates derived from spectral modelling. The present study characterises the impact of lineshape and SNR on GABA+ estimates from different modelling algorithms: FSL-MRS, Gannet, LCModel, Osprey, spant and Tarquin. Publicly available multi-site GABA-edited data (222 healthy subjects from 20 sites; conventional MEGA-PRESS editing; TE = 68 ms) were pre-processed with a standardised pipeline, then filtered to apply controlled levels of Lorentzian and Gaussian linebroadening and SNR reduction. Increased Lorentzian linewidth was associated with a 2-5% decrease in GABA+ estimates per Hz, observed consistently (albeit to varying degrees) across datasets and most algorithms. Weaker, often opposing effects were observed for Gaussian linebroadening. Variations are likely caused by differing baseline parametrization and lineshape constraints between models. Effects of linewidth on other metabolites (e.g., Glx and tCr) varied, suggesting that a linewidth confound may persist after scaling to an internal reference. These findings indicate a potentially significant confound for studies where linewidth may differ systematically between groups or experimental conditions, e.g. due to T2 differences between brain regions, age, or pathology, or varying T2* due to BOLD-related changes. We conclude that linewidth effects need to be rigorously considered during experimental design and data processing, for example by incorporating linewidth into statistical analysis of modelling outcomes or development of appropriate lineshape matching algorithms.

Multi-channel magnetic resonance spectroscopy graphical user interface (McMRSGUI)

This work introduces an open-sourced graphical user interface (GUI) software enabling the combination of multi-channel magnetic resonance spectroscopy data with different literature-based methods for the improvement of the quality and reliability of combined spectra. The multi-channel magnetic resonance spectroscopy graphical user interface (McMRSGUI) is a MATLAB-based spectroscopy processing GUI equipped to load multi-channel MRS data, pre-process, combine, and export combined data for evaluation with open-source quantification software (jMRUI). A literature-based, decision-tree process was incorporated into the combination type selection to serve as a guide to minimize spectral distortion in selecting between weighting methods. Multi-channel, simulated spectra were combined with the different combination techniques and evaluated for spectral distortion to validate the code. The incorporation of the combination methods into a single processing software enables multi-channel magnetic resonance spectroscopy (MRS) data to be combined and compared for improved spectral quality with little user knowledge of combination techniques. Through the spectral peak distortion simulation of the combination methods, combined signal-to-noise ratio (SNR) values from the literature were verified. The spectral peak distortion simulation provides a secondary tool for researchers to estimate the spectral SNR levels when spectral distortion could occur and use this knowledge to further guide the selection of their combination technique. The McMRSGUI provides a software toolkit for evaluating multi-channel MRS data and their combination. Simulations evaluating spectral distortion at different noise levels were performed for each combination method to validate the GUI and demonstrate a method for researchers to assess the combined SNR levels at which they could be introducing spectral distortion.

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 1 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.

In vivo cortical glutathione response to oral fumarate therapy in relapsing-remitting multiple sclerosis: A single-arm open-label phase IV trial using 7-Tesla 1H MRS

BACKGROUND

This is an open-label, single-arm, single-center pilot study using 7-Tesla in vivo proton magnetic resonance spectroscopy (1H MRS) to measure brain cortical glutathione concentration at baseline before and during the use of oral fumarates as a disease-modifying therapy for multiple sclerosis. The primary endpoint of this research was the change in prefrontal cortex glutathione concentration relative to a therapy-naïve baseline after one year of oral fumarate therapy.

METHODS

Brain glutathione concentrations were examined by 1H MRS in single prefrontal and occipital cortex cubic voxels (2.5 × 2.5 × 2.5 cm3) before and during initiation of oral fumarate therapy (120 mg b.i.d. for 7 days and 240 mg b.i.d. thereafter). Additional measurements of related metabolites glutamate, glutamine, myoinositol, total N-acetyl aspartate, and total choline were also acquired in voxels centered on the same regions. Seven relapsing-remitting multiple sclerosis patients (4 f / 3 m, age range 28-50 years, mean age 40 years) naïve to fumarate therapy were scanned at pre-therapy baseline and after 1, 3, 6 and 12 months of therapy. A group of 8 healthy volunteers (4 f / 4 m, age range 33-48 years, mean age 41 years) was also scanned at baseline and Month 6 to characterize 1H-MRS measurement reproducibility over a comparable time frame.

RESULTS

In the multiple sclerosis cohort, general linear models demonstrated a significant positive linear relationship between prefrontal glutathione and time either linearly across all time points (+0.05 ± 0.02 mM/month, t(27) = 2.6, p = 0.02) or specifically for factor variable Month 12 (+0.6 ± 0.3 mM/12 months, t(24) = 2.2, p = 0.04) relative to baseline. No such effects of time on glutathione concentration were demonstrated in the occipital cortex or in the healthy volunteer group. Changes in occipital total choline were further observed in the multiple sclerosis cohort as well as prefrontal total choline and occipital glutamine and myoinositol in the control cohort throughout the study duration.

CONCLUSIONS

While the open-label single-arm pilot study design and abbreviated control series cannot support firm conclusions about the influence of oral fumarate therapy independent of test-retest factors or normal biological variation in a state of either health or disease, these results do justify further investigation at a larger scale into the potential relationship between prefrontal cortex glutathione increases and oral fumarate therapy in relapsing-remitting multiple sclerosis.

Event-related functional magnetic resonance spectroscopy

Proton-Magnetic Resonance Spectroscopy (MRS) is a non-invasive brain imaging technique used to measure the concentration of different neurochemicals. "Single-voxel" MRS data is typically acquired across several minutes, before individual transients are averaged through time to give a measurement of neurochemical concentrations. However, this approach is not sensitive to more rapid temporal dynamics of neurochemicals, including those that reflect functional changes in neural computation relevant to perception, cognition, motor control and ultimately behaviour. In this review we discuss recent advances in functional MRS (fMRS) that now allow us to obtain event-related measures of neurochemicals. Event-related fMRS involves presenting different experimental conditions as a series of trials that are intermixed. Critically, this approach allows spectra to be acquired at a time resolution in the order of seconds. Here we provide a comprehensive user guide for event-related task designs, choice of MRS sequence, analysis pipelines, and appropriate interpretation of event-related fMRS data. We raise various technical considerations by examining protocols used to quantify dynamic changes in GABA, the primary inhibitory neurotransmitter in the brain. Overall, we propose that although more data is needed, event-related fMRS can be used to measure dynamic changes in neurochemicals at a temporal resolution relevant to computations that support human cognition and behaviour.

The future is 2D: spectral-temporal fitting of dynamic MRS data provides exponential gains in precision over conventional approaches

PURPOSE

Many MRS paradigms produce 2D spectral-temporal datasets, including diffusion-weighted, functional, and hyperpolarized and enriched (carbon-13, deuterium) experiments. Conventionally, temporal parameters-such as T₂ , T₁ , or diffusion constants-are assessed by first fitting each spectrum independently and subsequently fitting a temporal model (1D fitting). We investigated whether simultaneously fitting the entire dataset using a single spectral-temporal model (2D fitting) would improve the precision of the relevant temporal parameter.

METHODS

We derived a Cramer Rao lower bound for the temporal parameters for both 1D and 2D approaches for 2 experiments: a multi-echo experiment designed to estimate metabolite T₂ s, and a functional MRS experiment designed to estimate fractional change ( δ $$ \delta $$ ) in metabolite concentrations. We investigated the dependence of the relative standard deviation (SD) of T₂ in multi-echo and δ $$ \delta $$ in functional MRS.

RESULTS

When peaks were spectrally distant, 2D fitting improved precision by approximately 20% relative to 1D fitting, regardless of the experiment and other parameter values. These gains increased exponentially as peaks drew closer. Dependence on temporal model parameters was weak to negligible.

CONCLUSION

Our results strongly support a 2D approach to MRS fitting where applicable, and particularly in nuclei such as hydrogen and deuterium, which exhibit substantial spectral overlap.

Magnetic resonance imaging of cardiac metabolism in heart failure: how far have we come?

As one of the highest energy consumer organs in the body, the heart requires tremendous amount of adenosine triphosphate (ATP) to maintain its continuous mechanical work. Fatty acids, glucose, and ketone bodies are the primary fuel source of the heart to generate ATP with perturbations in ATP generation possibly leading to contractile dysfunction. Cardiac metabolic imaging with magnetic resonance imaging (MRI) plays a crucial role in understanding the dynamic metabolic changes occurring in the failing heart, where the cardiac metabolism is deranged. Also, targeting and quantifying metabolic changes in vivo noninvasively is a promising approach to facilitate diagnosis, determine prognosis, and evaluate therapeutic response. Here, we summarize novel MRI techniques used for detailed investigation of cardiac metabolism in heart failure including magnetic resonance spectroscopy (MRS), hyperpolarized MRS, and chemical exchange saturation transfer based on evidence from preclinical and clinical studies and to discuss the potential clinical application in heart failure.

Multiple sclerosis diagnosis and phenotype identification by multivariate classification of in vivo frontal cortex metabolite profiles

Multiple sclerosis (MS) is a heterogeneous autoimmune disease for which diagnosis continues to rely on subjective clinical judgment over a battery of tests. Proton magnetic resonance spectroscopy (¹H MRS) enables the noninvasive in vivo detection of multiple small-molecule metabolites and is therefore in principle a promising means of gathering information sufficient for multiple sclerosis diagnosis and subtype classification. Here we show that supervised classification using ¹H-MRS-visible normal-appearing frontal cortex small-molecule metabolites alone can indeed differentiate individuals with progressive MS from control (held-out validation sensitivity 79% and specificity 68%), as well as between relapsing and progressive MS phenotypes (held-out validation sensitivity 84% and specificity 74%). Post hoc assessment demonstrated the disproportionate contributions of glutamate and glutamine to identifying MS status and phenotype, respectively. Our finding establishes ¹H MRS as a viable means of characterizing progressive multiple sclerosis disease status and paves the way for continued refinement of this method as an auxiliary or mainstay of multiple sclerosis diagnostics.

Comparison of seven modelling algorithms for γ-aminobutyric acid-edited proton magnetic resonance spectroscopy

Edited MRS sequences are widely used for studying γ-aminobutyric acid (GABA) in the human brain. Several algorithms are available for modelling these data, deriving metabolite concentration estimates through peak fitting or a linear combination of basis spectra. The present study compares seven such algorithms, using data obtained in a large multisite study. GABA-edited (GABA+, TE = 68 ms MEGA-PRESS) data from 222 subjects at 20 sites were processed via a standardised pipeline, before modelling with FSL-MRS, Gannet, AMARES, QUEST, LCModel, Osprey and Tarquin, using standardised vendor-specific basis sets (for GE, Philips and Siemens) where appropriate. After referencing metabolite estimates (to water or creatine), systematic differences in scale were observed between datasets acquired on different vendors' hardware, presenting across algorithms. Scale differences across algorithms were also observed. Using the correlation between metabolite estimates and voxel tissue fraction as a benchmark, most algorithms were found to be similarly effective in detecting differences in GABA+. An interclass correlation across all algorithms showed single-rater consistency for GABA+ estimates of around 0.38, indicating moderate agreement. Upon inclusion of a basis set component explicitly modelling the macromolecule signal underlying the observed 3.0 ppm GABA peaks, single-rater consistency improved to 0.44. Correlation between discrete pairs of algorithms varied, and was concerningly weak in some cases. Our findings highlight the need for consensus on appropriate modelling parameters across different algorithms, and for detailed reporting of the parameters adopted in individual studies to ensure reproducibility and meaningful comparison of outcomes between different studies.

Proton magnetic resonance spectroscopy thermometry: Impact of separately acquired full water or partially suppressed water data on quantification and measurement error

In proton magnetic resonance spectroscopy (¹ H MRS) thermometry, separately acquired full water and partially suppressed water are commonly used for measuring temperature. This paper compares these two approaches. Single-voxel ¹ H MRS data were collected on a 3-T GE scanner from 26 human subjects. Every subject underwent five continuous MRS sessions, each separated by a 2-min phase. Each MRS session lasted 13 min and consisted of two free induction decays (FIDs) without water suppression (with full water [FW or w]) and 64 FIDs with partial water suppression (with partially suppressed water [PW or w']). Frequency differences between the two FWs, the first two PWs, the second FW and the first PW (FW₂ , PW₁ ), or between averaged water ( wav' ) and N-acetylaspartate (NAA), were measured. Intrasubject and intersubject variations of the frequency differences were used as a metric for the error in temperature measurement. The intrasubject variations of frequency differences between FW₂ and PW₁fw2-fw1' , calculated from the five MRS sessions for each subject, were larger than those between the two FWs or between the first two PWs (p = 1.54 x 10^-4 and p = 1.72 x 10^-4 , respectively). The mean values of intrasubject variations of fw2-fw1' for all subjects were 4.7 and 4.5 times those of fw2-fw1 and fw2'-fw1' , respectively. The intrasubject variations of the temperatures based on frequency differences, fw2-fNAA or ( fw1'-fNAA ), were about 2.5 times greater than those based on averaged water and NAA frequencies (fwav'-fNAA ). The mean temperature measured from (fwav'-fNAA ) (n = 26) was 0.29°C lower than that measured from fw2-fNAA and was 0.83°C higher than that from ( fw1'-fNAA ). It was concluded that the use of separately acquired unsuppressed or partially suppressed water signals may result in large errors in frequency and, consequently, temperature measurement.

In vivo brain MR spectroscopy in gliomas: clinical and pre-clinical chances

Purpose

Gliomas, the most common primary brain tumours, have recently been re-classified incorporating molecular aspects with important clinical, prognostic, and predictive implications. Concurrently, the reprogramming of metabolism, altering intracellular and extracellular metabolites affecting gene expression, differentiation, and the tumour microenvironment, is increasingly being studied, and alterations in metabolic pathways are becoming hallmarks of cancer. Magnetic resonance spectroscopy (MRS) is a complementary, non-invasive technique capable of quantifying multiple metabolites. The aim of this review focuses on the methodology and analysis techniques in proton MRS (1H MRS), including a brief look at X-nuclei MRS, and on its perspectives for diagnostic and prognostic biomarkers in gliomas in both clinical practice and preclinical research.

Methods

PubMed literature research was performed cross-linking the following key words: glioma, MRS, brain, in-vivo, human, animal model, clinical, pre-clinical, techniques, sequences, 1H, X-nuclei, Artificial Intelligence (AI), hyperpolarization.

Results

We selected clinical works (n = 51), preclinical studies (n = 35) and AI MRS application papers (n = 15) published within the last two decades. The methodological papers (n = 62) were taken into account since the technique first description.

Conclusions

Given the development of treatments targeting specific cancer metabolic pathways, MRS could play a key role in allowing non-invasive assessment for patient diagnosis and stratification, predicting and monitoring treatment responses and prognosis. The characterization of gliomas through MRS will benefit of a wide synergy among scientists and clinicians of different specialties within the context of new translational competences. Head coils, MRI hardware and post-processing analysis progress, advances in research, experts’ consensus recommendations and specific professionalizing programs will make the technique increasingly trustworthy, responsive, accessible.

In Vitro 31P MR Chemical Shifts of In Vivo-Detectable Metabolites at 3T as a Basis Set for a Pilot Evaluation of Skeletal Muscle and Liver 31P Spectra with LCModel Software

Most in vivo ³¹P MR studies are realized on 3T MR systems that provide sufficient signal intensity for prominent phosphorus metabolites. The identification of these metabolites in the in vivo spectra is performed by comparing their chemical shifts with the chemical shifts measured in vitro on high-field NMR spectrometers. To approach in vivo conditions at 3T, a set of phantoms with defined metabolite solutions were measured in a 3T whole-body MR system at 7.0 and 7.5 pH, at 37 °C. A free induction decay (FID) sequence with and without ¹H decoupling was used. Chemical shifts were obtained of phosphoenolpyruvate (PEP), phosphatidylcholine (PtdC), phosphocholine (PC), phosphoethanolamine (PE), glycerophosphocholine (GPC), glycerophosphoetanolamine (GPE), uridine diphosphoglucose (UDPG), glucose-6-phosphate (G6P), glucose-1-phosphate (G1P), 2,3-diphosphoglycerate (2,3-DPG), nicotinamide adenine dinucleotide (NADH and NAD⁺), phosphocreatine (PCr), adenosine triphosphate (ATP), adenosine diphosphate (ADP), and inorganic phosphate (Pi). The measured chemical shifts were used to construct a basis set of ³¹P MR spectra for the evaluation of ³¹P in vivo spectra of muscle and the liver using LCModel software (linear combination model). Prior knowledge was successfully employed in the analysis of previously acquired in vivo data.

In vivo evidence of differential frontal cortex metabolic abnormalities in progressive and relapsing-remitting multiple sclerosis

The pathophysiology of progressive multiple sclerosis remains elusive, significantly limiting available disease-modifying therapies. Proton MRS (¹ H-MRS) enables in vivo measurement of small molecules implicated in multiple sclerosis, but its application to key metabolites glutamate, γ-aminobutyric acid (GABA), and glutathione has been sparse. We employed, at 7 T, a previously validated ¹ H-MRS protocol to measure glutamate, GABA, and glutathione, as well as glutamine, N-acetyl aspartate, choline, and myoinositol, in the frontal cortex of individuals with relapsing-remitting (N = 26) or progressive (N = 21) multiple sclerosis or healthy control adults (N = 25) in a cross-sectional analysis. Only individuals with progressive multiple sclerosis demonstrated reduced glutamate (F_2,65  = 3.424, p = 0.04; 12.40 ± 0.62 mM versus control 13.17 ± 0.95 mM, p = 0.03) but not glutamine (F_2,65  = 0.352, p = 0.7; 4.71 ± 0.35 mM versus control 4.84 ± 0.42 mM), reduced GABA (F_2,65  = 3.89, p = 0.03; 1.29 ± 0.23 mM versus control 1.47 ± 0.25 mM, p = 0.05), and possibly reduced glutathione (F_2,65  = 0.352, p = 0.056; 2.23 ± 0.46 mM versus control 2.51 ± 0.48 mM, p < 0.1). As a group, multiple sclerosis patients demonstrated significant negative correlations between disease duration and glutamate or GABA (ρ = -0.4, p = 0.02) but not glutamine or glutathione. Alone, only relapsing-remitting multiple sclerosis patients exhibited a significant negative correlation between disease duration and GABA (ρ = -0.5, p = 0.03). Taken together, these results indicate that frontal cortex metabolism is differentially disturbed in progressive and relapsing-remitting multiple sclerosis.

Hippocampal single-voxel MR spectroscopy with a long echo time at 3 T using semi-LASER sequence

The hippocampus is one of the most challenging brain regions for proton MR spectroscopy (MRS) applications. Moreover, quantification of J-coupled species such as myo-inositol (m-Ins) and glutamate + glutamine (Glx) is affected by the presence of macromolecular background. While long echo time (TE) MRS eliminates the macromolecules, it also decreases the m-Ins and Glx signal and, as a result, these metabolites are studied mainly with short TE. Here, we investigate the feasibility of reproducibly measuring their concentrations at a long TE of 120 ms, using a semi-adiabatic localization by adiabatic selective refocusing (sLASER) sequence, as this sequence was recently recommended as a standard for clinical MRS. Gradient offset-independent adiabatic refocusing pulses were implemented, and an optimal long TE for the detection of m-Ins and Glx was determined using the T2 relaxation times of macromolecules. Metabolite concentrations and their coefficients of variation (CVs) were obtained for a 3.4-mL voxel centered on the left hippocampus on 3-T MR systems at two different sites with three healthy subjects (aged 32.5 ± 10.2 years [mean ± standard deviation]) per site, with each subject scanned over two sessions, and with each session comprising three scans. Concentrations of m-Ins, choline, creatine, Glx and N-acetyl-aspartate were 5.4 ± 1.5, 1.7 ± 0.2, 5.8 ± 0.3, 11.6 ± 1.2 and 5.9 ± 0.4 mM (mean ± standard deviation), respectively. Their respective mean within-session CVs were 14.5% ± 5.9%, 6.5% ± 5.3%, 6.0% ± 3.4%, 10.6% ± 6.2% and 3.5% ± 1.4%, and their mean within-subject CVs were 17.8% ± 18.2%, 7.5% ± 6.3%, 7.4% ± 6.4%, 12.4% ± 5.3% and 4.8% ± 3.0%. The between-subject CVs were 25.0%, 12.3%, 5.3%, 10.7% and 6.4%, respectively. Hippocampal long-TE sLASER single voxel spectroscopy can provide macromolecule-independent assessment of all major metabolites including Glx and m-Ins.

Are Cramér-Rao lower bounds an accurate estimate for standard deviations in in vivo magnetic resonance spectroscopy?

Due to inherent time constraints for in vivo experiments, it is infeasible to repeat multiple MRS scans to estimate standard deviations on the desired measured parameters. As such, the Cramér-Rao lower bounds (CRLBs) have become the routine method to approximate standard deviations for in vivo experiments. Cramér-Rao lower bounds, however, as the name suggests, are theoretically a lower bound on the standard deviation and it is not clear if and under what circumstances this approximation is valid. Realistic synthetic 3 T spectra were used to investigate the relationship between estimated CRLBs, true CRLBs and standard deviations. Here we demonstrate that, although the CRLBs are theoretically truly a lower bound on the standard deviation (not an equality) for the problem typically encountered in quantification, they are still an adequate approximation to standard deviation as long as the model perfectly characterizes the data. In the case when the macromolecule basis deviates from the measured macromolecules it was shown that the CRLBs can deviate from standard deviations by approximately 50% for N-acetylaspartic acid, creatine and glutamate and of the order of 100% or more for myo-inositol and γ-aminobutyric acid. In the case when the model perfectly reflects the data the CRLBs are within approximately 10% of standard deviations for all metabolites. The result of the CRLB being within 10% of standard deviations means that, for an accurate model, novel quantification methods such as machine learning or deep learning will not be able to obtain substantially more precise estimates for the desired parameters than traditional maximum-likelihood estimation.

Magnetic Resonance Spectroscopy of Hepatic Fat from Fundamental to Clinical Applications

The number of individuals suffering from fatty liver is increasing worldwide, leading to interest in the noninvasive study of liver fat. Magnetic resonance spectroscopy (MRS) is a powerful tool that allows direct quantification of metabolites in tissue or areas of interest. MRS has been applied in both research and clinical studies to assess liver fat noninvasively in vivo. MRS has also demonstrated excellent performance in liver fat assessment with high sensitivity and specificity compared to biopsy and other imaging modalities. Because of these qualities, MRS has been generally accepted as the reference standard for the noninvasive measurement of liver steatosis. MRS is an evolving technique with high potential as a diagnostic tool in the clinical setting. This review aims to provide a brief overview of the MRS principle for liver fat assessment and its application, and to summarize the current state of MRS study in comparison to other techniques.

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.