High-resolution dynamic 31 P-MRSI using a low-rank tensor model

Alternating low-rank tensor reconstruction for improved multiparametric mapping with cardiovascular MR Multitasking

PURPOSE

To develop a novel low-rank tensor reconstruction approach leveraging the complete acquired data set to improve precision and repeatability of multiparametric mapping within the cardiovascular MR Multitasking framework.

METHODS

A novel approach that alternated between estimation of temporal components and spatial components using the entire data set acquired (i.e., including navigator data and imaging data) was developed to improve reconstruction. The precision and repeatability of the proposed approach were evaluated on numerical simulations, 10 healthy subjects, and 10 cardiomyopathy patients at multiple scan times for 2D myocardial T1/T2 mapping with MR Multitasking and were compared with those of the previous navigator-derived fixed-basis approach.

RESULTS

In numerical simulations, the proposed approach outperformed the previous fixed-basis approach with lower T1 and T2 error against the ground truth at all scan times studied and showed better motion fidelity. In human subjects, the proposed approach showed no significantly different sharpness or T1/T2 measurement and significantly improved T1 precision by 20%-25%, T2 precision by 10%-15%, T1 repeatability by about 30%, and T2 repeatability by 25%-35% at 90-s and 50-s scan times The proposed approach at the 50-s scan time also showed comparable results with that of the previous fixed-basis approach at the 90-s scan time.

CONCLUSION

The proposed approach improved precision and repeatability for quantitative imaging with MR Multitasking while maintaining comparable motion fidelity, T1/T2 measurement, and septum sharpness and had the potential for further reducing scan time from 90 s to 50 s.

Learned Tensor Low-CP-Rank and Bloch Response Manifold Priors for Non-Cartesian MRF Reconstruction

Magnetic resonance fingerprinting (MRF) can rapidly perform simultaneous imaging of multiple tissue parameters. However, the rapid acquisition schemes used in MRF inevitably introduce aliasing artifacts in the recovered tissue fingerprints, reducing the accuracy of the predicted parameter maps. Current regularized reconstruction methods are based on iterative procedures which are usually time-consuming. In addition, most of the current deep learning-based methods for MRF often lack interpretability owing to the black-box nature, and most deep learning-based methods are not applicable for non-Cartesian scenarios, which limits the practical applications. In this paper, we propose a joint reconstruction model incorporating MRF-physics prior and the data correlation constraint for non-Cartesian MRF reconstruction. To avoid time-consuming iterative procedures, we unroll the reconstruction model into a deep neural network. Specifically, we propose a learned CANDECOMP/PARAFAC (CP) decomposition module to exploit the tensor low-rank priors of high-dimensional MRF data, which avoids computationally burdensome singular value decomposition. Inspired by the MRF-physics, we also propose a Bloch response manifold module to learn the mapping between reconstructed MRF data and the multiple parameter maps. Numerical experiments show that the proposed network can reconstruct high-quality MRF data and multiple parameter maps within significantly reduced computational time.

Deuterium echo-planar spectroscopic imaging (EPSI) in the human liver in vivo at 7 T

PURPOSE

To demonstrate the feasibility of deuterium echo-planar spectroscopic imaging (EPSI) to accelerate 3D deuterium metabolic imaging in the human liver at 7 T.

METHODS

A deuterium EPSI sequence, featuring a Hamming-weighted k-space acquisition pattern for the phase-encoding directions, was implemented. Three-dimensional deuterium EPSI and conventional MRSI were performed on a water/acetone phantom and in vivo in the human liver at natural abundance. Moreover, in vivo deuterium EPSI measurements were acquired after oral administration of deuterated glucose. The effect of acquisition time on SNR was evaluated by retrospectively reducing the number of averages.

RESULTS

The SNR of natural abundance deuterated water signal in deuterium EPSI was 6.5% and 5.9% lower than that of MRSI in the phantom and in vivo experiments, respectively. In return, the acquisition time of in vivo EPSI data could be reduced retrospectively to 2 min, beyond the minimal acquisition time of conventional MRSI (of 20 min in this case), while still leaving sufficient SNR. Three-dimensional deuterium EPSI, after administration of deuterated glucose, enabled monitoring of hepatic glucose dynamics with full liver coverage, a spatial resolution of 20 mm isotropic, and a temporal resolution of 9 min 50 s, which could retrospectively be shortened to 2 min.

CONCLUSION

In this work, we demonstrate the feasibility of accelerated 3D deuterium metabolic imaging of the human liver using deuterium EPSI. The acceleration obtained with EPSI can be used to increase temporal and/or spatial resolution, which will be valuable to study tissue metabolism of deuterated compounds over time.

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.

MR-Nucleomics: The study of pathological cellular processes with multinuclear magnetic resonance spectroscopy and imaging in vivo

Clinical medicine has experienced a rapid development in recent decades, during which therapies targeting specific cellular signaling pathways, or specific cell surface receptors, have been increasingly adopted. While these developments in clinical medicine call for improved precision in diagnosis and treatment monitoring, modern medical imaging methods are restricted mainly to anatomical imaging, lagging behind the requirements of precision medicine. Although positron emission tomography and single photon emission computed tomography have been used clinically for studies of metabolism, their applications have been limited by the exposure risk to ionizing radiation, the subsequent limitation in repeated and longitudinal studies, and the incapability in assessing downstream metabolism. Magnetic resonance spectroscopy (MRS) or spectroscopic imaging (MRSI) are, in theory, capable of assessing molecular activities in vivo, although they are often limited by sensitivity. Here, we review some recent developments in MRS and MRSI of multiple nuclei that have potential as molecular imaging tools in the clinic.

