Neurochemical Predictors of Generalized Learning Induced by Brain Stimulation and Training.

Methods of cognitive enhancement for humans are most impactful when they generalize across tasks. However, the extent to which such "transfer" is possible via interventions is widely debated. In addition, the contribution of excitatory and inhibitory processes to such transfer is unknown. Here, in a large-scale neuroimaging individual differences study with humans (both sexes), we paired multitasking training and noninvasive brain stimulation (transcranial direct current stimulation, tDCS) over multiple days and assessed performance across a range of paradigms. In addition, we varied tDCS dosage (1.0 and 2.0 mA), electrode montage (left or right prefrontal regions), and training task (multitasking vs a control task) and assessed GABA and glutamate concentrations via ultrahigh field 7T magnetic resonance spectroscopy. Generalized benefits were observed in spatial attention, indexed by visual search performance, when multitasking training was combined with 1.0 mA stimulation targeting either the left or right prefrontal cortex (PFC). This transfer effect persisted for ∼30 d post intervention. Critically, the transferred benefits associated with right prefrontal tDCS were predicted by pretraining concentrations of glutamate in the PFC. Thus, the effects of this combined stimulation and training protocol appear to be linked predominantly to excitatory brain processes.

Glutamate measurements using edited MRS.

PURPOSE: To demonstrate J-difference coediting of glutamate using Hadamard encoding and reconstruction of Mescher-Garwood-edited spectroscopy (HERMES). METHODS: Density-matrix simulations of HERMES (TE 80 ms) and 1D J-resolved (TE 31-229 ms) of glutamate (Glu), glutamine (Gln), γ-aminobutyric acid (GABA), and glutathione (GSH) were performed. HERMES comprised four sub-experiments with editing pulses applied as follows: (A) 1.9/4.56 ppm simultaneously (ONGABA /ONGSH ); (B) 1.9 ppm only (ONGABA /OFFGSH ); (C) 4.56 ppm only (OFFGABA /ONGSH ); and (D) 7.5 ppm (OFFGABA /OFFGSH ). Phantom HERMES and 1D J-resolved experiments of Glu were performed. Finally, in vivo HERMES (20-ms editing pulses) and 1D J-resolved (TE 31-229 ms) experiments were performed on 137 participants using 3 T MRI scanners. LCModel was used for quantification. RESULTS: HERMES simulation and phantom experiments show a Glu-edited signal at 2.34 ppm in the Hadamard sum combination A+B+C+D with no overlapping Gln signal. The J-resolved simulations and phantom experiments show substantial TE modulation of the Glu and Gln signals across the TEs, whose average yields a well-resolved Glu signal closely matching the Glu-edited signal from the HERMES sum spectrum. In vivo quantification of Glu show that the two methods are highly correlated (p < 0.001) with a bias of ∼10%, along with similar between-subject coefficients of variation (HERMES/TE-averaged: ∼7.3%/∼6.9%). Other Hadamard combinations produce the expected GABA-edited (A+B-C-D) or GSH-edited (A-B+C-D) signal. CONCLUSION: HERMES simulation and phantom experiments show the separation of Glu from Gln. In vivo HERMES experiments yield Glu (without Gln), GABA, and GSH in a single MRS scan.

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.

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

