Joint k-ω Space Image Reconstruction and Data Fitting for Chemical Exchange Saturation Transfer Magnetic Resonance Imaging

Chemical exchange saturation transfer (CEST) magnetic resonance imaging (MRI) is a novel MRI technology to image certain compounds at extremely low concentrations. Long acquisition time to measure signals at a set of offset frequencies of the Z-spectra and to repeat measurements to reduce noise pose significant challenges to its applications. This study explores correlations of CEST MR images along the spatial and Z-spectral dimensions to improve MR image quality and robustness of magnetization transfer ratio (MTR) asymmetry estimation via a joint k-ω reconstruction model. The model was formulated as an optimization problem with respect to MR images at all frequencies ω, while incorporating regularizations along the spatial and spectral dimensions. The solution was subject to a self-consistency condition that the Z-spectrum of each pixel follows a multi-peak data fitting model corresponding to different CEST pools. The optimization problem was solved using the alternating direction method of multipliers. The proposed joint reconstruction method was evaluated on a simulated CEST MRI phantom and semi-experimentally on choline and iopamidol phantoms with added Gaussian noise of various levels. Results demonstrated that the joint reconstruction method was more tolerable to noise and reduction in number of offset frequencies by improving signal-to-noise ratio (SNR) of the reconstructed images and reducing uncertainty in MTR asymmetry estimation. In the choline and iopamidol phantom cases with 10.5% noise in the measurement data, our method achieved an averaged SNR of 31.0 dB and 32.2 dB compared to the SNR of 24.7 dB and 24.4 dB in the conventional reconstruction approach. It reduced uncertainty of the MTR asymmetry estimation over all regions of interest by 54.4% and 43.7%, from 1.71 and 2.38 to 0.78 and 1.71, respectively.

Learning to deep learning: statistics and a paradigm test in selecting a UNet architecture to enhance MRI

OBJECTIVE

This study aims to assess the statistical significance of training parameters in 240 dense UNets (DUNets) used for enhancing low Signal-to-Noise Ratio (SNR) and undersampled MRI in various acquisition protocols. The objective is to determine the validity of differences between different DUNet configurations and their impact on image quality metrics.

MATERIALS AND METHODS

To achieve this, we trained all DUNets using the same learning rate and number of epochs, with variations in 5 acquisition protocols, 24 loss function weightings, and 2 ground truths. We calculated evaluation metrics for two metric regions of interest (ROI). We employed both Analysis of Variance (ANOVA) and Mixed Effects Model (MEM) to assess the statistical significance of the independent parameters, aiming to compare their efficacy in revealing differences and interactions among fixed parameters.

RESULTS

ANOVA analysis showed that, except for the acquisition protocol, fixed variables were statistically insignificant. In contrast, MEM analysis revealed that all fixed parameters and their interactions held statistical significance. This emphasizes the need for advanced statistical analysis in comparative studies, where MEM can uncover finer distinctions often overlooked by ANOVA.

DISCUSSION

These findings highlight the importance of utilizing appropriate statistical analysis when comparing different deep learning models. Additionally, the surprising effectiveness of the UNet architecture in enhancing various acquisition protocols underscores the potential for developing improved methods for characterizing and training deep learning models. This study serves as a stepping stone toward enhancing the transparency and comparability of deep learning techniques for medical imaging applications.

Sidebands in CEST MR-How to recognize and avoid them

PURPOSE

Clinical scanners require pulsed CEST sequences to maintain amplifier and specific absorption rate limits. During off-resonant RF irradiation and interpulse delay, the magnetization can accumulate specific relative phases within the pulse train. In this work, we show that these phases are important to consider, as they can lead to unexpected artifacts when no interpulse gradient spoiling is performed during the saturation train.

METHODS

We investigated sideband artifacts using a CEST-3D snapshot gradient-echo sequence at 3 T. Initially, Bloch-McConnell simulations were carried out with Pulseq-CEST, while measurements were performed in vitro and in vivo.

RESULTS

Sidebands can be hidden in Z-spectra, and their structure becomes clearly visible only at high sampling. Sidebands are further influenced by B0 inhomogeneities and the RF phase cycling within the pulse train. In vivo, sidebands are mostly visible in liquid compartments such as CSF. Multi-pulse sidebands can be suppressed by interpulse gradient spoiling.

CONCLUSION

We provide new insights into sidebands occurring in pulsed CEST experiments and show that, similar as in imaging sequences, gradient and RF spoiling play an important role. Gradient spoiling avoids misinterpretations of sidebands as CEST effects especially in liquid environments including pathological tissue or for CEST resonances close to water. It is recommended to simulate pulsed CEST sequences in advance to avoid artifacts.

Machine learning-based amide proton transfer imaging using partially synthetic training data

PURPOSE

Machine learning (ML) has been increasingly used to quantify CEST effect. ML models are typically trained using either measured data or fully simulated data. However, training with measured data often lacks sufficient training data, whereas training with fully simulated data may introduce bias because of limited simulations pools. This study introduces a new platform that combines simulated and measured components to generate partially synthetic CEST data, and to evaluate its feasibility for training ML models to predict amide proton transfer (APT) effect.

METHODS

Partially synthetic CEST signals were created using an inverse summation of APT effects from simulations and the other components from measurements. Training data were generated by varying APT simulation parameters and applying scaling factors to adjust the measured components, achieving a balance between simulation flexibility and fidelity. First, tissue-mimicking CEST signals along with ground truth information were created using multiple-pool model simulations to validate this method. Second, an ML model was trained individually on partially synthetic data, in vivo data, and fully simulated data, to predict APT effect in rat brains bearing 9 L tumors.

RESULTS

Experiments on tissue-mimicking data suggest that the ML method using the partially synthetic data is accurate in predicting APT. In vivo experiments suggest that our method provides more accurate and robust prediction than the training using in vivo data and fully synthetic data.

CONCLUSION

Partially synthetic CEST data can address the challenges in conventional ML methods.

Accelerating CEST imaging using a model-based deep neural network with synthetic training data

PURPOSE

To develop a model-based deep neural network for high-quality image reconstruction of undersampled multi-coil CEST data.

THEORY AND METHODS

Inspired by the variational network (VN), the CEST image reconstruction equation is unrolled into a deep neural network (CEST-VN) with a k-space data-sharing block that takes advantage of the inherent redundancy in adjacent CEST frames and 3D spatial-frequential convolution kernels that exploit correlations in the x-ω domain. Additionally, a new pipeline based on multiple-pool Bloch-McConnell simulations is devised to synthesize multi-coil CEST data from publicly available anatomical MRI data. The proposed network is trained on simulated data with a CEST-specific loss function that jointly measures the structural and CEST contrast. The performance of CEST-VN was evaluated on four healthy volunteers and five brain tumor patients using retrospectively or prospectively undersampled data with various acceleration factors, and then compared with other conventional and state-of-the-art reconstruction methods.

RESULTS

The proposed CEST-VN method generated high-quality CEST source images and amide proton transfer-weighted maps in healthy and brain tumor subjects, consistently outperforming GRAPPA, blind compressed sensing, and the original VN. With the acceleration factors increasing from 3 to 6, CEST-VN with the same hyperparameters yielded similar and accurate reconstruction without apparent loss of details or increase of artifacts. The ablation studies confirmed the effectiveness of the CEST-specific loss function and data-sharing block used.

