Apparent diffusion coefficients of the five major metabolites measured in the human brain in vivo at 3T

Age-Trajectories of Higher-Order Diffusion Properties of Major Brain Metabolites in Cerebral and Cerebellar Gray Matter Using In Vivo Diffusion-Weighted MR Spectroscopy at 3T

Healthy brain aging involves changes in both brain structure and function, including alterations in cellular composition and microstructure across brain regions. Unlike diffusion-weighted MRI (dMRI), diffusion-weighted MR spectroscopy (dMRS) can assess cell-type specific microstructural changes, providing indirect information on both cell composition and microstructure through the quantification and interpretation of metabolites' diffusion properties. This work investigates age-related changes in the higher-order diffusion properties of total N-Acetyl-aspartate (neuronal biomarker), total choline (glial biomarker), and total creatine (both neuronal and glial biomarker) beyond the classical apparent diffusion coefficient in cerebral and cerebellar gray matter of healthy human brain. Twenty-five subjects were recruited and scanned using a diffusion-weighted semi-LASER sequence in two brain regions-of-interest (ROI) at 3T: posterior-cingulate (PCC) and cerebellar cortices. Metabolites' diffusion was characterized by quantifying metrics from both Gaussian and non-Gaussian signal representations and biophysical models. All studied metabolites exhibited lower apparent diffusivities and higher apparent kurtosis values in the cerebellum compared to the PCC, likely stemming from the higher microstructural complexity of cellular composition in the cerebellum. Multivariate regression analysis (accounting for ROI tissue composition as a covariate) showed slight decrease (or no change) of all metabolites' diffusivities and slight increase of all metabolites' kurtosis with age, none of which statistically significant (p > 0.05). The proposed age-trajectories provide benchmarks for identifying anomalies in the diffusion properties of major brain metabolites which could be related to pathological mechanisms altering both the brain microstructure and cellular composition.

Single-shot diffusion trace spectroscopic imaging using radial echo planar trajectories

PURPOSE

Demonstrate the feasibility and evaluate the performance of single-shot diffusion trace-weighted radial echo planar spectroscopic imaging (Trace DW-REPSI) for quantifying the trace ADC in phantom and in vivo using a 3T clinical scanner.

THEORY AND METHODS

Trace DW-REPSI datasets were acquired in 10 phantom and 10 healthy volunteers, with a maximum b-value of 1601 s/mm2 and diffusion time of 10.75 ms. The self-navigation properties of radial acquisitions were used for corrections of shot-to-shot phase and frequency shift fluctuations of the raw data. In vivo trace ADCs of total NAA (tNAA), total creatine (tCr), and total choline (tCho) extrapolated to pure gray and white matter fractions were compared, as well as trace ADCs estimated in voxels within white or gray matter-dominant regions.

RESULTS

Trace ADCs in phantom show excellent agreement with reported values, and in vivo ADCs agree well with the expected differences between gray and white matter. For tNAA, tCr, and tCho, the trace ADCs extrapolated to pure gray and white matter ranged from 0.18-0.27 and 0.26-0.38 μm2/ms, respectively. In sets of gray and white matter-dominant voxels, the values ranged from 0.21 to 0.27 and 0.24 to 0.31 μm2/ms, respectively. The overestimated trace ADCs from this sequence can be attributed to the short diffusion time.

CONCLUSION

This study presents the first demonstration of the single-shot diffusion trace-weighted spectroscopic imaging sequence using radial echo planar trajectories. The Trace DW-REPSI sequence could provide an estimate of the trace ADC in a much shorter scan time compared to conventional approaches that require three separate measurements.

Chronic neuropathic pain components in whiplash-associated disorders correlate with metabolite concentrations in the anterior cingulate and dorsolateral prefrontal cortex: a consensus-driven MRS re-examination

Introduction

Whiplash injury (WHI) is characterised by a forced neck flexion/extension, which frequently occurs after motor vehicle collisions. Previous studies characterising differences in brain metabolite concentrations and correlations with neuropathic pain (NP) components with chronic whiplash-associated disorders (WAD) have been demonstrated in affective pain-processing areas such as the anterior cingulate cortex (ACC). However, the detection of a difference in metabolite concentrations within these cortical areas with chronic WAD pain has been elusive. In this study, single-voxel magnetic resonance spectroscopy (MRS), following the latest MRSinMRS consensus group guidelines, was performed in the anterior cingulate cortex (ACC), left dorsolateral prefrontal cortex (DLPFC), and occipital cortex (OCC) to quantify differences in metabolite concentrations in individuals with chronic WAD with or without neuropathic pain (NP) components.