PURPOSE: The J-difference edited γ-aminobutyric acid (GABA) signal is contaminated by other co-edited signals-the largest of which originates from co-edited macromolecules (MMs)-and is consequently often reported as "GABA+." MM signals are broader and less well-characterized than the metabolites, and are commonly approximated using a Gaussian model parameterization. Experimentally measured MM signals are a consensus-recommended alternative to parameterized modeling; however, they are relatively under-studied in the context of edited MRS. METHODS: To address this limitation in the literature, we have acquired GABA-edited MEGA-PRESS data with pre-inversion to null metabolite signals in 13 healthy controls. An experimental MM basis function was derived from the mean across subjects. We further derived a new parameterization of the MM signals from the experimental data, using multiple Gaussians to accurately represent their observed asymmetry. The previous single-Gaussian parameterization, mean experimental MM spectrum and new multi-Gaussian parameterization were compared in a three-way analysis of a public MEGA-PRESS dataset of 61 healthy participants. RESULTS: Both the experimental MMs and the multi-Gaussian parameterization exhibited reduced fit residuals compared to the single-Gaussian approach (p = 0.034 and p = 0.031, respectively), suggesting they better represent the underlying data than the single-Gaussian parameterization. Furthermore, both experimentally derived models estimated larger MM fractional contribution to the GABA+ signal for the experimental MMs (58%) and multi-Gaussian parameterization (58%), compared to the single-Gaussian approach (50%). CONCLUSIONS: Our results indicate that single-Gaussian parameterization of edited MM signals is insufficient and that both experimentally derived GABA+ spectra and their parameterized replicas improve the modeling of GABA+ spectra.

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.

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.

Image downsampling expedited adaptive least-squares (IDEAL) fitting improves intravoxel incoherent motion (IVIM) analysis in the human kidney.

PURPOSE: To improve the reliability of intravoxel incoherent motion (IVIM) model parameter estimation for the DWI in the kidney using a novel image downsampling expedited adaptive least-squares (IDEAL) approach. METHODS: The robustness of IDEAL was investigated using simulated DW-MRI data corrupted with different levels of Rician noise. Subsequently, the performance of the proposed method was tested by fitting bi- and triexponential IVIM model to in vivo renal DWI data acquired on a clinical 3 Tesla MRI scanner and compared to conventional approaches (fixed D* and segmented fitting). RESULTS: The numerical simulations demonstrated that the IDEAL algorithm provides robust estimates of the IVIM parameters in the presence of noise (SNR of 20) as indicated by relatively low absolute percentage bias (maximal sMdPB <20%) and normalized RMSE (maximal RMSE <28%). The analysis of the in vivo data showed that the IDEAL-based IVIM parameter maps were less noisy and more visually appealing than those obtained using the fixed D* and segmented methods. Further, coefficients of variation for nearly all IVIM parameters were significantly reduced in cortex and medulla for IDEAL-based biexponential (coefficients of variation: 4%-50%) and triexponential (coefficients of variation: 7.5%-75%) IVIM modelling compared to the segmented (coefficients of variation: 4%-120%) and fixed D* (coefficients of variation: 17%-174%) methods, reflecting greater accuracy of this method. CONCLUSION: The proposed fitting algorithm yields more robust IVIM parameter estimates and is less susceptible to poor SNR than the conventional fitting approaches. Thus, the IDEAL approach has the potential to improve the reliability of renal DW-MRI analysis for clinical applications.

Impact of gradient scheme and non-linear shimming on out-of-voxel echo artifacts in edited MRS.

Out-of-voxel (OOV) signals are common spurious echo artifacts in MRS. These signals often manifest in the spectrum as very strong "ripples," which interfere with spectral quantification by overlapping with targeted metabolite resonances. Dephasing optimization through coherence order pathway selection (DOTCOPS) gradient schemes are algorithmically optimized to suppress all potential alternative coherence transfer pathways (CTPs), and should suppress unwanted OOV echoes. In addition, second-order shimming uses non-linear gradient fields to maximize field homogeneity inside the voxel, which unfortunately increases the diversity of local gradient fields outside of the voxel. Given that strong local spatial B0 gradients can refocus unintended CTPs, it is possible that OOVs are less prevalent when only linear first-order shimming is applied. Here we compare the size of unwanted OOV signals in Hadamard-edited (HERMES) data acquired with either a local gradient scheme (which we refer to here as "Shared") or DOTCOPS, and with first- or second-order shimming. We collected data from 15 healthy volunteers in two brain regions (voxel size 30 × 26 × 26 mm3 ) from which it is challenging to acquire MRS data: medial prefrontal cortex and left temporal cortex. Characteristic OOV echoes were seen in both GABA- and GSH-edited spectra for both brain regions, gradient schemes, and shimming approaches. A linear mixed-effect model revealed a statistically significant difference in the average residual based on the gradient scheme in both GABA- (p < 0.001) and GSH-edited (p < 0.001) spectra: that is, the DOTCOPS gradient scheme resulted in smaller OOV artifacts compared with the Shared scheme. There were no significant differences in OOV artifacts associated with shimming method. Thus, these results suggest that the DOTCOPS gradient scheme for J-difference-edited PRESS acquisitions yields spectra with smaller OOV echo artifacts than the Shared gradient scheme implemented in a widely disseminated editing sequence.

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.

