Multiple-echo proton spectroscopic imaging using time domain parametric spectral analysis

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.

Multi-echo single-shot spectroscopy combined with simultaneous 2D model fitting for fast and accurate measurement of metabolite-specific concentrations and T2 relaxation times

The purpose of the current study was to develop a novel single-voxel MR spectroscopy acquisition scheme to simultaneously determine metabolite-specific concentrations and transverse relaxation times within realistic clinical scan times. Partly truncated multi-TE data are acquired as an echo train in a single acquisition (multi-echo single-shot [MESS]). A 2D multiparametric model fitting approach combines truncated, low-resolved short TE data with fully sampled, highly resolved, longer TE data to yield concentration and T2 estimates for major brain metabolites simultaneously. Cramer-Rao lower bounds (CRLB) are used as a measure of performance. The novel scheme was compared with traditional multi-echo multi-shot methods. In silico, in vitro, and in vivo experiments support the findings. MESS schemes, requiring only 2 min 12 s for the acquisition of three echo times, provide valid concentration and relaxation estimates for multiple metabolites and outperform traditional methods for simultaneous determinations of metabolite-specific T2 s and concentrations, with improvements ranging from 5% to 30% for T2 s and from 10% to 50% for concentrations. However, substantial unsuppressed residual water signals may hamper the method's reproducibility, as observed in an initial experiment setup that prioritizes short TEs with severely truncated acquisition for the benefit of signal-to-noise ratio (SNR). Nevertheless, CRLB have been confirmed to be well suited as design criteria, and within-session repeatability approaches CRLB when residual water is removed in postprocessing by exploiting longer and less truncated data recordings. MESS MRS combined with 2D model fitting promises comparable accuracy, increased precision, or inversely shorter experimental times compared with traditional approaches. However, the optimal design must be investigated as a trade-off between SNR, the truncation factor, and TE batch selections, all of which influence the robustness of estimations.

Rapid multi-echo measurement of brain metabolite T₂ values at 7 T using a single-shot spectroscopic Carr-Purcell-Meiboom-Gill sequence and prior information

We present a method for the robust and accurate estimation of brain metabolite transverse relaxation times (T2 ) from multiple spin-echo data acquired with a single-shot Carr-Purcell-Meiboom-Gill (CPMG) spectroscopic sequence. Each acquired echo consists of a small number of complex time-domain data points. The amplitudes of the spectral components in each echo are calculated by solving a set of linear equations in which previously estimated frequencies and linewidths serve as prior information. These priors are obtained from a short MRS experiment in which a large number of time-domain data points are acquired, and are subsequently estimated using linear prediction with singular value decomposition (LPSVD) processing. We show that this process can be used to accurately and rapidly measure the T2 values for the main singlet resonances in single-volume MRS measurements in the brain. The proposed method can be generalized to any set of MRS experiments comprising repeated measurements of amplitude changes, e.g. as a function of an experimental parameter, such as TE, inversion time or diffusion weighting.

Integrated data acquisition and processing to determine metabolite contents, relaxation times, and macromolecule baseline in single examinations of individual subjects

Absolute quantitation of clinical (1)H-MR spectra is virtually always incomplete for single subjects because the separate determination of spectrum, baseline, and transverse and longitudinal relaxation times in single subjects is prohibitively long. Integrated Processing and Acquisition of Data (IPAD) based on a combined 2-dimensional experimental and fitting strategy is suggested to substantially improve the information content from a given measurement time. A series of localized saturation-recovery spectra was recorded and combined with 2-dimensional prior-knowledge fitting to simultaneously determine metabolite T(1) (from analysis of the saturation-recovery time course), metabolite T(2) (from lineshape analysis based on metabolite and water peak shapes), macromolecular baseline (based on T(1) differences and analysis of the saturation-recovery time course), and metabolite concentrations (using prior knowledge fitting and conventional procedures of absolute standardization). The procedure was tested on metabolite solutions and applied in 25 subjects (15-78 years old). Metabolite content was comparable to previously found values. Interindividual variation was larger than intraindividual variation in repeated spectra for metabolite content as well as for some relaxation times. Relaxation times were different for various metabolite groups. Parts of the interindividual variation could be explained by significant age dependence of relaxation times.

[Fast 1H-NR spectroscopic imaging of the human brain]

Techniques allowing the reduction of the measurement time of MR spectroscopic imaging of protons (1H SI) were studied in a 1.5-T whole-body MR tomograph. The purpose was to facilitate the inclusion of the modality into an MRT routine protocol for examination of the human brain. The techniques were evaluated in a study with 6 healthy volunteers. Shortening of the measurement times can be achieved by reducing the number of phase-encoding steps with spherical k-space sampling, or by using delay times during TR periods for excitation and detection of several slices or echos. Since the tomograph enables echo-planar imaging, multi-gradiennt-echo techniques were implemented, thus sparing one phase-encoding gradient and shortening the measurement time by a factor of N. Due to signal loss, however, this technique is applicable only when the concentration of the metabolites of interest is sufficiently high to obtain a signal-to-noise ratio of at least 2 with a single acquisition.