Materials and methods

Healthy individuals (n = 29) and participants with chronic WAD (n = 29) were screened with the Douleur Neuropathique 4 Questionnaire (DN4) and divided into groups without (WAD-noNP, n = 15) or with NP components (WAD-NP, n = 14). Metabolites were quantified with LCModel following a single session in a 3 T MRI scanner within the ACC, DLPFC, and OCC.

Results

Participants with WAD-NP presented moderate pain intensity and interference compared with the WAD-noNP group. Single-voxel MRS analysis demonstrated a higher glutamate concentration in the ACC and lower total choline (tCho) in the DLPFC in the WAD-NP versus WAD-noNP group, with no intergroup metabolite difference detected in the OCC. Best fit and stepwise multiple regression revealed that the normalised ACC glutamate/total creatine (tCr) (p = 0.01), DLPFC n-acetyl-aspartate (NAA)/tCr (p = 0.001), and DLPFC tCho/tCr levels (p = 0.02) predicted NP components in the WAD-NP group (ACC r 2 = 0.26, α = 0.81; DLPFC r 2 = 0.62, α = 0.98). The normalised Glu/tCr concentration was higher in the healthy than the WAD-noNP group within the ACC (p < 0.05), but not in the DLPFC or OCC. Neither sex nor age affected key normalised metabolite concentrations related to WAD-NP components when compared to the WAD-noNP group.

Discussion

This study demonstrates that elevated glutamate concentrations within the ACC are related to chronic WAD-NP components, while higher NAA and lower tCho metabolite levels suggest a role for increased neuronal-glial signalling and cell membrane dysfunction in individuals with chronic WAD-NP components.

Single-shot diffusion trace-weighted MR spectroscopy: Comparison with unipolar and bipolar diffusion-weighted point-resolved spectroscopy

This study demonstrates the feasibility and performance of the point-resolved spectroscopy (PRESS)-based, single-shot diffusion trace-weighted sequence in quantifying the trace apparent diffusion coefficient (ADC) in phantom and in vivo using a 3-T MRI/MRS scanner. The single-shot diffusion trace-weighted PRESS sequence was implemented and compared with conventional diffusion-weighted (DW)-PRESS variants using bipolar and unipolar diffusion-sensitizing gradients. Nine phantom datasets were acquired using each sequence, and seven volunteers were scanned in three different brain regions to determine the range and variability of trace ADC values, and to allow a comparison of trace ADCs among the sequences. This sequence results in a comparatively stable range of trace ADC values that are statistically significantly higher than those produced from unipolar and bipolar DW-PRESS sequences. Only total n-acetylaspartate, total creatine, and total choline were reliably estimated in all sequences with Cramér-Rao lower bounds of, at most, 20%. The larger trace ADCs from the single-shot sequences are probably attributable to the shorter diffusion time relative to the other sequences. Overall, this study presents the first demonstration of the single-shot diffusion trace-weighted sequence in a clinical scanner at 3 T. The results show excellent agreement of phantom trace ADCs computed with all sequences, and in vivo ADCs agree well with the expected differences between gray and white matter. The diffusion trace-weighted sequence could provide an estimate of the trace ADC in a shorter scan time (by nearly a factor of 3) compared with conventional DW-PRESS approaches that require three separate orthogonal directions.

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

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

Scyllo-inositol: Transverse relaxation time constant at 3 T and concentration changes associated with aging and alcohol use

The goals of this study were to measure the apparent transverse relaxation time constant, T2 , of scyllo-inositol (sIns) in young and older healthy adults' brains and to investigate the effect of alcohol usage on sIns in young and older healthy adults' brains, using proton magnetic resonance spectroscopy (MRS) at 3 T. Twenty-nine young adults (age 21 ± 1 years) and 24 older adults (age 74 ± 3 years) participated in this study. MRS data were acquired from two brain regions (the occipital cortex and posterior cingulate cortex) at 3 T. The T2 of sIns was measured using a localization by adiabatic selective refocusing (LASER) sequence at various echo times, while the sIns concentrations were measured using a short-echo-time stimulated echo acquisition mode (STEAM) sequence. A trend towards lower T2 relaxation values of sIns in older adults was observed, although these were not significant. sIns concentration was higher with age in both brain regions and was significantly higher in the young when considering alcohol consumption of more than two drinks per week. This study shows that differences in sIns can be found in two distinct regions of the brain across two age groups, potentially reflecting normal aging. In addition, it is important to take into account alcohol consumption when reporting the sIns level in the brain.