Manifold Learning via Linear Tangent Space Alignment (LTSA) for Accelerated Dynamic MRI With Sparse Sampling

The spatial resolution and temporal frame-rate of dynamic magnetic resonance imaging (MRI) can be improved by reconstructing images from sparsely sampled k -space data with mathematical modeling of the underlying spatiotemporal signals. These models include sparsity models, linear subspace models, and non-linear manifold models. This work presents a novel linear tangent space alignment (LTSA) model-based framework that exploits the intrinsic low-dimensional manifold structure of dynamic images for accelerated dynamic MRI. The performance of the proposed method was evaluated and compared to state-of-the-art methods using numerical simulation studies as well as 2D and 3D in vivo cardiac imaging experiments. The proposed method achieved the best performance in image reconstruction among all the compared methods. The proposed method could prove useful for accelerating many MRI applications, including dynamic MRI, multi-parametric MRI, and MR spectroscopic imaging.

Emerging MRI techniques for molecular and functional phenotyping of the diseased heart

Recent advances in cardiac MRI (CMR) capabilities have truly transformed its potential for deep phenotyping of the diseased heart. Long known for its unparalleled soft tissue contrast and excellent depiction of three-dimensional (3D) structure, CMR now boasts a range of unique capabilities for probing disease at the tissue and molecular level. We can look beyond coronary vessel blockages and detect vessel disease not visible on a structural level. We can assess if early fibrotic tissue is being laid down in between viable cardiac muscle cells. We can measure deformation of the heart wall to determine early presentation of stiffening. We can even assess how cardiomyocytes are utilizing energy, where abnormalities are often precursors to overt structural and functional deficits. Finally, with artificial intelligence gaining traction due to the high computing power available today, deep learning has proven itself a viable contender with traditional acceleration techniques for real-time CMR. In this review, we will survey five key emerging MRI techniques that have the potential to transform the CMR clinic and permit early detection and intervention. The emerging areas are: (1) imaging microvascular dysfunction, (2) imaging fibrosis, (3) imaging strain, (4) imaging early metabolic changes, and (5) deep learning for acceleration. Through a concerted effort to develop and translate these areas into the CMR clinic, we are committing ourselves to actualizing early diagnostics for the most intractable heart disease phenotypes.

Metabolite-cycled echo-planar spectroscopic imaging of the human heart

Purpose

Spectroscopic imaging could provide insights into regional cardiac triglyceride variations, but is hampered by relatively long scan times. It is proposed to synergistically combine echo‐planar spectroscopic imaging (EPSI) with motion‐adapted gating, weighted acquisition and metabolite cycling to reduce scan times to less than 10 min while preserving spatial‐spectral quality. The method is compared to single‐voxel measurements and to metabolite‐cycled EPSI with conventional acquisition for assessing triglyceride‐to‐water (TG/W) ratios in the human heart.

Methods

Measurements were performed on 10 healthy volunteers using a clinical 1.5T system. EPSI data was acquired both without and with motion‐adapted gating in combination with weighted acquisition to assess TG/W ratios and relative Cramér‐Rao lower bounds (CRLB) of TG. For comparison, single‐voxel (PRESS) spectra were acquired in the interventricular septum.

Results

Bland–Altman analyses did not show a significant bias in TG/W when comparing both metabolite‐cycled EPSI methods to PRESS for any of the cardiac segments. Scan time was 8.05 ± 2.06 min and 17.91 ± 3.93 min for metabolite‐cycled EPSI with and without motion‐adapted gating and weighted acquisition, respectively, while relative CRLB of TG did not differ significantly between the two methods for any of the cardiac segments.

Conclusions

Metabolite‐cycled EPSI with motion‐adapted gating and weighted acquisition allows detecting TG/W ratios in different regions of the in vivo human heart. Scan time is reduced by more than 2‐fold to less than 10 min as compared to conventional acquisition, while keeping the quality of TG fitting constant.

Click here for author‐reader discussions

Arterial spin labeling MR image denoising and reconstruction using unsupervised deep learning

Arterial spin labeling (ASL) imaging is a powerful magnetic resonance imaging technique that allows to quantitatively measure blood perfusion non-invasively, which has great potential for assessing tissue viability in various clinical settings. However, the clinical applications of ASL are currently limited by its low signal-to-noise ratio (SNR), limited spatial resolution, and long imaging time. In this work, we propose an unsupervised deep learning-based image denoising and reconstruction framework to improve the SNR and accelerate the imaging speed of high resolution ASL imaging. The unique feature of the proposed framework is that it does not require any prior training pairs but only the subject's own anatomical prior, such as T1-weighted images, as network input. The neural network was trained from scratch in the denoising or reconstruction process, with noisy images or sparely sampled k-space data as training labels. Performance of the proposed method was evaluated using in vivo experiment data obtained from 3 healthy subjects on a 3T MR scanner, using ASL images acquired with 44-min acquisition time as the ground truth. Both qualitative and quantitative analyses demonstrate the superior performance of the proposed txtc framework over the reference methods. In summary, our proposed unsupervised deep learning-based denoising and reconstruction framework can improve the image quality and accelerate the imaging speed of ASL imaging.

Sensitivity schemes for dynamic glucose-enhanced magnetic resonance imaging to detect glucose uptake and clearance in mouse brain at 3 T

