High-resolution (1) H-MRSI of the brain using SPICE: Data acquisition and image reconstruction

Phase-constrained reconstruction of high-resolution multi-shot diffusion weighted image

Diffusion weighted imaging (DWI) is a unique examining method in tumor diagnosis, acute stroke evaluation. Single-shot echo planar imaging is currently conventional method for DWI. However, single-shot DWI suffers from image distortion, blurring and low spatial resolution. Although multi-shot DWI improves image resolution, it brings phase variations among different shots at the same time. In this paper, we introduce a smooth phase constraint of each shot image into multi-shot navigator-free DWI reconstruction by imposing the low-rankness of Hankel matrix constructed from the k-space data. Furthermore, we exploit the partial sum minimization of singular values to constrain the low-rankness of Hankel matrix. Results on brain imaging data show that the proposed method outperforms the state-of-the-art methods in terms of artifacts removal and our method potentially has the ability to reconstruct high number of shot of DWI.

Constrained Magnetic Resonance Spectroscopic Imaging by Learning Nonlinear Low-Dimensional Models

Magnetic resonance spectroscopic imaging (MRSI) is a powerful molecular imaging modality but has very limited speed, resolution, and SNR tradeoffs. Construction of a low-dimensional model to effectively reduce the dimensionality of the imaging problem has recently shown great promise in improving these tradeoffs. This paper presents a new approach to model and reconstruct the spectroscopic signals by learning a nonlinear low-dimensional representation of the general MR spectra. Specifically, we trained a deep neural network to capture the low-dimensional manifold, where the high-dimensional spectroscopic signals reside. A regularization formulation is proposed to effectively integrate the learned model and physics-based data acquisition model for MRSI reconstruction with the capability to incorporate additional spatiospectral constraints. An efficient numerical algorithm was developed to solve the associated optimization problem involving back-propagating the trained network. Simulation and experimental results were obtained to demonstrate the representation power of the learned model and the ability of the proposed formulation in producing SNR-enhancing reconstruction from the practical MRSI data.

Super-Resolution Whole-Brain 3D MR Spectroscopic Imaging for Mapping D-2-Hydroxyglutarate and Tumor Metabolism in Isocitrate Dehydrogenase 1-mutated Human Gliomas

Background Isocitrate dehydrogenase (IDH) mutations are highly frequent in glioma, producing high levels of the oncometabolite D-2-hydroxyglutarate (D-2HG). Hence, D-2HG represents a valuable imaging marker for IDH-mutated human glioma. Purpose To develop and evaluate a super-resolution three-dimensional (3D) MR spectroscopic imaging strategy to map D-2HG and tumor metabolism in IDH-mutated human glioma. Materials and Methods Between March and September 2018, participants with IDH1-mutated gliomas and healthy participants were prospectively scanned with a 3-T whole-brain 3D MR spectroscopic imaging protocol optimized for D-2HG. The acquired D-2HG maps with a voxel size of 5.2 × 5.2 × 12 mm were upsampled to a voxel size of 1.7 × 1.7 × 3 mm using a super-resolution method that combined weighted total variation, feature-based nonlocal means, and high-spatial-resolution anatomic imaging priors. Validation with simulated healthy and patient data and phantom measurements was also performed. The Mann-Whitney U test was used to check that the proposed super-resolution technique yields the highest peak signal-to-noise ratio and structural similarity index. Results Three participants with IDH1-mutated gliomas (mean age, 50 years ± 21 [standard deviation]; two men) and three healthy participants (mean age, 32 years ± 3; two men) were scanned. Twenty healthy participants (mean age, 33 years ± 5; 16 men) underwent a simulation of upsampled MR spectroscopic imaging. Super-resolution upsampling improved peak signal-to-noise ratio and structural similarity index by 62% (P < .05) and 7.3% (P < .05), respectively, for simulated data when compared with spline interpolation. Correspondingly, the proposed method significantly improved tissue contrast and structural information for the acquired 3D MR spectroscopic imaging data. Conclusion High-spatial-resolution whole-brain D-2-hydroxyglutarate imaging is possible in isocitrate dehydrogenase 1-mutated human glioma by using a super-resolution framework to upsample three-dimensional MR spectroscopic images acquired at lower resolution. © RSNA, 2020 Online supplemental material is available for this article. See also the editorial by Huang and Lin in this issue.

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.

Fast Imaging for Hyperpolarized MR Metabolic Imaging