Apparent diffusion coefficients of 31P metabolites in the human calf muscle at 7 T

PURPOSE

In this study, we aimed to measure the apparent diffusion coefficients (ADCs) of major phosphorous metabolites in the human calf muscle at 7 T with a diffusion-weighted (DW)-STEAM sequence.

METHODS

A DW-STEAM sequence with bipolar gradients was implemented at 7 T, and DW MR spectra were acquired in three orthogonal directions in the human calf muscle of six healthy volunteers (TE/TM/TR = 15 ms/750 ms/5 s) at three b-values (0, 800, and 1200 s/mm²). Frequency and phase alignments were applied prior to spectral averaging. Averaged DW MR spectra were analyzed with LCModel, and ADCs of ³¹P metabolites were estimated.

RESULTS

Four metabolites (phosphocreatine (PCr), adenosine triphosphate (ATP), inorganic phosphate (Pi) and glycerol phosphorylcholine (GPC)) were quantified at all b-values with mean CRLBs below 10%. The ADC values of PCr, ATP, Pi, and GPC were (0.24 ± 0.02, 0.15 ± 0.04, 0.43 ± 0.14, 0.40 ± 0.09) × 10^-3 mm²/s, respectively.

CONCLUSION

The ADCs of four ³¹P metabolites were successfully measured in the human calf muscle at 7 T, among which those of ATP, Pi and GPC were reported for the first time in humans. This study paves the way to investigate ³¹P metabolite diffusion properties in health and disease on the clinical MR scanner.

Macromolecular background signal and non-Gaussian metabolite diffusion determined in human brain using ultra-high diffusion weighting

Purpose

Definition of a macromolecular MR spectrum based on diffusion properties rather than relaxation time differences and characterization of non‐Gaussian diffusion of brain metabolites with strongly diffusion‐weighted MR spectroscopy.

Methods

Short echo time MRS with strong diffusion‐weighting with b‐values up to 25 ms/μm² at two diffusion times was implemented on a Connectom system and applied in combination with simultaneous spectral and diffusion decay modeling. Motion‐compensation was performed with a combined method based on the simultaneously acquired water and a macromolecular signal.

Results

The motion compensation scheme prevented spurious signal decay reflected in very small apparent diffusion constants for macromolecular signal. Macromolecular background signal patterns were determined using multiple fit strategies. Signal decay corresponding to non‐Gaussian metabolite diffusion was represented by biexponential fit models yielding parameter estimates for human gray matter that are in line with published rodent data. The optimal fit strategies used constraints for the signal decay of metabolites with limited signal contributions to the overall spectrum.

Conclusion

The determined macromolecular spectrum based on diffusion properties deviates from the conventional one derived from longitudinal relaxation time differences calling for further investigation before use as experimental basis spectrum when fitting clinical MR spectra. The biexponential characterization of metabolite signal decay is the basis for investigations into pathologic alterations of microstructure.

Click here for author‐reader discussions

White Matter Metabolite Relaxation and Diffusion Abnormalities in First-Episode Psychosis: A Longitudinal Study

Microstructural abnormalities in the white matter (WM) are implicated in the pathophysiology of psychosis. In vivo magnetic resonance spectroscopy (MRS) can probe the brain's intracellular microenvironment through the measurement of transverse relaxation and diffusion of neurometabolites and possibly provide cell-specific information. In our previous studies, we observed differential metabolite signal abnormalities in first episode and chronic stages of psychosis. In the present work, longitudinal data were presented for the first time on white matter cell-type specific abnormalities using a combination of diffusion tensor spectroscopy (DTS), T2 MRS, and diffusion tensor imaging (DTI) from a group of 25 first episode psychosis patients and nine matched controls scanned at baseline and one and two years of follow-up. We observed significantly reduced choline ADC in the year 1 of follow-up (0.194 µm2/ms) compared to baseline (0.229 µm2/ms), followed by a significant increase in NAA ADC in the year 2 follow-up (0.258 µm2/ms) from baseline (0.222 µm2/ms) and year 1 follow-up (0.217 µm2/ms). In contrast, NAA T2 relaxation, reflecting a related but different aspect of microenvironment from diffusion, was reduced at year 1 follow-up (257 ms) compared to baseline (278 ms). These abnormalities were observed in the absence of any abnormalities in water relaxation and diffusion at any timepoint. These findings indicate that abnormalities are seen in in glial-enriched (choline) signals in early stages of psychosis, followed by the subsequent emergence of neuronal-enriched (NAA) diffusion abnormalities, all in the absence of nonspecific water signal abnormalities.