CONCLUSIONS

The proposed CEST-VN method can offer high-quality CEST source images and amide proton transfer-weighted maps from highly undersampled multi-coil data by integrating the deep learning prior and multi-coil sensitivity encoding model.

An end-to-end LSTM-Attention based framework for quasi-steady-state CEST prediction

Chemical exchange saturation transfer (CEST)-magnetic resonance imaging (MRI) often takes prolonged saturation duration (Ts) and relaxation delay (Td) to reach the steady state, and yet the insufficiently long Ts and Td in actual experiments may underestimate the CEST measurement. In this study, we aimed to develop a deep learning-based model for quasi-steady-state (QUASS) prediction from non-steady-state CEST acquired in experiments, therefore overcoming the limitation of the CEST effect which needs prolonged saturation time to reach a steady state. To support network training, a multi-pool Bloch-McConnell equation was designed to derive wide-ranging simulated Z-spectra, so as to solve the problem of time and labor consumption in manual annotation work. Following this, we formulated a hybrid architecture of long short-term memory (LSTM)-Attention to improve the predictive ability. The multilayer perceptron, recurrent neural network, LSTM, gated recurrent unit, BiLSTM, and LSTM-Attention were included in comparative experiments of QUASS CEST prediction, and the best performance was obtained by the proposed LSTM-Attention model. In terms of the linear regression analysis, structural similarity index (SSIM), peak signal-to-noise ratio (PSNR), and mean-square error (MSE), the results of LSTM-Attention demonstrate that the coefficient of determination in the linear regression analysis was at least R2 = 0.9748 for six different representative frequency offsets, the mean values of prediction accuracies in terms of SSIM, PSNR and MSE were 0.9991, 49.6714, and 1.68 × 10-4 for all frequency offsets. It was concluded that the LSTM-Attention model enabled high-quality QUASS CEST prediction.

An uncertainty aided framework for learning based liverT1ρmapping and analysis

Objective. QuantitativeT1ρimaging has potential for assessment of biochemical alterations of liver pathologies. Deep learning methods have been employed to accelerate quantitativeT1ρimaging. To employ artificial intelligence-based quantitative imaging methods in complicated clinical environment, it is valuable to estimate the uncertainty of the predicatedT1ρvalues to provide the confidence level of the quantification results. The uncertainty should also be utilized to aid the post-hoc quantitative analysis and model learning tasks.Approach. To address this need, we propose a parametric map refinement approach for learning-basedT1ρmapping and train the model in a probabilistic way to model the uncertainty. We also propose to utilize the uncertainty map to spatially weight the training of an improvedT1ρmapping network to further improve the mapping performance and to remove pixels with unreliableT1ρvalues in the region of interest. The framework was tested on a dataset of 51 patients with different liver fibrosis stages.Main results. Our results indicate that the learning-based map refinement method leads to a relative mapping error of less than 3% and provides uncertainty estimation simultaneously. The estimated uncertainty reflects the actual error level, and it can be used to further reduce relativeT1ρmapping error to 2.60% as well as removing unreliable pixels in the region of interest effectively.Significance. Our studies demonstrate the proposed approach has potential to provide a learning-based quantitative MRI system for trustworthyT1ρmapping of the liver.

Unified Bayesian network for uncertainty quantification of physiological parameters in dynamic contrast enhanced (DCE) MRI of the liver

Objective. Physiological parameter estimation is affected by intrinsic ambiguity in the data such as noise and model inaccuracies. The aim of this work is to provide a deep learning framework for accurate parameter and uncertainty estimates for DCE-MRI in the liver.Approach. Concentration time curves are simulated to train a Bayesian neural network (BNN). Training of the BNN involves minimization of a loss function that jointly minimizes the aleatoric and epistemic uncertainties. Uncertainty estimation is evaluated for different noise levels and for different out of distribution (OD) cases, i.e. where the data during inference differs strongly to the data during training. The accuracy of parameter estimates are compared to a nonlinear least squares (NLLS) fitting in numerical simulations andin vivodata of a patient suffering from hepatic tumor lesions.Main results. BNN achieved lower root-mean-squared-errors (RMSE) than the NLLS for the simulated data. RMSE of BNN was on overage of all noise levels lower by 33% ± 1.9% forktrans, 22% ± 6% forveand 89% ± 5% forvpthan the NLLS. The aleatoric uncertainties of the parameters increased with increasing noise level, whereas the epistemic uncertainty increased when a BNN was evaluated with OD data. For thein vivodata, more robust parameter estimations were obtained by the BNN than the NLLS fit. In addition, the differences between estimated parameters for healthy and tumor regions-of-interest were significant (p< 0.0001).Significance. The proposed framework allowed for accurate parameter estimates for quantitative DCE-MRI. In addition, the BNN provided uncertainty estimates which highlighted cases of high noise and in which the training data did not match the data during inference. This is important for clinical application because it would indicate cases in which the trained model is inadequate and additional training with an adapted training data set is required.

Data-driven electrical conductivity brain imaging using 3 T MRI

Magnetic resonance electrical properties tomography (MR-EPT) is a non-invasive measurement technique that derives the electrical properties (EPs, e.g., conductivity or permittivity) of tissues in the radiofrequency range (64 MHz for 1.5 T and 128 MHz for 3 T MR systems). Clinical studies have shown the potential of tissue conductivity as a biomarker. To date, model-based conductivity reconstructions rely on numerical assumptions and approximations, leading to inaccuracies in the reconstructed maps. To address such limitations, we propose an artificial neural network (ANN)-based non-linear conductivity estimator trained on simulated data for conductivity brain imaging. Network training was performed on 201 synthesized T2-weighted spin-echo (SE) data obtained from the finite-difference time-domain (FDTD) electromagnetic (EM) simulation. The dataset was composed of an approximated T2-w SE magnitude and transceive phase information. The proposed method was tested three in-silico and in-vivo on two volunteers and three patients' data. For comparison purposes, various conventional phase-based EPT reconstruction methods were used that ignoreB 1 + magnitude information, such as Savitzky-Golay kernel combined with Gaussian filter (S-G Kernel), phase-based convection-reaction EPT (cr-EPT), magnitude-weighted polynomial-fitting phase-based EPT (Poly-Fit), and integral-based phase-based EPT (Integral-based). From the in-silico experiments, quantitative analysis showed that the proposed method provides more accurate and improved quality (e.g., high structural preservation) conductivity maps compared to conventional reconstruction methods. Representatively, in the healthy brain in-silico phantom experiment, the proposed method yielded mean conductivity values of 1.97 ± 0.20 S/m for CSF, 0.33 ± 0.04 S/m for WM, and 0.52 ± 0.08 S/m for GM, which were closer to the ground-truth conductivity (2.00, 0.30, 0.50 S/m) than the integral-based method (2.56 ± 2.31, 0.39 ± 0.12, 0.68 ± 0.33 S/m). In-vivo ANN-based conductivity reconstructions were also of improved quality compared to conventional reconstructions and demonstrated network generalizability and robustness to in-vivo data and pathologies. The reported in-vivo brain conductivity values were in agreement with literatures. In addition, the proposed method was observed for various SNR levels (SNR levels = 10, 20, 40, and 58) and repeatability conditions (the eight acquisitions with the number of signal averages = 1). The preliminary investigations on brain tumor patient datasets suggest that the network trained on simulated dataset can generalize to unforeseen in-vivo pathologies, thus demonstrating its potential for clinical applications.

