Metabolite T2 relaxation times decrease across the adult lifespan in a large multi-site cohort.

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 T2 changes with aging for multiple of the major metabolites. Here, we comprehensively investigate correlations between T2 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 T2. Results: Older age was correlated with shorter T2 for tNAA, tCr3.0, tCr3.9, tCho, Glx, and tissue water in CSO and PCC; rs = -0.21 to -0.65, all p<0.05, FDR-corrected for multiple comparisons. These associations remained statistically significant when controlling for cortical atrophy. T2 values did not differ across the adult lifespan for mI. By region, T2 values were longer in the CSO for tNAA, tCr3.0, tCr3.9, Glx, and tissue water and longer in the PCC for tCho and mI. Conclusion: These findings underscore the importance of considering metabolite T2 changes with aging in MRS quantification. We suggest that future 3T work utilize the equations presented here to estimate age-specific T2 values instead of relying on uniform default values.

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 on the default white matter and gray matter T2 reference values in Osprey and 2) shorter WM and GM T2 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). T2 correction based on shorter T2 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.

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.

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.

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.

Impact of acquisition and modeling parameters on the test-retest reproducibility of edited GABA.

Literature values vary widely for within-subject test-retest reproducibility of gamma-aminobutyric acid (GABA) measured with edited magnetic resonance spectroscopy (MRS). Reasons for this variation remain unclear. Here, we tested whether three acquisition parameters-(1) sequence complexity (two-experiment MEscher-GArwood Point RESolved Spectroscopy [MEGA-PRESS] vs. four-experiment Hadamard Encoding and Reconstruction of MEGA-Edited Spectroscopy [HERMES]); (2) editing pulse duration (14 vs. 20 ms); and (3) scanner frequency drift (interleaved water referencing [IWR] turned ON vs. OFF)-and two linear combination modeling variations-(1) three different coedited macromolecule models (called "1to1GABA", "1to1GABAsoft", and "3to2MM" in the Osprey software package); and (2) 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 (mean age: 30.8 ± 9.5 years; 10 males). Test and retest scans were separated by removing the participant from the scanner for 5-10 min. 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). We identified the best test-retest reproducibility following postprocessing with a composite model of the 0.9- and 3-ppm macromolecules ("3to2MM"); this model performed particularly well for the HERMES data. Furthermore, 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. However, reproducibility did not consistently differ for MEGA-PRESS compared with HERMES, for 14- compared with 20-ms editing pulses, or for IWR-ON versus IWR-OFF. These results highlight the importance of model selection for edited MRS studies of GABA+, particularly for clinical studies that focus on individual patient differences in GABA+ or changes following an intervention.

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.

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.

Application of a 1H Brain MRS Benchmark Dataset to Deep Learning for Out-of-Voxel Artifacts.

Neural networks are potentially valuable for many of the challenges associated with MRS data. The purpose of this manuscript is to describe the AGNOSTIC dataset, which contains 259,200 synthetic 1H MRS examples for training and testing neural networks. AGNOSTIC was created using 270 basis sets that were simulated across 18 field strengths and 15 echo times. The synthetic examples were produced to resemble in vivo brain data with combinations of metabolite, macromolecule, residual water signals, and noise. To demonstrate the utility, we apply AGNOSTIC to train two Convolutional Neural Networks (CNNs) to address out-of-voxel (OOV) echoes. A Detection Network was trained to identify the point-wise presence of OOV echoes, providing proof of concept for real-time detection. A Prediction Network was trained to reconstruct OOV echoes, allowing subtraction during post-processing. Complex OOV signals were mixed into 85% of synthetic examples to train two separate CNNs for the detection and prediction of OOV signals. AGNOSTIC is available through Dryad and all Python 3 code is available through GitHub. The Detection network was shown to perform well, identifying 95% of OOV echoes. Traditional modeling of these detected OOV signals was evaluated and may prove to be an effective method during linear-combination modeling. The Prediction Network greatly reduces OOV echoes within FIDs and achieved a median log10 normed-MSE of -1.79, an improvement of almost two orders of magnitude.

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.

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.

More than one-half of the variance in in vivo proton MR spectroscopy metabolite estimates is common to all metabolites.