We investigated three dynamic glucose-enhanced (DGE) MRI methods for sensitively monitoring glucose uptake and clearance in both brain parenchyma and cerebrospinal fluid (CSF) at clinical field strength (3 T). By comparing three sequences, namely, Carr-Purcell-Meiboom-Gill (CPMG), on-resonance variable delay multipulse (onVDMP), and on-resonance spin-lock (onSL), a high-sensitivity DGE MRI scheme with truncated multilinear singular value decomposition (MLSVD) denoising was proposed. The CPMG method showed the highest sensitivity in detecting the parenchymal DGE signal among the three methods, while both onVDMP and onSL were more robust for CSF DGE imaging. Here, onVDMP was applied for CSF imaging, as it displayed the best stability of the DGE results in this study. The truncated MLSVD denoising method was incorporated to further improve the sensitivity. The proposed DGE MRI scheme was examined in mouse brain with 50%/25%/12.5% w/w D-glucose injections. The results showed that this combination could detect DGE signal changes from the brain parenchyma and CSF with as low as a 12.5% w/w D-glucose injection. The proposed DGE MRI schemes could sensitively detect the glucose signal change from brain parenchyma and CSF after D-glucose injection at a clinically relevant concentration, demonstrating high potential for clinical translation.

Mapping an Extended Metabolic Profile of Gliomas Using High-Resolution 31P MRSI at 7T

Phosphorus magnetic resonance spectroscopic imaging (³¹P MRSI) is of particular interest for investigations of patients with brain tumors as it enables to non-invasively assess altered energy and phospholipid metabolism in vivo. However, the limited sensitivity of ³¹P MRSI hampers its broader application at clinical field strengths. This study aimed to identify the additional value of ³¹P MRSI in patients with glioma at ultra-high B₀ = 7T, where the increase in signal-to-noise ratio may foster its applicability for clinical research. High-quality, 3D ³¹P MRSI datasets with an effective voxel size of 5.7 ml were acquired from the brains of seven patients with newly diagnosed glioma. An optimized quantification model was implemented to reliably extract an extended metabolic profile, including low-concentrated metabolites such as extracellular inorganic phosphate, nicotinamide adenine dinucleotide [NAD(H)], and uridine diphosphoglucose (UDPG), which may act as novel tumor markers; a background signal was extracted as well, which affected measures of phosphomonoesters beneficially. Application of this model to the MRSI datasets yielded high-resolution maps of 12 different ³¹P metabolites, showing clear metabolic differences between white matter (WM) and gray matter, and between healthy and tumor tissues. Moreover, differences between tumor compartments in patients with high-grade glioma (HGG), i.e., gadolinium contrast-enhancing/necrotic regions (C+N) and peritumoral edema, could also be suggested from these maps. In the group of patients with HGG, the most significant changes in metabolite intensities were observed in C+N compared to WM, i.e., for phosphocholine +340%, UDPG +54%, glycerophosphoethanolamine −45%, and adenosine-5′-triphosphate −29%. Furthermore, a prominent signal from mobile phospholipids appeared in C+N. In the group of patients with low-grade glioma, only the NAD(H) intensity changed significantly by −28% in the tumor compared to WM. Besides the potential of ³¹P MRSI at 7T to provide novel insights into the biochemistry of gliomas in vivo, the attainable spatial resolutions improve the interpretability of ³¹P metabolite intensities obtained from malignant tissues, particularly when only subtle differences compared to healthy tissues are expected. In conclusion, this pilot study demonstrates that ³¹P MRSI at 7T has potential value for the clinical research of glioma.

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.

3D-printed, patient-specific DIEP flap templates for preoperative planning in breast reconstruction: a prospective case series

Background

Modern imaging technologies, such as computed tomographic angiography (CTA), can be useful for preoperative assessment in deep inferior epigastric artery perforator (DIEP) flap surgery. Planning perforator flap design can lead to improved surgical efficiency. However, current imaging modalities are limited by being displayed on a two-dimensional (2D) surface. In contrast, a 3D-printed model provides tactile feedback that facilitates superior understanding. Hence, we have 3D-printed patient-specific deep inferior epigastric artery perforator (DIEP) templates, in an affordable and convenient manner, for preoperative planning.

Methods

Twenty consecutive patients undergoing 25 immediate or delayed post-mastectomy autologous breast reconstruction with DIEP or muscle-sparing transverse rectus abdominis (MS-TRAM) flaps are recruited prospectively. Using free, open-source softwares (3D Slicer, Autodesk MeshMixer, and Cura) and desktop 3D printers (Ultimaker 3E and Moment), we created a template based on a patient's abdominal wall anatomy from CTA, with holes and lines indicating the position of perforators, their intramuscular course and the DIEA pedicle.

Results

The mean age of patients was 52 [38-67]. There were 15 immediate and 10 delayed reconstructions. 3D printing time took mean 18 hours and 123.7 g of plastic filament, which calculates to a mean material cost of AUD 8.25. DIEP templates accurately identified the perforators and reduced intraoperative perforator identification by 7.29 minutes (P=0.02). However, the intramuscular dissection time was not affected (P=0.34). Surgeons found the template useful for preoperative marking (8.6/10) and planning (7.9/10), but not for intramuscular dissection (5.9/10). There were no immediate flap-related complications.

Conclusions

Our 3D-printed, patient-specific DIEP template is accurate, significantly reduces intraoperative perforator identification time and, hence, may be a useful tool for preoperative planning in autologous breast reconstruction.

A multi-inversion multi-echo spin and gradient echo echo planar imaging sequence with low image distortion for rapid quantitative parameter mapping and synthetic image contrasts

PURPOSE