MRI with hyperpolarized carbon-13 agents has created a new type of noninvasive, in vivo metabolic imaging that can be applied in cell, animal, and human studies. The use of ¹³ C-labeled agents, primarily [1-¹³ C]pyruvate, enables monitoring of key metabolic pathways with the ability to image substrate and products based on their chemical shift. Over 10 sites worldwide are now performing human studies with this new approach for studies of cancer, heart disease, liver disease, and kidney disease. Hyperpolarized metabolic imaging studies must be performed within several minutes following creation of the hyperpolarized agent due to irreversible decay of the net magnetization back to equilibrium, so fast imaging methods are critical. The imaging methods must include multiple metabolites, separated based on their chemical shift, which are also undergoing rapid metabolic conversion (via label exchange), further exacerbating the challenges of fast imaging. This review describes the state-of-the-art in fast imaging methods for hyperpolarized metabolic imaging. This includes the approach and tradeoffs between three major categories of fast imaging methods-fast spectroscopic imaging, model-based strategies, and metabolite specific imaging-as well additional options of parallel imaging, compressed sensing, tailored RF flip angles, refocused imaging methods, and calibration methods that can improve the scan coverage, speed, signal-to-noise ratio (SNR), resolution, and/or robustness of these studies. To date, these approaches have produced extremely promising initial human imaging results. Improvements to fast hyperpolarized metabolic imaging methods will provide better coverage, SNR, resolution, and reproducibility for future human imaging studies. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY STAGE: 1.

Ultrafast magnetic resonance spectroscopic imaging using SPICE with learned subspaces

PURPOSE

To develop a subspace learning method for the recently proposed subspace-based MRSI approach known as SPICE, and achieve ultrafast 1 H-MRSI of the brain.

THEORY AND METHODS

A novel strategy is formulated to learn a low-dimensional subspace representation of MR spectra from specially acquired training data and use the learned subspace for general MRSI experiments. Specifically, the subspace learning problem is formulated as learning "empirical" distributions of molecule-specific spectral parameters (e.g., concentrations, lineshapes, and frequency shifts) by integrating physics-based model and the training data. The learned spectral parameters and quantum mechanical simulation basis can then be combined to construct acquisition-specific subspace for spatiospectral encoding and processing. High-resolution MRSI acquisitions combining ultrashort-TE/short-TR excitation, sparse sampling, and the elimination of water suppression have been performed to evaluate the feasibility of the proposed method.

RESULTS

The accuracy of the learned subspace and the capability of the proposed method in producing high-resolution 3D 1 H metabolite maps and high-quality spatially resolved spectra (with a nominal resolution of ∼2.4 × 2.4 × 3 mm3 in 5 minutes) were demonstrated using phantom and in vivo studies. By eliminating water suppression, we are also able to extract valuable information from the water signals for data processing ( B 0 map, frequency drift, and coil sensitivity) as well as for mapping tissue susceptibility and relaxation parameters.

CONCLUSIONS

The proposed method enables ultrafast 1 H-MRSI of the brain using a learned subspace, eliminating the need of acquiring subject-dependent navigator data (known as D 1 ) in the original SPICE technique. It represents a new way to perform MRSI experiments and an important step toward practical applications of high-resolution MRSI.

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.

Improved Low-Rank Filtering of MR Spectroscopic Imaging Data With Pre-Learnt Subspace and Spatial Constraints

OBJECTIVE

To investigate the use of pre-learnt subspace and spatial constraints for denoising magnetic resonance spectroscopic imaging (MRSI) data.

METHOD

We exploit the partial separability or subspace structures of high-dimensional MRSI data for denoising. More specifically, we incorporate a subspace model with pre-learnt spectral basis into the low-rank approximation (LORA) method. Spectral basis is determined based on empirical prior distributions of the spectral parameters variations learnt from auxiliary training data; spatial priors are also incorporated as is done in LORA to further improve denoising performance.

RESULTS

The effects of the explicit subspace and spatial constraints in reducing estimation bias and variance have been analyzed using Cramér-Rao Lower bound analysis, Monte-Carlo study, and experimental study.

CONCLUSION

The denoising effectiveness of LORA can be significantly improved by incorporating pre-learnt spectral basis and spatial priors into LORA.

SIGNIFICANCE

This study provides an effective method for denoising MRSI data along with comprehensive analyses of its performance. The proposed method is expected to be useful for a wide range of studies using MRSI.

Simultaneous metabolic and functional imaging of the brain using SPICE

PURPOSE

To enable simultaneous high-resolution mapping of brain function and metabolism.

METHODS

An encoding scheme was designed for interleaved acquisition of functional MRI (fMRI) data in echo volume imaging trajectories and MR spectroscopic imaging (MRSI) data in echo-planar spectroscopic imaging trajectories. The scheme eliminates water and lipid suppression and utilizes free induction decay signals to encode both functional and metabolic information with ultrashort TE, short TR, and sparse sampling of k , t -space. A subspace-based image reconstruction method was introduced for processing both the fMRI and MRSI data. The complementary information in the fMRI and MRSI data sets was also utilized to improve image reconstruction in the presence of intrascan head motion, field drift, and tissue susceptibility changes.

