Ultrafast magnetic resonance spectroscopic imaging using SPICE with learned subspaces

Learning-Assisted Fast Determination of Regularization Parameter in Constrained Image Reconstruction

OBJECTIVE

To leverage machine learning (ML) for fast selection of optimal regularization parameter in constrained image reconstruction.

METHODS

Constrained image reconstruction is often formulated as a regularization problem and selecting a good regularization parameter value is an essential step. We solved this problem using an ML-based approach by leveraging the finding that for a specific constrained reconstruction problem defined for a fixed class of image functions, the optimal regularization parameter value is weakly subject-dependent and the dependence can be captured using few experimental data. The proposed method has four key steps: a) solution of a given constrained reconstruction problem for a few (say, 3) pre-selected regularization parameter values, b) extraction of multiple approximated quality metrics from the initial reconstructions, c) predicting the true quality metrics values from the approximated values using pre-trained neural networks, and d) determination of the optimal regularization parameter by fusing the predicted quality metrics.

RESULTS

The effectiveness of the proposed method was demonstrated in two constrained reconstruction problems. Compared with L-curve-based method, the proposed method determined the regularization parameters much faster and produced substantially improved reconstructions. Our method also outperformed state-of-the-art learning-based methods when trained with limited experimental data.

CONCLUSION

This paper demonstrates the feasibility and improved reconstruction quality by using machine learning to determine the regularization parameter in constrained reconstruction.

SIGNIFICANCE

The proposed method substantially reduces the computational burden of the traditional methods (e.g., L-curve) or relaxes the requirement of large training data by modern learning-based methods, thus enhancing the practical utility of constrained reconstruction.

Mapping intracellular NAD content in entire human brain using phosphorus-31 MR spectroscopic imaging at 7 Tesla

Introduction

Nicotinamide adenine dinucleotide (NAD) is a crucial molecule in cellular metabolism and signaling. Mapping intracellular NAD content of human brain has long been of interest. However, the sub-millimolar level of cerebral NAD concentration poses significant challenges for in vivo measurement and imaging.

Methods

In this study, we demonstrated the feasibility of non-invasively mapping NAD contents in entire human brain by employing a phosphorus-31 magnetic resonance spectroscopic imaging (31P-MRSI)-based NAD assay at ultrahigh field (7 Tesla), in combination with a probabilistic subspace-based processing method.

Results

The processing method achieved about a 10-fold reduction in noise over raw measurements, resulting in remarkably reduced estimation errors of NAD. Quantified NAD levels, observed at approximately 0.4 mM, exhibited good reproducibility within repeated scans on the same subject and good consistency across subjects in group data (2.3 cc nominal resolution). One set of higher-resolution data (1.0 cc nominal resolution) unveiled potential for assessing tissue metabolic heterogeneity, showing similar NAD distributions in white and gray matter. Preliminary analysis of age dependence suggested that the NAD level decreases with age.

Discussion

These results illustrate favorable outcomes of our first attempt to use ultrahigh field 31P-MRSI and advanced processing techniques to generate a whole-brain map of low-concentration intracellular NAD content in the human brain.

High-Dimensional MR Reconstruction Integrating Subspace and Adaptive Generative Models

OBJECTIVE

To develop a new method that integrates subspace and generative image models for high-dimensional MR image reconstruction.

METHODS

We proposed a formulation that synergizes a low-dimensional subspace model of high-dimensional images, an adaptive generative image prior serving as spatial constraints on the sequence of "contrast-weighted" images or spatial coefficients of the subspace model, and a conventional sparsity regularization. A special pretraining plus subject-specific network adaptation strategy was proposed to construct an accurate generative-network-based representation for images with varying contrasts. An iterative algorithm was introduced to jointly update the subspace coefficients and the multi-resolution latent space of the generative image model that leveraged an recently proposed intermediate layer optimization technique for network inversion.

RESULTS

We evaluated the utility of the proposed method for two high-dimensional imaging applications: accelerated MR parameter mapping and high-resolution MR spectroscopic imaging. Improved performance over state-of-the-art subspace-based methods was demonstrated in both cases.

CONCLUSION

The proposed method provided a new way to address high-dimensional MR image reconstruction problems by incorporating an adaptive generative model as a data-driven spatial prior for constraining subspace reconstruction.

SIGNIFICANCE

Our work demonstrated the potential of integrating data-driven and adaptive generative priors with canonical low-dimensional modeling for high-dimensional imaging problems.

An Update on MR Spectroscopy in Cancer Management: Advances in Instrumentation, Acquisition, and Analysis