Brain imaging exams typically take 10-20 min and involve multiple sequential acquisitions. A low-distortion whole-brain echo planar imaging (EPI)-based approach was developed to efficiently encode multiple contrasts in one acquisition, allowing for calculation of quantitative parameter maps and synthetic contrast-weighted images.

METHODS

Inversion prepared spin- and gradient-echo EPI was developed with slice-order shuffling across measurements for efficient acquisition with T1 , T2 , and T 2 ∗ weighting. A dictionary-matching approach was used to fit the images to quantitative parameter maps, which in turn were used to create synthetic weighted images with typical clinical contrasts. Dynamic slice-optimized multi-coil shimming with a B0 shim array was used to reduce B0 inhomogeneity and, therefore, image distortion by >50%. Multi-shot EPI was also implemented to minimize distortion and blurring while enabling high in-plane resolution. A low-rank reconstruction approach was used to mitigate errors from shot-to-shot phase variation.

RESULTS

The slice-optimized shimming approach was combined with in-plane parallel-imaging acceleration of 4× to enable single-shot EPI with more than eight-fold distortion reduction. The proposed sequence efficiently obtained 40 contrasts across the whole-brain in just over 1 min at 1.2 × 1.2 × 3 mm resolution. The multi-shot variant of the sequence achieved higher in-plane resolution of 1 × 1 × 4 mm with good image quality in 4 min. Derived quantitative maps showed comparable values to conventional mapping methods.

CONCLUSION

The approach allows fast whole-brain imaging with quantitative parameter maps and synthetic weighted contrasts. The slice-optimized multi-coil shimming and multi-shot reconstruction approaches result in minimal EPI distortion, giving the sequence the potential to be used in rapid screening applications.

Denoising of hyperpolarized 13 C MR images of the human brain using patch-based higher-order singular value decomposition

PURPOSE

To improve hyperpolarized ¹³ C (HP-¹³ C) MRI by image denoising with a new approach, patch-based higher-order singular value decomposition (HOSVD).

METHODS

The benefit of using a patch-based HOSVD method to denoise dynamic HP-¹³ C MR imaging data was investigated. Image quality and the accuracy of quantitative analyses following denoising were evaluated first using simulated data of [1-¹³ C]pyruvate and its metabolic product, [1-¹³ C]lactate, and compared the results to a global HOSVD method. The patch-based HOSVD method was then applied to healthy volunteer HP [1-¹³ C]pyruvate EPI studies. Voxel-wise kinetic modeling was performed on both non-denoised and denoised data to compare the number of voxels quantifiable based on SNR criteria and fitting error.

RESULTS

Simulation results demonstrated an 8-fold increase in the calculated SNR of [1-¹³ C]pyruvate and [1-¹³ C]lactate with the patch-based HOSVD denoising. The voxel-wise quantification of k_PL (pyruvate-to-lactate conversion rate) showed a 9-fold decrease in standard errors for the fitted k_PL after denoising. The patch-based denoising performed superior to the global denoising in recovering k_PL information. In volunteer data sets, [1-¹³ C]lactate and [¹³ C]bicarbonate signals became distinguishable from noise across captured time points with over a 5-fold apparent SNR gain. This resulted in >3-fold increase in the number of voxels quantifiable for mapping k_PB (pyruvate-to-bicarbonate conversion rate) and whole brain coverage for mapping k_PL .

CONCLUSIONS

Sensitivity enhancement provided by this denoising significantly improved quantification of metabolite dynamics and could benefit future studies by improving image quality, enabling higher spatial resolution, and facilitating the extraction of metabolic information for clinical research.

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.

High-Resolution Oscillating Steady-State fMRI Using Patch-Tensor Low-Rank Reconstruction

The goals of fMRI acquisition include high spatial and temporal resolutions with a high signal to noise ratio (SNR). Oscillating Steady-State Imaging (OSSI) is a new fMRI acquisition method that provides large oscillating signals with the potential for high SNR, but does so at the expense of spatial and temporal resolutions. The unique oscillation pattern of OSSI images makes it well suited for high-dimensional modeling. We propose a patch-tensor low-rank model to exploit the local spatial-temporal low-rankness of OSSI images. We also develop a practical sparse sampling scheme with improved sampling incoherence for OSSI. With an alternating direction method of multipliers (ADMM) based algorithm, we improve OSSI spatial and temporal resolutions with a factor of 12 acquisition acceleration and 1.3 mm isotropic spatial resolution in prospectively undersampled experiments. The proposed model yields high temporal SNR with more activation than other low-rank methods. Compared to the standard grad- ient echo (GRE) imaging with the same spatial-temporal resolution, 3D OSSI tensor model reconstruction demonstrates 2 times higher temporal SNR with 2 times more functional activation.

Tensor image enhancement and optimal multichannel receiver combination analyses for human hyperpolarized 13 C MRSI

PURPOSE

With the initiation of human hyperpolarized ¹³ C (HP-¹³ C) trials at multiple sites and the development of improved acquisition methods, there is an imminent need to maximally extract diagnostic information to facilitate clinical interpretation. This study aims to improve human HP-¹³ C MR spectroscopic imaging through means of Tensor Rank truncation-Image enhancement (TRI) and optimal receiver combination (ORC).

METHODS

A data-driven processing framework for dynamic HP ¹³ C MR spectroscopic imaging (MRSI) was developed. Using patient data sets acquired with both multichannel arrays and single-element receivers from the brain, abdomen, and pelvis, we examined the theory and application of TRI, as well as 2 ORC techniques: whitened singular value decomposition (WSVD) and first-point phasing. Optimal conditions for TRI were derived based on bias-variance trade-off.

RESULTS

