4722

Accelerated in-vivo 23Na Multi-Quantum Coherences MRI by utilizing Low-Rank Matrix Completion
Christian Licht1, Stanislas Rapacchi2, and Lothar R. Schad1
1Computer Assisted Clinical Medicine, Heidelberg University, Mannheim, Germany, 2CRMBM, Aix-Marseille Université, Marseille, France

Synopsis

Keywords: Non-Proton, Multi-Contrast, Sodium

Sodium (23Na) MRI has received increased attention as a potential biomarker for disease states thanks to the advent of ultra-high field MRI. One interesting asset from 23Na MRI is the distinction between single and triple quantum signal to further characterize tissues. However, 3D 23Na multi-quantum coherences (MQC) imaging requires multiple radiofrequency phase-cycling, which is inherently time-consuming and therefore difficult to include in protocols. In this work, we propose to accelerate 23Na MQC MRI by leveraging multi-dimensional under-sampling coupled with a dedicated Low-Rank matrix completion image reconstruction.

Introduction

Conventional sodium (23Na) MRI is a promising tool to probe tissue ionic homeostasis, but often remains limited to the analysis of the sodium signal intensity1. Thanks to its 3/2 spin, 23Na multi-quantum coherences MRI on the other hand enables to disentangle the underlying multi-quantum coherences (MQC) and therefore holds potentially richer sodium tissue characterization2 when compared to conventional 23Na signal. Especially the triple quantum (TQ) signal has shown to be sensitive to the intracellular sodium compartment and the molecular environment surrounding sodium atoms3. Unfortunately, 23Na MQC MRI requires a minimum of 12 RF phase cycling steps, which prolongs acquisition time up to one hour and is thus impractical for large scale in-vivo studies. Alternatively, 23Na MQC MRI exhibits multi-dimensional data, which is therefore intrinsically redundant and can be exploited to reconstruct undersampled 23Na MQC MRI data. Thus, for the first time we propose to utilize low-rank matrix completion to reliably reconstruct phase-cycled under-sampled phantom, in-vivo knee and brain 23Na MQC MRI.

Methods

All measurements were performed on 3T MRI (Trio, Siemens Heathineers, Erlangen,Germany) with a bird-cage dual-tuned 23Na / 1H head coil (RapidBiomedical, Rimpar, Germany). 23Na MQC MR images were obtained by utilizing the CRISTINA sequence4, with the following parameters.
In phantom: FoV 260x260x240mm3, matrix size 32x32x12, $$$\tau$$$ =8.1ms, BW=200Hz/px,TE/$$$\Delta$$$TE/NTE= 1.74/6.4ms/10, TR=140ms, 7 averages resulting in TA=2x30 minutes.
For knee in-vivo, one volunteer (42, m), FoV 270x270x200mm3,matrix size 26x26x10, $$$\tau$$$=5.7ms, BW=220Hz/px, TE/ $$$\Delta$$$TE/NTE=1.71/6.2ms/10, TR=150ms, TA=2x31.
For brain in-vivo, in three volunteers (28+/-2, f): FoV 230x230x160mm, matrix size 24x24x8, $$$\tau$$$=12.9ms, BW=210Hz/px, TE/$$$\Delta$$$TE/NTE= 1.62/6.2ms/10, TR=150ms resulting in TA=2x31min.
Image reconstruction: 23Na MQC raw data were processed either fully-sampled or with a retrospective variable density 3D Poisson-disk undersampling5 to simulate various acceleration rates. Sampling pattern was varied along the phase-cycle dimension to create complementary patterns. Finally, the images were reconstructed utilizing the SAKE6 framework, replacing the coil dimension with the phase-cycled dimension, and solving the following optimization problem:
$$ \min_u || \Phi_F(u) - f ||^2_2 < \sigma^2 s.t. u = H^*(A), rank(A) = k' $$
with $$$u$$$ being the image to reconstruct, $$$\Phi_F$$$ being the Fourier transform operator, $$$f$$$ the acquired data and $$$\sigma$$$ the noise signal variance . $$$k'$$$ is a prior estimate of $$$rank(A)$$$ of matrix $$$A$$$, enforcing subspace SVD and $$$H^*$$$ being the inverse structured Hankel matrix operator6. By reshaping multi-echo and phase-cycle dimensions, rank deficiency of the data matrix is thus enforced. Reconstruction performance was evaluated on phantom and in-vivo brain data by structural similarity (SSIM), root mean squared error (RMSE) and sodium quantification.

Results

For 3-fold undersampled in-vivo knee 23Na MQC MRI, SAKE improved SSIM for the SQ image significantly by 0.07, when compared to 3D CS(Fig.3). Especially for higher acceleration factors, SAKE outperformed compressed sensing by reducing RMSE in SQ and TQ/SQ ratio. In addition, SAKE retained anatomical structures in the TQ images, due to shared information in 23Na MQC's multi-dimensional space, which resulted in degraded RMSE and SSIM. These findings were spatially consistent across several slices (Fig.4). SAKE images exhibited superior reconstruction performance on in-vivo brain data by reducing the RMSE for the SQ images by a factor of 2 but not for TQ images when compared to 3D CS( Fig.5).

Discussion & Conclusion