RESULTS

In-vivo experimental results were obtained from healthy human subjects in resting-state fMRI/MRSI experiments. In these experiments, the proposed method was able to simultaneously acquire metabolic and functional information from the brain in high resolution. For scans of 6.5 minutes, we achieved 3.0 × 3.0 × 1.8 mm³ spatial resolution for fMRI, 1.9 × 2.5 × 3.0 mm³ nominal spatial resolution for MRSI, and 1.9 × 1.9 × 1.8 mm³ nominal spatial resolution for quantitative susceptibility maps.

CONCLUSION

This work demonstrates the feasibility of simultaneous high-resolution mapping of brain function and metabolism with improved spatial resolution and synergistic image reconstruction.

Echo planar time-resolved imaging (EPTI)

PURPOSE

To develop an efficient distortion- and blurring-free multi-shot EPI technique for time-resolved multiple-contrast and/or quantitative imaging.

METHODS

EPI is a commonly used sequence but suffers from geometric distortions and blurring. Here, we introduce a new multi-shot EPI technique termed echo planar time-resolved imaging (EPTI), which has the ability to rapidly acquire distortion- and blurring-free multi-contrast data set. The EPTI approach performs encoding in ky -t space and uses a new highly accelerated spatio-temporal CAIPI sampling trajectory to take advantage of signal correlation along these dimensions. Through this acquisition and a B0 -informed parallel imaging reconstruction, hundreds of "time-resolved" distortion- and blurring-free images at different TEs across the EPI readout window can be created at sub-millisecond temporal increments using a small number of EPTI shots. Moreover, a method for self-estimation and correction of shot-to-shot B0 variations was developed. Simultaneous multi-slice acquisition was also incorporated to further improve the acquisition efficiency.

RESULTS

We evaluated EPTI under varying simulated acceleration factors, B0 -inhomogeneity, and shot-to-shot B0 variations to demonstrate its ability to provide distortion- and blurring-free images at multiple TEs. Two variants of EPTI were demonstrated in vivo at 3T: (1) a combined gradient- and spin-echo EPTI for quantitative mapping of T2 , T2* , proton density, and susceptibility at 1.1 × 1.1 × 3 mm3 whole-brain in 28 s (0.8 s/slice), and (2) a gradient-echo EPTI, for multi-echo and quantitative T2* fMRI at 2 × 2 × 3 mm3 whole-brain at a 3.3 s temporal resolution.

CONCLUSION

EPTI is a new approach for multi-contrast and/or quantitative imaging that can provide fast acquisition of distortion- and blurring-free images at multiple TEs.

Fast data acquisition techniques in magnetic resonance spectroscopic imaging

Magnetic resonance spectroscopic imaging (MRSI) is an important technique for assessing the spatial variation of metabolites in vivo. The long scan times in MRSI limit clinical applicability due to patient discomfort, increased costs, motion artifacts, and limited protocol flexibility. Faster acquisition strategies can address these limitations and could potentially facilitate increased adoption of MRSI into routine clinical protocols with minimal addition to the current anatomical and functional acquisition protocols in terms of imaging time. Not surprisingly, a lot of effort has been devoted to the development of faster MRSI techniques that aim to capture the same underlying metabolic information (relative metabolite peak areas and spatial distribution) as obtained by conventional MRSI, in greatly reduced time. The gain in imaging time results, in some cases, in a loss of signal-to-noise ratio and/or in spatial and spectral blurring. This review examines the current techniques and advances in fast MRSI in two and three spatial dimensions and their applications. This review categorizes the acceleration techniques according to their strategy for acquisition of the k-space. Techniques such as fast/turbo-spin echo MRSI, echo-planar spectroscopic imaging, and non-Cartesian MRSI effectively cover the full k-space in a more efficient manner per T_R . On the other hand, techniques such as parallel imaging and compressed sensing acquire fewer k-space points and employ advanced reconstruction algorithms to recreate the spatial-spectral information, which maintains statistical fidelity in test conditions (ie no statistically significant differences on voxel-wise comparisions) with the fully sampled data. The advantages and limitations of each state-of-the-art technique are reviewed in detail, concluding with a note on future directions and challenges in the field of fast spectroscopic imaging.

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.

Pros and cons of ultra-high-field MRI/MRS for human application

