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

Diffusion-weighted MR spectroscopy of the prostate

PURPOSE

Prostate tissue has a complex microstructure, mainly composed of epithelial and stromal cells, and of extracellular (acinar-luminal) spaces. Diffusion-weighted MR spectroscopy (DW-MRS) is ideally suited to explore complex microstructure in vivo with metabolites selectively distributed in different subspaces. To date, this technique has been applied to brain and muscle. This study presents the development and pioneering utilization of 1H-DW-MRS in the prostate, accompanied by in vitro studies to support interpretations of in vivo findings.

METHODS

Nine healthy volunteers underwent a prostate MR examination (mean age, 56 years; range, 31-66). Metabolic complexation was studied in vitro using solutions with major compounds found in prostatic fluid of the lumen. DW-MRS was performed at 3 T with a non-water-suppressed single-voxel sequence with metabolite-cycling to concurrently measure metabolite and water signals. The water signal was used in postprocessing as a reference in a motion-compensation scheme. The spectra were fitted simultaneously in the spectral and diffusion-weighting dimensions. Apparent diffusion coefficients (ADCs) were derived by fitting signal decays that were assumed to be mono-exponential for metabolites and biexponential for water.

RESULTS

DW-MRS of the prostate revealed relatively low ADCs for Cho and Cr compounds, aligning with their intracellular location and higher ADCs for citrate and spermine supporting their luminal origin. In vitro assessments of the ADCs of citrate and spermine demonstrated their complex formation and protein binding. Tissue concentrations of MRS-detectable metabolites were as expected for the voxel location.

CONCLUSIONS

This work successfully demonstrates the feasibility of 1H-DW-MRS of the prostate and its potential for providing valuable microstructural information.

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.

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.

Diffusion-weighted SPECIAL improves the detection of J-coupled metabolites at ultrahigh magnetic field

PURPOSE

To improve the detection and subsequent estimation of the diffusion properties of strongly J-coupled metabolites in diffusion-weighted MRS (DWS).

METHODS

A new sequence for single-voxel diffusion-weighted 1 H MR spectroscopy, named DW-SPECIAL, is proposed. It combines the semi-adiabatic SPECIAL sequence with a stimulated echo diffusion block. Acquisitions with DW-SPECIAL and STE-LASER, the current gold standard for rodent DWS experiments at high fields, were performed at 14.1T on phantoms and in vivo on the rat brain. The apparent diffusion coefficient and intra-stick diffusivity (Callaghan's model, randomly-oriented sticks) were fitted and compared between the sequences for glutamate, glutamine, myo-inositol, taurine, total NAA, total Cho, total Cr, and the macromolecules.

RESULTS

The shorter TE achieved with DW-SPECIAL (18 ms against 33 ms with STE-LASER) substantially limited the metabolites' signal loss caused by J-evolution. In addition, DW-SPECIAL preserved the main advantages of STE-LASER: absence of cross-terms, diffusion time during a stimulated echo, and limited sensitivity to B1 inhomogeneities. In vivo, compared to STE-LASER, DW-SPECIAL yielded the same spectral quality and reduced the Cramer Rao Lower Bounds for J-coupled metabolites, irrespective of the b-value. DW-SPECIAL also reduced the SD of the metabolites' diffusion estimates based on individual animal fitting without loss of accuracy compared to the fit on the averaged decay.

CONCLUSION

We conclude that due to its reduced TE, DW-SPECIAL can serve as an alternative to STE-LASER when strongly J-coupled metabolites like glutamine are investigated, thereby extending the range of accessible metabolites in the context of DWS acquisitions.

Practical considerations of diffusion-weighted MRS with ultra-strong diffusion gradients

Introduction

Diffusion-weighted magnetic resonance spectroscopy (DW-MRS) offers improved cellular specificity to microstructure-compared to water-based methods alone-but spatial resolution and SNR is severely reduced and slow-diffusing metabolites necessitate higher b-values to accurately characterize their diffusion properties. Ultra-strong gradients allow access to higher b-values per-unit time, higher SNR for a given b-value, and shorter diffusion times, but introduce additional challenges such as eddy-current artefacts, gradient non-uniformity, and mechanical vibrations.

Methods

In this work, we present initial DW-MRS data acquired on a 3T Siemens Connectom scanner equipped with ultra-strong (300 mT/m) gradients. We explore the practical issues associated with this manner of acquisition, the steps that may be taken to mitigate their impact on the data, and the potential benefits of ultra-strong gradients for DW-MRS. An in-house DW-PRESS sequence and data processing pipeline were developed to mitigate the impact of these confounds. The interaction of TE, b-value, and maximum gradient amplitude was investigated using simulations and pilot data, whereby maximum gradient amplitude was restricted. Furthermore, two DW-MRS voxels in grey and white matter were acquired using ultra-strong gradients and high b-values.

Results

Simulations suggest T2-based SNR gains that are experimentally confirmed. Ultra-strong gradient acquisitions exhibit similar artefact profiles to those of lower gradient amplitude, suggesting adequate performance of artefact mitigation strategies. Gradient field non-uniformity influenced ADC estimates by up to 4% when left uncorrected. ADC and Kurtosis estimates for tNAA, tCho, and tCr align with previously published literature.

Discussion

In conclusion, we successfully implemented acquisition and data processing strategies for ultra-strong gradient DW-MRS and results indicate that confounding effects of the strong gradient system can be ameliorated, while achieving shorter diffusion times and improved metabolite SNR.

Quantification of short echo time MRS signals with improved version of QUantitation based on quantum ESTimation algorithm

Magnetic resonance spectroscopy offers information about metabolite changes in the organism, which can be used in diagnosis. While short echo time proton spectra exhibit more distinguishable metabolites compared with proton spectra acquired with long echo times, their quantification (and providing estimates of metabolite concentrations) is more challenging. They are hampered by a background signal, which originates mainly from macromolecules (MM) and mobile lipids. An improved version of the quantification algorithm QUantitation based on quantum ESTimation (QUEST), with MM prior knowledge (QUEST-MM), dedicated to proton signals and invoking appropriate prior knowledge on MM, is proposed and tested. From a single acquisition, it enables better metabolite quantification, automatic estimation of the background, and additional automatic quantification of MM components, thus improving its applicability in the clinic. The proposed algorithm may facilitate studies that involve patients with pathological MM in the brain. QUEST-MM and three QUEST-based strategies for quantifying short echo time signals are compared in terms of bias-variance trade-off and Cramér-Rao lower bound estimates. The performances of the methods are evaluated through extensive Monte Carlo studies. In particular, the histograms of the metabolite and MM amplitude distributions demonstrate the performances of the estimators. They showed that QUEST-MM works better than QUEST (Subtract approach) and is a good alternative to QUEST when measured MM signal is unavailable or unsuitable. Quantification with QUEST-MM is shown for 1 H in vivo rat brain signals obtained with the SPECIAL pulse sequence at 9.4 T, and human brain signals obtained, respectively, with STEAM at 4 T and PRESS at 3 T. QUEST-MM is implemented in jMRUI and will be available for public use from version 7.1.

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.