Search

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 2D 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 SD 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 FPC 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, for example, long TEs or strong diffusion weighting., (© 2024 International Society for Magnetic Resonance in Medicine.)

Purpose: Relaxation correction is crucial for accurately estimating metabolite concentrations measured using in vivo MRS. However, the majority of MRS quantification routines assume that relaxation values remain constant across the lifespan, despite prior evidence of T 2 changes with aging for multiple of the major metabolites. Here, we comprehensively investigate correlations between T 2 and age in a large, multi-site cohort., Methods: We recruited approximately 10 male and 10 female participants from each decade of life: 18-29, 30-39, 40-49, 50-59, and 60+ y old (n = 101 total). We collected PRESS data at eight TEs (30, 50, 74, 101, 135, 179, 241, and 350 ms) from voxels placed in white-matter-rich centrum semiovale (CSO) and gray-matter-rich posterior cingulate cortex (PCC). We quantified metabolite amplitudes using Osprey and fit exponential decay curves to estimate T 2 ., Results: Older age was correlated with shorter T 2 for tNAA 2.0 , tCr 3.0 , tCr 3.9 , tCho, and tissue water (CSO and PCC), as well as mI and Glx (PCC only); r s  = -0.22 to -0.63, all p < 0.05, false discovery rate (FDR)-corrected. These associations largely remained statistically significant when controlling for cortical atrophy. By region, T 2 values were longer in the CSO for tNAA 2.0 , tCr 3.9 , Glx, and tissue water and longer in the PCC for tCho and mI. T 2 did not differ by region for tCr 3.0 ., Conclusion: These findings underscore the importance of considering metabolite T 2 differences with aging in MRS quantification. We suggest that future 3T work utilize the equations presented here to estimate age-specific T 2 values instead of relying on uniform default values., (© 2024 International Society for Magnetic Resonance in Medicine.)

Positron emission tomography (PET) and magnetic resonance spectroscopy ( 1 H-MRS) are complementary techniques that can be applied to study how proteinopathy and neurometabolism relate to cognitive deficits in preclinical stages of Alzheimer's disease (AD)-mild cognitive impairment (MCI) and late-life depression (LLD). We acquired beta-amyloid (Aβ) PET and 7 T 1 H-MRS measures of GABA, glutamate, glutathione, N-acetylaspartate, N-acetylaspartylglutamate, myo-inositol, choline, and lactate in the anterior and posterior cingulate cortices (ACC, PCC) in 13 MCI and 9 LLD patients, and 13 controls. We used linear regression to examine associations between metabolites, Aβ, and cognitive scores, and whether metabolites and Aβ explained cognitive scores better than Aβ alone. In the ACC, higher Aβ was associated with lower GABA in controls but not MCI or LLD patients, but results depended upon MRS data quality control criteria. Greater variance in California Verbal Learning Test scores was better explained by a model that combined ACC glutamate and Aβ deposition than by models that only included one of these variables. These findings identify preliminary associations between Aβ, neurometabolites, and cognition., (Copyright © 2024 Elsevier Inc. All rights reserved.)

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., (© 2024 International Society for Magnetic Resonance in Medicine.)