Deep learning-based Lorentzian fitting of water saturation shift referencing spectra in MRI

PURPOSE

Water saturation shift referencing (WASSR) Z-spectra are used commonly for field referencing in chemical exchange saturation transfer (CEST) MRI. However, their analysis using least-squares (LS) Lorentzian fitting is time-consuming and prone to errors because of the unavoidable noise in vivo. A deep learning-based single Lorentzian Fitting Network (sLoFNet) is proposed to overcome these shortcomings.

METHODS

A neural network architecture was constructed and its hyperparameters optimized. Training was conducted on a simulated and in vivo-paired data sets of discrete signal values and their corresponding Lorentzian shape parameters. The sLoFNet performance was compared with LS on several WASSR data sets (both simulated and in vivo 3T brain scans). Prediction errors, robustness against noise, effects of sampling density, and time consumption were compared.

RESULTS

LS and sLoFNet performed comparably in terms of RMS error and mean absolute error on all in vivo data with no statistically significant difference. Although the LS method fitted well on samples with low noise, its error increased rapidly when increasing sample noise up to 4.5%, whereas the error of sLoFNet increased only marginally. With the reduction of Z-spectral sampling density, prediction errors increased for both methods, but the increase occurred earlier (at 25 vs. 15 frequency points) and was more pronounced for LS. Furthermore, sLoFNet performed, on average, 70 times faster than the LS-method.

CONCLUSION

Comparisons between LS and sLoFNet on simulated and in vivo WASSR MRI Z-spectra in terms of robustness against noise and decreased sample resolution, as well as time consumption, showed significant advantages for sLoFNet.

Acquisition sequences and reconstruction methods for fast chemical exchange saturation transfer imaging

Chemical exchange saturation transfer (CEST) imaging is an emerging molecular magnetic resonance imaging (MRI) technique that has been developed and employed in numerous diseases. Based on the unique saturation transfer principle, a family of CEST-detectable biomolecules in vivo have been found capable of providing valuable diagnostic information. However, CEST MRI needs a relatively long scan time due to the common long saturation labeling module and typical acquisition of multiple frequency offsets and signal averages, limiting its widespread clinical applications. So far, a plethora of imaging schemes and techniques has been developed to accelerate CEST MRI. In this review, the key acquisition and reconstruction methods for fast CEST imaging are summarized from a practical and systematic point of view. The first acquisition sequence section describes the major development of saturation schemes, readout patterns, ultrafast z-spectroscopy, and saturation-editing techniques for rapid CEST imaging. The second reconstruction method section lists the important advances of parallel imaging, compressed sensing, sparsity in the z-spectrum, and algorithms beyond the Fourier transform for speeding up CEST MRI.

Frequency importance analysis for chemical exchange saturation transfer magnetic resonance imaging using permuted random forest

Chemical exchange saturation transfer magnetic resonance imaging (CEST MRI) is a promising molecular imaging tool that allows sensitive detection of endogenous metabolic changes. However, because the CEST spectrum does not display a clear peak like MR spectroscopy, its signal interpretation is challenging, especially under 3-T field strength or with a large saturation B₁ . Herein, as an alternative to conventional Z-spectral fitting approaches, a permuted random forest (PRF) method is developed to determine featured saturation frequencies for lesion identification, so-called CEST frequency importance analysis. Briefly, voxels in the CEST dataset were labeled as lesion and control according to multicontrast MR images. Then, by considering each voxel's saturation signal series as a sample, a permutation importance algorithm was employed to rank the contribution of saturation frequency offsets in the differentiation of lesion and normal tissue. Simulations demonstrated that PRF could correctly determine the frequency offsets (3.5 or -3.5 ppm) for classifying two groups of Z-spectra, under a range of B₀ , B₁ conditions and sample sizes. For ischemic rat brains, PRF only displayed high feature importance around amide frequency at 2 h postischemia, reflecting that the pH changes occurred at an early stage. By contrast, the data acquired at 24 h postischemia exhibited high feature importance at multiple frequencies (amide, water, and lipids), which suggested the complex tissue changes that occur during the later stages. Finally, PRF was assessed using 3-T CEST data from four brain tumor patients. By defining the tumor region on amide proton transfer-weighted images, PRF analysis identified different CEST frequency importance for two types of tumors (glioblastoma and metastatic tumor) (p < 0.05, with each image slice as a subject). In conclusion, the PRF method was able to rank and interpret the contribution of all acquired saturation offsets to lesion identification; this may facilitate CEST analysis in clinical applications, and open up new doors for comprehensive CEST analysis tools other than model-based approaches.

Chemical exchange saturation transfer imaging of creatine, phosphocreatine, and protein arginine residue in tissues

Chemical exchange saturation transfer (CEST) MRI has become a promising technique to assay target proteins and metabolites through their exchangeable protons, noninvasively. The ubiquity of creatine (Cr) and phosphocreatine (PCr) due to their pivotal roles in energy homeostasis through the creatine phosphate pathway has made them prime targets for CEST in the diagnosis and monitoring of disease pathologies, particularly in tissues heavily dependent on the maintenance of rich energy reserves. Guanidinium CEST from protein arginine residues (i.e. arginine CEST) can also provide information about the protein profile in tissue. However, numerous obfuscating factors stand as obstacles to the specificity of arginine, Cr, and PCr imaging through CEST, such as semisolid magnetization transfer, fast chemical exchanges such as primary amines, and the effects of nuclear Overhauser enhancement from aromatic and amide protons. In this review, the specific exchange properties of protein arginine residues, Cr, and PCr, along with their validation, are discussed, including the considerations necessary to target and tune their signal effects through CEST imaging. Additionally, strategies that have been employed to enhance the specificity of these exchanges in CEST imaging are described, along with how they have opened up possible applications of protein arginine residues, Cr and PCr CEST imaging in the study and diagnosis of pathology. A clear understanding of the capabilities and caveats of using CEST to image these vital metabolites and mitigation strategies is crucial to expanding the possibilities of this promising technology.

Glutamate-weighted CEST (gluCEST) imaging for mapping neurometabolism: An update on the state of the art and emerging findings from in vivo applications

Glutamate is the primary excitatory neurotransmitter in the mammalian central nervous system. As such, its proper regulation is essential to the healthy function of the human brain, and dysregulation of glutamate metabolism and compartmentalization underlies numerous neurological and neuropsychiatric pathologies. Glutamate-weighted chemical exchange saturation transfer (gluCEST) MRI is one of the only ways to non-invasively observe the relative concentration and spatial distribution of glutamate in the human brain. In the past 10 years, gluCEST has developed from a proof-of-concept experiment carried out in imaging phantoms and model systems to an increasingly sophisticated technique applied to reveal deviations from baseline neural metabolism in human beings, most notably in patients experiencing seizures of various origins or those on the psychosis spectrum. This article traces that progress, including in-depth discussion of the technical specifics of gluCEST and potential challenges to performing these experiments rigorously. We discuss the neurobiological context of glutamate, including the widely accepted hypotheses and models in the literature regarding its involvement in neurodegenerative diseases and other pathology. We then review the state of the art of in vivo glutamate detection by magnetic resonance imaging and the limitations on this front of in vivo MR spectroscopy. The gluCEST experiment is introduced and its advantages, challenges and limitations are thoroughly explored, beginning with the phantom experiment results demonstrated in the initial publication, through the latest approaches to correcting human brain images for B₁ inhomogeneity. We then give a comprehensive overview of preclinical applications demonstrated to date, including Alzheimer's disease, Parkinson's disease, Huntington's disease, Traumatic brain injury and cancer, followed by a similar discussion of human studies. Finally, we highlight emerging applications, and discuss technical improvements on the horizon that hold promise for improving the robustness and versatility of gluCEST and its increasing presence in the arena of translational and precision medicine.