MR spectroscopy (MRS) is a noninvasive imaging method enabling chemical and molecular profiling of tissues in a localized, multiplexed, and nonionizing manner. As metabolic reprogramming is a hallmark of cancer, MRS provides valuable metabolic and molecular information for cancer diagnosis, prognosis, treatment monitoring, and patient management. This review provides an update on the use of MRS for clinical cancer management. The first section includes an overview of the principles of MRS, current methods, and conventional metabolites of interest. The remainder of the review is focused on three key areas: advances in instrumentation, specifically ultrahigh-field-strength MRI scanners and hybrid systems; emerging methods for acquisition, including deuterium imaging, hyperpolarized carbon 13 MRI and MRS, chemical exchange saturation transfer, diffusion-weighted MRS, MR fingerprinting, and fast acquisition; and analysis aided by artificial intelligence. The review concludes with future recommendations to facilitate routine use of MRS in cancer management. Keywords: MR Spectroscopy, Spectroscopic Imaging, Molecular Imaging in Oncology, Metabolic Reprogramming, Clinical Cancer Management © RSNA, 2024.

Boosting quantification accuracy of chemical exchange saturation transfer MRI with a spatial-spectral redundancy-based denoising method

Chemical exchange saturation transfer (CEST) is a versatile technique that enables noninvasive detections of endogenous metabolites present in low concentrations in living tissue. However, CEST imaging suffers from an inherently low signal-to-noise ratio (SNR) due to the decreased water signal caused by the transfer of saturated spins. This limitation challenges the accuracy and reliability of quantification in CEST imaging. In this study, a novel spatial-spectral denoising method, called BOOST (suBspace denoising with nOnlocal lOw-rank constraint and Spectral local-smooThness regularization), was proposed to enhance the SNR of CEST images and boost quantification accuracy. More precisely, our method initially decomposes the noisy CEST images into a low-dimensional subspace by leveraging the global spectral low-rank prior. Subsequently, a spatial nonlocal self-similarity prior is applied to the subspace-based images. Simultaneously, the spectral local-smoothness property of Z-spectra is incorporated by imposing a weighted spectral total variation constraint. The efficiency and robustness of BOOST were validated in various scenarios, including numerical simulations and preclinical and clinical conditions, spanning magnetic field strengths from 3.0 to 11.7 T. The results demonstrated that BOOST outperforms state-of-the-art algorithms in terms of noise elimination. As a cost-effective and widely available post-processing method, BOOST can be easily integrated into existing CEST protocols, consequently promoting accuracy and reliability in detecting subtle CEST effects.

Subspace Model-Assisted Deep Learning for Improved Image Reconstruction

Image reconstruction from limited and/or sparse data is known to be an ill-posed problem and a priori information/constraints have played an important role in solving the problem. Early constrained image reconstruction methods utilize image priors based on general image properties such as sparsity, low-rank structures, spatial support bound, etc. Recent deep learning-based reconstruction methods promise to produce even higher quality reconstructions by utilizing more specific image priors learned from training data. However, learning high-dimensional image priors requires huge amounts of training data that are currently not available in medical imaging applications. As a result, deep learning-based reconstructions often suffer from two known practical issues: a) sensitivity to data perturbations (e.g., changes in data sampling scheme), and b) limited generalization capability (e.g., biased reconstruction of lesions). This paper proposes a new method to address these issues. The proposed method synergistically integrates model-based and data-driven learning in three key components. The first component uses the linear vector space framework to capture global dependence of image features; the second exploits a deep network to learn the mapping from a linear vector space to a nonlinear manifold; the third is an unrolling-based deep network that captures local residual features with the aid of a sparsity model. The proposed method has been evaluated with magnetic resonance imaging data, demonstrating improved reconstruction in the presence of data perturbation and/or novel image features. The method may enhance the practical utility of deep learning-based image reconstruction.

Denoising magnetic resonance spectroscopy (MRS) data using stacked autoencoder for improving signal-to-noise ratio and speed of MRS

BACKGROUND

While magnetic resonance imaging (MRI) provides high resolution anatomical images with sharp soft tissue contrast, magnetic resonance spectroscopy (MRS) enables non-invasive detection and measurement of biochemicals and metabolites. However, MRS has low signal-to-noise ratio (SNR) when concentrations of metabolites are in the range of millimolar. Standard approach of using a high number of signal averaging (NSA) to achieve sufficient SNR comes at the cost of a long acquisition time.

PURPOSE

We propose to use deep-learning approaches to denoise MRS data without increasing NSA. This method has potential to reduce the acquisition time as well as improve SNR and quality of spectra, which could enhance the diagnostic value and broaden the clinical applications of MRS.

METHODS

The study was conducted using data collected from the brain spectroscopy phantom and human subjects. We utilized a stack auto-encoder (SAE) network to train deep learning models for denoising low NSA data (NSA = 1, 2, 4, 8, and 16) randomly truncated from high SNR data collected with high NSA (NSA = 192), which were also used to obtain the ground truth. We applied both self-supervised and fully-supervised training approaches and compared their performance of denoising low NSA data based on improvement in SNR. To prevent overfitting, the SAE network was trained in a patch-based manner. We then tested the denoising methods on noise-containing data collected from the phantom and human subjects, including data from brain tumor patients. We evaluated their performance by comparing the SNR levels and mean squared errors (MSEs) calculated for the whole spectra against high SNR "ground truth", as well as the value of chemical shift of N-acetyl-aspartate (NAA) before and after denoising.