TRI and ORC techniques together provided a 63-fold mean apparent signal-to-noise ratio (aSNR) gain for receiver arrays and a 31-fold gain for single-element configurations, which particularly improved quantification of the lower-SNR-[¹³ C]bicarbonate and [1-¹³ C]alanine signals that were otherwise not detectable in many cases. Substantial SNR enhancements were observed for data sets that were acquired even with suboptimal experimental conditions, including delayed (114 s) injection (8× aSNR gain solely by TRI), or from challenging anatomy or geometry, as in the case of a pediatric patient with brainstem tumor (597× using combined TRI and WSVD). Improved correlation between elevated pyruvate-to-lactate conversion, biopsy-confirmed cancer, and mp-MRI lesions demonstrated that TRI recovered quantitative diagnostic information.

CONCLUSION

Overall, this combined approach was effective across imaging targets and receiver configurations and could greatly benefit ongoing and future HP ¹³ C MRI research through major aSNR improvements.

High-sensitivity CEST mapping using a spatiotemporal correlation-enhanced method

PURPOSE

To obtain high-sensitivity CEST maps by exploiting the spatiotemporal correlation between CEST images.

METHODS

A postprocessing method accomplished by multilinear singular value decomposition (MLSVD) was used to enhance the CEST SNR by exploiting the correlation between the Z-spectrum for each voxel and the low-rank property of the overall CEST data. The performance of this method was evaluated using CrCEST in ischemic mouse brain at 11.7 tesla. Then, MLSVD CEST was applied to obtain Cr, amide, and amine CEST maps of the ischemic mouse brain to demonstrate its general applications.

RESULTS

Complex-valued Gaussian noise was added to CEST k-space data to mimic a low SNR situation. MLSVD CEST analysis was able to suppress the noise, recover the degraded CEST peak, and provide better CrCEST quality compared to the smoothing and singular value decomposition (SVD)-based denoising methods. High-resolution Cr, amide, and amine CEST maps of an ischemic stroke using MLSVD CEST suggest that CrCEST is also a sensitive pH mapping method, and a wide range of pH changes can be detected by combing CrCEST with amine CEST at high magnetic fields.

CONCLUSION

MLSVD CEST provides a simple and efficient way to improve the SNR of CEST images.

High-Resolution Dynamic 31P-MR Spectroscopic Imaging for Mapping Mitochondrial Function

OBJECTIVE

To enable non-invasive dynamic metabolic mapping in rodent model studies of mitochondrial function using 31P-MR spectroscopic imaging (MRSI).

METHODS

We developed a novel method for high-resolution dynamic 31P-MRSI. The method synergistically integrates physics-based models of spectral structures, biochemical modeling of molecular dynamics, and subspace learning to capture spatiospectral variations. Fast data acquisition was achieved using rapid spiral trajectories and sparse sampling of (k, t, T)-space; image reconstruction was accomplished using a low-rank tensor-based framework.

RESULTS

The proposed method provided high-resolution dynamic metabolic mapping in rat hindlimb at spatial and temporal resolutions of 4[Formula: see text]2 mm3 and 1.28 s, respectively. This allowed for in vivo mapping of the time-constant of phosphocreatine resynthesis, a well established index of mitochondrial oxidative capacity. Multiple rounds of in vivo experiments were performed to demonstrate reproducibility, and in vitro experiments were used to validate the accuracy of the estimated metabolite maps.

CONCLUSIONS

A new model-based method is proposed to achieve high-resolution dynamic 31P-MRSI. The proposed method's ability to delineate metabolic heterogeneity was demonstrated in rat hindlimb.

SIGNIFICANCE

Abnormal mitochondrial metabolism is a key cellular dysfunction in many prevalent diseases such as diabetes and heart disease; however, current understanding of mitochondrial function is mostly gained from studies on isolated mitochondria under nonphysiological conditions. The proposed method has the potential to open new avenues of research by allowing in vivo and longitudinal studies of mitochondrial dysfunction in disease development and progression.

Dynamic hyperpolarized 13 C MR spectroscopic imaging using SPICE in mouse kidney at 9.4 T

This study aims to investigate the feasibility of dynamic hyperpolarized ¹³ C MR spectroscopic imaging (MRSI) using the SPectroscopic Imaging by exploiting spatiospectral CorrElation (SPICE) technique and an estimation of the spatially resolved conversion constant rate (k_pl ). An acquisition scheme comprising a single training dataset and several imaging datasets was proposed considering hyperpolarized ¹³ C circumstances. The feasibility and advantage of the scheme were investigated in two parts: (a) consistency of spectral basis over time and (b) accuracy of the estimated k_pl . The simulations and in vivo experiments support accurate k_pl estimation with consistent spectral bases. The proposed method was implemented in an enzyme phantom and via in vivo experiments. In the enzyme phantom experiments, spatially resolved homogeneous k_pl maps were observed. In the in vivo experiments, normal diet (ND) mice and high-fat diet (HFD) mice had k_pl (s^-1 ) values of medullar (ND: 0.0119 ± 0.0022, HFD: 0.0195 ± 0.0005) and cortical (ND: 0.0148 ±0.0023, HFD: 0.0224 ±0.0054) regions which were higher than vascular (ND: 0.0087 ±0.0013, HFD: 0.0132 ±0.0050) regions. In particular, the k_pl value in the medullar region exhibited a significant difference between the two diet groups. In summary, the feasibility of using modified SPICE for dynamic hyperpolarized ¹³ C MRSI was demonstrated via simulations and in vivo experiments. The consistency of spectral bases over time and the accuracy of the estimated k_pl values validate the proposed acquisition scheme, which comprises only a single training dataset. The proposed method improved the spatial resolution of dynamic hyperpolarized ¹³ C MRSI, which could be used for k_pl estimation using high signal-to-noise ratio spectral bases.

