Frequency drift in MR spectroscopy at 3T

Real-time radiation beam imaging on an MR linear accelerator using quantitative T1 mapping

BACKGROUND

Direct three-dimensional imaging of radiation beams could enable more accurate radiation dosimetry. It has been previously reported that changes in T1-weighted magnetic resonance imaging (MRI) intensity could be observed during radiation due to radiochemical oxygen depletion. Quantitative T1 mapping could increase sensitivity for dosimetry applications.

PURPOSE

We use an MRI linear accelerator (MR-Linac) to visualize radiation delivery through the real-time effects of dose on the spin-lattice magnetic relaxation time (T1) of water. We quantify the relationships between dose, spin-lattice relaxation rates (R1) and dissolved oxygen concentration to further investigate the mechanisms of T1 change.

METHODS

An ultrapure water phantom and a 1% agarose gel phantom were irradiated and imaged on a 1.5 T Elekta Unity MR-Linac. Radiation plans were created using the Monaco treatment planning system. Images were acquired before, during and after radiation. A dual-echo Look-Locker inversion recovery pulse sequence was used for simultaneous dynamic T1/B0 mapping. The change in R1 with respect to dose (∆R1/∆Dose) and the radiochemical oxygen depletion (ROD = ∆O2/∆Dose) were measured. The relaxivity of oxygen (r1,O2 = ∆R1/∆O2) in water was also measured in a separate experiment with samples of various dissolved oxygen concentrations. The minimum measurable dose over a 20-min period was estimated using a single-tailed 99th quantile Student's t-distribution.

RESULTS

Changes to R1 were found to be spatiotemporally correlated to the predicted delivered radiation dose and persisted for at least 1 h after radiation. A complex dose plan could be imaged in the 1% agarose gel phantom, as the gel limits diffusion and convective mixing. In water, the ∆R1/∆Dose was found to be -1.0 × 10-4 s-1/Gy, the r1,O2 was found to be 5.4 × 10-3 s-1/(mg/L), and the ROD was found to be -0.010 (mg/L)/Gy. Both r1,O2 and ROD agree with published values. However, combining these two values yields a predicted ∆R1/∆Dose of -5.4 × 10-5 s-1/Gy, indicating that radiochemical oxygen depletion alone under-predicts the MRI effect. The detection limit of R1 was 1.1 × 10-3 s-1 which corresponded to a single-voxel minimum detectable dose of 11.1 Gy for this specific sequence.

CONCLUSION

Quantitative T1 mapping was used to image radiation dose patterns in real-time in water and agarose gel. Radiochemical oxygen depletion only partially explains the T1 changes measured. Agarose gel could be used as a simple system for three-dimensional patient-specific quality assurance. Future applications may include in vivo dosimetry for FLASH radiotherapy, though improvements in acquisition methods and hardware are likely needed.

Deep learning and its applications in nuclear magnetic resonance spectroscopy

Nuclear Magnetic Resonance (NMR), as an advanced technology, has widespread applications in various fields like chemistry, biology, and medicine. However, issues such as long acquisition times for multidimensional spectra and low sensitivity limit the broader application of NMR. Traditional algorithms aim to address these issues but have limitations in speed and accuracy. Deep Learning (DL), a branch of Artificial Intelligence (AI) technology, has shown remarkable success in many fields including NMR. This paper presents an overview of the basics of DL and current applications of DL in NMR, highlights existing challenges, and suggests potential directions for improvement.

Test-retest precision of brain metabolites in healthy participants using 31P-MRS and 1H MEGA-PRESS on a 3T multi-nuclear MRI system

Background

Magnetic resonance spectroscopy (MRS) enables the non-invasive quantification of brain metabolites, and its reliability is crucial for accurate interpretation of disease state. This study assessed the test-retest precision of phosphorus-31 (31P)-MRS and hydrogen (1H)-MEscher-GArwood Point RESolved Spectroscopy (MEGA-PRESS) in measuring 31P metabolites, γ-aminobutyric acid (GABA), and glutathione (GSH) using a 3T multi-nucleus magnetic resonance imaging (MRI) system.

Methods