In vivo 1 H MR spectroscopy with J-refocusing

PURPOSE

The goal of this study was to propose a novel localized proton MR spectroscopy (MRS) sequence that reduces signal loss due to J-modulation in the rat brain in vivo.

METHODS

Sprague-Dawley rats were studied at 9.4 T. A semi-LASER sequence with evenly distributed echo-time (T_E ) was used, and a 90° J-refocusing pulse was inserted at T_E /2. Proton spectra were acquired at two T_E s (30 and 68 ms), with and without the J-refocused pulse. Data were processed in MATLAB and quantified with LCModel.

RESULTS

The J-refocused spectrum acquired at T_E = 30 ms did not show any signal losses due to J-modulation and had comparable spectral pattern to the one acquired with semi-LASER using the minimum achievable T_E . Higher signal amplitudes for glutamine, γ-aminobutyric acid and glutathione led to more reliable quantification precision for these metabolites. The refocused signal intensities at T_E = 68 ms were also unaffected by J-modulation but were smaller than the signals at T_E = 30 ms mainly due to transverse T₂ relaxation of metabolites.

CONCLUSION

The proposed localized MRS sequence will be beneficial in both animal and human MRS studies when using ultra-short T_E is not possible while also providing more reliable quantification precision for J-coupled metabolites.

Mapping complex cell morphology in the grey matter with double diffusion encoding MR: A simulation study

This paper investigates the impact of cell body (namely soma) size and branching of cellular projections on diffusion MR imaging (dMRI) and spectroscopy (dMRS) signals for both standard single diffusion encoding (SDE) and more advanced double diffusion encoding (DDE) measurements using numerical simulations. The aim is to investigate the ability of dMRI/dMRS to characterize the complex morphology of brain cells focusing on these two distinctive features of brain grey matter. To this end, we employ a recently developed computational framework to create three dimensional meshes of neuron-like structures for Monte Carlo simulations, using diffusion coefficients typical of water and brain metabolites. Modelling the cellular structure as realistically connected spherical soma and cylindrical cellular projections, we cover a wide range of combinations of sphere radii and branching order of cellular projections, characteristic of various grey matter cells. We assess the impact of spherical soma size and branching order on the b-value dependence of the SDE signal as well as the time dependence of the mean diffusivity (MD) and mean kurtosis (MK). Moreover, we also assess the impact of spherical soma size and branching order on the angular modulation of DDE signal at different mixing times, together with the mixing time dependence of the apparent microscopic anisotropy (μA), a promising contrast derived from DDE measurements. The SDE results show that spherical soma size has a measurable impact on both the b-value dependence of the SDE signal and the MD and MK diffusion time dependence for both water and metabolites. On the other hand, we show that branching order has little impact on either, especially for water. In contrast, the DDE results show that spherical soma size has a measurable impact on the DDE signal's angular modulation at short mixing times and the branching order of cellular projections significantly impacts the mixing time dependence of the DDE signal's angular modulation as well as of the derived μA, for both water and metabolites. Our results confirm that SDE based techniques may be sensitive to spherical soma size, and most importantly, show for the first time that DDE measurements may be more sensitive to the dendritic tree complexity (as parametrized by the branching order of cellular projections), paving the way for new ways of characterizing grey matter morphology, non-invasively using dMRS and potentially dMRI.

Rapid acquisition diffusion MR spectroscopy of metabolites in human brain

Few studies have focused on metabolite diffusion in the human brain using ¹ H-MRS due to significant technical challenges. Moreover, such studies have required lengthy acquisition times and are therefore impractical to implement clinically. By first characterizing and then minimizing the effects of linear and oscillating eddy currents, which arise from the diffusion gradients, and by implementing phase-cycle and slice-order strategies, as well as introducing a new phase-alignment methodology, we report a method that allows data acquisition requiring 20 seconds per spectrum. This remained feasible, even for b-values >8000 s/mm² , with a rapid acquisition diffusion MRS methodology. It has allowed the nonlinear characterization of signal intensity with multiple b-values, and has improved the measurement of rotationally invariant diffusion parameters via six-direction, six b-value diffusion tensor spectroscopy (DTS) in 12 minutes at 4.7 T. The shorter DTS acquisition will enable its application to white matter regions not aligned with the gradients and permit clinical studies in a feasible time.

In vivo diffusion-weighted MRS using semi-LASER in the human brain at 3 T: Methodological aspects and clinical feasibility