Conclusively, we have demonstrated that low-rank matrix completion is an interesting framework to accelerate 23Na MQC MRI that can be used to reconstruct 3-fold undersampled in-vivo knee data. SAKE is especially well-suited to reconstruct in-vivo brain and knee SQ images, but fails at reconstructing TQ brain signals. This issue is attributed to shared information along the dimensions and the exploitation of subspace singular value decomposition. For instance, SQ signal exhibits relatively high signal amplitude when compared to TQ. Hence, SQ signal is coherent along temporal dimension, potentially biasing TQ image reconstruction. Further integration combining low-rank and MQC signal modeling as well as prospective undersampling to increase resolution will be investigated in the future to fully demonstrate the potential of low-rank matrix completion for 23Na MQC MRI.

Acknowledgements

The author received financial support for the research by ISMRM research exchange grant program and PROCOPE Mobility 2022.

References

1. Huhn K, Engelhorn T, Linker RA, Nagel AM. Potential of Sodium MRI as a Biomarker for Neurodegeneration and Neuroinflammation in Multiple Sclerosis. Front Neurol. 2019 Feb.

2. Hoesl, M.A., Kleimaier, D., Hu, R., Malzacher, M., Nies, C., Gottwald, E. and Schad, L.R. (2019), 23Na Triple-quantum signal of in vitro human liver cells, liposomes, and nanoparticles: Cell viability assessment vs. separation of intra- and extracellular signal. J. Magn. Reson. Imaging, 50: 435-444.

3. LaVerde, G., Nemoto, E., Jungreis, C., Tanase, C. and Boada, F. (2007), Serial triple quantum sodium MRI during non-human primate focal brain ischemia. Magn. Reson. Med., 57: 201-205.

4. Hoesl, M.A.U., Schad, L.R., Rapacchi, S.,Volumetric 23Na Single and Triple-Quantum Imaging at 7T: 3D-CRISTINA, Zeitschrift für Medizinische Physik, Volume 32, Issue 2, 2022, Pages 199-208.

5. Kustner, T., Schwartz, M., Würslin, C., Martirosian, P., Schwenzer, N., Yang, B., & Schmidt, H. (2016). Compressed Sensing LAB: An MR acquisition and reconstruction system. In Proceedings of the ISMRM Workshop on Data Sampling and Reconstruction.

6. Shin PJ, Larson PE, Ohliger MA, Elad M, Pauly JM, Vigneron DB, Lustig M. Calibrationless parallel imaging reconstruction based on structured low-rank matrix completion. Magn Reson Med. 2014 Oct;72(4):959-70.

Figures

Fig.1: Construction of data matrix. 5D 23Na MQC k-space data is reshaped into 4D array [x,y,z,TExPC] with the last dimension being TE= multi-echo and PC= phase-cycle. These 3D volumes are vectorized into a structured 2D matrix , $$$A=H(u)$$$, by sliding a kernel (green box) across k-space. Finally, a structured 2D matrix is formed. Exploiting coherent information along TE and PC dimensions simultaneously, by enforcing structured low-rankness, missing data are synthesized. Finally, the low-rank matrix is transformed back to k-space, $$$u=H^*(A_{LR})$$$, for data consistency

Fig.2: Phantom study, (A) used phantom with different NaCl, but consistent agar gel concentration. (B) Image reconstruction comparison of original, 3D CS and SAKE with undersampling factor (USF) of 3. Note the significant SNR improvement and reduced blurring extent in the SAKE reconstruction. (C) Linear regression for SQ intensity, SI, and (D) TQ SI against known NaCl concentration. Both graphs exhibit the same linear pattern, however, for both reconstructions SAKE provides increased accuracy for SQ and TQ signal intensity reconstruction, especially for higher SQ signal.

Fig.3: One slice in-vivo knee. (A) Anatomical PD 1H image of the shown slice. (B) SQ, TQ and TQ/SQ images for original and 3-fold undersampled reconstructions. Note the blurred SQ image when 3D CS is used, while 5D CS offers improved SQ signal reconstruction. (C) RMSE evaluation and (D) SSIM evaluation in the knee for different undersampling factors reveal that 5D CS significantly improves SQ image reconstruction and allows for increased TQ/SQ ratio accuracy. *SQ images were used as reference for SSIM in TQ images . Note the superior structural consistency in the SAKE TQ image.

Fig.4: In-vivo 23Na MQC knee images. The bright vertical lines in the first slice only appear in the SQ, but not TQ images, indicating it to be blood vessels. All SQ as well as TQ images exhibit high signal intensities in the articular cartilage, with the TQ signal intensity being higher in lateral and medial compartments. Note the well preserved anatomical structure in the SAKE reconstructed TQ images, due to exploitation of coherent information along 23Na MQC's multi-dimensional space. The TQ/SQ ratio shows slight intensity drops in the lateral compartment for the undersampled images.

Fig.5: One slice in-vivo brain. (A) Anatomical T2w 1H image of the shown slice. (B) SQ, TQ and TQ/SQ images for original and 2-fold undersampled reconstructions. (C) RMSE evaluation and (D) SSIM evaluation in the brain for different undersampling factors reveal that SAKE enables improved SQ image reconstruction. However, it fails to reconstruct the TQ signal accurately leading to blurred TQ/SQ maps and thus, increased RMSE.

Proc. Intl. Soc. Mag. Reson. Med. 31 (2023)
4722
DOI: https://doi.org/10.58530/2023/4722