Linear projection-based chemical exchange saturation transfer parameter estimation

Isolated evaluation of multiparametric in vivo chemical exchange saturation transfer (CEST) MRI often requires complex computational processing for both correction of B0 and B1 inhomogeneity and contrast generation. For that, sufficiently densely sampled Z-spectra need to be acquired. The list of acquired frequency offsets largely determines the total CEST acquisition time, while potentially representing redundant information. In this work, a linear projection-based multiparametric CEST evaluation method is introduced that offers fast B0 and B1 inhomogeneity correction, contrast generation and feature selection for CEST data, enabling reduction of the overall measurement time. To that end, CEST data acquired at 7 T in six healthy subjects and in one brain tumor patient were conventionally evaluated by interpolation-based inhomogeneity correction and Lorentzian curve fitting. Linear regression was used to obtain coefficient vectors that directly map uncorrected data to corrected Lorentzian target parameters. L1-regularization was applied to find subsets of the originally acquired CEST measurements that still allow for such a linear projection mapping. The linear projection method allows fast and interpretable mapping from acquired raw data to contrast parameters of interest, generalizing from healthy subject training data to unseen healthy test data and to the tumor patient dataset. The L1-regularization method shows that a fraction of the acquired CEST measurements is sufficient to preserve tissue contrasts, offering up to a 2.8-fold reduction of scan time. Similar observations as for the 7-T data can be made for data from a clinical 3-T scanner. Being a fast and interpretable computation step, the proposed method is complementary to neural networks that have recently been employed for similar purposes. The scan time acceleration offered by the L1-regularization ("CEST-LASSO") constitutes a step towards better applicability of multiparametric CEST protocols in a clinical context.

Accelerated and quantitative three-dimensional molecular MRI using a generative adversarial network

PURPOSE

To substantially shorten the acquisition time required for quantitative three-dimensional (3D) chemical exchange saturation transfer (CEST) and semisolid magnetization transfer (MT) imaging and allow for rapid chemical exchange parameter map reconstruction.

METHODS

Three-dimensional CEST and MT magnetic resonance fingerprinting (MRF) datasets of L-arginine phantoms, whole-brains, and calf muscles from healthy volunteers, cancer patients, and cardiac patients were acquired using 3T clinical scanners at three different sites, using three different scanner models and coils. A saturation transfer-oriented generative adversarial network (GAN-ST) supervised framework was then designed and trained to learn the mapping from a reduced input data space to the quantitative exchange parameter space, while preserving perceptual and quantitative content.

RESULTS