RESULTS

With the SAE model, the SNR of low NSA data (NSA = 1) obtained from the phantom increased by 28.5% and the MSE decreased by 42.9%. For low NSA data of the human parietal and temporal lobes, the SNR increased by 32.9% and the MSE decreased by 63.1%. In all cases, the chemical shift of NAA in the denoised spectra closely matched with the high SNR spectra without significant distortion to the spectra after denoising. Furthermore, the denoising performance of the SAE model was more effective in denoising spectra with higher noise levels.

CONCLUSIONS

The reported SAE denoising method is a model-free approach to enhance the SNR of MRS data collected with low NSA. With the denoising capability, it is possible to acquire MRS data with a few NSA, shortening the scan time while maintaining adequate spectroscopic information for detecting and quantifying the metabolites of interest. This approach has the potential to improve the efficiency and effectiveness of clinical MRS data acquisition by reducing the scan time and increasing the quality of spectroscopic data.

Prediction of Stroke Onset Time With Combined Fast High-Resolution Magnetic Resonance Spectroscopic and Quantitative T2 Mapping

OBJECTIVE

The purpose of this work is to develop a multispectral imaging approach that combines fast high-resolution 3D magnetic resonance spectroscopic imaging (MRSI) and fast quantitative T2 mapping to capture the multifactorial biochemical changes within stroke lesions and evaluate its potentials for stroke onset time prediction.

METHODS

Special imaging sequences combining fast trajectories and sparse sampling were used to obtain whole-brain maps of both neurometabolites (2.0 × 3.0 × 3.0 mm3) and quantitative T2 values (1.9 × 1.9 × 3.0 mm3) within a 9-minute scan. Participants with ischemic stroke at hyperacute (0-24 h, n = 23) or acute (24 h-7d, n = 33) phase were recruited in this study. Lesion N-acetylaspartate (NAA), lactate, choline, creatine, and T2 signals were compared between groups and correlated with patient symptomatic duration. Bayesian regression analyses were employed to compare the predictive models of symptomatic duration using multispectral signals.

RESULTS

In both groups, increased T2 and lactate levels, as well as decreased NAA and choline levels were detected within the lesion (all p < 0.001). Changes in T2, NAA, choline, and creatine signals were correlated with symptomatic duration for all patients (all p < 0.005). Predictive models of stroke onset time combining signals from MRSI and T2 mapping achieved the best performance (hyperacute: R2 = 0.438; all: R2 = 0.548).

CONCLUSION

The proposed multispectral imaging approach provides a combination of biomarkers that index early pathological changes after stroke in a clinical-feasible time and improves the assessment of the duration of cerebral infarction.

SIGNIFICANCE

Developing accurate and efficient neuroimaging techniques to provide sensitive biomarkers for prediction of stroke onset time is of great importance for maximizing the proportion of patients eligible for therapeutic intervention. The proposed method provides a clinically feasible tool for the assessment of symptom onset time post ischemic stroke, which will help guide time-sensitive clinical management.

B1 mapping using pre-learned subspaces for quantitative brain imaging

PURPOSE

To develop a machine learning-based method for estimation of both transmitter and receiver B1 fields desired for correction of the B1 inhomogeneity effects in quantitative brain imaging.

THEORY AND METHODS

A subspace model-based machine learning method was proposed for estimation of B1t and B1r fields. Probabilistic subspace models were used to capture scan-dependent variations in the B1 fields; the subspace basis and coefficient distributions were learned from pre-scanned training data. Estimation of the B1 fields for new experimental data was achieved by solving a linear optimization problem with prior distribution constraints. We evaluated the performance of the proposed method for B1 inhomogeneity correction in quantitative brain imaging scenarios, including T1 and proton density (PD) mapping from variable-flip-angle spoiled gradient-echo (SPGR) data as well as neurometabolic mapping from MRSI data, using phantom, healthy subject and brain tumor patient data.

RESULTS

In both phantom and healthy subject data, the proposed method produced high-quality B1 maps. B1 correction on SPGR data using the estimated B1 maps produced significantly improved T1 and PD maps. In brain tumor patients, the proposed method produced more accurate and robust B1 estimation and correction results than conventional methods. The B1 maps were also applied to MRSI data from tumor patients and produced improved neurometabolite maps, with better separation between pathological and normal tissues.

CONCLUSION

This work presents a novel method to estimate B1 variations using probabilistic subspace models and machine learning. The proposed method may make correction of B1 inhomogeneity effects more robust in practical applications.

Sodium Hydrosulfide Treatment During Porcine Kidney Ex Vivo Perfusion and Transplantation

Background

In rodents, hydrogen sulfide (H2S) reduces ischemia-reperfusion injury and improves renal graft function after transplantation. Here, we hypothesized that the benefits of H2S are conserved in pigs, a more clinically relevant model.

Methods