In total, 32 participants, who underwent two scanning sessions within three days, using two dimensional (2D)-chemical shift imaging (CSI)-31P-MRS and 1H-MEGA-PRESS sequences, were enrolled in the study. γ-aminobutyric acid and macromolecules (GABA+), glutamate and glutamine (Glx), GSH, and 12 31P metabolites were analyzed using the MATLAB-based tool Gannet and jMRUI software. Precision was assessed based on the coefficients of variation (CVs) and Bland-Altman plots.

Results

The results revealed that potential of hydrogen (pH) and phosphocreatine (PCr) showed the greatest stability as evidenced by low CVs, suggesting reliable measurements across sessions. The adenosine triphosphates (ATPs) showed considerable stability. Conversely, metabolites, such as phosphomonoesters (PMEs) and phosphodiesters (PDEs), located to the left of PCr, showed reduced stability, while glycerophosphatidylcholine (GPTC) had the highest CV, indicating significant variability in clinical practice. Among the various brain regions, intermediate areas such as the temporal lobe and thalamus exhibited greater stability than peripheral regions such as the frontal and occipital lobes. Single-voxel MEGA-PRESS measurements showed that Glx and GABA+ had higher precision than GSH.

Conclusions

Both the 31P-MRS and 1H-MEGA-PRESS sequences showed high precision in measuring brain metabolites, but some metabolites showed higher stability than others. These results are crucial for exploring the clinical and research applications of these methods, and provide a solid foundation for subsequent investigations.

Simultaneous frequency and phase corrections of single-shot MRS data using cross-correlation

PURPOSE

The objective of this study was to propose a novel preprocessing approach to simultaneously correct for the frequency and phase drifts in MRS data using cross-correlation technique.

METHODS

The performance of the proposed method was first investigated at different SNR levels using simulation. Random frequency and phase offsets were added to a previously acquired STEAM human data at 7 T, simulating two different noise levels with and without baseline artifacts. Alongside the proposed spectral cross-correlation (SC) method, three other simultaneous alignment approaches were evaluated. Validation was performed on human brain data at 3 T and mouse brain data at 16.4 T.

RESULTS

The results showed that the SC technique effectively corrects for both small and large frequency and phase drifts, even at low SNR levels. Furthermore, the mean square measurement error of the SC algorithm was comparable to the other three methods used, with much faster processing time. The efficacy of the proposed technique was successfully demonstrated in both human brain MRS data and in a noisy MRS dataset acquired from a small volume-of-interest in the mouse brain.

CONCLUSION

The study demonstrated the availability of a fast and robust technique that accurately corrects for both small and large frequency and phase shifts in MRS.

Dixon-based B0 self-navigation in radial stack-of-stars multi-echo gradient echo imaging

PURPOSE

To develop a Dixon-basedB 0 $$ {\mathrm{B}}_0 $$ self-navigation approach to estimate and correct temporalB 0 $$ {\mathrm{B}}_0 $$ variations in radial stack-of-stars gradient echo imaging for quantitative body MRI.

METHODS

The proposed method estimates temporalB 0 $$ {\mathrm{B}}_0 $$ variations using aB 0 $$ {\mathrm{B}}_0 $$ self-navigator estimated by a graph-cut-based water-fat separation algorithm on the oversampled k-space center. TheB 0 $$ {\mathrm{B}}_0 $$ self-navigator was employed to correct for phase differences between radial spokes (one-dimensional [1D] correction) and to perform a motion-resolved reconstruction to correct spatiotemporal pseudo-periodicB 0 $$ {\mathrm{B}}_0 $$ variations (three-dimensional [3D] correction). Numerical simulations, phantom experiments and in vivo neck scans were performed to evaluate the effects of temporalB 0 $$ {\mathrm{B}}_0 $$ variations on the field-map, proton density fat fraction (PDFF) andT 2 ∗ $$ {\mathrm{T}}_2^{\ast } $$ map, and to validate the proposed method.

RESULTS

TemporalB 0 $$ {\mathrm{B}}_0 $$ variations were found to cause signal loss and phase shifts on the multi-echo images that lead to an underestimation ofT 2 ∗ $$ {\mathrm{T}}_2^{\ast } $$ , while PDFF mapping was less affected. TheB 0 $$ {\mathrm{B}}_0 $$ self-navigator captured slowly varying temporalB 0 $$ {\mathrm{B}}_0 $$ drifts and temporal variations caused by respiratory motion. While the 1D correction effectively correctedB 0 $$ {\mathrm{B}}_0 $$ drifts in phantom studies, it was insufficient in vivo due to 3D spatially varying temporalB 0 $$ {\mathrm{B}}_0 $$ variations with amplitudes of up to 25 Hz at 3 T near the lungs. The proposed 3D correction locally improved the correction of field-map andT 2 ∗ $$ {\mathrm{T}}_2^{\ast } $$ and reduced image artifacts.