Diffusion-weighted (DW-) MRS investigates non-invasively microstructural properties of tissue by probing metabolite diffusion in vivo. Despite the growing interest in DW-MRS for clinical applications, little has been published on the reproducibility of this technique. In this study, we explored the optimization of a single-voxel DW-semi-LASER sequence for clinical applications at 3 T, and evaluated the reproducibility of the method under different experimental conditions. DW-MRS measurements were carried out in 10 healthy participants and repeated across three sessions. Metabolite apparent diffusion coefficients (ADCs) were calculated from mono-exponential fits (ADC_exp ) up to b = 3300 s/mm² , and from the diffusional kurtosis approach (ADC_K ) up to b = 7300 s/mm² . The inter-subject variabilities of ADCs of N-acetylaspartate + N-acetylaspartylglutamate (tNAA), creatine + phosphocreatine, choline containing compounds, and myo-inositol were calculated in the posterior cingulate cortex (PCC) and in the corona radiata (CR). We explored the effect of physiological motion on the DW-MRS signal and the importance of cardiac gating and peak thresholding to account for signal amplitude fluctuations. Additionally, we investigated the dependence of the intra-subject variability on the acquisition scheme using a bootstrapping resampling method. Coefficients of variation were lower in PCC than CR, likely due to the different sensitivities to motion artifacts of the two regions. Finally, we computed coefficients of repeatability for ADC_exp and performed power calculations needed for designing clinical studies. The power calculation for ADC_exp of tNAA showed that in the PCC seven subjects per group are sufficient to detect a difference of 5% between two groups with an acquisition time of 4 min, suggesting that ADC_exp of tNAA is a suitable marker for disease-related intracellular alteration even in small case-control studies. In the CR, further work is needed to evaluate the voxel size and location that minimize the motion artifacts and variability of the ADC measurements.

Measuring transverse relaxation rates of the major brain metabolites from single-voxel PRESS acquisitions at a single TE

PURPOSE

To compare transverse relaxation rates of brain metabolites estimated from single-TE PRESS acquisitions with more conventionally derived rates estimated from multiple-TE PRESS acquisitions.

METHODS

Single-voxel (8 mL) PRESS data within white matter from 6 subjects were acquired at five different TEs. Transverse relaxation rates R₂ of N-acetylaspartate, creatine, and choline were estimated from a single TE using full versus right-side-only sampling of the echo. These R₂ values were compared with R_2Hahn values obtained from the multiple-TE PRESS acquisitions.

RESULTS

Following baseline subtraction and RMS weighting, interindividual mean R₂ values from TE = 288 ms magnitude spectra for choline, creatine, and N-acetylaspartate were highly correlated with respective R_2Hahn values (r² = 0.99). Paired individual measurements at this TE showed less correlation (r² = 0.48), primarily due to the N-acetylaspartate resonance. Using TE = 360 ms data for N-acetylaspartate and 288 ms for choline and creatine resulted in an improved correlation coefficient (r² = 0.80). The average absolute intra-individual differences in the estimated R₂ s between single-TE and Hahn method was 9.6 ± 7.7%.

CONCLUSION

For the major brain metabolite singlets, R_2Hahn values showed correlations with more fragile measurements of R₂ from a single TE that are worthy of interest. Because the left side of long-TE spin echoes is available "for free" from an acquisition perspective, and although the single-TE method for estimating R₂ values is associated with lower precision, the reduction in scan time may be clinically helpful.

Changes in the intracellular microenvironment in the aging human brain

Normal brain aging is associated with changes occurring at all levels. This study investigates age-related differences in the brain intracellular microenvironment by comparing the apparent diffusion coefficients (ADC) and apparent transverse relaxation time constants (T₂) of 5 neurochemicals (i.e., total N-acetyl-aspartate, total creatine, total choline, glutamate, and myo-inositol) between young and older adults. Thirty-two young healthy adults (18-22 years) and 26 older healthy adults (70-83 years) were recruited. Three brain regions were studied at 3 T: prefrontal, posterior cingulate and occipital cortices. ADC and T₂ were measured using stimulated echo acquisition mode and localization by adiabatic selective refocusing sequences, respectively. This study shows that the diffusivities of several neurochemicals are higher in older than in younger adults. In contrast, shorter apparent T₂ values for several metabolites were measured in older adults. Age-related difference in ADC and apparent T₂ of metabolites seem to be region-specific. Furthermore, this study shows that it is feasible to observe age-related differences in the cellular microenvironment of neurochemicals in the normal aging brain.