Adult porcine kidneys retrieved immediately or after 60 min of warm ischemia (WI) were exposed to 100 µM sodium hydrosulfide (NaHS) (1) during the hypothermic ex vivo perfusion only, (2) during WI only, and (3) during both WI and ex vivo perfusion. Kidney perfusion was evaluated with dynamic contrast-enhanced MRI. MRI spectroscopy was further employed to assess energy metabolites including ATP. Renal biopsies were collected at various time points for histopathological analysis.

Results

Perfusion for 4 h pig kidneys with Belzer MPS UW + NaHS resulted in similar renal perfusion and ATP levels than perfusion with UW alone. Similarly, no difference was observed when NaHS was administered in the renal artery before ischemia. After autotransplantation, no improvement in histologic lesions or cortical/medullary kidney perfusion was observed upon H2S administration. In addition, AMP and ATP levels were identical in both groups.

Conclusions

In conclusion, treatment of porcine kidney grafts using NaHS did not result in a significant reduction of ischemia-reperfusion injury or improvement of kidney metabolism. Future studies will need to define the benefits of H2S in human, possibly using other molecules as H2S donors.

High-Dimensional MR Spatiospectral Imaging by Integrating Physics-Based Modeling and Data-Driven Machine Learning: Current progress and future directions

Magnetic resonance spectroscopic imaging (MRSI) offers a unique molecular window into the physiological and pathological processes in the human body. However, the applications of MRSI have been limited by a number of long-standing technical challenges due to high dimensionality and low signal-to-noise ratio (SNR). Recent technological developments integrating physics-based modeling and data-driven machine learning that exploit unique physical and mathematical properties of MRSI signals have demonstrated impressive performance in addressing these challenges for rapid, high-resolution, quantitative MRSI. This paper provides a systematic review of these progresses in the context of MRSI physics and offers perspectives on promising future directions.

Subnormothermic Ex Vivo Porcine Kidney Perfusion Improves Energy Metabolism: Analysis Using 31P Magnetic Resonance Spectroscopic Imaging

The ideal preservation temperature for donation after circulatory death kidney grafts is unknown. We investigated whether subnormothermic (22 °C) ex vivo kidney machine perfusion could improve kidney metabolism and reduce ischemia-reperfusion injury.

T 2 ' mapping of the brain from water-unsuppressed 1 H-MRSI and turbo spin-echo data

PURPOSE

To obtain high-quality T 2 ' $$ {\mathrm{T}}_2^{\prime } $$ maps of brain tissues from water-unsuppressed magnetic resonance spectroscopic imaging (MRSI) and turbo spin-echo (TSE) data.

METHODS

T 2 ' $$ {\mathrm{T}}_2^{\prime } $$ mapping can be achieved using T 2 * $$ {\mathrm{T}}_2^{\ast } $$ mapping from water-unsuppressed MRSI data and T 2 $$ {\mathrm{T}}_2 $$ mapping from TSE data. However, T 2 * $$ {\mathrm{T}}_2^{\ast } $$ mapping often suffers from signal dephasing and distortions caused by B 0 $$ {\mathrm{B}}_0 $$ field inhomogeneity; T 2 $$ {\mathrm{T}}_2 $$ measurements may be biased due to system imperfections, especially for T 2 $$ {\mathrm{T}}_2 $$ -weighted image with small number of TEs. In this work, we corrected the B 0 $$ {\mathrm{B}}_0 $$ field inhomogeneity effect on T 2 * $$ {\mathrm{T}}_2^{\ast } $$ mapping using a subspace model-based method, incorporating pre-learned spectral basis functions of the water signals. T 2 $$ {\mathrm{T}}_2 $$ estimation bias was corrected using a TE-adjustment method, which modeled the deviation between measured and reference T 2 $$ {\mathrm{T}}_2 $$ decays as TE shifts.

RESULTS

In vivo experiments were performed to evaluate the performance of the proposed method. High-quality T 2 * $$ {\mathrm{T}}_2^{\ast } $$ maps were obtained in the presence of large field inhomogeneity in the prefrontal cortex. Bias in T 2 $$ {\mathrm{T}}_2 $$ measurements obtained from TSE data was effectively reduced. Based on the T 2 * $$ {\mathrm{T}}_2^{\ast } $$ and T 2 $$ {\mathrm{T}}_2 $$ measurements produced by the proposed method, high-quality T 2 ' $$ {\mathrm{T}}_2^{\prime } $$ maps were obtained, along with neurometabolite maps, from MRSI and TSE data that were acquired in about 9 min. The results obtained from acute stroke and glioma patients demonstrated the feasibility of the proposed method in the clinical setting.

CONCLUSIONS

High-quality T 2 ' $$ {\mathrm{T}}_2^{\prime } $$ maps can be obtained from water-unsuppressed ¹ H-MRSI and TSE data using the proposed method. With further development, this method may lay a foundation for simultaneously imaging oxygenation and neurometabolic alterations of brain disorders.

SNR Enhancement for Multi-TE MRSI Using Joint Low-Dimensional Model and Spatial Constraints

