Model-Based Chemical Exchange Saturation Transfer MRI for Robust z-Spectrum Analysis

Prospective acceleration of whole-brain CEST imaging by golden-angle view ordering in Cartesian coordinates and joint k-space and image-space parallel imaging (KIPI)

PURPOSE

To prospectively accelerate whole-brain CEST acquisition by joint k-space and image-space parallel imaging (KIPI) with a proposed golden-angle view ordering technique (GAVOT) in Cartesian coordinates.

THEORY AND METHODS

The T2-decay effect will vary across frames with variable acceleration factors (AF) in the prospective acquisition using sequences with long echo trains. The GAVOT method uses a subset strategy to eliminate the T2-decay inconsistency, where all frames use a subset of shots from the calibration frame to form their k-space view ordering. The golden-angle rule is adapted to ensure uniform k-space coverage for arbitrary AFs in Cartesian coordinates. Phantom and in vivo studies were conducted on a 3 T scanner.

RESULTS

The GAVOT view ordering yielded a higher g-factor than conventional uniformly centric ordering, whereas the noise propagation in amide proton transfer (APT) weighted images was similar between different view ordering. Compared to centric ordering, GAVOT successfully eliminated the T2-decay inconsistency across all frames, resulting in fewer image artifacts for both KIPI and conventional parallel imaging methods. The synergy of GAVOT and KIPI mitigated strong aliasing artifacts and achieved high-quality reconstruction of prospective variable-AF datasets. GAVOT-KIPI reduced the scan time to 2.1 min for whole-brain APT weighted imaging and 4.7 min for quantitative APT signal (APT#) mapping.

CONCLUSION

GAVOT makes the prospective variable AF strategy flexible and practical, and, in conjunction with KIPI, ensures high-quality reconstruction from highly undersampled data, facilitating the clinical translation of whole-brain CEST imaging.

Nonlinear parameter estimation with physics-constrained spectral-spatial priors for highly accelerated chemical exchange saturation transfer MRI

Objective.To develop a nonlinear, model-based parameter estimation method directly from incomplete measurements ink - wspace for robust spectral analysis in highly accelerated chemical exchange saturation transfer (CEST) magnetic resonance imaging (MRI).Approach. A CEST-specific, separable nonlinear model, which describes spectral decomposition using multi-pool Lorentzian functions (conventional magnetization transfer (MT), direct saturation of water signals (DS), amide, amine, and nuclear Overhauser effect) derived from the steady-state Bloch McConnel equation, is incorporated into a measurement model in CEST MRI. Furthermore, signal drop in saturation transfer experiments is formulated by an additional, separable nonlinear spectral prior indicating that the symmetric z-spectra synthesized using conventional MT and DS always remain higher or equal to the whole z-spectra with all pools. Given the above considerations, linear and nonlinear parameters in the proposed method are estimated in an alternating fashion directly from highly incomplete measurements ink - wspace by solving a constrained optimization problem with the physics-constrained spectral priors while imposing additional sparsity priors on spatial parameter maps.Main results.Compared with conventional methods, the proposed method yields clearer delineation of tumor-specific CEST maps without apparent artifact and noise.Significance.We successfully demonstrated the feasibility of the proposed method for CEST MRI with highly incomplete measurements thus enabling high-resolution whole brain CEST MRI in clinically reasonable imaging time.

Highly-accelerated CEST MRI using frequency-offset-dependent k-space sampling and deep-learning reconstruction

PURPOSE

To develop a highly accelerated CEST Z-spectral acquisition method using a specifically-designed k-space sampling pattern and corresponding deep-learning-based reconstruction.

METHODS

For k-space down-sampling, a customized pattern was proposed for CEST, with the randomized probability following a frequency-offset-dependent (FOD) function in the direction of saturation offset. For reconstruction, the convolution network (CNN) was enhanced with a Partially Separable (PS) function to optimize the spatial domain and frequency domain separately. Retrospective experiments on a self-acquired human brain dataset (13 healthy adults and 15 brain tumor patients) were conducted using k-space resampling. The prospective performance was also assessed on six healthy subjects.

RESULTS

In retrospective experiments, the combination of FOD sampling and PS network (FOD + PSN) showed the best quantitative metrics for reconstruction, outperforming three other combinations of conventional sampling with varying density and a regular CNN (nMSE and SSIM, p < 0.001 for healthy subjects). Across all acceleration factors from 4 to 14, the FOD + PSN approach consistently outperformed the comparative methods in four contrast maps including MTRasym, MTRrex, as well as the Lorentzian Difference maps of amide and nuclear Overhauser effect (NOE). In the subspace replacement experiment, the error distribution demonstrated the denoising benefits achieved in the spatial subspace. Finally, our prospective results obtained from healthy adults and brain tumor patients (14×) exhibited the initial feasibility of our method, albeit with less accurate reconstruction than retrospective ones.

CONCLUSION

The combination of FOD sampling and PSN reconstruction enabled highly accelerated CEST MRI acquisition, which may facilitate CEST metabolic MRI for brain tumor patients.

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.

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.

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

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

3D CEST MRI with an unevenly segmented RF irradiation scheme: A feasibility study in brain tumor imaging

PURPOSE

To integrate 3D CEST EPI with an unevenly segmented RF irradiation module and preliminarily demonstrate it in the clinical setting.

METHODS

A CEST MRI with unevenly segmented RF saturation was implemented, including a long primary RF saturation to induce the steady-state CEST effect, maintained with repetitive short secondary RF irradiation between readouts. This configuration reduces relaxation-induced blur artifacts during acquisition, allowing fast 3D spatial coverage. Numerical simulations were performed to select parameters such as flip angle (FA), short RF saturation duration (Ts2), and the number of readout segments. The sequence was validated experimentally with data from a phantom, healthy volunteers, and a brain tumor patient.

RESULTS

Based on the numerical simulation and l-carnosine gel phantom experiment, FA, Ts2, and the number of segments were set to 20°, 0.3 s, and the range from 4 to 8, respectively. The proposed method minimized signal modulation in the human brain images in the kz direction during the acquisition and provided the blur artifacts-free CEST contrast over the whole volume. Additionally, the CEST contrast in the tumor tissue region is higher than in the contralateral normal tissue region.

CONCLUSIONS

It is feasible to implement a highly accelerated 3D EPI CEST imaging with unevenly segmented RF irradiation.

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.

Joint K-space and Image-space Parallel Imaging (KIPI) for accelerated chemical exchange saturation transfer acquisition

PURPOSE

To develop an auto-calibrated technique by joint K-space and Image-space Parallel Imaging (KIPI) for accelerated CEST acquisition.

THEORY AND METHODS

The KIPI method selects a calibration frame with a low acceleration factor (AF) and auto-calibration signals (ACS) acquired, from which the coil sensitivity profiles and artifact correction maps are calculated after restoring the k-space by GRAPPA. Then the other frames with high AF and without ACS can be reconstructed by SENSE and artifact suppression. The signal leakage due to the T₂ -decay filtering in k-space compromises the SENSE reconstruction, which can be corrected by the artifact suppression algorithm of KIPI. The 2D and 3D imaging experiments were done on the phantom, healthy volunteer, and brain tumor patient with a 3T scanner.

RESULTS

The proposed KIPI method was evaluated by retrospectively undersampled data with variable AFs and compared against existing parallel imaging methods (SENSE/auto, GRAPPA, and ESPIRiT). KIPI enabled CEST frames with random AFs to achieve similar image quality, eliminated the strong aliasing artifacts, and generated significantly smaller errors than the other methods (p < 0.01). The KIPI method permitted an AF up to 12-fold in both phase-encoding and slice-encoding directions for 3D CEST source images, achieving an overall 8.2-fold speedup in scan time.

CONCLUSION

KIPI is a novel auto-calibrated parallel imaging method that enables variable AFs for different CEST frames, achieves a significant reduction in scan time, and does not compromise the accuracy of CEST maps.

B0 Correction for 3T Amide Proton Transfer (APT) MRI Using a Simplified Two-Pool Lorentzian Model of Symmetric Water and Asymmetric Solutes

Amide proton transfer (APT)-weighted MRI is a promising molecular imaging technique that has been employed in clinic for detection and grading of brain tumors. MTRasym, the quantification method of APT, is easily influenced by B0 inhomogeneity and causes artifacts. Current model-free interpolation methods have enabled moderate B0 correction for middle offsets, but have performed poorly at limbic offsets. To address this shortcoming, we proposed a practical B0 correction approach that is suitable under time-limited sparse acquisition scenarios and for B1 ≥ 1 μT under 3T. In this study, this approach employed a simplified Lorentzian model containing only two pools of symmetric water and asymmetric solutes, to describe the Z-spectral shape with wide and ‘invisible’ CEST peaks. The B0 correction was then performed on the basis of the fitted two-pool Lorentzian lines, instead of using conventional model-free interpolation. The approach was firstly evaluated on densely sampled Z-spectra data by using the spline interpolation of all acquired 16 offsets as the gold standard. When only six offsets were available for B0 correction, our method outperformed conventional methods. In particular, the errors at limbic offsets were significantly reduced (n = 8, p < 0.01). Secondly, our method was assessed on the six-offset APT data of nine brain tumor patients. Our MTRasym (3.5 ppm), using the two-pool model, displayed a similar contrast to the vendor-provided B0-orrected MTRasym (3.5 ppm). While the vendor failed in correcting B0 at 4.3 and 2.7 ppm for a large portion of voxels, our method enabled well differentiation of B0 artifacts from tumors. In conclusion, the proposed approach could alleviate analysis errors caused by B0 inhomogeneity, which is useful for facilitating the comprehensive metabolic analysis of brain tumors.

Rapid three-dimensional steady-state chemical exchange saturation transfer magnetic resonance imaging

PURPOSE

To make clinically feasible whole-brain chemical exchange saturation transfer (CEST) magnetic resonance imaging (MRI) by enhancing imaging efficiency.

METHODS

A novel, whole-brain three-dimensional (3D) steady-state CEST MRI method was introduced by utilizing a time-efficient, fat-suppressed excitation followed by rapid, segmented 3D echo-planar-imaging with incoherent undersampling in k-ω space. Missing signals and CEST-specific spectral images were then jointly estimated directly from incomplete measurements using model-based reconstruction and robust spectral analysis. In vivo studies were performed at 3T both retrospectively (using a fully sampled reference) and prospectively to validate the effectiveness of the proposed method in patients with brain cancer.

RESULTS

In retrospective studies, the proposed method exhibits superior accuracies to existing methods in estimating images, z-spectra, and APTw relative to the reference. In prospective patient studies, compared with existing methods, the proposed method is statistically significantly different in contrast-to-noise ratio of the APTw contrast between tumor and normal appearing white matter (NAWM) and amide proton transfer weighted contrast (p < 0.05) while not being significantly different in signal-to-noise ratio in an NAWM region.

CONCLUSIONS

We successfully demonstrated that it is feasible to perform whole-brain CEST MRI roughly within 4 minutes for patients with brain cancer. It is expected that the proposed method widens clinical utilities of CEST MRI.