What is the optimal schedule for multiparametric MRS? A magnetic resonance fingerprinting perspective

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.

Inhibitory mechanisms in the prefrontal-cortex differentially mediate Putamen activity during valence-based learning

Learning from appetitive and aversive stimuli is important for survival. It involves interactions between the prefrontal cortex and subcortical structures, with inhibition playing a crucial role. However, direct evidence for this in humans is limited. Here, we overcome the difficulty of measuring inhibition in the human brain and find that GABA, the main inhibitory neurotransmitter, affects how the dACC interacts with subcortical structures during appetitive and aversive learning differently. We used 7T magnetic resonance spectroscopy (MRS) to track GABA levels in the dACC alongside whole-brain fMRI scans while participants engaged in appetitive and aversive learning tasks. During appetitive learning, dACC GABA levels were negatively correlated with learning performance and BOLD activity measured from the dACC and the Putamen. While under aversive learning, dACC GABA concentration negatively correlated with the functional connectivity between the dACC and the Putamen. Our results show that inhibition in the dACC mediates appetitive and aversive learning in humans through distinct mechanisms.

Diffusion-weighted MR spectroscopy: Consensus, recommendations, and resources from acquisition to modeling

Brain cell structure and function reflect neurodevelopment, plasticity, and aging; and changes can help flag pathological processes such as neurodegeneration and neuroinflammation. Accurate and quantitative methods to noninvasively disentangle cellular structural features are needed and are a substantial focus of brain research. Diffusion-weighted MRS (dMRS) gives access to diffusion properties of endogenous intracellular brain metabolites that are preferentially located inside specific brain cell populations. Despite its great potential, dMRS remains a challenging technique on all levels: from the data acquisition to the analysis, quantification, modeling, and interpretation of results. These challenges were the motivation behind the organization of the Lorentz Center workshop on "Best Practices & Tools for Diffusion MR Spectroscopy" held in Leiden, the Netherlands, in September 2021. During the workshop, the dMRS community established a set of recommendations to execute robust dMRS studies. This paper provides a description of the steps needed for acquiring, processing, fitting, and modeling dMRS data, and provides links to useful resources.

Multi-echo single-shot spectroscopy combined with simultaneous 2D model fitting for fast and accurate measurement of metabolite-specific concentrations and T2 relaxation times

The purpose of the current study was to develop a novel single-voxel MR spectroscopy acquisition scheme to simultaneously determine metabolite-specific concentrations and transverse relaxation times within realistic clinical scan times. Partly truncated multi-TE data are acquired as an echo train in a single acquisition (multi-echo single-shot [MESS]). A 2D multiparametric model fitting approach combines truncated, low-resolved short TE data with fully sampled, highly resolved, longer TE data to yield concentration and T2 estimates for major brain metabolites simultaneously. Cramer-Rao lower bounds (CRLB) are used as a measure of performance. The novel scheme was compared with traditional multi-echo multi-shot methods. In silico, in vitro, and in vivo experiments support the findings. MESS schemes, requiring only 2 min 12 s for the acquisition of three echo times, provide valid concentration and relaxation estimates for multiple metabolites and outperform traditional methods for simultaneous determinations of metabolite-specific T2 s and concentrations, with improvements ranging from 5% to 30% for T2 s and from 10% to 50% for concentrations. However, substantial unsuppressed residual water signals may hamper the method's reproducibility, as observed in an initial experiment setup that prioritizes short TEs with severely truncated acquisition for the benefit of signal-to-noise ratio (SNR). Nevertheless, CRLB have been confirmed to be well suited as design criteria, and within-session repeatability approaches CRLB when residual water is removed in postprocessing by exploiting longer and less truncated data recordings. MESS MRS combined with 2D model fitting promises comparable accuracy, increased precision, or inversely shorter experimental times compared with traditional approaches. However, the optimal design must be investigated as a trade-off between SNR, the truncation factor, and TE batch selections, all of which influence the robustness of estimations.

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.

FSL-MRS: An end-to-end spectroscopy analysis package

PURPOSE

We introduce FSL-MRS, an end-to-end, modular, open-source MRS analysis toolbox. It provides spectroscopic data conversion, preprocessing, spectral simulation, fitting, quantitation, and visualization.