We present a novel method to enhance the SNR for multi-TE MR spectroscopic imaging (MRSI) data by integrating learned nonlinear low-dimensional model and spatial constraints. A deep complex convolutional autoen-coder (DCCAE) was developed to learn a nonlinear low-dimensional representation of the high-dimensional multi-TE 1 H spectroscopy signals. The learned model significantly reduces the data dimension thus serving as an effective constraint for noise reduction. A reconstruction formulation was proposed to integrate the spatiospectral encoding model, the learned model, and a spatial constraint for an SNR-enhancing reconstruction from multi-TE data. The proposed method has been evaluated using both numerical simulations and in vivo brain MRSI experiments. The superior denoising performance of the proposed over alternative methods was demonstrated, both qualitatively and quantitatively. In vivo multi-TE data was used to assess the improved metabolite quantification reproducibility and accuracy achieved by the proposed method. We expect the proposed SNR-enhancing reconstruction to enable faster and/or higher-resolution multi-TE 1 H-MRSI of the brain, potentially useful for various clinical applications.

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.

Whole-brain high-resolution metabolite mapping with 3D compressed-sensing SENSE low-rank 1 H FID-MRSI

There is a growing interest in the neuroscience community to map the distribution of brain metabolites in vivo. Magnetic resonance spectroscopic imaging (MRSI) is often limited by either a poor spatial resolution and/or a long acquisition time, which severely restricts its applications for clinical and research purposes. Building on a recently developed technique of acquisition-reconstruction for 2D MRSI, we combined a fast Cartesian ¹ H-FID-MRSI acquisition sequence, compressed-sensing acceleration, and low-rank total-generalized-variation constrained reconstruction to produce 3D high-resolution whole-brain MRSI with a significant acquisition time reduction. We first evaluated the acceleration performance using retrospective undersampling of a fully sampled dataset. Second, a 20 min accelerated MRSI acquisition was performed on three healthy volunteers, resulting in metabolite maps with 5 mm isotropic resolution. The metabolite maps exhibited the detailed neurochemical composition of all brain regions and revealed parts of the underlying brain anatomy. The latter assessment used previous reported knowledge and a atlas-based analysis to show consistency of the concentration contrasts and ratio across all brain regions. These results acquired on a clinical 3 T MRI scanner successfully combined 3D ¹ H-FID-MRSI with a constrained reconstruction to produce detailed mapping of metabolite concentrations at high resolution over the whole brain, with an acquisition time suitable for clinical or research settings.

Machine Learning-Enabled High-Resolution Dynamic Deuterium MR Spectroscopic Imaging

Deuterium magnetic resonance spectroscopic imaging (DMRSI) has recently been recognized as a potentially powerful tool for noninvasive imaging of brain energy metabolism and tumor. However, the low sensitivity of DMRSI has significantly limited its utility for both research and clinical applications. This work presents a novel machine learning-based method to address this limitation. The proposed method synergistically integrates physics-based subspace modeling and data-driven deep learning for effective denoising, making high-resolution dynamic DMRSI possible. Specifically, a novel subspace model was used to represent the dynamic DMRSI signals; deep neural networks were trained to capture the low-dimensional manifolds of the spectral and temporal distributions of practical dynamic DMRSI data. The learned subspace and manifold structures were integrated via a regularization formulation to remove measurement noise. Theoretical analysis, computer simulations, and in vivo experiments have been conducted to demonstrate the denoising efficacy of the proposed method which enabled high-resolution imaging capability. The translational potential was demonstrated in tumor-bearing rats, where the Warburg effect associated with cancer metabolism and tumor heterogeneity were successfully captured. The new method may not only provide an effective tool to enhance the sensitivity of DMRSI for basic research and clinical applications but also provide a framework for denoising other spatiospectral data.

High-resolution, 3D multi-TE 1 H MRSI using fast spatiospectral encoding and subspace imaging

PURPOSE

To develop a novel method to achieve fast, high-resolution, 3D multi-TE ¹ H-MRSI of the brain.

METHODS

A new multi-TE MRSI acquisition strategy was developed that integrates slab selective excitation with adiabatic refocusing for better volume coverage, rapid spatiospectral encoding, sparse multi-TE sampling, and interleaved water navigators for field mapping and calibration. Special data processing strategies were developed to interpolate the sparsely sampled data, remove nuisance signals, and reconstruct multi-TE spatiospectral distributions with high SNR. Phantom and in vivo experiments have been carried out to demonstrate the capability of the proposed method.

RESULTS

The proposed acquisition can produce multi-TE ¹ H-MRSI data with three TEs at a nominal spatial resolution of 3.4 × 3.4 × 5.3 mm³ in around 20 min. High-SNR brain metabolite spatiospectral reconstructions can be obtained from both a metabolite phantom and in vivo experiments by the proposed method.

CONCLUSION

High-resolution, 3D multi-TE ¹ H-MRSI of the brain can be achieved within clinically feasible time. This capability, with further optimizations, could be translated to clinical applications and neuroscience studies where simultaneously mapping metabolites and neurotransmitters and TE-dependent molecular spectral changes are of interest.

High-Resolution Label-Free Molecular Imaging of Brain Tumor