Free-Breathing Three-Dimensional T1 Mapping of the Heart Using Subspace-Based Data Acquisition and Image Reconstruction

Mapping the longitudinal relaxation time constant (T1) of the myocardium using Magnetic Resonance Imaging (MRI) is an emerging technique for quantitative assessment of the morphology and viability of the myocardium. However, three-dimensional (3D) T1 mapping of the heart is challenging due to the high dimensionality of the signal and the presence of cardiac and respiratory motions. We propose a subspace-based method for free-breathing 3D T1 mapping of the heart without respiratory gating. The image function is represented as a high-order partially separable (PS) function to explore the inherent spatiotemporal correlations of the underlying signal. A special data acquisition scheme enabled by the high-order PS model is used for sparse sampling of the (k,t)-space, where complementary sparse datasets are acquired, each covering only a small portion of the (k,t)-space to characterize a single subspace (spatial or temporal). High-resolution dynamic MR images are reconstructed from the highly undersampled (k,t)-space using low-rank tensor and sparsity constraints. We demonstrate the feasibility of our proposed method using in vivo data obtained from healthy subjects on a 3T MR scanner. The proposed method can enable new clinical applications of T1 mapping in cardiac MR.

Super-Resolution 1H Magnetic Resonance Spectroscopic Imaging Utilizing Deep Learning

Magnetic resonance spectroscopic imaging (SI) is a unique imaging technique that provides biochemical information from in vivo tissues. The 1H spectra acquired from several spatial regions are quantified to yield metabolite concentrations reflective of tissue metabolism. However, since these metabolites are found in tissues at very low concentrations, SI is often acquired with limited spatial resolution. In this work, we test the hypothesis that deep learning is able to upscale low resolution SI, together with the T1-weighted (T1w) image, to reconstruct high resolution SI. We report on a novel densely connected UNet (D-UNet) architecture capable of producing super-resolution spectroscopic images. The inputs for the D-UNet are the T1w image and the low resolution SI image while the output is the high resolution SI. The results of the D-UNet are compared both qualitatively and quantitatively to simulated and in vivo high resolution SI. It is found that this deep learning approach can produce high quality spectroscopic images and reconstruct entire 1H spectra from low resolution acquisitions, which can greatly advance the current SI workflow.

Low-Rank Tensor Models for Improved Multi-Dimensional MRI: Application to Dynamic Cardiac T 1 Mapping

Multi-dimensional, multi-contrast magnetic resonance imaging (MRI) has become increasingly available for comprehensive and time-efficient evaluation of various pathologies, providing large amounts of data and offering new opportunities for improved image reconstructions. Recently, a cardiac phase-resolved myocardial T 1 mapping method has been introduced to provide dynamic information on tissue viability. Improved spatio-temporal resolution in clinically acceptable scan times is highly desirable but requires high acceleration factors. Tensors are well-suited to describe inter-dimensional hidden structures in such multi-dimensional datasets. In this study, we sought to utilize and compare different tensor decomposition methods, without the use of auxiliary navigator data. We explored multiple processing approaches in order to enable high-resolution cardiac phase-resolved myocardial T 1 mapping. Eight different low-rank tensor approximation and processing approaches were evaluated using quantitative analysis of accuracy and precision in T 1 maps acquired in six healthy volunteers. All methods provided comparable T 1 values. However, the precision was significantly improved using local processing, as well as a direct tensor rank approximation. Low-rank tensor approximation approaches are well-suited to enable dynamic T 1 mapping at high spatio-temporal resolutions.

GRASP-Pro: imProving GRASP DCE-MRI through self-calibrating subspace-modeling and contrast phase automation

PURPOSE

To propose a highly accelerated, high-resolution dynamic contrast-enhanced MRI (DCE-MRI) technique called GRASP-Pro (golden-angle radial sparse parallel imaging with imProved performance) through a joint sparsity and self-calibrating subspace constraint with automated selection of contrast phases.

METHODS

GRASP-Pro reconstruction enforces a combination of an explicit low-rank subspace-constraint and a temporal sparsity constraint. The temporal basis used to construct the subspace is learned from an intermediate reconstruction step using the low-resolution portion of radial k-space, which eliminates the need for generating the basis using auxiliary data or a physical signal model. A convolutional neural network was trained to generate the contrast enhancement curve in the artery, from which clinically relevant contrast phases are automatically selected for evaluation. The performance of GRASP-Pro was demonstrated for high spatiotemporal resolution DCE-MRI of the prostate and was compared against standard GRASP in terms of overall image quality, image sharpness, and residual streaks and/or noise level.

RESULTS

Compared to GRASP, GRASP-Pro reconstructed dynamic images with enhanced sharpness, less residual streaks and/or noise, and finer delineation of the prostate without prolonging reconstruction time. The image quality improvement reached statistical significance (P < 0.05) in all the assessment categories. The neural network successfully generated contrast enhancement curves in the artery, and corresponding peak enhancement indexes correlated well with that from the manual selection.

CONCLUSION

GRASP-Pro is a promising method for rapid and continuous DCE-MRI. It enables superior reconstruction performance over standard GRASP and allows reliable generation of artery enhancement curve to guide the selection of desired contrast phases for improving the efficiency of GRASP MRI workflow.

Comparison of compressed sensing reconstruction algorithms for 31P magnetic resonance spectroscopic imaging