CONCLUSION

TemporalB 0 $$ {\mathrm{B}}_0 $$ variations particularly affectT 2 ∗ $$ {\mathrm{T}}_2^{\ast } $$ mapping in radial stack-of-stars imaging. The self-navigation approach can be applied without modifying the MR acquisition to correct forB 0 $$ {\mathrm{B}}_0 $$ drift and physiological motion-inducedB 0 $$ {\mathrm{B}}_0 $$ variations, especially in the presence of fat.

Quantifying brain development in the HEALthy Brain and Child Development (HBCD) Study: The magnetic resonance imaging and spectroscopy protocol

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.

Model-based frequency-and-phase correction of 1H MRS data with 2D linear-combination modeling

PURPOSE

Retrospective frequency-and-phase correction (FPC) methods attempt to remove frequency-and-phase variations between transients to improve the quality of the averaged MR spectrum. However, traditional FPC methods like spectral registration struggle at low SNR. Here, we propose a method that directly integrates FPC into a 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.

Accelerated 3D metabolite T1 mapping of the brain using variable-flip-angle SPICE

PURPOSE

To develop a practical method to enable 3D T1 mapping of brain metabolites.

THEORY AND METHODS

Due to the high dimensionality of the imaging problem underlying metabolite T1 mapping, measurement of metabolite T1 values has been currently limited to a single voxel or slice. This work achieved 3D metabolite T1 mapping by leveraging a recent ultrafast MRSI technique called SPICE (spectroscopic imaging by exploiting spatiospectral correlation). The Ernst-angle FID MRSI data acquisition used in SPICE was extended to variable flip angles, with variable-density sparse sampling for efficient encoding of metabolite T1 information. In data processing, a novel generalized series model was used to remove water and subcutaneous lipid signals; a low-rank tensor model with prelearned subspaces was used to reconstruct the variable-flip-angle metabolite signals jointly from the noisy data.

RESULTS

The proposed method was evaluated using both phantom and healthy subject data. Phantom experimental results demonstrated that high-quality 3D metabolite T1 maps could be obtained and used for correction of T1 saturation effects. In vivo experimental results showed metabolite T1 maps with a large spatial coverage of 240 × 240 × 72 mm3 and good reproducibility coefficients (< 11%) in a 14.5-min scan. The metabolite T1 times obtained ranged from 0.99 to 1.44 s in gray matter and from 1.00 to 1.35 s in white matter.

CONCLUSION

We successfully demonstrated the feasibility of 3D metabolite T1 mapping within a clinically acceptable scan time. The proposed method may prove useful for both T1 mapping of brain metabolites and correcting the T1-weighting effects in quantitative metabolic imaging.

Simultaneous multi-transient linear-combination modeling of MRS data improves uncertainty estimation

PURPOSE

The interest in applying and modeling dynamic MRS has recently grown. Two-dimensional modeling yields advantages for the precision of metabolite estimation in interrelated MRS data. However, it is unknown whether including all transients simultaneously in a 2D model without averaging (presuming a stable signal) performs similarly to one-dimensional (1D) modeling of the averaged spectrum. Therefore, we systematically investigated the accuracy, precision, and uncertainty estimation of both described model approaches.

METHODS

Monte Carlo simulations of synthetic MRS data were used to compare the accuracy and uncertainty estimation of simultaneous 2D multitransient linear-combination modeling (LCM) with 1D-LCM of the average. A total of 2,500 data sets per condition with different noise representations of a 64-transient MRS experiment at six signal-to-noise levels for two separate spin systems (scyllo-inositol and gamma-aminobutyric acid) were analyzed. Additional data sets with different levels of noise correlation were also analyzed. Modeling accuracy was assessed by determining the relative bias of the estimated amplitudes against the ground truth, and modeling precision was determined by SDs and Cramér-Rao lower bounds (CRLBs).

RESULTS