Magnetic resonance imaging and spectroscopic techniques are widely used in humans both for clinical diagnostic applications and in basic research areas such as cognitive neuroimaging. In recent years, new human MR systems have become available operating at static magnetic fields of 7 T or higher (≥300 MHz proton frequency). Imaging human-sized objects at such high frequencies presents several challenges including non-uniform radiofrequency fields, enhanced susceptibility artifacts, and higher radiofrequency energy deposition in the tissue. On the other side of the scale are gains in signal-to-noise or contrast-to-noise ratio that allow finer structures to be visualized and smaller physiological effects to be detected. This review presents an overview of some of the latest methodological developments in human ultra-high field MRI/MRS as well as associated clinical and scientific applications. Emphasis is given to techniques that particularly benefit from the changing physical characteristics at high magnetic fields, including susceptibility-weighted imaging and phase-contrast techniques, imaging with X-nuclei, MR spectroscopy, CEST imaging, as well as functional MRI. In addition, more general methodological developments such as parallel transmission and motion correction will be discussed that are required to leverage the full potential of higher magnetic fields, and an overview of relevant physiological considerations of human high magnetic field exposure is provided.

Reduction of Acquisition time using Partition of the sIgnal Decay in Spectroscopic Imaging technique (RAPID-SI)

To overcome long acquisition times of Chemical Shift Imaging (CSI), a new Magnetic Resonance Spectroscopic Imaging (MRSI) technique called Reduction of Acquisition time by Partition of the sIgnal Decay in Spectroscopic Imaging (RAPID-SI) using blipped phase encoding gradients inserted during signal acquisition was developed. To validate the results using RAPID-SI and to demonstrate its usefulness in terms of acquisition time and data quantification; simulations, phantom and in vivo studies were conducted, and the results were compared to standard CSI. The method was based upon the partition of a magnetic resonance spectroscopy (MRS) signal into sequential sub-signals encoded using blipped phase encoding gradients inserted during signal acquisition at a constant time interval. The RAPID-SI technique was implemented on a clinical 3 T Siemens scanner to demonstrate its clinical utility. Acceleration of data collection was performed by inserting R (R = acceleration factor) blipped gradients along a given spatial direction during data acquisition. Compared to CSI, RAPID-SI reduced acquisition time by the acceleration factor R. For example, a 2D 16x16 data set acquired in about 17 min with CSI, was reduced to approximately 2 min with the RAPID-SI (R = 8). While the SNR of the acquired RAPID-SI signal was lower compared to CSI by approximately the factor √R, it can be improved after data pre-processing and reconstruction. Compared to CSI, RAPID-SI reduces acquisition time, while preserving metabolites information. Furthermore, the method is flexible and could be combined with other acceleration methods such as Parallel Imaging.

The application of magnetic resonance fingerprinting to single voxel proton spectroscopy

Magnetic resonance fingerprinting has been proposed as a method for undersampling k-space while simultaneously yielding multiparametric tissue maps. In the context of single voxel spectroscopy, fingerprinting can provide a unified framework for parameter estimation. We demonstrate the utility of such a magnetic resonance spectroscopic fingerprinting (MRSF) framework for simultaneously quantifying metabolite concentrations, T₁ and T₂ relaxation times and transmit inhomogeneity for major singlets of N-acetylaspartate, creatine and choline. This is achieved by varying T_R , T_E and the flip angle of the first pulse in a PRESS sequence between successive excitations (i.e. successive T_R values). The need for multiparametric schemes such as MRSF for accurate medical diagnostics is demonstrated with the aid of realistic in vivo simulations; these show that certain schemes lead to substantial increases to the area under receiver operating characteristics of metabolite concentrations, when viewed as classifiers of pathologies. Numerical simulations and phantom and in vivo experiments using several different schedules of variable length demonstrate superior precision and accuracy for metabolite concentrations and longitudinal relaxation, and similar performance for the quantification of transverse relaxation.

Miniature pig magnetic resonance spectroscopy model of normal adolescent brain development

BACKGROUND

We are developing the miniature pig (Sus scrofa domestica), an in-vivo translational, gyrencephalic model for brain development, as an alternative to laboratory rodents/non-human primates. We analyzed longitudinal changes in adolescent pigs using proton magnetic resonance spectroscopy (¹H-MRS) and examined the relationship with white matter (WM) integrity derived from diffusion weighted imaging (DWI).

NEW METHOD

Twelve female Sinclair^™ pigs underwent three imaging/spectroscopy sessions every 23.95 ± 3.73 days beginning at three months of age using a clinical 3 T scanner. ¹H-MRS data were collected using 1.2 × 1.0 × 3.0 cm voxels placed in left and right hemisphere WM using a Point Resolved Spectroscopy sequence (TR = 2000 ms, TE = 30 ms). Concentrations of N-acetylaspartate, myo-inositol (MI), glutamate + glutamine, choline, creatine, and macromolecules (MM) 09 and 14 were averaged from both hemispheres. DWI data were collected using 15 shells of b-values (b = 0-3500 s/mm2) with 32 directions/shell and fit using the WM Tract Integrity model to calculate fractional anisotropy (FA), kurtosis anisotropy (KA) and permeability-diffusivity index.