Phosphorus MR spectroscopy and spectroscopic imaging (³¹P-MRS/MRSI) provide information about energy metabolism, membrane degradation and pH in vivo. In spite of their proven utility, ³¹P-MRS/MRSI are not often used primarily because of the challenges imposed by the low sensitivity and low concentration of metabolites leading to low signal to noise ratio (SNR), coarse spatial resolution and prolonged acquisition time. More recently there has been considerable interest in compressed sensing as an acceleration method for MR signal acquisition. This approach takes advantage of the intrinsic sparsity of the spectral data. In this work, we present a ³¹P-MRSI sequence that combines a flyback EPSI trajectory and compressed sensing, and we compared two different reconstruction methods, L1 norm minimization and low rank Hankel matrix completion. Our phantom results showed good preservation of spectral quality for both ×2.0 and ×3.0 acceleration factors, using both CS reconstruction methods. However, in vivo ³¹P-MRS brain data showed the low rank reconstruction approach was most suitable. Overall, this study shows the feasibility of combining a flyback EPSI trajectory and compressed sensing in the acquisition of ³¹P-MRSI as well as the better suitability of a low rank reconstruction approach.

Dynamic Imaging of Glucose and Lactate Metabolism by 13C-MRS without Hyperpolarization

Metabolic reprogramming is one of the defining features of cancer and abnormal metabolism is associated with many other pathologies. Molecular imaging techniques capable of detecting such changes have become essential for cancer diagnosis, treatment planning, and surveillance. In particular, ¹⁸F-FDG (fluorodeoxyglucose) PET has emerged as an essential imaging modality for cancer because of its unique ability to detect a disturbed molecular pathway through measurements of glucose uptake. However, FDG-PET has limitations that restrict its usefulness in certain situations and the information gained is limited to glucose uptake only.¹³C magnetic resonance spectroscopy theoretically has certain advantages over FDG-PET, but its inherent low sensitivity has restricted its use mostly to single voxel measurements unless dissolution dynamic nuclear polarization (dDNP) is used to increase the signal, which brings additional complications for clinical use. We show here a new method of imaging glucose metabolism in vivo by MRI chemical shift imaging (CSI) experiments that relies on a simple, but robust and efficient, post-processing procedure by the higher dimensional analog of singular value decomposition, tensor decomposition. Using this procedure, we achieve an order of magnitude increase in signal to noise in both dDNP and non-hyperpolarized non-localized experiments without sacrificing accuracy. In CSI experiments an approximately 30-fold increase was observed, enough that the glucose to lactate conversion indicative of the Warburg effect can be imaged without hyper-polarization with a time resolution of 12s and an overall spatial resolution that compares favorably to ¹⁸F-FDG PET.

Dynamic 31 P spectroscopic imaging of skeletal muscles combining flyback echo-planar spectroscopic imaging and compressed sensing

PURPOSE

Dynamic phosphorus MR spectroscopic imaging (³¹ P-MRSI) experiments require temporal resolution on the order of seconds to concurrently assess different muscle groups. A highly accelerated pulse sequence combining flyback echo-planar spectroscopic imaging (EPSI) and compressed sensing was developed and tested in a phantom and healthy humans during an exercise-recovery challenge of the lower leg muscles, using a clinical 3T MRI.

METHODS

A flyback EPSI readout designed to achieve 2.25 × 2.25 cm² resolution over a 18 × 18 cm² field of view (i.e., 8 × 8 matrix) was combined with compressed sensing through the inclusion of pseudorandom gradient blips to sub-sample the ky-kt dimensions by a factor of 2.7×, achieving a temporal resolution of 9 s. The sequence was first tested in a phantom to assess performance compared to fully sampled EPSI (fidEPSI) and phase encoded chemical shift imaging (fidCSI). Then, tests were performed in 11 healthy volunteers during an exercise-recovery challenge of the lower leg muscles. Voxels containing tissue from different muscle groups were evaluated measuring percentage phosphocreatine (%PCr) depletion, time constant of PCr recovery (τPCr) and intracellular pH at rest and following exercise.

RESULTS

The sequence was capable to track the dynamic PCr response of multiple muscles simultaneously. No statistical differences were found in the metabolite ratio, pH or linewidth when compared with fidEPSI and fidCSI in the phantom study. Dynamic experiments showed differences in PCr depletion when comparing soleus with gastrocnemius muscles. Intracellular pH, τPCr and %PCr decrease were consistent with reported values.

CONCLUSION

Highly accelerated ³¹ P-MRSI combining flyback EPSI and compressed sensing is capable of assessing concurrent energy metabolism in multiple muscle groups using a clinical 3T MR system.

Fast high-resolution brain metabolite mapping on a clinical 3T MRI by accelerated 1 H-FID-MRSI and low-rank constrained reconstruction

PURPOSE

Epitomizing the advantages of ultra short echo time and no chemical shift displacement error, high-resolution-free induction decay magnetic resonance spectroscopic imaging (FID-MRSI) sequences have proven to be highly effective in providing unbiased characterizations of metabolite distributions. However, its merits are often overshadowed in high-resolution settings by reduced signal-to-noise ratios resulting from the smaller voxel volumes procured by extensive phase encoding and the related acquisition times.

METHODS