The GAN-ST 3D acquisition time was 42-52 s, 70% shorter than CEST-MRF. The quantitative reconstruction of the entire brain took 0.8 s. An excellent agreement was observed between the ground truth and GAN-based L-arginine concentration and pH values (Pearson's r > 0.95, ICC > 0.88, NRMSE < 3%). GAN-ST images from a brain-tumor subject yielded a semi-solid volume fraction and exchange rate NRMSE of3 . 8 ± 1 . 3 % $$ 3.8\pm 1.3\% $$ and4 . 6 ± 1 . 3 % $$ 4.6\pm 1.3\% $$ , respectively, and SSIM of96 . 3 ± 1 . 6 % $$ 96.3\pm 1.6\% $$ and95 . 0 ± 2 . 4 % $$ 95.0\pm 2.4\% $$ , respectively. The mapping of the calf-muscle exchange parameters in a cardiac patient, yielded NRMSE < 7% and SSIM > 94% for the semi-solid exchange parameters. In regions with large susceptibility artifacts, GAN-ST has demonstrated improved performance and reduced noise compared to MRF.

CONCLUSION

GAN-ST can substantially reduce the acquisition time for quantitative semi-solid MT/CEST mapping, while retaining performance even when facing pathologies and scanner models that were not available during training.

DeepCEST 7 T: Fast and homogeneous mapping of 7 T CEST MRI parameters and their uncertainty quantification

PURPOSE

In this work, we investigated the ability of neural networks to rapidly and robustly predict Lorentzian parameters of multi-pool CEST MRI spectra at 7 T with corresponding uncertainty maps to make them quickly and easily available for routine clinical use.

METHODS

We developed a deepCEST 7 T approach that generates CEST contrasts from just 1 scan with robustness against B₁ inhomogeneities. The input data for a neural feed-forward network consisted of 7 T in vivo uncorrected Z-spectra of a single B₁ level, and a B₁ map. The 7 T raw data were acquired using a 3D snapshot gradient echo multiple interleaved mode saturation CEST sequence. These inputs were mapped voxel-wise to target data consisting of Lorentzian amplitudes generated conventionally by 5-pool Lorentzian fitting of normalized, denoised, B₀ - and B₁ -corrected Z-spectra. The deepCEST network was trained with Gaussian negative log-likelihood loss, providing an uncertainty quantification in addition to the Lorentzian amplitudes.

RESULTS

The deepCEST 7 T network provides fast and accurate prediction of all Lorentzian parameters also when only a single B₁ level is used. The prediction was highly accurate with respect to the Lorentzian fit amplitudes, and both healthy tissues and hyperintensities in tumor areas are predicted with a low uncertainty. In corrupted cases, high uncertainty indicated wrong predictions reliably.

CONCLUSION

The proposed deepCEST 7 T approach reduces scan time by 50% to now 6:42 min, but still delivers both B₀ - and B₁ -corrected homogeneous CEST contrasts along with an uncertainty map, which can increase diagnostic confidence. Multiple accurate 7 T CEST contrasts are delivered within seconds.

calf - Software for CEST Analysis with Lorentzian Fitting

Analysis of chemical exchange saturation transfer (CEST) MRI data requires sophisticated methods to obtain reliable results about metabolites in the tissue under study. CEST generates z-spectra with multiple components, each originating from individual molecular groups. The individual lines with Lorentzian line shape are mostly overlapping and disturbed by various effects. We present an elaborate method based on an adaptive nonlinear least squares algorithm that provides robust quantification of z-spectra and incorporates prior knowledge in the fitting process. To disseminate CEST to the research community, we developed software as part of this study that runs on the Microsoft Windows operating system and will be made freely available to the community. Special attention has been paid to establish a low entrance threshold and high usability, so that even less experienced users can successfully analyze CEST data.

Reassembled saturation transfer (REST) MR images at 2 B1 values for in vivo exchange-dependent imaging of amide and nuclear Overhauser enhancement

PURPOSE

Design an efficient CEST scheme for exchange-dependent images with high contrast-to-noise ratio.

THEORY

Reassembled saturation transfer (REST) signals were defined as Δ $$ \Delta $$ r.Z = r.Z_ref - r.Z_CEST and the reassembled exchange-dependen magnetization transfer ratio r.MTR_Rex  = r.1/Z_ref - r.1/Z_CEST , utilizing the averages over loosely sampled reference frequency offsets as Z_ref and over densely sampled target offsets as Z_CEST . Using r.MTR_Rex measured under 2 B_1,sat values, exchange rate could be estimated.

METHODS

The REST approach was optimized and assessed quantitatively by simulations for various exchange rates, pool concentration, and water T₁ . In vivo evaluation was performed on ischemic rat brains at 7 Tesla and human brains at 3 Tesla, in comparison with conventional asymmetrical analysis, Lorentzian difference (LD), an MTR_{Rex_} LD.

RESULTS

For a broad choice of Δ ω ref $$ \Delta {\omega}_{ref} $$ ranges and numbers, Δr.Z and r.MTR_Rex exhibited comparable quantification features with conventional LD and MTR_Rex _LD, respectively, when B_1,sat  ≤ 1 μT. The subtraction of 2 REST values under distinct B_1,sat values showed linear relationships with exchange rate and obtained immunity to field inhomogeneity and variation in MT and water T₁ . For both rat and human studies, REST images exhibited similar contrast distribution to MTR_Rex _LD, with superiority in contrast-to-noise ratio and acquisition efficiency. Compared with MTR_Rex _LD, 2-B_1,sat subtraction REST images displayed better resistance to B₁ inhomogeneity, with more specific enhanced regions. They also showed higher signals for amide than for nuclear Overhauser enhancement effect in human brain, presumably reflecting the higher increment from faster-exchanging species as B_1,sat increased.

CONCLUSION

Featuring high contrast-to-noise ratio efficiency, REST could be a practical exchange-dependent approach readily applicable to either retrospective Z-spectra analysis or perspective 6-offset acquisition.

Chemical Exchange Saturation Transfer Magnetic Resonance Imaging for Longitudinal Assessment of Intracerebral Hemorrhage and Deferoxamine Treatment at 3T in a Mouse Model

BACKGROUND

Noninvasive imaging of molecular alterations after intracerebral hemorrhage (ICH) could provide valuable information to guide and monitor treatments. Chemical exchange saturation transfer (CEST) magnetic resonance imaging has demonstrated promises in identifying proliferation, necrosis, and changes in cellularity in brain tumors. Here, we applied CEST magnetic resonance imaging to monitor molecular changes in hematoma without and with treatment noninvasively over 2 weeks at 3T using endogenous contrast.

METHODS

CEST contrast related to proteins at 3.5 ppm (amide proton transfer) and proteins/lipids at -3.5 ppm (relayed nuclear overhauser effect [rNOE]) were examined over 14 days in a collagenase-induced ICH mouse model. Imaging findings were validated with immunohistochemistry based on the ICH neuropathology. We also examined iron-containing phantoms that mimicked iron concentrations in hematoma to ensure the iron will not attenuate the CEST contrast during disease progression. Based on the validity of the CEST contrast of hematoma, we further examined related molecular alterations under iron-chelation treatment with deferoxamine.

RESULTS

We observed the temporal and spatial differences of CEST contrasts between rNOE at -3.5 ppm and amide proton transfer at 3.5 ppm, in which the core and perihematoma could be identified by rNOE on day 3 and day 14, and amide proton transfer on day 1, day 7, and day 14. Moreover, we observed a 25.7% significant reduction (P<0.05) of rNOE contrast after deferoxamine treatment to the ICH mice on day 3, which was not observable in amide proton transfer contrast. Our histology data indicated that rNOE primarily correlated with the myelin pathology, and amide proton transfer could reflect the cellularity increase at hematoma up to day 7.

CONCLUSIONS

Significant rNOE changes correlated well with histologic findings, especially myelin lipids, and regional characteristics in hematoma indicate the uniqueness of CEST magnetic resonance imaging in monitoring molecular changes during ICH and treatment.

Uncertainty-aware self-supervised neural network for liverT1ρmapping with relaxation constraint

Objective.T1ρmapping is a promising quantitative MRI technique for the non-invasive assessment of tissue properties. Learning-based approaches can mapT1ρfrom a reduced number ofT1ρweighted images but requires significant amounts of high-quality training data. Moreover, existing methods do not provide the confidence level of theT1ρestimation. We aim to develop a learning-based liverT1ρmapping approach that can mapT1ρwith a reduced number of images and provide uncertainty estimation.Approach. We proposed a self-supervised neural network that learns aT1ρmapping using the relaxation constraint in the learning process. Epistemic uncertainty and aleatoric uncertainty are modelled for theT1ρquantification network to provide a Bayesian confidence estimation of theT1ρmapping. The uncertainty estimation can also regularize the model to prevent it from learning imperfect data. Main results. We conducted experiments onT1ρdata collected from 52 patients with non-alcoholic fatty liver disease. The results showed that when only collecting twoT1ρ-weighted images, our method outperformed the existing methods forT1ρquantification of the liver. Our uncertainty estimation can further regularize the model to improve the performance of the model and it is consistent with the confidence level of liverT1ρvalues.Significance. Our method demonstrates the potential for accelerating theT1ρmapping of the liver by using a reduced number of images. It simultaneously provides uncertainty ofT1ρquantification which is desirable in clinical applications.

IMPULSED model based cytological feature estimation with U-Net: Application to human brain tumor at 3T

PURPOSE

This work introduces and validates a deep-learning-based fitting method, which can rapidly provide accurate and robust estimation of cytological features of brain tumor based on the IMPULSED (imaging microstructural parameters using limited spectrally edited diffusion) model fitting with diffusion-weighted MRI data.

METHODS

The U-Net was applied to rapidly quantify extracellular diffusion coefficient (D_ex ), cell size (d), and intracellular volume fraction (v_in ) of brain tumor. At the training stage, the image-based training data, synthesized by randomizing quantifiable microstructural parameters within specific ranges, was used to train U-Net. At the test stage, the pre-trained U-Net was applied to estimate the microstructural parameters from simulated data and the in vivo data acquired on patients at 3T. The U-Net was compared with conventional non-linear least-squares (NLLS) fitting in simulations in terms of estimation accuracy and precision.

RESULTS

Our results confirm that the proposed method yields better fidelity in simulations and is more robust to noise than the NLLS fitting. For in vivo data, the U-Net yields obvious quality improvement in parameter maps, and the estimations of all parameters are in good agreement with the NLLS fitting. Moreover, our method is several orders of magnitude faster than the NLLS fitting (from about 5 min to <1 s).

CONCLUSION

The image-based training scheme proposed herein helps to improve the quality of the estimated parameters. Our deep-learning-based fitting method can estimate the cell microstructural parameters fast and accurately.

High-resolution neural network-driven mapping of multiple diffusion metrics leveraging asymmetries in the balanced steady-state free precession frequency profile

We propose to utilize the rich information content about microstructural tissue properties entangled in asymmetric balanced steady-state free precession (bSSFP) profiles to estimate multiple diffusion metrics simultaneously by neural network (NN) parameter quantification. A 12-point bSSFP phase-cycling scheme with high-resolution whole-brain coverage is employed at 3 and 9.4 T for NN input. Low-resolution target diffusion data are derived based on diffusion-weighted spin-echo echo-planar-imaging (SE-EPI) scans, that is, mean, axial, and radial diffusivity (MD, AD, and RD), fractional anisotropy (FA), as well as the spherical coordinates (azimuth Φ and inclination ϴ) of the principal diffusion eigenvector. A feedforward NN is trained with incorporated probabilistic uncertainty estimation. The NN predictions yielded highly reliable results in white matter (WM) and gray matter structures for MD. The quantification of FA, AD, and RD was overall in good agreement with the reference but the dependence of these parameters on WM anisotropy was somewhat biased (e.g. in corpus callosum). The inclination ϴ was well predicted for anisotropic WM structures, while the azimuth Φ was overall poorly predicted. The findings were highly consistent across both field strengths. Application of the optimized NN to high-resolution input data provided whole-brain maps with rich structural details. In conclusion, the proposed NN-driven approach showed potential to provide distortion-free high-resolution whole-brain maps of multiple diffusion metrics at high to ultrahigh field strengths in clinically relevant scan times.

Quantitative imaging of apoptosis following oncolytic virotherapy by magnetic resonance fingerprinting aided by deep learning

Non-invasive imaging methods for detecting intratumoural viral spread and host responses to oncolytic virotherapy are either slow, lack specificity or require the use of radioactive or metal-based contrast agents. Here we show that in mice with glioblastoma multiforme, the early apoptotic responses to oncolytic virotherapy (characterized by decreased cytosolic pH and reduced protein synthesis) can be rapidly detected via chemical-exchange-saturation-transfer magnetic resonance fingerprinting (CEST-MRF) aided by deep learning. By leveraging a deep neural network trained with simulated magnetic resonance fingerprints, CEST-MRF can generate quantitative maps of intratumoural pH and of protein and lipid concentrations by selectively labelling the exchangeable amide protons of endogenous proteins and the exchangeable macromolecule protons of lipids, without requiring exogenous contrast agents. We also show that in a healthy volunteer, CEST-MRF yielded molecular parameters that are in good agreement with values from the literature. Deep-learning-aided CEST-MRF may also be amenable to the characterization of host responses to other cancer therapies and to the detection of cardiac and neurological pathologies.

Learning-based optimization of acquisition schedule for magnetization transfer contrast MR fingerprinting

Magnetization transfer contrast MR fingerprinting (MTC-MRF) is a novel quantitative imaging method that simultaneously quantifies free bulk water and semisolid macromolecule parameters using pseudo-randomized scan parameters. To improve acquisition efficiency and reconstruction accuracy, the optimization of MRF sequence design has been of recent interest in the MRF field, but has been challenging due to the large number of degrees of freedom to be optimized in the sequence. Herein, we propose a framework for learning-based optimization of the acquisition schedule (LOAS), which optimizes RF saturation-encoded MRF acquisitions with a minimal number of scan parameters for tissue parameter determination. In a supervised learning framework, scan parameters were subsequently updated to minimize a predefined loss function that can directly represent tissue quantification errors. We evaluated the performance of the proposed approach with a numerical phantom and in in vivo experiments. For validation, MRF images were synthesized using the tissue parameters estimated from a fully connected neural network framework and compared with references. Our results showed that LOAS outperformed existing indirect optimization methods with regard to quantification accuracy and acquisition efficiency. The proposed LOAS method could be a powerful optimization tool in the design of MRF pulse sequences.

Bayesian deep learning-based 1 H-MRS of the brain: Metabolite quantification with uncertainty estimation using Monte Carlo dropout

PURPOSE

To develop a Bayesian convolutional neural network (BCNN) with Monte Carlo dropout sampling for metabolite quantification with simultaneous uncertainty estimation in deep learning-based proton MRS of the brain.

METHODS

Human brain spectra were simulated using basis spectra for 17 metabolites and macromolecules (N = 100 000) at 3.0 Tesla. In addition, actual in vivo spectra (N = 5) were modified by adjusting SNR and linewidth with increasing severity of spectral degradation (N = 50). A BCNN was trained on the simulated spectra to generate a noise-free, line-narrowed, macromolecule signal-removed, metabolite-only spectrum from a typical human brain spectrum. At inference, each input spectrum was Monte Carlo dropout sampled (50 times), and the resulting mean spectrum and variance spectrum were used for metabolite quantification and uncertainty estimation, respectively.

RESULTS

Using the simulated spectra, the mean absolute percent errors of the BCNN-predicted metabolite content were < 10% for Cr, Glu, Gln, mI, NAA, and Tau (< 5% for Glu, NAA, and mI). For all metabolites, the correlations (r's) between the ground-truth error and BCNN-predicted uncertainty ranged 0.72-0.94 (0.83 ± 0.06; p < 0.001). Using the modified in vivo spectra, the extent of variation in the estimated metabolite content against the increasing severity of spectral degradation tended to be smaller with BCNN than with linear combination of model spectra (LCModel). Overall, the variation in metabolite content tended to be more highly correlated with the uncertainty from BCNN than with the Cramér-Rao lower-bounds from LCModel (0.938 ± 0.019 vs. 0.881 ± 0.057 [p = 0.115]).

CONCLUSION

The BCNN with Monte Carlo dropout sampling may be used in deep learning-based MRS for the estimation of uncertainty in the machine-predicted metabolite content, which is important in the clinical application of deep learning-based MRS.

Molecular Imaging of Brain Tumors and Drug Delivery Using CEST MRI: Promises and Challenges

Chemical exchange saturation transfer (CEST) magnetic resonance imaging (MRI) detects molecules in their natural forms in a sensitive and non-invasive manner. This makes it a robust approach to assess brain tumors and related molecular alterations using endogenous molecules, such as proteins/peptides, and drugs approved for clinical use. In this review, we will discuss the promises of CEST MRI in the identification of tumors, tumor grading, detecting molecular alterations related to isocitrate dehydrogenase (IDH) and O-6-methylguanine-DNA methyltransferase (MGMT), assessment of treatment effects, and using multiple contrasts of CEST to develop theranostic approaches for cancer treatments. Promising applications include (i) using the CEST contrast of amide protons of proteins/peptides to detect brain tumors, such as glioblastoma multiforme (GBM) and low-grade gliomas; (ii) using multiple CEST contrasts for tumor stratification, and (iii) evaluation of the efficacy of drug delivery without the need of metallic or radioactive labels. These promising applications have raised enthusiasm, however, the use of CEST MRI is not trivial. CEST contrast depends on the pulse sequences, saturation parameters, methods used to analyze the CEST spectrum (i.e., Z-spectrum), and, importantly, how to interpret changes in CEST contrast and related molecular alterations in the brain. Emerging pulse sequence designs and data analysis approaches, including those assisted with deep learning, have enhanced the capability of CEST MRI in detecting molecules in brain tumors. CEST has become a specific marker for tumor grading and has the potential for prognosis and theranostics in brain tumors. With increasing understanding of the technical aspects and associated molecular alterations detected by CEST MRI, this young field is expected to have wide clinical applications in the near future.

MR fingerprinting for semisolid magnetization transfer and chemical exchange saturation transfer quantification

Chemical exchange saturation transfer (CEST) MRI has positioned itself as a promising contrast mechanism, capable of providing molecular information at sufficient resolution and amplified sensitivity. However, it has not yet become a routinely employed clinical technique, due to a variety of confounding factors affecting its contrast-weighted image interpretation and the inherently long scan time. CEST MR fingerprinting (MRF) is a novel approach for addressing these challenges, allowing simultaneous quantitation of several proton exchange parameters using rapid acquisition schemes. Recently, a number of deep-learning algorithms have been developed to further boost the performance and speed of CEST and semi-solid macromolecule magnetization transfer (MT) MRF. This review article describes the fundamental theory behind semisolid MT/CEST-MRF and its main applications. It then details supervised and unsupervised learning approaches for MRF image reconstruction and describes artificial intelligence (AI)-based pipelines for protocol optimization. Finally, practical considerations are discussed, and future perspectives are given, accompanied by basic demonstration code and data.

Deep learning-based classification of preclinical breast cancer tumor models using chemical exchange saturation transfer magnetic resonance imaging

Chemical exchange saturation transfer (CEST) magnetic resonance imaging has shown promise for classifying tumors based on their aggressiveness, but CEST contrast is complicated by multiple signal sources and thus prolonged acquisition times are often required to extract the signal of interest. We investigated whether deep learning could help identify pertinent Z-spectral features for distinguishing tumor aggressiveness as well as the possibility of acquiring only the pertinent spectral regions for more efficient CEST acquisition. Human breast cancer cells, MDA-MB-231 and MCF-7, were used to establish bi-lateral tumor xenografts in mice to represent higher and lower aggressive tumors, respectively. A convolutional neural network (CNN)-based classification model, trained on simulated data, utilized Z-spectral features as input to predict labels of different tissue types, including MDA-MB-231, MCF-7, and muscle tissue. Saliency maps reported the influence of Z-spectral regions on classifying tissue types. The model was robust to noise with an accuracy of more than 91.5% for low and moderate noise levels in simulated testing data (SD of noise less than 2.0%). For in vivo CEST data acquired with a saturation pulse amplitude of 2.0 μT, the model had a superior ability to delineate tissue types compared with Lorentzian difference (LD) and magnetization transfer ratio asymmetry (MTR_asym ) analysis, classifying tissues to the correct types with a mean accuracy of 85.7%, sensitivity of 81.1%, and specificity of 94.0%. The model's performance did not improve substantially when using data acquired at multiple saturation pulse amplitudes or when adding LD or MTR_asym spectral features, and did not change when using saliency map-based partial or downsampled Z-spectra. This study demonstrates the potential of CNN-based classification to distinguish between different tumor types and muscle tissue, and speed up CEST acquisition protocols.

A Fully Convolutional Network (FCN) based Automated Ischemic Stroke Segment Method using Chemical Exchange Saturation Transfer Imaging

BACKGROUND

Chemical exchange saturation transfer (CEST) MRI is a promising imaging modality in ischemic stroke detection for its sensitivity in sensing post-ischemic pH alteration. However, the accurate segmentation of pH-altered regions remains difficult due to the complicated sources in water signal changes of CEST MRI. Meanwhile, manual localization and quantification of stroke lesions are laborious and time-consuming, which cannot meet the urgent need for timely therapeutic interventions.

PURPOSE

The goal of this study was to develop an automatic lesion segmentation approach of ischemic region based on CEST MR images. A novel segmentation framework based on fully convolutional neural network was investigated for our task.

METHODS

Z-spectra from 10 rats were manually labeled as ground truth and split into two datasets, where the training dataset including 3 rats was used to generate a segmentation model, and the remaining rats were used as test datasets to evaluate the model's performance. Then a 1-D fully convolutional neural network equipped with bottleneck structures was set up, and a Grad-CAM approach was used to produce a coarse localization map, which can reflect the relevancy to the 'ischemia' class of each pixel.

RESULTS

As compared with the ground truth, the proposed network model achieved satisfying segmentation results with high values of evaluation metrics including specificity (SPE), sensitivity (SEN), accuracy (ACC), and Dice similarity coefficient (DSC), especially in some intractable situations where conventional MRI modalities and CEST quantitative method failed to distinguish between ischemic and normal tissues, and the model with augmentation was robust to input perturbations. The Grad-CAM maps performed clear tissue change distributions and interpreted the segmentations, and showed a strong correlation with the quantitative method, gave extended thinking to the function of networks.

CONCLUSIONS

The proposed method can segment ischemia region from CEST images, with the Grad-CAM maps give access to interpretative information about the segmentations, which demonstrates great potential in clinical routines. This article is protected by copyright. All rights reserved.

Deep neural network based CEST and AREX processing: Application in imaging a model of Alzheimer's disease at 3 T

PURPOSE

To optimize and apply deep neural network based CEST (deepCEST) and apparent exchange dependent-relaxation (deepAREX) for imaging the mouse brain with Alzheimer's disease (AD) at 3T MRI.

METHODS

CEST and T₁ data of central and anterior brain slices of 10 AD mice and 10 age-matched wild type (WT) mice were acquired at a 3T animal MRI scanner. The networks of deepCEST/deepAREX were optimized and trained on the WT data. The CEST/AREX contrasts of AD and WT mice predicted by the networks were analyzed and further validated by immunohistochemistry.

RESULTS

After optimization and training on CEST data of WT mice, deepCEST/deepAREX could rapidly (~1 s) generate precise CEST and AREX results for unseen CEST data of AD mice, indicating the accuracy and generalization of the networks. Significant lower amide weighted (3.5 ppm) signal related to amyloid β-peptide (Aβ) plaque depositions, which was validated by immunohistochemistry results, was detected in both central and anterior brain slices of AD mice compared to WT mice. Decreased magnetization transfer (MT) signal was also found in AD mice especially in the anterior slice.

CONCLUSION

DeepCEST/deepAREX could rapidly generate accurate CEST/AREX contrasts in animal study. The well-optimized deepCEST/deepAREX have potential for AD differentiation at 3T MRI.

Pulseq-CEST: Towards multi-site multi-vendor compatibility and reproducibility of CEST experiments using an open-source sequence standard

PURPOSE

As the field of CEST grows, various novel preparation periods using different parameters are being introduced. At the same time, large, multisite clinical studies require clearly defined protocols, especially across different vendors. Here, we propose a CEST definition standard using the open Pulseq format for a shareable, simple, and exact definition of CEST protocols.

METHODS

We present the benefits of such a standard in three ways: (1) an open database on GitHub, where fully defined, human-readable CEST protocols can be shared; (2) an open-source Bloch-McConnell simulation to test and optimize CEST preparation periods in silico; and (3) a hybrid MR sequence that plays out the CEST preparation period and can be combined with any existing readout module.

RESULTS

The exact definition of the CEST preparation period, in combination with the flexible simulation, leads to a good match between simulations and measurements. The standard allowed finding consensus on three amide proton transfer-weighted protocols that could be compared in healthy subjects and a tumor patient. In addition, we could show coherent multisite results for a sophisticated CEST method, highlighting the benefits regarding protocol sharing and reproducibility.

CONCLUSION

With Pulseq-CEST, we provide a straightforward approach to standardize, share, simulate, and measure different CEST preparation schemes, which are inherently completely defined.

Hyperpolarized MRI, functional MRI, MR spectroscopy and CEST to provide metabolic information in vivo

Access to metabolic information in vivo using magnetic resonance (MR) technologies has generally been the niche of MR spectroscopy (MRS) and spectroscopic imaging (MRSI). Metabolic fluxes can be studied using the infusion of substrates labeled with magnetic isotopes, with the use of hyperpolarization especially powerful. Unfortunately, these promising methods are not yet accepted clinically, where fast, simple, and reliable measurement and diagnosis are key. Recent advances in functional MRI and chemical exchange saturation transfer (CEST) MRI allow the use of water imaging to study oxygen metabolism and tissue metabolite levels. These, together with the use of novel data analysis approaches such as machine learning for all of these metabolic MR approaches, are increasing the likelihood of their clinical translation.

Accelerated multi-target chemical exchange saturation transfer magnetic resonance imaging of the mouse heart

Cardiac chemical exchange saturation transfer-magnetic resonance imaging (CEST-MRI) has been used to probe levels of various metabolites that provide insight into myocardial structure and function. However, imaging of the heart using CEST-MRI is prolonged by the need to repeatedly acquire multiple images for a full Z-spectrum and to perform saturation and acquisition around cardiac and respiratory cycles. Compressed sensing (CS) reconstruction of sparse data enables accelerated acquisition, but reconstruction artifacts may bias subsequently derived measures of CEST contrast. In this study, we examine the impact of CS reconstruction of increasingly under-sampled cardiac CEST-MRI data on subsequent CEST contrasts of amine-containing metabolites and amide-containing proteins. Cardiac CEST-MRI data sets were acquired in six mice using low and high RF saturation for single and dual contrast generation, respectively. CEST-weighted images were reconstructed using CS methods at 2-5× levels of under-sampling. CEST contrasts were derived from corresponding Z-spectra and the impact of accelerated imaging on accuracy was assessed via analysis of variance. CS reconstruction preserved myocardial signal to noise ratio as compared to conventional reconstruction. However, greater absolute error and distribution of derived contrasts was observed with increasing acceleration factors. The results from this study indicate that acquisition of radial cardiac CEST-MRI data can be modestly, but meaningfully, accelerated via CS reconstructions with little error in CEST contrast quantification.

Uncertainty quantification implementations in human hemodynamic flows

BACKGROUND AND OBJECTIVE

Human hemodynamic modeling is usually influenced by uncertainties occurring from a considerable unavailability of information linked to the boundary conditions and the physical properties used in the numerical models. Calculating the effect of these uncertainties on the numerical findings along the cardiovascular system is a demanding process due to the complexity of the morphology of the body and the area dynamics. To cope with all these difficulties, Uncertainty Quantification (UQ) methods seem to be an ideal tool.

RESULTS

This study focuses on analyzing and summarizing some of the recent research efforts and directions of implementing UQ in human hemodynamic flows by analyzing 139 research papers. Initially, the suitability of applying this approach is analyzed and demonstrated. Then, an overview of the most significant research work in various fields of biomedical hemodynamic engineering is presented. Finally, it is attempted to identify any possible forthcoming directions for research and methodological progress of UQ in biomedical sciences.

CONCLUSION

This review concludes that by finding the best statistical methods and parameters to represent the propagated uncertainties, while achieving a good interpretation of the interaction between input-output, is crucial for implementing UQ in biomedical sciences.

Brain tissues have single-voxel signatures in multi-spectral MRI

Since the seminal works by Brodmann and contemporaries, it is well-known that different brain regions exhibit unique cytoarchitectonic and myeloarchitectonic features. Transferring the approach of classifying brain tissues - and other tissues - based on their intrinsic features to the realm of magnetic resonance (MR) is a longstanding endeavor. In the 1990s, atlas-based segmentation replaced earlier multi-spectral classification approaches because of the large overlap between the class distributions. Here, we explored the feasibility of performing global brain classification based on intrinsic MR features, and used several technological advances: ultra-high field MRI, q-space trajectory diffusion imaging revealing voxel-intrinsic diffusion properties, chemical exchange saturation transfer and semi-solid magnetization transfer imaging as a marker of myelination and neurochemistry, and current neural network architectures to analyze the data. In particular, we used the raw image data as well to increase the number of input features. We found that a global brain classification of roughly 97 brain regions was feasible with gross classification accuracy of 60%; and that mapping from voxel-intrinsic MR data to the brain region to which the data belongs is possible. This indicates the presence of unique MR signals of different brain regions, similar to their cytoarchitectonic and myeloarchitectonic fingerprints.

Fast chemical exchange saturation transfer imaging based on PROPELLER acquisition and deep neural network reconstruction

PURPOSE

To develop a method for fast chemical exchange saturation transfer (CEST) imaging.

METHODS

The periodically rotated overlapping parallel lines enhanced reconstruction (PROPELLER) sampling scheme was introduced to shorten the acquisition time. Deep neural network was employed to reconstruct CEST contrast images. Numerical simulation and experiments on a creatine phantom, hen egg, and in vivo tumor rat brain were performed to test the feasibility of this method.

RESULTS

The results from numerical simulation and experiments show that there is no significant difference between reference images and CEST-PROPELLER reconstructed images under an acceleration factor of 8.

CONCLUSION

Although the deep neural network is trained entirely on synthesized data, it works well on reconstructing experimental data. The proof of concept study demonstrates that the combination of the PROPELLER sampling scheme and the deep neural network enables considerable acceleration of saturated image acquisition and may find applications in CEST MRI.

Multi-parametric artificial neural network fitting of phase-cycled balanced steady-state free precession data

PURPOSE

Standard relaxation time quantification using phase-cycled balanced steady-state free precession (bSSFP), eg, motion-insensitive rapid configuration relaxometry (MIRACLE), is subject to a considerable underestimation of tissue T₁ and T₂ due to asymmetric intra-voxel frequency distributions. In this work, an artificial neural network (ANN) fitting approach is proposed to simultaneously extract accurate reference relaxation times (T₁ , T₂ ) and robust field map estimates ( B 1 + , ΔB₀ ) from the bSSFP profile.

METHODS

Whole-brain bSSFP data acquired at 3T were used for the training of a feedforward ANN with N = 12, 6, and 4 phase-cycles. The magnitude and phase of the Fourier transformed complex bSSFP frequency response served as input and the multi-parametric reference set [T₁ , T₂ , B 1 + , ∆B₀ ] as target. The ANN predicted relaxation times were validated against the target and MIRACLE.

RESULTS

The ANN prediction of T₁ and T₂ for trained and untrained data agreed well with the reference, even for only four acquired phase-cycles. In contrast, relaxometry based on 4-point MIRACLE was prone to severe off-resonance-related artifacts. ANN predicted B 1 + and ∆B₀ maps showed the expected spatial inhomogeneity patterns in high agreement with the reference measurements for 12-point, 6-point, and 4-point bSSFP phase-cycling schemes.

CONCLUSION

ANNs show promise to provide accurate brain tissue T₁ and T₂ values as well as reliable field map estimates. Moreover, the bSSFP acquisition can be accelerated by reducing the number of phase-cycles while still delivering robust T₁ , T₂ , B 1 + , and ∆B₀ estimates.