Water removal in MR spectroscopic imaging with L2 regularization

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

PURPOSE

The J-difference edited γ-aminobutyric acid (GABA) signal is contaminated by other co-edited signals-the largest of which originates from co-edited macromolecules (MMs)-and is consequently often reported as "GABA+." MM signals are broader and less well-characterized than the metabolites, and are commonly approximated using a Gaussian model parameterization. Experimentally measured MM signals are a consensus-recommended alternative to parameterized modeling; however, they are relatively under-studied in the context of edited MRS.

METHODS

To address this limitation in the literature, we have acquired GABA-edited MEGA-PRESS data with pre-inversion to null metabolite signals in 13 healthy controls. An experimental MM basis function was derived from the mean across subjects. We further derived a new parameterization of the MM signals from the experimental data, using multiple Gaussians to accurately represent their observed asymmetry. The previous single-Gaussian parameterization, mean experimental MM spectrum and new multi-Gaussian parameterization were compared in a three-way analysis of a public MEGA-PRESS dataset of 61 healthy participants.

RESULTS

Both the experimental MMs and the multi-Gaussian parameterization exhibited reduced fit residuals compared to the single-Gaussian approach (p = 0.034 and p = 0.031, respectively), suggesting they better represent the underlying data than the single-Gaussian parameterization. Furthermore, both experimentally derived models estimated larger MM fractional contribution to the GABA+ signal for the experimental MMs (58%) and multi-Gaussian parameterization (58%), compared to the single-Gaussian approach (50%).

CONCLUSIONS

Our results indicate that single-Gaussian parameterization of edited MM signals is insufficient and that both experimentally derived GABA+ spectra and their parameterized replicas improve the modeling of GABA+ spectra.

Water removal in MR spectroscopic imaging with Casorati singular value decomposition

PURPOSE

Water removal is one of the computational bottlenecks in the processing of high-resolution MRSI data. The purpose of this work is to propose an approach to reduce the computing time required for water removal in large MRS data.

METHODS

In this work, we describe a singular value decomposition-based approach that uses the partial position-time separability and the time-domain linear predictability of MRSI data to reduce the computational time required for water removal. Our approach arranges MRS signals in a Casorati matrix form, applies low-rank approximations utilizing singular value decomposition, removes residual water from the most prominent left-singular vectors, and finally reconstructs the water-free matrix using the processed left-singular vectors.

RESULTS