In vivo spectral editing of phosphorylethanolamine.

PURPOSE: To demonstrate J-difference editing of phosphorylethanolamine (PE) with chemical shifts at 3.22 (PE3.22 ) and 3.98 (PE3.98 ) ppm, and compare the merits of two editing strategies. METHODS: Density-matrix simulations of MEGA-PRESS (Mescher-Garwood PRESS) for PE were performed at TEs ranging from 80 to 200 ms in steps of 2 ms, applying 20-ms editing pulses (ON/OFF) at (1) 3.98/7.5 ppm to detect PE3.22 and (2) 3.22/7.5 ppm to detect PE3.98 . Phantom experiments were performed using a PE phantom to validate simulation results. Ten subjects were scanned using a Philips 3T MRI scanner at TEs of 90 ms and 110 ms to edit PE3.22 and PE3.98 . Osprey was used for data processing, modeling, and quantification. RESULTS: Simulations show substantial TE modulation of the intensity and shape of the edited signals due to coupling evolution. Simulated and phantom integrals suggest that TEs of 110 ms and 90 ms were optimal for the edited detection of PE3.22 and PE3.98 , respectively. Phantom results indicated strong agreement with the simulated spectra and integrals. In vivo quantification of the PE3.22 /total creatine and PE3.98 /total creatine concentration ratio yielded values of 0.26 ± 0.04 (between-subject coefficient of variation [CV]: 15.4%) and 0.18 ± 0.04 (CV: 22.8%), respectively, at TE = 90 ms, and 0.24 ± 0.02 (CV: 8.2%) and 0.23 ± 0.04 (CV: 18.0%), respectively, at TE = 110 ms. CONCLUSION: Simulations and in vivo MEGA-PRESS of PE demonstrate that both PE3.22 and PE3.98 are potential candidates for editing, but PE3.22 at TE = 110 ms yields lower variation across TEs.

Influence of editing pulse flip angle on J-difference MR spectroscopy.

PURPOSE: To investigate the editing-pulse flip angle (FA) dependence of editing efficiency and ultimately to maximize the edited signal of commonly edited MR spectroscopy (MRS) signals, such as gamma-aminobutyric acid (GABA) and lactate. METHODS: Density-matrix simulations were performed for a range of spin systems to find the editing-pulse FA for maximal editing efficiency. Simulations were confirmed by phantom experiments and in vivo measurements in 10 healthy participants using a 3T Philips scanner. Four MEGA-PRESS in vivo measurements targeting GABA+ and lactate were performed, comparing the conventional editing-pulse FA (FA = 180°) to the optimal one suggested by simulations (FA = 210°). RESULTS: Simulations and phantom experiments show that edited GABA and lactate signals are maximal at FA = 210°. Compared to conventional editing (FA = 180°), in vivo signals from GABA+ and lactate signals increase on average by 8.5% and 9.3%, respectively. CONCLUSION: Increasing the FA of editing-pulses in the MEGA-PRESS experiment from 180° to 210° increases the edited signals from GABA+ and lactate by about 9% in vivo.

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