RESULTS

MI and MM09 significantly declined with age. Increased FA and KA significantly correlated with decline in MI and MM09. Correlations lost significance once corrected for age.

COMPARISON WITH EXISTING METHODS

MRI scanners/protocols can be used to collect 1H-MRS and DWI data in pigs. Pigs have a larger, more complex, gyrencephalic brain than laboratory rodents but are less complex than non-human primates, thus satisfying the "replacement" principle of animal research.

CONCLUSIONS

Longitudinal effects in MRS measurements were similar to those reported in adolescent humans. MRS changes correlated with diffusion measurements indicating ongoing WM myelination/maturation.

Compressed sensing for high-resolution nonlipid suppressed 1 H FID MRSI of the human brain at 9.4T

PURPOSE

The aim of this study was to apply compressed sensing to accelerate the acquisition of high resolution metabolite maps of the human brain using a nonlipid suppressed ultra-short TR and TE ¹ H FID MRSI sequence at 9.4T.

METHODS

X-t sparse compressed sensing reconstruction was optimized for nonlipid suppressed ¹ H FID MRSI data. Coil-by-coil x-t sparse reconstruction was compared with SENSE x-t sparse and low rank reconstruction. The effect of matrix size and spatial resolution on the achievable acceleration factor was studied. Finally, in vivo metabolite maps with different acceleration factors of 2, 4, 5, and 10 were acquired and compared.

RESULTS

Coil-by-coil x-t sparse compressed sensing reconstruction was not able to reliably recover the nonlipid suppressed data, rather a combination of parallel and sparse reconstruction was necessary (SENSE x-t sparse). For acceleration factors of up to 5, both the low-rank and the compressed sensing methods were able to reconstruct the data comparably well (root mean squared errors [RMSEs] ≤ 10.5% for Cre). However, the reconstruction time of the low rank algorithm was drastically longer than compressed sensing. Using the optimized compressed sensing reconstruction, acceleration factors of 4 or 5 could be reached for the MRSI data with a matrix size of 64 × 64. For lower spatial resolutions, an acceleration factor of up to R∼4 was successfully achieved.

CONCLUSION

By tailoring the reconstruction scheme to the nonlipid suppressed data through parameter optimization and performance evaluation, we present high resolution (97 µL voxel size) accelerated in vivo metabolite maps of the human brain acquired at 9.4T within scan times of 3 to 3.75 min.

Pulse sequences and parallel imaging for high spatiotemporal resolution MRI at ultra-high field

The SNR and CNR benefits of ultra-high field (UHF) have helped push the envelope of achievable spatial resolution in MRI. For applications based on susceptibility contrast where there is a large CNR gain, high quality sub-millimeter resolution imaging is now being routinely performed, particularly in fMRI and phase imaging/QSM. This has enabled the study of structure and function of very fine-scale structures in the brain. UHF has also helped push the spatial resolution of many other MRI applications as will be outlined in this review. However, this push in resolution comes at a cost of a large encoding burden leading to very lengthy scans. Developments in parallel imaging with controlled aliasing and the move away from 2D slice-by-slice imaging to much more SNR-efficient simultaneous multi-slice (SMS) and 3D acquisitions have helped address this issue. In particular, these developments have revolutionized the efficiency of UHF MRI to enable high spatiotemporal resolution imaging at an order of magnitude faster acquisition. In addition to describing the main approaches to these techniques, this review will also outline important key practical considerations in using these methods in practice. Furthermore, new RF pulse design to tackle the B1+ and SAR issues of UHF and the increased SAR and power requirement of SMS RF pulses will also be touched upon. Finally, an outlook into new developments of smart encoding in more dimensions, particularly through using better temporal/across-contrast encoding and reconstruction will be described. Just as controlled aliasing fully exploits spatial encoding in parallel imaging to provide large multiplicative gains in accelerations, the complimentary use of these new approaches in temporal and across-contrast encoding are expected to provide exciting opportunities for further large gains in efficiency to further push the spatiotemporal resolution of MRI.

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.

A minimum-phase Shinnar-Le Roux spectral-spatial excitation RF pulse for simultaneous water and lipid suppression in 1H-MRSI of body extremities

PURPOSE

To develop a spectral-spatial (SPSP) excitation RF pulse for simultaneous water and lipid suppression in proton (1H) magnetic resonance spectroscopic imaging (MRSI) of body extremities.

METHODS

An SPSP excitation pulse is designed to excite Creatine (Cr) and Choline (Cho) metabolite signals while suppressing the overwhelming water and lipid signals. The SPSP pulse is designed using a recently proposed multidimensional Shinnar-Le Roux (SLR) RF pulse design method. A minimum-phase spectral selectivity profile is used to minimize signal loss from T2⁎ decay.