Molecular imaging has long been recognized as an important tool for diagnosis, characterization, and monitoring of treatment responses of brain tumors. Magnetic resonance spectroscopic imaging (MRSI) is a label-free molecular imaging technique capable of mapping metabolite distributions non-invasively. Several metabolites detectable by MRSI, including Choline, Lactate and N-Acetyl Aspartate, have been proved useful biomarkers for brain tumor characterization. However, clinical application of MRSI has been limited by poor resolution, small spatial coverage, low signal-to-noise ratio and long scan time. This work presents a novel MRSI method for fast, high-resolution metabolic imaging of brain tumor. This method synergistically integrates fast acquisition sequence, sparse sampling, subspace modeling and machine learning to enable 3D mapping of brain metabolites with a spatial resolution of 2.0×3.0×3.0 mm³ in a 7-minute scan. Experimental results obtained from patients with diagnosed brain tumor showed great promise for capturing small-size tumors and revealing intra-tumor metabolic heterogeneities.Clinical Relevance - This paper presents a novel technique for label-free molecular imaging of brain tumor. With further development, this technology may enable many potential clinical applications, from tumor detection, characterization, to assessment of treatment efficacy.

Achieving high-resolution 1H-MRSI of the human brain with compressed-sensing and low-rank reconstruction at 7 Tesla

Low sensitivity MR techniques such as magnetic resonance spectroscopic imaging (MRSI) greatly benefit from the gain in signal-to-noise provided by ultra-high field MR. High-resolution and whole-slab brain MRSI remains however very challenging due to lengthy acquisition, low signal, lipid contamination and field inhomogeneity. In this study, we propose an acquisition-reconstruction scheme that combines 1H free-induction-decay (FID)-MRSI sequence, short TR acquisition, compressed sensing acceleration and low-rank modeling with total-generalized-variation constraint to achieve metabolite imaging in two and three dimensions at 7 Tesla. The resulting images and volumes reveal highly detailed distributions that are specific to each metabolite and follow the underlying brain anatomy. The MRSI method was validated in a high-resolution phantom containing fine metabolite structures, and in five healthy volunteers. This new application of compressed sensing acceleration paves the way for high-resolution MRSI in clinical setting with acquisition times of 5 min for 2D MRSI at 2.5 mm and of 20 min for 3D MRSI at 3.3 mm isotropic.

Absolute choline tissue concentration mapping for prostate cancer localization and characterization using 3D 1 H MRSI without water-signal suppression

PURPOSE

Until now, ¹ H MRSI of the prostate has been performed with suppression of the large water signal to avoid distortions of metabolite signals. However, this signal can be used for absolute quantification and spectral corrections. We investigated the feasibility of water-unsuppressed MRSI in patients with prostate cancer for water signal-mediated spectral quality improvement and determination of absolute tissue levels of choline.

METHODS

Eight prostate cancer patients scheduled for radical prostatectomy underwent multi-parametric MRI at 3 T, including 3D water-unsuppressed semi-LASER MRSI. A postprocessing algorithm was developed to remove the water signal and its artifacts and use the extracted water signal as intravoxel reference for phase and frequency correction of metabolite signals and for absolute metabolite quantification.

RESULTS

Water-unsuppressed MRSI with dedicated postprocessing produced water signal and artifact-free MR spectra throughout the prostate. In all patients, the absolute choline tissue concentration was significantly higher in tumorous than in benign tissue areas (mean ± SD: 7.2 ± 1.4 vs 3.8 ± 0.7 mM), facilitating tumor localization by choline mapping. Tumor tissue levels of choline correlated better with the commonly used (choline + spermine + creatine)/citrate ratio (r = 0.78 ± 0.1) than that of citrate (r = 0.21 ± 0.06). The highest maximum choline concentrations occurred in high-risk cancer foci.

CONCLUSION

This report presents the first successful water-unsuppressed MRSI of the whole prostate. The water signal enabled amelioration of spectral quality and absolute metabolite quantification. In this way, choline tissue levels were identified as tumor biomarker. Choline mapping may serve as a tool in prostate cancer localization and risk scoring in multi-parametric MRI for diagnosis and biopsy procedures.

Uncertainty in denoising of MRSI using low-rank methods

PURPOSE

Low-rank denoising of MRSI data results in an apparent increase in spectral SNR. However, it is not clear if this translates to a lower uncertainty in metabolite concentrations after spectroscopic fitting. Estimation of the true uncertainty after denoising is desirable for downstream analysis in spectroscopy. In this work, the uncertainty reduction from low-rank denoising methods based on spatiotemporal separability and linear predictability in MRSI are assessed. A new method for estimating metabolite concentration uncertainty after denoising is proposed. Automatic rank threshold selection methods are also assessed in simulated low SNR regimes.

METHODS

Assessment of denoising methods is conducted using Monte Carlo simulation of proton MRSI data and by reproducibility of repeated in vivo acquisitions in 5 subjects.

RESULTS