To address these limitations, we here propose an acquisition and reconstruction scheme that offers both implicit dataset denoising and acquisition acceleration. Specifically, a slice selective high-resolution FID-MRSI sequence was implemented. Spectroscopic datasets were processed to remove fat contamination, and then reconstructed using a total generalized variation (TGV) regularized low-rank model. We further measured reconstruction performance for random undersampled data to assess feasibility of a compressed-sensing SENSE acceleration scheme. Performance of the lipid suppression was assessed using an ad hoc phantom, while that of the low-rank TGV reconstruction model was benchmarked using simulated MRSI data. To assess real-world performance, 2D FID-MRSI acquisitions of the brain in healthy volunteers were reconstructed using the proposed framework.

RESULTS

Results from the phantom and simulated data demonstrate that skull lipid contamination is effectively removed and that data reconstruction quality is improved with the low-rank TGV model. Also, we demonstrated that the presented acquisition and reconstruction methods are compatible with a compressed-sensing SENSE acceleration scheme.

CONCLUSIONS

An original reconstruction pipeline for 2D 1 H-FID-MRSI datasets was presented that places high-resolution metabolite mapping on 3T MR scanners within clinically feasible limits.

Quantitative imaging of brain energy metabolisms and neuroenergetics using in vivo X-nuclear 2H, 17O and 31P MRS at ultra-high field

Brain energy metabolism relies predominantly on glucose and oxygen utilization to generate biochemical energy in the form of adenosine triphosphate (ATP). ATP is essential for maintaining basal electrophysiological activities in a resting brain and supporting evoked neuronal activity under an activated state. Studying complex neuroenergetic processes in the brain requires sophisticated neuroimaging techniques enabling noninvasive and quantitative assessment of cerebral energy metabolisms and quantification of metabolic rates. Recent state-of-the-art in vivo X-nuclear MRS techniques, including 2H, 17O and 31P MRS have shown promise, especially at ultra-high fields, in the quest for understanding neuroenergetics and brain function using preclinical models and in human subjects under healthy and diseased conditions.

Magnetic resonance multitasking for motion-resolved quantitative cardiovascular imaging

Quantitative cardiovascular magnetic resonance (CMR) imaging can be used to characterize fibrosis, oedema, ischaemia, inflammation and other disease conditions. However, the need to reduce artefacts arising from body motion through a combination of electrocardiography (ECG) control, respiration control, and contrast-weighting selection makes CMR exams lengthy. Here, we show that physiological motions and other dynamic processes can be conceptualized as multiple time dimensions that can be resolved via low-rank tensor imaging, allowing for motion-resolved quantitative imaging with up to four time dimensions. This continuous-acquisition approach, which we name cardiovascular MR multitasking, captures — rather than avoids — motion, relaxation and other dynamics to efficiently perform quantitative CMR without the use of ECG triggering or breath holds. We demonstrate that CMR multitasking allows for T₁ mapping, T₁-T₂ mapping and time-resolved T₁ mapping of myocardial perfusion without ECG information and/or in free-breathing conditions. CMR multitasking may provide a foundation for the development of setup-free CMR imaging for the quantitative evaluation of cardiovascular health.

High resolution hyperpolarized 13 C MRSI using SPICE at 9.4T

PURPOSE

To test the feasibility of using the SPICE (SPectroscopic Imaging by exploiting spatiospectral CorrElation) technique, which uses the partial separability of spectroscopic data, for high resolution hyperpolarized (HP) ¹³ C spectroscopic imaging.

METHODS

Numerical simulations were performed to investigate the impact of transient HP signals on SPICE reconstruction. Furthermore, spectroscopic imaging exams from SPICE and conventional EPSI (echo-planar spectroscopic imaging) were simulated for comparison. For in vivo experiments, HP ¹³ C SPICE was performed in a mouse kidney by means of the injection of HP [1-¹³ C] pyruvate at 9.4T.

RESULTS

The variation of lactate/pyruvate from the simulated SPICE was less than 4% under various factors that affect the transient HP signal, suggesting that the impact is negligible. We found that while HP ¹³ C EPSI was limited to the low signal-to-noise ratio (SNR) of lactate, these limitations were mitigated through HP ¹³ C SPICE, facilitating the improved SNR of lactate and the distinction of tissues. Acquisition of a high resolution HP ¹³ C spectroscopic image was possible for the in vivo experiments. With the fine structural information, the acquired image showed higher signal of pyruvate and lactate in the renal cortices than in the medullas, which is known to be attributed to higher activity of lactate dehydrogenase.

CONCLUSION

The feasibility of HP ¹³ C SPICE was investigated. Simulation studies were conducted and in vivo experiments were performed in the mouse kidney at 9.4T. Results confirmed that a high resolution HP ¹³ C spectroscopic image with adequate spectral resolution can be obtained. Magn Reson Med 80:703-710, 2018. © 2018 International Society for Magnetic Resonance in Medicine.

Assessing tissue metabolism by phosphorous-31 magnetic resonance spectroscopy and imaging: a methodology review

Many human diseases are caused by an imbalance between energy production and demand. Magnetic resonance spectroscopy (MRS) and magnetic resonance imaging (MRI) provide the unique opportunity for in vivo assessment of several fundamental events in tissue metabolism without the use of ionizing radiation. Of particular interest, phosphate metabolites that are involved in ATP generation and utilization can be quantified noninvasively by phosphorous-31 (31P) MRS/MRI. Furthermore, 31P magnetization transfer (MT) techniques allow in vivo measurement of metabolic fluxes via creatine kinase (CK) and ATP synthase. However, a major impediment for the clinical applications of 31P-MRS/MRI is the prohibitively long acquisition time and/or the low spatial resolution that are necessary to achieve adequate signal-to-noise ratio. In this review, current 31P-MRS/MRI techniques used in basic science and clinical research are presented. Recent advances in the development of fast 31P-MRS/MRI methods are also discussed.