RESULTS

The performance of the SPSP pulse is evaluated via Bloch equation simulations and phantom experiments. The feasibility of the proposed method is demonstrated using three-dimensional, short repetition-time, free induction decay-based 1H-MRSI in the thigh muscle at 3T.

CONCLUSION

The proposed SPSP excitation pulse is useful for simultaneous water and lipid suppression. The proposed method enables new applications of high-resolution 1H-MRSI in body extremities.

Density-weighted concentric rings k-space trajectory for 1 H magnetic resonance spectroscopic imaging at 7 T

It has been shown that density-weighted (DW) k-space sampling with spiral and conventional phase encoding trajectories reduces spatial side lobes in magnetic resonance spectroscopic imaging (MRSI). In this study, we propose a new concentric ring trajectory (CRT) for DW-MRSI that samples k-space with a density that is proportional to a spatial, isotropic Hanning window. The properties of two different DW-CRTs were compared against a radially equidistant (RE) CRT and an echo-planar spectroscopic imaging (EPSI) trajectory in simulations, phantoms and in vivo experiments. These experiments, conducted at 7 T with a fixed nominal voxel size and matched acquisition times, revealed that the two DW-CRT designs improved the shape of the spatial response function by suppressing side lobes, also resulting in improved signal-to-noise ratio (SNR). High-quality spectra were acquired for all trajectories from a specific region of interest in the motor cortex with an in-plane resolution of 7.5 × 7.5 mm2 in 8 min 3 s. Due to hardware limitations, high-spatial-resolution spectra with an in-plane resolution of 5 × 5 mm2 and an acquisition time of 12 min 48 s were acquired only for the RE and one of the DW-CRT trajectories and not for EPSI. For all phantom and in vivo experiments, DW-CRTs resulted in the highest SNR. The achieved in vivo spectral quality of the DW-CRT method allowed for reliable metabolic mapping of eight metabolites including N-acetylaspartylglutamate, γ-aminobutyric acid and glutathione with Cramér-Rao lower bounds below 50%, using an LCModel analysis. Finally, high-quality metabolic mapping of a whole brain slice using DW-CRT was achieved with a high in-plane resolution of 5 × 5 mm2 in a healthy subject. These findings demonstrate that our DW-CRT MRSI technique can perform robustly on MRI systems and within a clinically feasible acquisition time.

Simultaneous QSM and metabolic imaging of the brain using SPICE

PURPOSE

To map brain metabolites and tissue magnetic susceptibility simultaneously using a single three-dimensional 1 H-MRSI acquisition without water suppression.

METHODS

The proposed technique builds on a subspace imaging method called spectroscopic imaging by exploiting spatiospectral correlation (SPICE), which enables ultrashort echo time (TE)/short pulse repetition time (TR) acquisitions for 1 H-MRSI without water suppression. This data acquisition scheme simultaneously captures both the spectral information of brain metabolites and the phase information of the water signals that is directly related to tissue magnetic susceptibility variations. In extending this scheme for simultaneous QSM and metabolic imaging, we increase k-space coverage by using dual density sparse sampling and ramp sampling to achieve spatial resolution often required by QSM, while maintaining a reasonable signal-to-noise ratio (SNR) for the spatiospectral data used for metabolite mapping. In data processing, we obtain high-quality QSM from the unsuppressed water signals by taking advantage of the larger number of echoes acquired and any available anatomical priors; metabolite spatiospectral distributions are reconstructed using a union-of-subspaces model.

RESULTS

In vivo experimental results demonstrate that the proposed method can produce susceptibility maps at a resolution higher than 1.8 × 1.8 × 2.4 mm3 along with metabolite spatiospectral distributions at a nominal spatial resolution of 2.4 × 2.4 × 2.4 mm3 from a single 7-min MRSI scan. The estimated susceptibility values are consistent with those obtained using the conventional QSM method with 3D multi-echo gradient echo acquisitions.

CONCLUSION

This article reports a new capability for simultaneous susceptibility mapping and metabolic imaging of the brain from a single 1 H-MRSI scan, which has potential for a wide range of applications. Magn Reson Med 79:13-21, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Predictive Systems Toxicology