The present study characterized associations among brain metabolite levels, applying bivariate and multivariate (i.e., factor analysis) statistical methods to total creatine (tCr)-referenced estimates of the major Point RESolved Spectroscopy (PRESS) proton MR spectroscopy (1 H-MRS) metabolites (i.e., total NAA/tCr, total choline/tCr, myo-inositol/tCr, glutamate + glutamine/tCr) acquired at 3 T from medial parietal lobe in a large (n = 299), well-characterized international cohort of healthy volunteers. Results supported the hypothesis that 1 H-MRS-measured metabolite estimates are moderately intercorrelated (Mr = 0.42, SDr = 0.11, ps < 0.001), with more than one-half (i.e., 57%) of the total variability in metabolite estimates explained by a single common factor. Older age was significantly associated with lower levels of the identified common metabolite variance (CMV) factor (ß = -0.09, p = 0.048), despite not being associated with levels of any individual metabolite. Holding CMV factor levels constant, females had significantly lower levels of total choline (i.e., unique metabolite variance; ß = -0.19, p < 0.001), mirroring significant bivariate correlations between sex and total choline reported previously. Supplementary analysis of water-referenced metabolite estimates (i.e., including tCr/water) demonstrated lower, although still substantial, intercorrelations among metabolites, with 37% of total metabolite variance explained by a single common factor. If replicated, these results would suggest that applied 1 H-MRS researchers shift their analytical framework from examining bivariate associations between individual metabolites and specialty-dependent (e.g., clinical, research) variables of interest (e.g., using t-tests) to examining multivariable (i.e., covariate) associations between multiple metabolites and specialty-dependent variables of interest (e.g., using multiple regression).

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

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

Neurometabolic timecourse of healthy aging.

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.

Quantitative T1-relaxation corrected metabolite mapping of 12 metabolites in the human brain at 9.4 T.

Magnetic resonance spectroscopic imaging (MRSI) is a non-invasive imaging modality that enables observation of metabolites. Applications of MRSI for neuroimaging have shown promise for monitoring and detecting various diseases. This study builds off previously developed techniques of short TR, 1H FID MRSI by correcting for T1-weighting of the metabolites and utilizing an internal water reference to produce quantitative (mmol kg-1) metabolite maps. This work reports and shows quantitative metabolite maps for 12 metabolites for a single slice. Voxel-specific T1-corrections for water are common in MRSI studies; however, most studies use either averaged T1-relaxation times to correct for T1-weighting of metabolites or omit this correction step entirely. This work employs the use of voxel-specific T1-corrections for metabolites in addition to water. Utilizing averaged T1-relaxation times for metabolites can bias metabolite maps for metabolites that have strong differences between T1-relaxation for GM and WM (i.e. Glu). This work systematically compares quantitative metabolite maps to single voxel quantitative results and qualitatively compares metabolite maps to previous works.

Relaxation-corrected macromolecular model enables determination of 1 H longitudinal T1 -relaxation times and concentrations of human brain metabolites at 9.4T.

PURPOSE: Ultrahigh field MRS has improved characterization of the neurochemical profile. To compare results obtained at 9.4T to those from lower field strengths, it is of interest to quantify the concentrations of metabolites measured. Thus, measuring T1 -relaxation times is necessary to correct for T1 -weighting that occurs in acquisitions for single-voxel spectroscopy and spectroscopic imaging. A macromolecule (MM) simulation model was developed to fit MM contributions to the short TE inversion series used to measure T1 -relaxation times. METHODS: An inversion series with seven time points was acquired with metabolite-cycled STEAM to estimate T1 -relaxation times of metabolites. A short TE was employed in this study to retain signals from metabolites with short T2 -relaxation times and J-couplings. The underlying macromolecule spectrum was corrected by developing a sequence-specific, relaxation-corrected simulated MM model. Quantification of metabolite peaks was performed using internal water referencing and relaxation corrections. RESULTS: T1 -relaxation times for metabolites range from approximately 750 to approximately 2000 ms and approximately 1000 to approximately 2400 ms in gray matter (GM)- and white matter (WM)- rich voxels, respectively. Quantification of metabolites was compared between GM and WM voxels, as well as between results that used a simulated MM spectrum against those that used an experimentally acquired MM spectrum. Metabolite concentrations are reported in mmol/kg quantities. CONCLUSION: T1 -relaxation times are reported for nonsinglet resonances for the first time at 9.4T by use of a MM simulation model to account for contributions from the MM spectrum. In addition to T1 -relaxation times, quantification results of metabolites from GM- and WM-rich voxels are reported.

The macromolecular MR spectrum does not change with healthy aging.