The HEALthy Brain and Child Development (HBCD) Study, a multi-site prospective longitudinal cohort study, will examine human brain, cognitive, behavioral, social, and emotional development beginning prenatally and planned through early childhood. The acquisition of multimodal magnetic resonance-based brain development data is central to the study's core protocol. However, application of Magnetic Resonance Imaging (MRI) methods in this population is complicated by technical challenges and difficulties of imaging in early life. Overcoming these challenges requires an innovative and harmonized approach, combining age-appropriate acquisition protocols together with specialized pediatric neuroimaging strategies. The HBCD MRI Working Group aimed to establish a core acquisition protocol for all 27 HBCD Study recruitment sites to measure brain structure, function, microstructure, and metabolites. Acquisition parameters of individual modalities have been matched across MRI scanner platforms for harmonized acquisitions and state-of-the-art technologies are employed to enable faster and motion-robust imaging. Here, we provide an overview of the HBCD MRI protocol, including decisions of individual modalities and preliminary data. The result will be an unparalleled resource for examining early neurodevelopment which enables the larger scientific community to assess normative trajectories from birth through childhood and to examine the genetic, biological, and environmental factors that help shape the developing brain., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. Tobias Kober and Tom Hilbert are employees of Siemens Healthineers International AG, Switzerland. Yulin Chang is an employee of Siemens Medical Solutions USA Inc. Dan Rettmann and Ralph Noeske are employed by GE HealthCare. Guillaume Gilbert, Yansong Zhao, Sandeep Ganji, and Maarten Versluis are employed by Philips Healthcare. Carina Lucena, Lucky Heisler-Roman, and Dhruman Goradia are employed by PrimeNeuro Inc. Under a license agreement between Philips and the Johns Hopkins University, Dr. van Zijl and the University are entitled to fees related to an imaging device used in the study discussed for publication. Dr. van Zijl also is a paid lecturer for Philips and receives research support from Philips. This arrangement has been reviewed and approved by the Johns Hopkins University in accordance with its conflict of interest policies. Damien Fair is a patent holder on the Framwise Integrated Real-Time Motion Monitoring (FIRMM) software. He is also a co-founder of Turing Medical Technologies, Inc. The nature of this financial interest and the design of the study have been reviewed by two committees at the University of Minnesota. They have put in place a plan to help ensure that this research is not affected by the financial interest. All other authors report no biomedical financial interests or potential conflicts of interest., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Purpose: Metabolite amplitude estimates derived from linear combination modeling of MR spectra depend upon the precise list of constituent metabolite basis functions used (the "basis set"). The absence of clear consensus on the "ideal" composition or objective criteria to determine the suitability of a particular basis set contributes to the poor reproducibility of MRS. In this proof-of-concept study, we demonstrate a novel, data-driven approach for deciding the basis-set composition using Bayesian information criteria (BIC)., Methods: We have developed an algorithm that iteratively adds metabolites to the basis set using iterative modeling, informed by BIC scores. We investigated two quantitative "stopping conditions", referred to as max-BIC and zero-amplitude, and whether to optimize the selection of basis set on a per-spectrum basis or at the group level. The algorithm was tested using two groups of synthetic in-vivo-like spectra representing healthy brain and tumor spectra, respectively, and the derived basis sets (and metabolite amplitude estimates) were compared to the ground truth., Results: All derived basis sets correctly identified high-concentration metabolites and provided reasonable fits of the spectra. At the single-spectrum level, the two stopping conditions derived the underlying basis set with 77-87% accuracy. When optimizing across a group, basis set determination accuracy improved to 84-92%., Conclusion: Data-driven determination of the basis set composition is feasible. With refinement, this approach could provide a valuable data-driven way to derive or refine basis sets, reducing the operator bias of MRS analyses, enhancing the objectivity of quantitative analyses, and increasing the clinical viability of 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. T 2 correction based on the dual-TE water integrals were compared with: 1) T 2 correction based on the default white matter and gray matter T 2 reference values in Osprey and 2) shorter WM and GM T 2 values from recent literature., 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). T 2 correction based on shorter T 2 values showed better agreement to the dual-TE water integral ratio., Conclusions: ISTHMUS facilitated data acquisition and post-processing and reduced operator workload to eliminate potential human error., Competing Interests: Declaration of Competing Interest None., (Copyright © 2024 Elsevier B.V. All rights reserved.)

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 T 1 relaxation across age, it is possible that aging alters T 1 relaxation rates, as is seen for T 2 relaxation. Here, we investigate the age dependence of metabolite T 1 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 T 1 relaxation times for total N-acetyl aspartate at 2.0 ppm (tNAA 2.0 ) and total creatine at 3.0 ppm (tCr 3.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 T 1 and subject age were evaluated. Spearman correlations revealed that estimated T 1 relaxation times of tNAA 2.0 (r s  = -0.27; p < 0.006) and tCr 3.0 (r s  = -0.40; p < 0.001) decreased significantly with age in white-matter-rich CSO, and less steeply for tNAA 2.0 (r s  = -0.228; p = 0.005) and (not significantly for) tCr 3.0 (r s  = -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 T 1 relaxation times change with age. This preliminary study stresses the importance of further work to measure age-normed metabolite T 1 relaxation times for accurate quantification of metabolite levels in studies of aging., (© 2024 John Wiley & Sons Ltd.)

Purpose: Relaxometry, specifically T 1 and T 2 mapping, has become an essential technique for assessing the properties of biological tissues related to various physiological and pathological conditions. Many techniques are being used to estimate T 1 and T 2 relaxation times, ranging from the traditional inversion or saturation recovery and spin-echo sequences to more advanced methods. Choosing the appropriate method for a specific application is critical since the precision and accuracy of T 1 and T 2 measurements are influenced by a variety of factors including the pulse sequence and its parameters, the inherent properties of the tissue being examined, the MRI hardware, and the image reconstruction. The aim of this study is to evaluate and compare the test-retest reproducibility of two advanced MRI relaxometry techniques (Driven Equilibrium Single Pulse Observation of T 1 and T 2 , DESPOT, and 3D Quantification using an interleaved Look-Locker acquisition Sequence with a T 2 preparation pulse, QALAS), for T 1 and T 2 mapping in a healthy volunteer cohort., Methods: 10 healthy volunteers underwent brain MRI at 1.3 mm 3 isotropic resolution, acquiring DESPOT and QALAS data (~11.8 and ~5 minutes duration, including field maps, respectively), test-retest with subject repositioning, on a 3.0 Tesla Philips Ingenia Elition scanner. To reconstruct the T 1 and T 2 maps, we used an equation-based algorithm for DESPOT and a dictionary-based algorithm that incorporates inversion efficiency and B 1 -field inhomogeneity for QALAS. The test-retest reproducibility was assessed using the coefficient of variation (CoV), intraclass correlation coefficient (ICC) and Bland-Altman plots., Results: Our results indicate that both the DESPOT and QALAS techniques demonstrate good levels of test-retest reproducibility for T 1 and T 2 mapping across the brain. Higher whole-brain voxel-to-voxel ICCs are observed in QALAS for T 1 (0.84 ± 0.039) and in DESPOT for T 2 (0.897 ± 0.029). The Bland-Altman plots show smaller bias and variability of T 1 estimates for QALAS (mean of -0.02 s, and upper and lower limits of -0.14 and 0.11 s, 95% CI) than for DESPOT (mean of -0.02 s, and limits of -0.31 and 0.27 s). QALAS also showed less variability (mean 1.08 ms, limits -1.88 to 4.04 ms) for T 2 compared to DESPOT (mean of 2.56 ms, and limits -17.29 to 22.41 ms). The within-subject CoVs for QALAS range from 0.6% ( T 2 in CSF) to 5.8% ( T 2 in GM), while for DESPOT they range from 2.1% ( T 2 in CSF) to 6.7% ( T 2 in GM). The between-subject CoVs for QALAS range from 2.5% ( T 2 in GM) to 12% ( T 2 in CSF), and for DESPOT they range from 3.7% ( T 2 in WM) to 9.3% ( T 2 in CSF)., Conclusion: Overall, QALAS demonstrated better reproducibility for T 1 and T 2 measurements than DESPOT, in addition to reduced acquisition time., Competing Interests: Declaration of Competing Interests No Competing Interests

Purpose: Relaxation correction is crucial for accurately estimating metabolite concentrations measured using in vivo magnetic resonance spectroscopy (MRS). However, the majority of MRS quantification routines assume that relaxation values remain constant across the lifespan, despite prior evidence of T 2 changes with aging for multiple of the major metabolites. Here, we comprehensively investigate correlations between T 2 and age in a large, multi-site cohort., Methods: We recruited approximately 10 male and 10 female participants from each decade of life: 18-29, 30-39, 40-49, 50-59, and 60+ years old ( n =101 total). We collected PRESS data at 8 TEs (30, 50, 74, 101, 135, 179, 241, and 350 ms) from voxels placed in white-matter-rich centrum semiovale (CSO) and gray-matter-rich posterior cingulate cortex (PCC). We quantified metabolite amplitudes using Osprey and fit exponential decay curves to estimate T 2 ., Results: Older age was correlated with shorter T 2 for tNAA, tCr 3.0 , tCr 3.9 , tCho, Glx, and tissue water in CSO and PCC; r s = -0.21 to -0.65, all p <0.05, FDR-corrected for multiple comparisons. These associations remained statistically significant when controlling for cortical atrophy. T 2 values did not differ across the adult lifespan for mI. By region, T 2 values were longer in the CSO for tNAA, tCr 3.0 , tCr 3.9 , Glx, and tissue water and longer in the PCC for tCho and mI., Conclusion: These findings underscore the importance of considering metabolite T 2 changes with aging in MRS quantification. We suggest that future 3T work utilize the equations presented here to estimate age-specific T 2 values instead of relying on uniform default values., Competing Interests: Competing Interests All authors declare that they have no competing interests.

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 3 Tesla scanner., Results: The AN group reported greater anxiety to eat 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.

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., (© 2024 International Society for Autism Research and Wiley Periodicals LLC.)

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., (© 2023 International Society for Magnetic Resonance in Medicine.)

Introduction: Diffusion-weighted magnetic resonance spectroscopy (DW-MRS) offers improved cellular specificity to microstructure-compared to water-based methods alone-but spatial resolution and SNR is severely reduced and slow-diffusing metabolites necessitate higher b -values to accurately characterize their diffusion properties. Ultra-strong gradients allow access to higher b -values per-unit time, higher SNR for a given b -value, and shorter diffusion times, but introduce additional challenges such as eddy-current artefacts, gradient non-uniformity, and mechanical vibrations., Methods: In this work, we present initial DW-MRS data acquired on a 3T Siemens Connectom scanner equipped with ultra-strong (300 mT/m) gradients. We explore the practical issues associated with this manner of acquisition, the steps that may be taken to mitigate their impact on the data, and the potential benefits of ultra-strong gradients for DW-MRS. An in-house DW-PRESS sequence and data processing pipeline were developed to mitigate the impact of these confounds. The interaction of TE, b -value, and maximum gradient amplitude was investigated using simulations and pilot data, whereby maximum gradient amplitude was restricted. Furthermore, two DW-MRS voxels in grey and white matter were acquired using ultra-strong gradients and high b -values., Results: Simulations suggest T 2 -based SNR gains that are experimentally confirmed. Ultra-strong gradient acquisitions exhibit similar artefact profiles to those of lower gradient amplitude, suggesting adequate performance of artefact mitigation strategies. Gradient field non-uniformity influenced ADC estimates by up to 4% when left uncorrected. ADC and Kurtosis estimates for tNAA, tCho, and tCr align with previously published literature., Discussion: In conclusion, we successfully implemented acquisition and data processing strategies for ultra-strong gradient DW-MRS and results indicate that confounding effects of the strong gradient system can be ameliorated, while achieving shorter diffusion times and improved metabolite SNR., Competing Interests: FF was employed by company Siemens Healthcare Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2023 Davies-Jenkins, Döring, Fasano, Kleban, Mueller, Evans, Afzali, Jones, Ronen, Branzoli and Tax.)

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., (© 2023. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.)

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

Purpose: The neurometabolic timecourse of healthy aging is not well-established, in part due to diversity of quantification methodology. In this study, a large structured cross-sectional cohort of male and female subjects throughout adulthood was recruited to investigate neurometabolic changes as a function of age, using consensus-recommended magnetic resonance spectroscopy quantification methods., Methods: 102 healthy volunteers, with approximately equal numbers of male and female participants in each decade of age from the 20s, 30s, 40s, 50s, and 60s, were recruited with IRB approval. MR spectroscopic data were acquired on a 3T MRI scanner. Metabolite spectra were acquired using PRESS localization (TE=30 ms; 96 transients) in the centrum semiovale (CSO) and posterior cingulate cortex (PCC). Water-suppressed spectra were modeled using the Osprey algorithm, employing a basis set of 18 simulated metabolite basis functions and a cohort-mean measured macromolecular spectrum. Pearson correlations were conducted to assess relationships between metabolite concentrations and age for each voxel; Spearman correlations were conducted where metabolite distributions were non-normal. Paired t-tests were run to determine whether metabolite concentrations differed between the PCC and CSO. Finally, robust linear regressions were conducted to assess both age and sex as predictors of metabolite concentrations in the PCC and CSO and separately, to assess age, signal-noise ratio, and full width half maximum (FWHM) linewidth as predictors of metabolite concentrations., Results: Data from four voxels were excluded (2 ethanol; 2 unacceptably large lipid signal). Statistically-significant age*metabolite Pearson correlations were observed for tCho (r(98)=0.33, p<0.001), tCr (r(98)=0.60, p<0.001), and mI (r(98)=0.32, p=0.001) in the CSO and for NAAG (r(98)=0.26, p=0.008), tCho(r(98)=0.33, p<0.001), tCr (r(98)=0.39, p<0.001), and Gln (r(98)=0.21, p=0.034) in the PCC. Spearman correlations for non-normal variables revealed a statistically significant correlation between sI and age in the CSO (r(86)=0.26, p=0.013). No significant correlations were seen between age and tNAA, NAA, Glx, Glu, GSH, PE, Lac, or Asp in either region (all p>0.20). Age associations for tCho, tCr, mI and sI in the CSO and for NAAG, tCho, and tCr in the PCC remained when controlling for sex in robust regressions. CSO NAAG and Asp, as well as PCC tNAA, sI, and Lac were higher in women; PCC Gln was higher in men. When including an age*sex interaction term in robust regression models, a significant age*sex interaction was seen for tCho (F(1,96)=11.53, p=0.001) and GSH (F(1,96)=7.15, p=0.009) in the CSO and tCho (F(1,96)=9.17, p=0.003), tCr (F(1,96)=9.59, p=0.003), mI (F(1,96)=6.48, p=0.012), and Lac (F(1,78)=6.50, p=0.016) in the PCC. In all significant interactions, metabolite levels increased with age in females, but not males. There was a significant positive correlation between linewidth and age. Age relationships with tCho, tCr, and mI in the CSO and tCho, tCr, mI, and sI in the PCC were significant after controlling for linewidth and FWHM in robust regressions., Conclusion: The primary (correlation) results indicated age relationships for tCho, tCr, mI, and sI in the CSO and for NAAG, tCho, tCr, and Gln in the PCC, while no age correlations were found for tNAA, NAA, Glx, Glu, GSH, PE, Lac, or Asp in either region. Our results provide a normative foundation for future work investigating the neurometabolic time course of healthy aging using MRS., Competing Interests: Declaration of Competing Interest All authors declare no conflicts of interest., (Copyright © 2022. Published by Elsevier Inc.)