Amplitude estimates for 1D- and 2D-LCM agreed well and showed a similar level of bias compared with the ground truth. Estimated CRLBs agreed well between both models and with ground-truth CRLBs. For correlated noise, the estimated CRLBs increased with the correlation strength for the 1D-LCM but remained stable for the 2D-LCM.

CONCLUSION

Our results indicate that the model performance of 2D multitransient LCM is similar to averaged 1D-LCM. This validation on a simplified scenario serves as a necessary basis for further applications of 2D modeling.

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.

Intra-scan RF power amplifier drift correction

PURPOSE

The drift in radiofrequency (RF) power amplifiers (RFPAs) is assessed and several contributing factors are investigated. Two approaches for prospective correction of drift are proposed and their effectiveness is evaluated.

METHODS

RFPA drift assessment encompasses both intra-pulse and inter-pulse drift analyses. Scan protocols with varying flip angle (FA), RF length, and pulse repetition time (TR) are used to gauge the influence of these parameters on drift. Directional couplers (DICOs) monitor the forward waveforms of the RFPA outputs. DICOs data is stored for evaluation, allowing calculation of correction factors to adjust RFPAs' transmit voltage. Two correction methods, predictive and run-time, are employed: predictive correction necessitates a calibration scan, while run-time correction calculates factors during the ongoing scan.

RESULTS

RFPA drift is indeed influenced by the RF duty-cycle, and in the cases examined with a maximum duty-cycle of 66%, the potential drift is approximately 41% or 15%, depending on the specific RFPA revision. Notably, in low transmit voltage scenarios, FA has minimal impact on RFPA drift. The application of predictive and run-time drift correction techniques effectively reduces the average drift from 10.0% to less than 1%, resulting in enhanced MR signal stability.

CONCLUSION

Utilizing DICO recordings and implementing a feedback mechanism enable the prospective correction of RFPA drift. Having a calibration scan, predictive correction can be utilized with fewer complexity; for enhanced performance, a run-time approach can be employed.

Neurometabolite differences in Autism as assessed with Magnetic Resonance Spectroscopy: A systematic review and meta-analysis

1H-Magnetic Resonance Spectroscopy (MRS) is a non-invasive technique that can be used to quantify the concentrations of metabolites in the brain in vivo. MRS findings in the context of autism are inconsistent and conflicting. We performed a systematic review and meta-analysis of MRS studies measuring glutamate and gamma-aminobutyric acid (GABA), as well as brain metabolites involved in energy metabolism (glutamine, creatine), neural and glial integrity (e.g. n-acetyl aspartate (NAA), choline, myo-inositol) and oxidative stress (glutathione) in autism cohorts. Data were extracted and grouped by metabolite, brain region and several other factors before calculation of standardised effect sizes. Overall, we find significantly lower concentrations of GABA and NAA in autism, indicative of disruptions to the balance between excitation/inhibition within brain circuits, as well as neural integrity. Further analysis found these alterations are most pronounced in autistic children and in limbic brain regions relevant to autism phenotypes. Additionally, we show how study outcome varies due to demographic and methodological factors , emphasising the importance of conforming with standardised consensus study designs and transparent reporting.

Universal dynamic fitting of magnetic resonance spectroscopy

PURPOSE

Dynamic (2D) MRS is a collection of techniques where acquisitions of spectra are repeated under varying experimental or physiological conditions. Dynamic MRS comprises a rich set of contrasts, including diffusion-weighted, relaxation-weighted, functional, edited, or hyperpolarized spectroscopy, leading to quantitative insights into multiple physiological or microstructural processes. Conventional approaches to dynamic MRS analysis ignore the shared information between spectra, and instead proceed by independently fitting noisy individual spectra before modeling temporal changes in the parameters. Here, we propose a universal dynamic MRS toolbox which allows simultaneous fitting of dynamic spectra of arbitrary type.

METHODS

A simple user-interface allows information to be shared and precisely modeled across spectra to make inferences on both spectral and dynamic processes. We demonstrate and thoroughly evaluate our approach in three types of dynamic MRS techniques. Simulations of functional and edited MRS are used to demonstrate the advantages of dynamic fitting.

RESULTS

Analysis of synthetic functional 1H-MRS data shows a marked decrease in parameter uncertainty as predicted by prior work. Analysis with our tool replicates the results of two previously published studies using the original in vivo functional and diffusion-weighted data. Finally, joint spectral fitting with diffusion orientation models is demonstrated in synthetic data.