PURPOSE: To acquire the mobile macromolecule (MM) spectrum from healthy participants, and to investigate changes in the signals with age and sex. METHODS: 102 volunteers (49 M/53 F) between 20 and 69 years were recruited for in vivo data acquisition in the centrum semiovale (CSO) and posterior cingulate cortex (PCC). Spectral data were acquired at 3T using PRESS localization with a voxel size of 30 × 26 × 26 mm3 , pre-inversion (TR/TI 2000/600 ms) and CHESS water suppression. Metabolite-nulled spectra were modeled to eliminate residual metabolite signals, which were then subtracted out to yield a "clean" MM spectrum using the Osprey software. Pearson's correlation coefficient was calculated between integrals and age for the 14 MM signals. One-way ANOVA was performed to determine differences between age groups. An independent t-test was carried out to determine differences between sexes. RESULTS: MM spectra were successfully acquired in 99 (CSO) and 96 (PCC) of 102 subjects. No significant correlations were seen between age and MM signals. One-way ANOVA also suggested no age-group differences for any MM peak (all p > .004). No differences were observed between sex groups. WM and GM voxel fractions showed a significant (p < .05) negative linear association with age in the WM-predominant CSO (R = -0.29) and GM-predominant PCC regions (R = -0.57) respectively while CSF increased significantly with age in both regions. CONCLUSION: Our findings suggest that a pre-defined MM basis function can be used for linear combination modeling of metabolite data from different age and sex groups.

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

PURPOSE: Accurate and precise MRS fitting is crucial for metabolite concentration quantification of 1 H-MRS spectra. LCModel, a spectral fitting software, has shown to have certain limitations to perform advanced spectral fitting by previous literature. Herein, we propose an open-source spectral fitting algorithm with adaptive spectral baseline determination and more complex cost functions. THEORY: The MRS spectra are characterized by several parameters, which reflect the environment of the contributing metabolites, properties of the acquisition sequence, or additional disturbances. Fitting parameters should accurately describe these parameters. Baselines are also a major contributor to MRS spectra, in which smoothness of the spline baselines used for fitting can be adjusted based on the properties of the spectra. Three different cost functions used for the minimization problem were also investigated. METHODS: The newly developed ProFit-1D fitting algorithm is systematically evaluated for simulations of several types of possible in vivo parameter variations. Although accuracy and precision are tested with simulated spectra, spectra measured in vivo at 9.4 T are used for testing precision using subsets of averages. ProFit-1D fitting results are also compared with LCModel. RESULTS: Both ProFit-1D and LCModel fitted the spectra well with induced parameter and baseline variations. ProFit-1D proved to be more accurate than LCModel for simulated spectra. However, LCModel showed a somewhat increased precision for some spectral simulations and for in vivo data. CONCLUSION: The open-source ProFit-1D fitting algorithm demonstrated high accuracy while maintaining precise metabolite concentration quantification. Finally, through the newly proposed cost functions, new ways to improve fitting were shown.

In vivo characterization of downfield peaks at 9.4 T: T2 relaxation times, quantification, pH estimation, and assignments.

PURPOSE: Relaxation times are a valuable asset when determining spectral assignments. In this study, apparent T2 relaxation times ( T2app ) of downfield peaks are reported in the human brain at 9.4 T and are used to guide spectral assignments of some downfield metabolite peaks. METHODS: Echo time series of downfield metabolite spectra were acquired at 9.4 T using a metabolite-cycled semi-LASER sequence. Metabolite spectral fitting was performed using LCModel V6.3-1L while fitting a pH sweep to estimate the pH of the homocarnosine (hCs) imidazole ring. T2app were calculated by fitting the resulting relative amplitudes of the peaks to a mono-exponential decay across the TE series. Furthermore, estimated tissue concentrations of molecules were calculated using the relaxation times and internal water as a reference. RESULTS: T2app of downfield metabolites are reported within a range from 16 to 32 ms except for homocarnosine with T2app of 50 ms. Correcting T2app for exchange rates ( T2corr ) resulted in relaxation times between 20 and 33 ms. The estimated pH values based on hCs imidazole range from 7.07 to 7.12 between subjects. Furthermore, analyzing the linewidths of the downfield peaks and their T2app contribution led to possible peak assignments. CONCLUSION: T2app relaxation times were longer for the assigned metabolite peaks compared to the unassigned peaks. Tissue pH estimation in vivo with proton MRS and simultaneous quantification of amide protons at 8.30 ± 0.15 ppm is likely possible. Based on concentration, linewidth, and exchange rates measurements, tentative peak assignments are discussed for adenosine triphosphate (ATP), N-acetylaspartylglutamate (NAAG), and urea.