PURPOSE: The purpose of this study is to present a cloud-based spectral simulation tool "MRSCloud," which allows MRS users to simulate a vendor-specific and sequence-specific basis set online in a convenient and time-efficient manner. This tool can simulate basis sets for GE, Philips, and Siemens MR scanners, including conventional acquisitions and spectral editing schemes with PRESS and semi-LASER localization at 3 T. METHODS: The MRSCloud tool was built on the spectral simulation functionality in the FID-A software package. We added three extensions to accelerate computation (ie, one-dimensional projection method, coherence pathways filters, and precalculation of propagators). The RF waveforms were generated based on vendors' generic pulse shapes and timings. Simulations were compared within MRSCloud using different numbers of spatial resolution (21 × 21, 41 × 41, and 101 × 101). Simulated metabolite basis functions from MRSCloud were compared with those generated by the generic FID-A and MARSS, and a phantom-acquired basis set from LCModel. Intraclass correlation coefficients were calculated to measure the agreement between individual metabolite basis functions. Statistical analysis was performed using R in RStudio. RESULTS: Simulation time for a full PRESS basis set is approximately 11 min on the server. The interclass correlation coefficients ICCs were at least 0.98 between MRSCloud and FID-A and were at least 0.96 between MRSCloud and MARSS. The interclass correlation coefficients between simulated MRSCloud basis spectra and acquired LCModel basis spectra were lowest for glutamine at 0.68 and highest for N-acetylaspartate at 0.96. CONCLUSIONS: Substantial reductions in runtime have been achieved. High ICC values indicated that the accelerating features are running correctly and produce comparable and accurate basis sets.

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.

Comparison of linear combination modeling strategies for edited magnetic resonance spectroscopy at 3 T.

J-difference-edited spectroscopy is a valuable approach for the in vivo detection of γ-aminobutyric-acid (GABA) with magnetic resonance spectroscopy (MRS). A recent expert consensus article recommends linear combination modeling (LCM) of edited MRS but does not give specific details regarding implementation. This study explores different modeling strategies to adapt LCM for GABA-edited MRS. Sixty-one medial parietal lobe GABA-edited MEGA-PRESS spectra from a recent 3-T multisite study were modeled using 102 different strategies combining six different approaches to account for co-edited macromolecules (MMs), three modeling ranges, three baseline knot spacings, and the use of basis sets with or without homocarnosine. The resulting GABA and GABA+ estimates (quantified relative to total creatine), the residuals at different ranges, standard deviations and coefficients of variation (CVs), and Akaike information criteria, were used to evaluate the models' performance. Significantly different GABA+ and GABA estimates were found when a well-parameterized MM3co basis function was included in the model. The mean GABA estimates were significantly lower when modeling MM3co , while the CVs were similar. A sparser spline knot spacing led to lower variation in the GABA and GABA+ estimates, and a narrower modeling range-only including the signals of interest-did not substantially improve or degrade modeling performance. Additionally, the results suggest that LCM can separate GABA and the underlying co-edited MM3co . Incorporating homocarnosine into the modeling did not significantly improve variance in GABA+ estimates. In conclusion, GABA-edited MRS is most appropriately quantified by LCM with a well-parameterized co-edited MM3co basis function with a constraint to the nonoverlapped MM0.93 , in combination with a sparse spline knot spacing (0.55 ppm) and a modeling range of 0.5-4 ppm.

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.

Edited magnetic resonance spectroscopy in the neonatal brain.

J-difference-edited spectroscopy is a valuable approach for the detection of low-concentration metabolites with magnetic resonance spectroscopy (MRS). Currently, few edited MRS studies are performed in neonates due to suboptimal signal-to-noise ratio, relatively long acquisition times, and vulnerability to motion artifacts. Nonetheless, the technique presents an exciting opportunity in pediatric imaging research to study rapid maturational changes of neurotransmitter systems and other metabolic systems in early postnatal life. Studying these metabolic processes is vital to understanding the widespread and rapid structural and functional changes that occur in the first years of life. The overarching goal of this review is to provide an introduction to edited MRS for neonates, including the current state-of-the-art in editing methods and editable metabolites, as well as to review the current literature applying edited MRS to the neonatal brain. Existing challenges and future opportunities, including the lack of age-specific reference data, are also discussed.

Frequency drift in MR spectroscopy at 3T.