We have demonstrated the effectiveness of our proposed algorithm for water removal using both simulated and in vivo data. The proposed algorithm encompasses a pip-installable tool ( https://pypi.org/project/CSVD/), available on GitHub ( https://github.com/amirshamaei/CSVD), empowering researchers to use it in future studies. Additionally, to further promote transparency and reproducibility, we provide comprehensive code for result replication.

CONCLUSIONS

The findings of this study suggest that the proposed method is a promising alternative to existing water removal methods due to its low processing time and good performance in removing water signals.

Application of a 1H Brain MRS Benchmark Dataset to Deep Learning for Out-of-Voxel Artifacts

Neural networks are potentially valuable for many of the challenges associated with MRS data. The purpose of this manuscript is to describe the AGNOSTIC dataset, which contains 259,200 synthetic 1H MRS examples for training and testing neural networks. AGNOSTIC was created using 270 basis sets that were simulated across 18 field strengths and 15 echo times. The synthetic examples were produced to resemble in vivo brain data with combinations of metabolite, macromolecule, residual water signals, and noise. To demonstrate the utility, we apply AGNOSTIC to train two Convolutional Neural Networks (CNNs) to address out-of-voxel (OOV) echoes. A Detection Network was trained to identify the point-wise presence of OOV echoes, providing proof of concept for real-time detection. A Prediction Network was trained to reconstruct OOV echoes, allowing subtraction during post-processing. Complex OOV signals were mixed into 85% of synthetic examples to train two separate CNNs for the detection and prediction of OOV signals. AGNOSTIC is available through Dryad and all Python 3 code is available through GitHub. The Detection network was shown to perform well, identifying 95% of OOV echoes. Traditional modeling of these detected OOV signals was evaluated and may prove to be an effective method during linear-combination modeling. The Prediction Network greatly reduces OOV echoes within FIDs and achieved a median log10 normed-MSE of -1.79, an improvement of almost two orders of magnitude.

Test-retest reproducibility of human brain multi-slice 1 H FID-MRSI data at 9.4T after optimization of lipid regularization, macromolecular model, and spline baseline stiffness

PURPOSE

This study analyzes the effects of retrospective lipid suppression, a simulated macromolecular prior knowledge and different spline baseline stiffness values on 9.4T multi-slice proton FID-MRSI data spanning the whole cerebrum of human brain and the reproducibility of respective metabolite ratio to total creatine (/tCr) maps for 10 brain metabolites.

METHODS

Measurements were performed twice on 5 volunteers using a short TR and TE FID MRSI 2D sequence at 9.4T. The effects of retrospective lipid L2-regularization, macromolecular spectrum and different LCModel baseline flexibilities on SNR, FWHM, fitting residual, Cramér-Rao lower bound, and metabolite ratio maps were investigated. Intra-subject, inter-session coefficient of variation and the test-retest reproducibility of the mean metabolite ratios (/tCr) of each slice was calculated.

RESULTS

Transversal, sagittal, and coronal slices of many metabolite ratio maps correspond to the anatomically expected concentration relations in gray and white matter for the majority of the cerebrum when using a flexible baseline in LCModel fit. Results from the second measurements of the same subjects show that slice positioning and data quality correlate significantly to the first measurement. L2-regularization provided effective suppression of lipid-artifacts, but should be avoided if no artifacts are detected.

CONCLUSION

Reproducible concentration ratio maps (/tCr) for 4 metabolites (total choline, N-acetylaspartate, glutamate, and myoinositol) spanning the majority of the cerebrum and 6 metabolites (N-acetylaspartylglutamate, γ-aminobutyric acid, glutathione, taurine, glutamine, and aspartate) covering 32 mm in the upper part of the brain were acquired at 9.4T using multi-slice FID MRSI with retrospective lipid suppression, a macromolecular spectrum and a flexible LCModel baseline.

A novel method for Magnetic Resonance Spectroscopy lipid signal suppression using Semi-classical signal analysis and Bidirectional Long short-term memory

Magnetic resonance spectroscopy (MRS) is a non-invasive method that enables the analysis and quantification of brain metabolites, which provide useful information about the neuro-biological substrates of brain function. Lactate plays a pivotal role in the diagnosis of various brain diseases. However, accurate lactate quantification is generally difficult to achieve due to the presence of large lipid peaks resonating at a similar spectral position. To overcome this problem several techniques have been proposed. However, most of them suffer from lactate signal loss or poor lipid peak removal. In this paper, a novel method for lipid suppression for MRS signal is proposed. The method combines a semi-classical signal analysis method and a bidirectional long short term memory technique. The method is validated using simulated data that mimics real MRS data.

Emerging methods and applications of ultra-high field MR spectroscopic imaging in the human brain

Magnetic Resonance Spectroscopic Imaging (MRSI) of the brain enables insights into the metabolic changes and fluxes in diseases such as tumors, multiple sclerosis, epilepsy, or hepatic encephalopathy, as well as insights into general brain functionality. However, the routine application of MRSI is mostly hampered by very low signal-to-noise ratios (SNR) due to the low concentrations of metabolites, about 10000 times lower than water. Furthermore, MRSI spectra have a dense information content with many overlapping metabolite resonances, especially for proton MRSI. MRI scanners at ultra-high field strengths, like 7 T or above, offer the opportunity to increase SNR, as well as the separation between resonances, thus promising to solve both challenges. Yet, MRSI at ultra-high field strengths is challenged by decreased B0- and B1-homogeneity, shorter T2 relaxation times, stronger chemical shift displacement errors, and aggravated lipid contamination. Therefore, to capitalize on the advantages of ultra-high field strengths, these challenges must be overcome. This review focuses on the challenges MRSI of the human brain faces at ultra-high field strength, as well as the possible applications to this date.

Magnetic resonance spectroscopic imaging of downfield proton resonances in the human brain at 3 T

PURPOSE

To develop an MRSI technique capable of mapping downfield proton resonances in the human brain.

METHODS

A spectral-spatial excitation and frequency-selective refocusing scheme, in combination with 2D phase encoding, was developed for mapping of downfield resonances without any perturbation of the water magnetization. An alternative scheme using spectral-spatial refocusing was also investigated for simultaneous detection of both downfield and upfield resonances. The method was tested in 5 healthy human volunteers.

RESULTS

Downfield metabolite maps with a nominal spatial resolution of 1.5 cm³ were recorded at 3 T in a scan time of 12 minutes. Cramer-Rao lower bounds for nine different downfield peaks were 20% or less over a single supraventricular slice. Downfield spectral profiles were similar to those in the literature recorded previously using single-voxel localization methods. The same approach was also used for upfield MRSI, and simultaneous upfield and downfield acquisitions.

CONCLUSION

The developed MRSI pulse sequence was shown to be an efficient way of rapidly mapping downfield resonances in the human brain at 3 T, maximizing sensitivity through the relaxation enhancement effect. Because the MRSI approach is efficient in terms of data collection and can be readily implemented at short TE, somewhat higher spatial resolution can be achieved than has been reported in previous single-voxel downfield MRS studies. With this approach, nine downfield resonances could be mapped in a single slice for the first time using MRSI at 3 T.

Accelerated J-resolved 1 H-MRSI with limited and sparse sampling of ( k , t 1 , t 2 -space

PURPOSE

To accelerate the acquisition of J-resolved proton magnetic resonance spectroscopic imaging (¹ H-MRSI) data for high-resolution mapping of brain metabolites and neurotransmitters.

METHODS

The proposed method used a subspace model to represent multidimensional spatiospectral functions, which significantly reduced the number of parameters to be determined from J-resolved ¹ H-MRSI data. A semi-LASER-based (Localization by Adiabatic SElective Refocusing) echo-planar spectroscopic imaging (EPSI) sequence was used for data acquisition. The proposed data acquisition scheme sampled k , t 1 , t 2 -space in variable density, where t₁ and t₂ specify the J-coupling and chemical-shift encoding times, respectively. Selection of the J-coupling encoding times (or, echo time values) was based on a Cramer-Rao lower bound analysis, which were optimized for gamma-aminobutyric acid (GABA) detection. In image reconstruction, parameters of the subspace-based spatiospectral model were determined by solving a constrained optimization problem.

RESULTS

Feasibility of the proposed method was evaluated using both simulated and experimental data from a spectroscopic phantom. The phantom experimental results showed that the proposed method, with a factor of 12 acceleration in data acquisition, could determine the distribution of J-coupled molecules with expected accuracy. In vivo study with healthy human subjects also showed that 3D maps of brain metabolites and neurotransmitters can be obtained with a nominal spatial resolution of 3.0 × 3.0 × 4.8 mm³ from J-resolved ¹ H-MRSI data acquired in 19.4 min.

CONCLUSIONS

This work demonstrated the feasibility of highly accelerated J-resolved ¹ H-MRSI using limited and sparse sampling of k , t 1 , t 2 -space and subspace modeling. With further development, the proposed method may enable high-resolution mapping of brain metabolites and neurotransmitters in clinical applications.