CONCLUSION

A toolbox for generalized and universal fitting of dynamic, interrelated MR spectra has been released and validated. The toolbox is shared as a fully open-source software with comprehensive documentation, example data, and tutorials.

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.

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.

GABA, glutamatergic dynamics and BOLD contrast assessed concurrently using functional MRS during a cognitive task

A recurring issue in functional neuroimaging is how to link task-driven haemodynamic blood oxygen level dependent functional MRI (BOLD-fMRI) responses to underlying neurochemistry at the synaptic level. Glutamate and γ-aminobutyric acid (GABA), the major excitatory and inhibitory neurotransmitters respectively, are typically measured with MRS sequences separately from fMRI, in the absence of a task. The present study aims to resolve this disconnect, developing acquisition and processing techniques to simultaneously assess GABA, glutamate and glutamine (Glx) and BOLD in relation to a cognitive task, at 3 T. Healthy subjects (N = 81) performed a cognitive task (Eriksen flanker), which was presented visually in a task-OFF, task-ON block design, with individual event onset timing jittered with respect to the MRS readout. fMRS data were acquired from the medial anterior cingulate cortex during task performance, using an adapted MEGA-PRESS implementation incorporating unsuppressed water-reference signals at a regular interval. These allowed for continuous assessment of BOLD activation, through T2 *-related changes in water linewidth. BOLD-fMRI data were additionally acquired. A novel linear model was used to extract modelled metabolite spectra associated with discrete functional stimuli, building on well established processing and quantification tools. Behavioural outcomes from the flanker task, and activation patterns from the BOLD-fMRI sequence, were as expected from the literature. BOLD response assessed through fMRS showed a significant correlation with fMRI, specific to the fMRS-targeted region of interest; fMRS-assessed BOLD additionally correlated with lengthening of response time in the incongruent flanker condition. While no significant task-related changes were observed for GABA+, a significant increase in measured Glx levels (~8.8%) was found between task-OFF and task-ON periods. These findings verify the efficacy of our protocol and analysis pipelines for the simultaneous assessment of metabolite dynamics and BOLD. As well as establishing a robust basis for further work using these techniques, we also identify a number of clear directions for further refinement in future studies.

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

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

Simultaneous multi-transient linear-combination modeling of MRS data improves uncertainty estimation

Purpose

The interest in applying and modeling dynamic MRS has recently grown. 2D modeling yields advantages for the precision of metabolite estimation in interrelated MRS data. However, it is unknown whether including all transients simultaneously in a 2D model without averaging (presuming a stable signal) performs similarly to 1D modeling of the averaged spectrum. Therefore, we systematically investigated the accuracy, precision, and uncertainty estimation of both described model approaches.

Methods

Monte Carlo simulations of synthetic MRS data were used to compare the accuracy and uncertainty estimation of simultaneous 2D multi-transient LCM with 1D-LCM of the average. 2,500 datasets per condition with different noise representations of a 64-transient MRS experiment at 6 signal-to-noise levels for two separate spin systems (scyllo-inositol and GABA) were analyzed. Additional datasets with different levels of noise correlation were also analyzed. Modeling accuracy was assessed by determining the relative bias of the estimated amplitudes against the ground truth, and modeling precision was determined by standard deviations and Cramér-Rao Lower Bounds (CRLB).

Results

Amplitude estimates for 1D- and 2D-LCM agreed well and showed similar level of bias compared to the ground truth. Estimated CRLBs agreed well between both models and with ground truth CRLBs. For correlated noise the estimated CRLBs increased with the correlation strength for the 1D-LCM but remained stable for the 2D-LCM.

Conclusion

Our results indicate that the model performance of 2D multi-transient LCM is similar to averaged 1D-LCM. This validation on a simplified scenario serves as necessary basis for further applications of 2D modeling.

Real-Time Semi-Automated and Automated Voxel Placement using fMRI Targets for Repeated Acquisition Magnetic Resonance Spectroscopy

BACKGROUND

Currently, magnetic resonance spectroscopy (MRS) is dependent on the investigative team to manually prescribe, or demarcate, the desired tissue volume-of-interest. The need for a new method to automate precise voxel placements is warranted to improve the utility and interpretability of MRS data.