In simulated and in vivo data, spatiotemporal based denoising is shown to reduce the concentration uncertainty, but linear prediction denoising increases uncertainty. Uncertainty estimates provided by fitting algorithms after denoising consistently underestimate actual metabolite uncertainty. However, the proposed uncertainty estimation, based on an analytical expression for entry-wise variance after denoising, is more accurate. It is also shown automated rank threshold selection using Marchenko-Pastur distribution can bias the data in low SNR conditions. An alternative soft-thresholding function is proposed.

CONCLUSION

Low-rank denoising methods based on spatiotemporal separability do reduce uncertainty in MRS(I) data. However, thorough assessment is needed as assessment by SNR measured from residual baseline noise is insufficient given the presence of non-uniform variance. It is also important to select the right rank thresholding method in low SNR cases.

Spectral Wavelet-feature Analysis and Classification Assisted Denoising for enhancing magnetic resonance spectroscopy

Magnetic resonance spectroscopy (MRS) is capable of revealing important biochemical and metabolic information of tissues noninvasively. However, the low concentrations of metabolites often lead to poor signal-to-noise ratio (SNR) and a long acquisition time. Therefore, the applications of MRS in detection and quantitative measurements of metabolites in vivo remain limited. Reducing or even eliminating noise can improve SNR sufficiently to obtain high quality spectra in addition to increasing the number of signal averaging (NSA) or the field strength, both of which are limited in clinical applications. We present a Spectral Wavelet-feature ANalysis and Classification Assisted Denoising (SWANCAD) approach to differentiate signal and noise peaks in magnetic resonance spectra based on their respective wavelet features, followed by removing the identified noise components to improve SNR. The performance of this new denoising approach was evaluated by measuring and comparing SNRs and quantified metabolite levels of low NSA spectra (e.g. NSA = 8) before and after denoising using the SWANCAD approach or by conventional spectral fitting and denoising methods, such as LCModel and wavelet threshold methods, as well as the high NSA spectra (e.g. NSA = 192) recorded in the same sampling volumes. The results demonstrated that SWANCAD offers a more effective way to detect the signals and improve SNR by removing noise from the noisy spectra collected with low NSA or in the subminute scan time (e.g. NSA = 8 or 16 s). The potential applications of SWANCAD include using low NSA to accelerate MRS acquisition while maintaining adequate spectroscopic information for detection and quantification of the metabolites of interest when a limited time is available for an MRS examination in the clinical setting.

Accelerated MR spectroscopic imaging-a review of current and emerging techniques

Over more than 30 years in vivo MR spectroscopic imaging (MRSI) has undergone an enormous evolution from theoretical concepts in the early 1980s to the robust imaging technique that it is today. The development of both fast and efficient sampling and reconstruction techniques has played a fundamental role in this process. State-of-the-art MRSI has grown from a slow purely phase-encoded acquisition technique to a method that today combines the benefits of different acceleration techniques. These include shortening of repetition times, spatial-spectral encoding, undersampling of k-space and time domain, and use of spatial-spectral prior knowledge in the reconstruction. In this way in vivo MRSI has considerably advanced in terms of spatial coverage, spatial resolution, acquisition speed, artifact suppression, number of detectable metabolites and quantification precision. Acceleration not only has been the enabling factor in high-resolution whole-brain 1 H-MRSI, but today is also common in non-proton MRSI (31 P, 2 H and 13 C) and applied in many different organs. In this process, MRSI techniques had to constantly adapt, but have also benefitted from the significant increase of magnetic field strength boosting the signal-to-noise ratio along with high gradient fidelity and high-density receive arrays. In combination with recent trends in image reconstruction and much improved computation power, these advances led to a number of novel developments with respect to MRSI acceleration. Today MRSI allows for non-invasive and non-ionizing mapping of the spatial distribution of various metabolites' tissue concentrations in animals or humans, is applied for clinical diagnostics and has been established as an important tool for neuro-scientific and metabolism research. This review highlights the developments of the last five years and puts them into the context of earlier MRSI acceleration techniques. In addition to 1 H-MRSI it also includes other relevant nuclei and is not limited to certain body regions or specific applications.

Separation of Metabolites and Macromolecules for Short-TE 1H-MRSI Using Learned Component-Specific Representations

Short-echo-time (TE) proton magnetic resonance spectroscopic imaging (MRSI) allows for simultaneously mapping a number of molecules in the brain, and has been recognized as an important tool for studying in vivo biochemistry in various neuroscience and disease applications. However, separation of the metabolite and macromolecule (MM) signals present in the short-TE data with significant spectral overlaps remains a major technical challenge. This work introduces a new approach to solve this problem by integrating imaging physics and representation learning. Specifically, a mixed unsupervised and supervised learning-based strategy was developed to learn the metabolite and MM-specific low-dimensional representations using deep autoencoders. A constrained reconstruction formulation is proposed to integrate the MRSI spatiospectral encoding model and the learned representations as effective constraints for signal separation. An efficient algorithm was developed to solve the resulting optimization problem with provable convergence. Simulation and experimental results have been obtained to demonstrate the component-specific representation power of the learned models and the capability of the proposed method in separating metabolite and MM signals for practical short-TE [Formula: see text]-MRSI data.