METHODS

The FSL-MRS package is modular. Its programs operate on data in a standard format (Neuroimaging Informatics Technology Initiative [NIfTI]) capable of storing single-voxel and multivoxel spectroscopy, including spatial orientation information. The FSL-MRS toolbox includes tools for preprocessing of raw spectroscopy data, including coil combination, frequency and phase alignment, and filtering. A density matrix simulation program is supplied for generation of basis spectra from simple text-based descriptions of pulse sequences. Fitting is based on linear combination of basis spectra and implements Markov chain Monte Carlo optimization for the estimation of the full posterior distribution of metabolite concentrations. Validation of the fitting is carried out on independently created simulated data, phantom data, and three in vivo human data sets (257 single-voxel spectroscopy and 8 MRSI data sets) at 3 T and 7 T. Interactive HTML reports are automatically generated by processing and fitting stages of the toolbox. The FSL-MRS package can be used on the command line or interactively in the Python language.

RESULTS

Validation of the fitting shows low error in simulation (median error of 11.9%) and in phantom (3.4%). Average correlation between a third-party toolbox (LCModel) and FSL-MRS was high (0.53-0.81) in all three in vivo data sets.

CONCLUSION

The FSL-MRS toolbox is designed to be flexible and extensible to new forms of spectroscopic acquisitions. Custom fitting models can be specified within the framework for dynamic or multivoxel spectroscopy. It is available as part of the FMRIB Software Library.

Quantification of creatine kinase reaction rate in mouse hindlimb using phosphorus-31 magnetic resonance spectroscopic fingerprinting

The goal of this study was to evaluate the accuracy, reproducibility, and efficiency of a 31 P magnetic resonance spectroscopic fingerprinting (31 P-MRSF) method for fast quantification of the forward rate constant of creatine kinase (CK) in mouse hindlimb. The 31 P-MRSF method acquired spectroscopic fingerprints using interleaved acquisition of phosphocreatine (PCr) and γATP with ramped flip angles and a saturation scheme sensitive to chemical exchange between PCr and γATP. Parameter estimation was performed by matching the acquired fingerprints to a dictionary of simulated fingerprints generated from the Bloch-McConnell model. The accuracy of 31 P-MRSF measurements was compared with the magnetization transfer (MT-MRS) method in mouse hindlimb at 9.4 T (n = 8). The reproducibility of 31 P-MRSF was also assessed by repeated measurements. Estimation of the CK rate constant using 31 P-MRSF (0.39 ± 0.03 s-1 ) showed a strong agreement with that using MT-MRS measurements (0.40 ± 0.05 s-1 ). Variations less than 10% were achieved with 2 min acquisition of 31 P-MRSF data. Application of the 31 P-MRSF method to mice subjected to an electrical stimulation protocol detected an increase in CK rate constant in response to stimulation-induced muscle contraction. These results demonstrated the potential of the 31 P-MRSF framework for rapid, accurate, and reproducible quantification of the chemical exchange rate of CK in vivo.

An information theory model for optimizing quantitative magnetic resonance imaging acquisitions

Acquisition parameter selection is currently performed empirically for many quantitative MRI (qMRI) acquisitions. Tuning parameters for different scan times, tissues, and resolutions requires some amount of trial and error. There is an opportunity to quantitatively optimize these acquisition parameters in order to minimize variability of quantitative maps and post-processing techniques such as synthetic image generation. The objective of this work is to introduce and evaluate a quantitative method for selecting parameters that minimize image variability. An information theory framework was developed for this purpose and applied to a 3D-quantification using an interleaved Look-Locker acquisition sequence with T2 preparation pulse (3D-QALAS) signal model for qMRI. In this framework, mutual information is used to measure the information gained by a measurement as a function of acquisition parameters, quantifying the information content of potential acquisitions and allowing informed parameter selection. The information theory framework was tested on artificial data generated from a representative mathematical phantom, measurements acquired on a qMRI multiparametric imaging standard phantom, and in vivo measurements in a human brain. The phantom measurements showed that higher mutual information calculated by the model correlated with smaller coefficient of variation in the reconstructed parametric maps, and in vivo measurements demonstrated that information-based calibration of acquisition parameters resulted in a decrease in parametric map variability consistent with model predictions.