NEW METHOD

We propose and validate robust and real-time methods to automate MRS voxel placement using functionally defined coordinates within the prefrontal cortex. Data were collected and analyzed using two independent prospective studies: 1) two independent imaging days with each consisting of a multi-session sandwich design (MRS data only collected on one of the days determined based on scan time) and 2) a longitudinal design. Participants with fibromyalgia syndrome (N = 50) and major depressive disorder (N = 35) underwent neuroimaging. MRS acquisitions were acquired at 3-tesla. Evaluation of the reproducibility of spatial location and tissue segmentation was assessed for: 1) manual, 2) semi-automated, and 3) automated voxel prescription approaches RESULTS: Variability of voxel grey and white matter tissue composition was reduced using automated placement protocols. Spatially, post- to pre-voxel center-of-gravity distance was reduced and voxel overlap increased significantly across datasets using automated compared to manual procedures COMPARISON WITH EXISTING METHODS: Manual prescription, the current standard in the field, can produce inconsistent data across repeated acquisitions. Using automated voxel placement, we found reduced variability and more consistent voxel placement across multiple acquisitions CONCLUSIONS: These results demonstrate the within subject reliability and reproducibility of a method for reducing variability introduced by spatial inconsistencies during MRS acquisitions. The proposed method is a meaningful advance toward improved consistency of MRS data in neuroscience and can be utilized for multi-session and longitudinal studies.

Model-informed unsupervised deep learning approaches to frequency and phase correction of MRS signals

PURPOSE

A supervised deep learning (DL) approach for frequency and phase correction (FPC) of MRS data recently showed encouraging results, but obtaining transients with labels for supervised learning is challenging. This work investigates the feasibility and efficiency of unsupervised deep learning-based FPC.

METHODS

Two novel deep learning-based FPC methods (deep learning-based Cr referencing and deep learning-based spectral registration), which use a priori physics domain knowledge, are presented. The proposed networks were trained, validated, and evaluated using simulated, phantom, and publicly accessible in vivo MEGA-edited MRS data. The performance of our proposed FPC methods was compared with other generally used FPC methods, in terms of precision and time efficiency. A new measure was proposed in this study to evaluate the FPC method performance. The ability of each of our methods to carry out FPC at varying SNR levels was evaluated. A Monte Carlo study was carried out to investigate the performance of our proposed methods.

RESULTS

The validation using low-SNR manipulated simulated data demonstrated that the proposed methods could perform FPC comparably with other methods. The evaluation showed that the deep learning-based spectral registration over a limited frequency range method achieved the highest performance in phantom data. The applicability of the proposed method for FPC of GABA-edited in vivo MRS data was demonstrated. Our proposed networks have the potential to reduce computation time significantly.

CONCLUSIONS

The proposed physics-informed deep neural networks trained in an unsupervised manner with complex data can offer efficient FPC of large MRS data in a shorter time.

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.

Harmonization of multi-scanner in vivo magnetic resonance spectroscopy: ENIGMA consortium task group considerations

Magnetic resonance spectroscopy is a powerful, non-invasive, quantitative imaging technique that allows for the measurement of brain metabolites that has demonstrated utility in diagnosing and characterizing a broad range of neurological diseases. Its impact, however, has been limited due to small sample sizes and methodological variability in addition to intrinsic limitations of the method itself such as its sensitivity to motion. The lack of standardization from a data acquisition and data processing perspective makes it difficult to pool multiple studies and/or conduct multisite studies that are necessary for supporting clinically relevant findings. Based on the experience of the ENIGMA MRS work group and a review of the literature, this manuscript provides an overview of the current state of MRS data harmonization. Key factors that need to be taken into consideration when conducting both retrospective and prospective studies are described. These include (1) MRS acquisition issues such as pulse sequence, RF and B0 calibrations, echo time, and SNR; (2) data processing issues such as pre-processing steps, modeling, and quantitation; and (3) biological factors such as voxel location, age, sex, and pathology. Various approaches to MRS data harmonization are then described including meta-analysis, mega-analysis, linear modeling, ComBat and artificial intelligence approaches. The goal is to provide both novice and experienced readers with the necessary knowledge for conducting MRS data harmonization studies.

Functional MRS studies of GABA and glutamate/Glx - A systematic review and meta-analysis