PURPOSE: Heating of gradient coils and passive shim components is a common cause of instability in the B0 field, especially when gradient intensive sequences are used. The aim of the study was to set a benchmark for typical drift encountered during MR spectroscopy (MRS) to assess the need for real-time field-frequency locking on MRI scanners by comparing field drift data from a large number of sites. METHOD: A standardized protocol was developed for 80 participating sites using 99 3T MR scanners from 3 major vendors. Phantom water signals were acquired before and after an EPI sequence. The protocol consisted of: minimal preparatory imaging; a short pre-fMRI PRESS; a ten-minute fMRI acquisition; and a long post-fMRI PRESS acquisition. Both pre- and post-fMRI PRESS were non-water suppressed. Real-time frequency stabilization/adjustment was switched off when appropriate. Sixty scanners repeated the protocol for a second dataset. In addition, a three-hour post-fMRI MRS acquisition was performed at one site to observe change of gradient temperature and drift rate. Spectral analysis was performed using MATLAB. Frequency drift in pre-fMRI PRESS data were compared with the first 5:20 minutes and the full 30:00 minutes of data after fMRI. Median (interquartile range) drifts were measured and showed in violin plot. Paired t-tests were performed to compare frequency drift pre- and post-fMRI. A simulated in vivo spectrum was generated using FID-A to visualize the effect of the observed frequency drifts. The simulated spectrum was convolved with the frequency trace for the most extreme cases. Impacts of frequency drifts on NAA and GABA were also simulated as a function of linear drift. Data from the repeated protocol were compared with the corresponding first dataset using Pearson's and intraclass correlation coefficients (ICC). RESULTS: Of the data collected from 99 scanners, 4 were excluded due to various reasons. Thus, data from 95 scanners were ultimately analyzed. For the first 5:20 min (64 transients), median (interquartile range) drift was 0.44 (1.29) Hz before fMRI and 0.83 (1.29) Hz after. This increased to 3.15 (4.02) Hz for the full 30 min (360 transients) run. Average drift rates were 0.29 Hz/min before fMRI and 0.43 Hz/min after. Paired t-tests indicated that drift increased after fMRI, as expected (p < 0.05). Simulated spectra convolved with the frequency drift showed that the intensity of the NAA singlet was reduced by up to 26%, 44 % and 18% for GE, Philips and Siemens scanners after fMRI, respectively. ICCs indicated good agreement between datasets acquired on separate days. The single site long acquisition showed drift rate was reduced to 0.03 Hz/min approximately three hours after fMRI. DISCUSSION: This study analyzed frequency drift data from 95 3T MRI scanners. Median levels of drift were relatively low (5-min average under 1 Hz), but the most extreme cases suffered from higher levels of drift. The extent of drift varied across scanners which both linear and nonlinear drifts were observed.

Frequency and phase correction of J-difference edited MR spectra using deep learning.

PURPOSE: To investigate whether a deep learning-based (DL) approach can be used for frequency-and-phase correction (FPC) of MEGA-edited MRS data. METHODS: Two neural networks (1 for frequency, 1 for phase) consisting of fully connected layers were trained and validated using simulated MEGA-edited MRS data. This DL-FPC was subsequently tested and compared to a conventional approach (spectral registration [SR]) and to a model-based SR implementation (mSR) using in vivo MEGA-edited MRS datasets. Additional artificial offsets were added to these datasets to further investigate performance. RESULTS: The validation showed that DL-based FPC was capable of correcting within 0.03 Hz of frequency and 0.4°of phase offset for unseen simulated data. DL-based FPC performed similarly to SR for the unmanipulated in vivo test datasets. When additional offsets were added to these datasets, the networks still performed well. However, although SR accurately corrected for smaller offsets, it often failed for larger offsets. The mSR algorithm performed well for larger offsets, which was because the model was generated from the in vivo datasets. In addition, the computation times were much shorter using DL-based FPC or mSR compared to SR for heavily distorted spectra. CONCLUSION: These results represent a proof of principle for the use of DL for preprocessing MRS data.