In this review we address to what extent computational techniques can augment our ability to predict toxicity. The first section provides a brief history of empirical observations on toxicity dating back to the dawn of Sumerian civilization. Interestingly, the concept of dose emerged very early on, leading up to the modern emphasis on kinetic properties, which in turn encodes the insight that toxicity is not solely a property of a compound but instead depends on the interaction with the host organism. The next logical step is the current conception of evaluating drugs from a personalized medicine point of view. We review recent work on integrating what could be referred to as classical pharmacokinetic analysis with emerging systems biology approaches incorporating multiple omics data. These systems approaches employ advanced statistical analytical data processing complemented with machine learning techniques and use both pharmacokinetic and omics data. We find that such integrated approaches not only provide improved predictions of toxicity but also enable mechanistic interpretations of the molecular mechanisms underpinning toxicity and drug resistance. We conclude the chapter by discussing some of the main challenges, such as how to balance the inherent tension between the predicitive capacity of models, which in practice amounts to constraining the number of features in the models versus allowing for rich mechanistic interpretability, i.e., equipping models with numerous molecular features. This challenge also requires patient-specific predictions on toxicity, which in turn requires proper stratification of patients as regards how they respond, with or without adverse toxic effects. In summary, the transformation of the ancient concept of dose is currently successfully operationalized using rich integrative data encoded in patient-specific models.

A Subspace Approach to Spectral Quantification for MR Spectroscopic Imaging

OBJECTIVE

To provide a new approach to spectral quantification for magnetic resonance spectroscopic imaging (MRSI), incorporating both spatial and spectral priors.

METHODS

A novel signal model is proposed, which represents the spectral distributions of each molecule as a subspace and the entire spectrum as a union of subspaces. Based on this model, the spectral quantification can be solved in two steps: 1) subspace estimation based on the empirical distributions of the spectral parameters estimated using spectral priors; and 2) parameter estimation for the union-of-subspaces model incorporating spatial priors.

RESULTS

The proposed method has been evaluated using both simulated and experimental data, producing impressive results.

CONCLUSION

The proposed union-of-subspaces representation of spatiospectral functions provides an effective computational framework for solving the MRSI spectral quantification problem with spatiospectral constraints.

SIGNIFICANCE

The proposed approach transforms how the MRSI spectral quantification problem is solved and enables efficient and effective use of spatiospectral priors to improve parameter estimation. The resulting algorithm is expected to be useful for a wide range of quantitative metabolic imaging studies using MRSI.

Non-water-suppressed short-echo-time magnetic resonance spectroscopic imaging using a concentric ring k-space trajectory

Water-suppressed MRS acquisition techniques have been the standard MRS approach used in research and for clinical scanning to date. The acquisition of a non-water-suppressed MRS spectrum is used for artefact correction, reconstruction of phased-array coil data and metabolite quantification. Here, a two-scan metabolite-cycling magnetic resonance spectroscopic imaging (MRSI) scheme that does not use water suppression is demonstrated and evaluated. Specifically, the feasibility of acquiring and quantifying short-echo (TE  = 14 ms), two-dimensional stimulated echo acquisition mode (STEAM) MRSI spectra in the motor cortex is demonstrated on a 3 T MRI system. The increase in measurement time from the metabolite-cycling is counterbalanced by a time-efficient concentric ring k-space trajectory. To validate the technique, water-suppressed MRSI acquisitions were also performed for comparison. The proposed non-water-suppressed metabolite-cycling MRSI technique was tested for detection and correction of resonance frequency drifts due to subject motion and/or hardware instability, and the feasibility of high-resolution metabolic mapping over a whole brain slice was assessed. Our results show that the metabolite spectra and estimated concentrations are in agreement between non-water-suppressed and water-suppressed techniques. The achieved spectral quality, signal-to-noise ratio (SNR) > 20 and linewidth <7 Hz allowed reliable metabolic mapping of five major brain metabolites in the motor cortex with an in-plane resolution of 10 × 10 mm2 in 8 min and with a Cramér-Rao lower bound of less than 20% using LCModel analysis. In addition, the high SNR of the water peak of the non-water-suppressed technique enabled voxel-wise single-scan frequency, phase and eddy current correction. These findings demonstrate that our non-water-suppressed metabolite-cycling MRSI technique can perform robustly on 3 T MRI systems and within a clinically feasible acquisition time.

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

PURPOSE

To develop a rapid ³¹ P-MRSI method with high spatiospectral resolution using low-rank tensor-based data acquisition and image reconstruction.

METHODS

The multidimensional image function of ³¹ P-MRSI is represented by a low-rank tensor to capture the spatial-spectral-temporal correlations of data. A hybrid data acquisition scheme is used for sparse sampling, which consists of a set of "training" data with limited k-space coverage to capture the subspace structure of the image function, and a set of sparsely sampled "imaging" data for high-resolution image reconstruction. An explicit subspace pursuit approach is used for image reconstruction, which estimates the bases of the subspace from the "training" data and then reconstructs a high-resolution image function from the "imaging" data.

RESULTS