Functional magnetic resonance spectroscopy (fMRS) can be used to investigate neurometabolic responses to external stimuli in-vivo, but findings are inconsistent. We performed a systematic review and meta-analysis on fMRS studies of the primary neurotransmitters Glutamate (Glu), Glx (Glutamate + Glutamine), and GABA. Data were extracted, grouped by metabolite, stimulus domain, and brain region, and analysed by determining standardized effect sizes. The quality of individual studies was rated. When results were analysed by metabolite type small to moderate effect sizes of 0.29-0.47 (p < 0.05) were observed for changes in Glu and Glx regardless of stimulus domain and brain region, but no significant effects were observed for GABA. Further analysis suggests that Glu, Glx and GABA responses differ by stimulus domain or task and vary depending on the time course of stimulation and data acquisition. Here, we establish effect sizes and directionality of GABA, Glu and Glx response in fMRS. This work highlights the importance of standardised reporting and minimal best practice for fMRS research.

Improving the reproducibility of proton magnetic resonance spectroscopy brain thermometry: Theoretical and empirical approaches

In proton magnetic resonance spectroscopy (¹ H MRS)-based thermometry of brain, averaging temperatures measured from more than one reference peak offers several advantages, including improving the reproducibility (i.e., precision) of the measurement. This paper proposes theoretically and empirically optimal weighting factors to improve the weighted average of temperatures measured from three references. We first proposed concepts of equivalent noise and equivalent signal-to-noise ratio in terms of frequency measurement and a concept of relative frequency that allows the combination of different peaks in a spectrum for improving the precision of frequency measurement. Based on these, we then derived a theoretically optimal weighting factor and proposed an empirical weighting factor, both involving equivalent noise levels, for a weighted average of temperatures measured from three references (i.e., the singlets of NAA, Cr, and Ch in the ¹ H MR spectrum). We assessed these two weighting factors by comparing their errors in measurement of temperatures with the errors of temperatures measured from individual references; we also compared these two new weighting factors with two previously proposed weighting factors. These errors were defined as the standard deviations in repeated measurements or in Monte Carlo studies. Both the proposed theoretical and empirical weighting factors outperformed the two previously proposed weighting factors as well as the three individual references in all phantom and in vivo experiments. In phantom experiments with 4- or 10-Hz line broadening, the theoretical weighting factor outperformed the empirical one, but the latter was superior in all other repeated and Monte Carlo tests performed on phantom and in vivo data. The proposed weighting factors are superior to the two previously proposed weighting factors and can improve the reproducibility of temperature measurement using ¹ H MRS-based thermometry.

On the relationship between GABA+ and glutamate across the brain

Equilibrium between excitation and inhibition (E/I balance) is key to healthy brain function. Conversely, disruption of normal E/I balance has been implicated in a range of central neurological pathologies. Magnetic resonance spectroscopy (MRS) provides a non-invasive means of quantifying in vivo concentrations of excitatory and inhibitory neurotransmitters, which could be used as diagnostic biomarkers. Using the ratio of excitatory and inhibitory neurotransmitters as an index of E/I balance is common practice in MRS work, but recent studies have shown inconsistent evidence for the validity of this proxy. This is underscored by the fact that different measures are often used in calculating E/I balance such as glutamate and Glx (glutamate and glutamine). Here we used a large MRS dataset obtained at ultra-high field (7 T) measured from 193 healthy young adults and focused on two brain regions - prefrontal and occipital cortex - to resolve this inconsistency. We find evidence that there is an inter-individual common ratio between GABA+ (γ-aminobutyric acid and macromolecules) and Glx in the occipital, but not prefrontal cortex. We further replicate the prefrontal result in a legacy dataset (n = 78) measured at high-field (3 T) strength. By contrast, with ultra-high field MRS data, we find extreme evidence that there is a common ratio between GABA+ and glutamate in both prefrontal and occipital cortices, which cannot be explained by participant demographics, signal quality, fractional tissue volume, or other metabolite concentrations. These results are consistent with previous electrophysiological and theoretical work supporting E/I balance. Our findings indicate that MRS-detected GABA+ and glutamate (but not Glx), are a reliable measure of E/I balance .

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

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

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 MM_3co basis function was included in the model. The mean GABA estimates were significantly lower when modeling MM_3co , 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 MM_3co . 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 MM_3co basis function with a constraint to the nonoverlapped MM_0.93 , in combination with a sparse spline knot spacing (0.55 ppm) and a modeling range of 0.5-4 ppm.