Novel proton exchange rate MRI presents unique contrast in brains of ischemic stroke patients

BACKGROUND

To map and quantify the proton exchange rate (k_ex) of brain tissues using improved omega plots in ischemic stroke patients and to investigate whether k_ex can serve as a potential endogenous surrogate imaging biomarker for detecting the metabolic state and the pathologic changes due to ischemic stroke.

NEW METHOD

Three sets of Z-spectra were acquired from seventeen ischemic stroke patients using a spin echo-echo planar imaging sequence with pre-saturation chemical exchange saturation transfer (CEST) pulse at B₁ of 1.5, 2.5, and 3.5 μT, respectively. Pixel-wise k_ex was calculated from improved omega plot of water direct saturation (DS)-removed Z-spectral signals.

RESULTS

The derived k_ex maps can differentiate infarcts from contralateral normal brain tissues with significantly increased signal (893 ± 52 s^-1vs. 739 ± 34 s^-1, P < 0.001).

COMPARISON WITH EXISTING METHOD(S)

The k_ex maps were found to be different from conventional contrasts from diffusion-weighted imaging (DWI), CEST, and semi-solid magnetization transfer (MT) MRI. In brief, k_ex MRI showed larger lesion areas than DWI with different degrees and different lesion contrast compared to CEST and MT.

CONCLUSIONS

In this preliminary translational research, the k_ex MRI based on DS-removed omega plots has been demonstrated for in vivo imaging of clinical ischemic stroke patients. As a noninvasive and unique MRI contrast, k_ex MRI at 3 T may serve as a potential surrogate imaging biomarker for the metabolic changes of stroke and help for monitoring the evolution and the treatment of stroke.

Simultaneous QSM and metabolic imaging of the brain using SPICE: Further improvements in data acquisition and processing

PURPOSE

To achieve high-resolution mapping of brain tissue susceptibility in simultaneous QSM and metabolic imaging.

METHODS

Simultaneous QSM and metabolic imaging was first achieved using SPICE (spectroscopic imaging by exploiting spatiospectral correlation), but the QSM maps thus obtained were at relatively low-resolution (2.0 × 3.0 × 3.0 mm³ ). We overcome this limitation using an improved SPICE data acquisition method with the following novel features: 1) sampling (k, t)-space in dual densities, 2) sampling central k-space fully to achieve nominal spatial resolution of 3.0 × 3.0 × 3.0 mm³ for metabolic imaging, and 3) sampling outer k-space sparsely to achieve spatial resolution of 1.0 × 1.0 × 1.9 mm³ for QSM. To keep the scan time short, we acquired spatiospectral encodings in echo-planar spectroscopic imaging trajectories in central k-space but in CAIPIRINHA (controlled aliasing in parallel imaging results in higher acceleration) trajectories in outer k-space using blipped phase encodings. For data processing and image reconstruction, a union-of-subspaces model was used, effectively incorporating sensitivity encoding, spatial priors, and spectral priors of individual molecules.

RESULTS

In vivo experiments were carried out to evaluate the feasibility and potential of the proposed method. In a 6-min scan, QSM maps at 1.0 × 1.0 × 1.9 mm³ resolution and metabolic maps at 3.0 × 3.0 × 3.0 mm³ nominal resolution were obtained simultaneously. Compared with the original method, the QSM maps obtained using the new method reveal fine-scale brain structures more clearly.

CONCLUSION

We demonstrated the feasibility of achieving high-resolution QSM simultaneously with metabolic imaging using a modified SPICE acquisition method. The improved capability of SPICE may further enhance its practical utility in brain mapping.

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.

High-Resolution 3D Spin-Echo MRSI Using Interleaved Water Navigators, Sparse Sampling and Subspace-Based Processing

This work presents a new method to achieve accelerated, high-resolution magnetic resonance spectroscopic imaging (MRSI) with spin-echo excitations. A new data acquisition strategy is proposed that integrates adiabatic refocusing, elimination of lipid suppression, rapid spatiospectral encoding with sparse (k,t)-space sampling, and interleaved water navigators. This integration leads to a significantly improved combination of volume coverage, spatial resolution (approximately 3 × 3.4 × 4 mm³) and speed (< 10 minutes), while eliminating additional scans for field mapping and coil sensitivity estimation. A data processing strategy that integrates parallel imaging reconstruction and subspace-based processing is devised to produce high-SNR spatiospectral reconstruction from the sparsely sampled, noisy and highresolution MRSI data. Promising in vivo results have been obtained to demonstrate the potential of the proposed method.Clinical relevance- The proposed method enabled volumetric MRSI with a nominal resolution of 3 × 3.4 × 4 mm³ in less than 10 minutes. With further developments and optimizations, the proposed method is expected to be useful for providing molecular-level information of brain functions and diseases, and has the potential to provide new biomarkers for disease diagnosis and treatment monitoring.