We have validated the feasibility of the proposed method using phantom and in vivo studies on a 3T whole-body scanner and a 9.4T preclinical scanner. The proposed method produced high-resolution static ³¹ P-MRSI images (i.e., 6.9 × 6.9 × 10 mm³ nominal resolution in a 15-min acquisition at 3T) and high-resolution, high-frame-rate dynamic ³¹ P-MRSI images (i.e., 1.5 × 1.5 × 1.6 mm³ nominal resolution, 30 s/frame at 9.4T).

CONCLUSIONS

Dynamic spatiospectral variations of ³¹ P-MRSI signals can be efficiently represented by a low-rank tensor. Exploiting this mathematical structure for data acquisition and image reconstruction can lead to fast ³¹ P-MRSI with high resolution, frame-rate, and SNR. Magn Reson Med 78:419-428, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

High-resolution 1 H-MRSI of the brain using short-TE SPICE

PURPOSE

To improve signal-to-noise ratio (SNR) for high-resolution spectroscopic imaging using a subspace-based technique known as SPectroscopic Imaging by exploiting spatiospectral CorrElation (SPICE).

METHODS

The proposed method is based on a union-of-subspaces model of MRSI signals, which exploits the partial separability properties of water, lipid, baseline and metabolite signals. Enabled by this model, a special scheme is used for accelerated data acquisition, which includes a double-echo CSI component used to collect a "training" dataset (for determination of the basis functions) and a short-TE EPSI component used to collect a sparse "imaging" dataset (for determination of the overall spatiospectral distributions). A set of signal processing algorithms are developed to remove the water and lipid signals and jointly reconstruct the metabolite and baseline signals.

RESULTS

In vivo 1 H-MRSI results show that the proposed method can effectively remove the remaining water and lipid signals from sparse MRSI data acquired at 20 ms TE. Spatiospectral distributions of metabolite signals at 2 mm in-plane resolution with good SNR were obtained in a 15.5 min scan.

CONCLUSIONS

The proposed method can effectively remove nuisance signals and reconstruct high-resolution spatiospectral functions from sparse data to make short-TE SPICE possible. The method should prove useful for high-resolution 1 H-MRSI of the brain. Magn Reson Med 77:467-479, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Imaging based magnetic resonance spectroscopy (MRS) localization for quantitative neurochemical analysis and cerebral metabolism studies

Accurate quantitative metabolic imaging of the brain presents significant challenges due to the complexity and heterogeneity of its structures and compositions with distinct compartmentations of brain tissue types (e.g., gray and white matter). The brain is compartmentalized into various regions based on their unique functions and locations. In vivo magnetic resonance spectroscopy (MRS) techniques allow non-invasive measurements of neurochemicals in either single voxel or multiple voxels, yet the spatial resolution and detection sensitivity of MRS are significantly lower compared with MRI. A fundamentally different approach, namely spectral localization by imaging (SLIM) provides a new framework that overcomes major limitations of conventional MRS techniques. Conventional MRS allows only rectangular voxel shapes that do not conform to the shapes of brain structures or lesions, while SLIM allows compartments with arbitrary shapes. However, the restrictive assumption proposed in the original concept of SLIM, i.e., compartmental homogeneity, led to spectral localization errors, which have limited its broad applications. This review focuses on the recent technical frontiers of image-based MRS localization techniques that overcome the limitations of SLIM through the development and implementation of various new strategies, including incorporation of magnetic field inhomogeneity corrections, the use of multiple receiver coils, and prospective optimization of data acquisition.

Compartmentalized low-rank recovery for high-resolution lipid unsuppressed MRSI

PURPOSE

To introduce a novel algorithm for the recovery of high-resolution magnetic resonance spectroscopic imaging (MRSI) data with minimal lipid leakage artifacts, from dual-density spiral acquisition.

METHODS

The reconstruction of MRSI data from dual-density spiral data is formulated as a compartmental low-rank recovery problem. The MRSI dataset is modeled as the sum of metabolite and lipid signals, each of which is support limited to the brain and extracranial regions, respectively, in addition to being orthogonal to each other. The reconstruction problem is formulated as an optimization problem, which is solved using iterative reweighted nuclear norm minimization.

RESULTS

The comparisons of the scheme against dual-resolution reconstruction algorithm on numerical phantom and in vivo datasets demonstrate the ability of the scheme to provide higher spatial resolution and lower lipid leakage artifacts. The experiments demonstrate the ability of the scheme to recover the metabolite maps, from lipid unsuppressed datasets with echo time (TE) = 55 ms.

CONCLUSION

The proposed reconstruction method and data acquisition strategy provide an efficient way to achieve high-resolution metabolite maps without lipid suppression. This algorithm would be beneficial for fast metabolic mapping and extension to multislice acquisitions. Magn Reson Med 78:1267-1280, 2017. © 2016 International Society for Magnetic Resonance in Medicine.