Quantifying the excitatory-inhibitory balance: A comparison of SemiLASER and MEGA-SemiLASER for simultaneously measuring GABA and glutamate at 7T

The importance of the excitatory-inhibitory (E/I) balance in a wide range of cognitive and behavioral processes has prompted a commensurate interest in methods for reliably quantifying it. Proton Magnetic Resonance Spectroscopy (¹H-MRS) remains the only method capable of safely and non-invasively measuring the concentrations of the brain's major excitatory (glutamate) and inhibitory (γ-aminobutyric-acid, GABA) neurotransmitters in-vivo. MRS relies on spectral Mescher-Garwood (MEGA) editing techniques at 3T to distinguish GABA from its overlapping resonances. However, with the increased spectral resolution at ultrahigh field strengths of 7T and above, non-edited spectroscopic techniques become potential viable alternatives to MEGA based approaches, and also address some of their shortcomings, such as signal loss, sensitivity to transmitter inhomogeneities and temporal resolution. We present a comprehensive comparison of both edited and non-edited strategies at 7T for simultaneously quantifying glutamate and GABA from the dorsal anterior cingulate cortex (dACC), and evaluate their reproducibility and relative bias. The combined root-mean-square test-retest reproducibility of Glu and GABA (CV_E/I) was as low as 13.3% for unedited MRS at TE=80 ms using SemiLASER localization, while edited MRS at TE=80 ms yielded CV_E/I=20% and 21% for asymmetric and symmetric MEGA editing, respectively. An unedited SemiLASER acquisition using a shorter echo time of TE=42 ms yielded CV_E/I as low as 24.9%. Our results show that non-edited sequences at an echo time of 80 ms provide better reproducibility than either edited sequences at the same TE, or non-edited sequences at a shorter TE of 42 ms. This is supported by numerical simulations and is driven in part by a pseudo-singlet appearance of the GABA multiplets at TE=80 ms, and the excellent spectral resolution at 7T. Our results uphold a transition to non-edited MRS for monitoring the E/I balance at ultrahigh fields, and stress the importance of using a properly-optimized echo time.

Spatially resolved free-induction decay spectroscopy using a 3D ultra-short echo time multi-echo imaging sequence with systematic echo shifting and compensation of B0 field drifts

PURPOSE

Biologically interesting signals can exhibit fast transverse relaxation and frequency shifts compared to free water. For spectral assignment, a ultra-short echo time (UTE) imaging sequence was modified to provide pixel-wise free-induction decay (FID) acquisition.

METHODS

The UTE-FID approach presented relies on a multi-echo 3D spiral UTE sequence with six echoes per radiofrequency (RF) excitation (TE_min 0.05 ms, echo spacing 3 ms). A complex pixel-wise raw data set for FID spectroscopy is obtained by several multi-echo UTE measurements with systematic shifting of the readout by 0.25 or 0.5 ms, until the time domain is filled for 18 or 45 ms. B₀ drifts are compensated by mapping and according phase correction. Autoregressive extrapolation of the signal is performed before Gaussian filtering. This method was applied to a phantom containing collagen-water solutions of different concentrations. To calculate the collagen content, a 19-peak collagen model was extracted from a non-selective FID spectrum (50% collagen solution). Proton-density-collagen-fraction (PDCF) was calculated for 10 collagen solutions (2%-50%). Furthermore, an in vivo UTE-FID spectrum of adipose tissue was recorded.

RESULTS

UTE-FID signal patterns agreed well with the non-spatially selective pulse-acquire FID spectrum from a sphere filled with 50% collagen. Differentiation of collagen solution from distilled water in the PDCF map was possible from 4% collagen concentration for a UTE-FID sequence with 128 × 128 × 64 matrix (voxel size 1 × 1 × 2.85 mm³ ). The mean values of the PDCF correlate linearly with collagen concentration.

CONCLUSION

The presented UTE-FID approach allows pixel-wise raw data acquisition similar to non-spatially selective pulse-acquire spectroscopy. Spatially resolved applications for assessment of spectra of rapidly decaying signals seem feasible.