Accelerated multiple-quantum-filtered sodium magnetic resonance imaging using compressed sensing at 7 T

PURPOSE

Multiple-quantum-filtered (MQF) sodium magnetic resonance imaging (MRI), such as enhanced single-quantum and triple-quantum-filtered imaging of 23Na (eSISTINA), enables images to be weighted towards restricted sodium, a promising biomarker in clinical practice, but often suffers from clinically infeasible acquisition times and low image quality. This study aims to mitigate the above limitation by implementing a novel eSISTINA sequence at 7 T with the application of compressed sensing (CS) to accelerate eSISTINA acquisitions without a noticeable loss of information.

METHODS

A novel eSISTINA sequence with a 3D spiral-based sampling scheme was implemented at 7 T for the application of CS. Fully sampled datasets were obtained from one phantom and ten healthy subjects, and were then retrospectively undersampled by various undersampling factors. CS undersampled reconstructions were compared to fully sampled and undersampled nonuniform fast Fourier transform (NUFFT) reconstructions. Reconstruction performance was evaluated based on structural similarity (SSIM), signal-to-noise ratio (SNR), weightings towards total and compartmental sodium, and in vivo quantitative estimates.

RESULTS

CS-based phantom and in vivo images have less noise and better structural delineation while maintaining the weightings towards total, non-restricted (predominantly extracellular), and restricted (primarily intracellular) sodium. CS generally outperforms NUFFT with a higher SNR and a better SSIM, except for the SSIM in TQ brain images, which is likely due to substantial noise contamination. CS enables in vivo quantitative estimates with <15% errors at an undersampling factor of up to two.

CONCLUSIONS

Successful implementation of an eSISTINA sequence with an incoherent sampling scheme at 7 T was demonstrated. CS can accelerate eSISTINA by up to twofold at 7 T with reduced noise levels compared to NUFFT, while maintaining major structural information, reasonable weightings towards total and compartmental sodium, and relatively reliable in vivo quantification. The associated reduction in acquisition time has the potential to facilitate the clinical applicability of MQF sodium MRI.

Recent technical developments and clinical research applications of sodium (23Na) MRI

Sodium is an essential ion that plays a central role in many physiological processes including the transmembrane electrochemical gradient and the maintenance of the body's homeostasis. Due to the crucial role of sodium in the human body, the sodium nucleus is a promising candidate for non-invasively assessing (patho-)physiological changes. Almost 10 years ago, Madelin et al. provided a comprehensive review of methods and applications of sodium (23Na) MRI (Madelin et al., 2014) [1]. More recent review articles have focused mainly on specific applications of 23Na MRI. For example, several articles covered 23Na MRI applications for diseases such as osteoarthritis (Zbyn et al., 2016, Zaric et al., 2020) [2,3], multiple sclerosis (Petracca et al., 2016, Huhn et al., 2019) [4,5] and brain tumors (Schepkin, 2016) [6], or for imaging certain organs such as the kidneys (Zollner et al., 2016) [7], the brain (Shah et al., 2016, Thulborn et al., 2018) [8,9], and the heart (Bottomley, 2016) [10]. Other articles have reviewed technical developments such as radiofrequency (RF) coils for 23Na MRI (Wiggins et al., 2016, Bangerter et al., 2016) [11,12], pulse sequences (Konstandin et al., 2014) [13], image reconstruction methods (Chen et al., 2021) [14], and interleaved/simultaneous imaging techniques (Lopez Kolkovsky et al., 2022) [15]. In addition, 23Na MRI topics have been covered in review articles with broader topics such as multinuclear MRI or ultra-high-field MRI (Niesporek et al., 2019, Hu et al., 2019, Ladd et al., 2018) [16-18]. During the past decade, various research groups have continued working on technical improvements to sodium MRI and have investigated its potential to serve as a diagnostic and prognostic tool. Clinical research applications of 23Na MRI have covered a broad spectrum of diseases, mainly focusing on the brain, cartilage, and skeletal muscle (see Fig. 1). In this article, we aim to provide a comprehensive summary of methodological and hardware developments, as well as a review of various clinical research applications of sodium (23Na) MRI in the last decade (i.e., published from the beginning of 2013 to the end of 2022).

Compressed Sensing in Sodium Magnetic Resonance Imaging: Techniques, Applications, and Future Prospects

Sodium (²³ Na) yields the second strongest nuclear magnetic resonance (NMR) signal in biological tissues and plays a vital role in cell physiology. Sodium magnetic resonance imaging (MRI) can provide insights into cell integrity and tissue viability relative to pathologies without significant anatomical alternations, and thus it is considered to be a potential surrogate biomarker that provides complementary information for standard hydrogen (¹ H) MRI in a noninvasive and quantitative manner. However, sodium MRI suffers from a relatively low signal-to-noise ratio and long acquisition times due to its relatively low NMR sensitivity. Compressed sensing-based (CS-based) methods have been shown to accelerate sodium imaging and/or improve sodium image quality significantly. In this manuscript, the basic concepts of CS and how CS might be applied to improve sodium MRI are described, and the historical milestones of CS-based sodium MRI are briefly presented. Representative advanced techniques and evaluation methods are discussed in detail, followed by an expose of clinical applications in multiple anatomical regions and diseases as well as thoughts and suggestions on potential future research prospects of CS in sodium MRI. EVIDENCE LEVEL: 5 TECHNICAL EFFICACY: Stage 1.

Frontiers of Sodium MRI Revisited: From Cartilage to Brain Imaging

Sodium magnetic resonance imaging (²³ Na-MRI) is a highly promising imaging modality that offers the possibility to noninvasively quantify sodium content in the tissue, one of the most relevant parameters for biochemical investigations. Despite its great potential, due to the intrinsically low signal-to-noise ratio (SNR) of sodium imaging generated by low in vivo sodium concentrations, low gyromagnetic ratio, and substantially shorter relaxation times than for proton (¹ H) imaging, ²³ Na-MRI is extremely challenging. In this article, we aim to provide a comprehensive overview of the literature that has been published in the last 10-15 years and which has demonstrated different technical designs for a range of ²³ Na-MRI methods applicable for disease diagnoses and treatment efficacy evaluations. Currently, a wider use of 3.0T and 7.0T systems provide imaging with the expected increase in SNR and, consequently, an increased image resolution and a reduced scanning time. A great interest in translational research has enlarged the field of sodium MRI applications to almost all parts of the body: articular cartilage tendons, spine, heart, breast, muscle, kidney, and brain, etc., and several pathological conditions, such as tumors, neurological and degenerative diseases, and others. The quantitative parameter, tissue sodium concentration, which reflects changes in intracellular sodium concentration, extracellular sodium concentration, and intra-/extracellular volume fractions is becoming acknowledged as a reliable biomarker. Although the great potential of this technique is evident, there must be steady technical development for ²³ Na-MRI to become a standard imaging tool. The future role of sodium imaging is not to be considered as an alternative to ¹ H MRI, but to provide early, diagnostically valuable information about altered metabolism or tissue function associated with disease genesis and progression. LEVEL OF EVIDENCE: 1 TECHNICAL EFFICACY STAGE: 1.

23 Na-MRI as a Noninvasive Biomarker for Cancer Diagnosis and Prognosis

The influx of sodium (Na⁺) ions into a resting cell is regulated by Na⁺ channels and by Na⁺/H⁺ and Na⁺/Ca²⁺ exchangers, whereas Na⁺ ion efflux is mediated by the activity of Na⁺/K⁺‐ATPase to maintain a high transmembrane Na⁺ ion gradient. Dysfunction of this system leads to changes in the intracellular sodium concentration that promotes cancer metastasis by mediating invasion and migration. In addition, the accumulation of extracellular Na⁺ ions in cancer due to inflammation contributes to tumor immunogenicity. Thus, alterations in the Na⁺ ion concentration may potentially be used as a biomarker for malignant tumor diagnosis and prognosis. However, current limitations in detection technology and a complex tumor microenvironment present significant challenges for the in vivo assessment of Na⁺ concentration in tumor. ²³Na‐magnetic resonance imaging (²³Na‐MRI) offers a unique opportunity to study the effects of Na⁺ ion concentration changes in cancer. Although challenged by a low signal‐to‐noise ratio, the development of ultrahigh magnetic field scanners and specialized sodium acquisition sequences has significantly advanced ²³Na‐MRI. ²³Na‐MRI provides biochemical information that reflects cell viability, structural integrity, and energy metabolism, and has been shown to reveal rapid treatment response at the molecular level before morphological changes occur. Here we review the basis of ²³Na‐MRI technology and discuss its potential as a direct noninvasive in vivo diagnostic and prognostic biomarker for cancer therapy, particularly in cancer immunotherapy. We propose that ²³Na‐MRI is a promising method with a wide range of applications in the tumor immuno‐microenvironment research field and in cancer immunotherapy monitoring.

Level of Evidence

2

Technical Efficacy

Stage 2

Multinuclear MRI at Ultrahigh Fields

In this article, an overview of the current developments and research applications for non-proton magnetic resonance imaging (MRI) at ultrahigh magnetic fields (UHFs) is given. Due to technical and methodical advances, efficient MRI of physiologically relevant nuclei, such as Na, Cl, Cl, K, O, or P has become feasible and is of interest to obtain spatially and temporally resolved information that can be used for biomedical and diagnostic applications. Sodium (Na) MRI is the most widespread multinuclear imaging method with applications ranging over all regions of the human body. Na MRI yields the second largest in vivo NMR signal after the clinically used proton signal (H). However, other nuclei such as O and P (energy metabolism) or Cl and K (cell viability) are used in an increasing number of MRI studies at UHF. One major advancement has been the increased availability of whole-body MR scanners with UHFs (B0 ≥7T) expanding the range of detectable nuclei. Nevertheless, efforts in terms of pulse sequence and post-processing developments as well as hardware designs must be made to obtain valuable information in clinically feasible measurement times. This review summarizes the available methods in the field of non-proton UHF MRI, especially for Na MRI, as well as introduces potential applications in clinical research.

Potential of Sodium MRI as a Biomarker for Neurodegeneration and Neuroinflammation in Multiple Sclerosis

In multiple sclerosis (MS), experimental and ex vivo studies indicate that pathologic intra- and extracellular sodium accumulation may play a pivotal role in inflammatory as well as neurodegenerative processes. Yet, in vivo assessment of sodium in the microenvironment is hard to achieve. Here, sodium magnetic resonance imaging (²³NaMRI) with its non-invasive properties offers a unique opportunity to further elucidate the effects of sodium disequilibrium in MS pathology in vivo in addition to regular proton based MRI. However, unfavorable physical properties and low in vivo concentrations of sodium ions resulting in low signal-to-noise-ratio (SNR) as well as low spatial resolution resulting in partial volume effects limited the application of ²³NaMRI. With the recent advent of high-field MRI scanners and more sophisticated sodium MRI acquisition techniques enabling better resolution and higher SNR, ²³NaMRI revived. These studies revealed pathologic total sodium concentrations in MS brains now even allowing for the (partial) differentiation of intra- and extracellular sodium accumulation. Within this review we (1) demonstrate the physical basis and imaging techniques of ²³NaMRI and (2) analyze the present and future clinical application of ²³NaMRI focusing on the field of MS thus highlighting its potential as biomarker for neuroinflammation and -degeneration.

Double quantum filtered 23 Na MRI with magic angle excitation of human skeletal muscle in the presence of B0 and B1 inhomogeneities

Double quantum filtered ²³ Na MRI with magic angle excitation (DQF-MA) can be used to selectively detect sodium ions located within anisotropic structures such as muscle fibers. It might therefore be a promising tool to analyze the microscopic environment of sodium ions, for example in the context of osmotically neutral sodium retention. However, DQF-MA imaging is challenging due to various signal dependences, on both measurement parameters and external influences. The aim of this work was to examine how B₀ in combination with B₁ inhomogeneities alter the DQF-MA signal intensity. We showed that, in the presence of B₀ inhomogeneities, flip angle schemes with only one 54.7° pulse can be favorable compared with the classical 90°-54.7°-54.7° scheme. DQF-MA images of the human lower leg were acquired at B₀  = 3 T with a nominal spatial resolution of 12 × 12 × 36 mm³ within an acquisition time of T_Acq  < 10 min, and compared with spin density weighted (DW), as well as triple quantum filtration (TQF) ²³ Na images. We found mean normalized signal-to-noise ratios of TQF/DW = 13.7 ± 2.3% (tibialis anterior), 11.9 ± 2.3% (soleus) and 11.4 ± 2.2% (gastrocnemius medialis), as well as DQF-MA/DW = 4.7 ± 1.1% (tibialis anterior), 3.3 ± 0.73% (soleus) and 3.4 ± 0.6% (gastrocnemius medialis). These ratios might serve as additional measures in future clinical studies of sodium retention within human skeletal muscle. However, the influence of B₀ and B₁ inhomogeneities should be considered when interpreting DQF-MA images.

Acute changes in extracellular volume fraction in skeletal muscle monitored by 23Na NMR spectroscopy

In this article, we induced acute changes in extracellular volume fraction in skeletal muscle tissue and compared the sensitivity of a standard ¹H T₂ imaging method with different ²³Na‐NMR spectroscopy parameters within acquisition times compatible with clinical investigations. First, we analyzed the effect of a short ischemia on the sodium distribution in the skeletal muscle. Then, the lower leg of 21 healthy volunteers was scanned under different vascular filling conditions (vascular draining, filling, and normal condition) expected to modify exclusively the extracellular volume. The first experiment showed no change in the total sodium content during a 15 min ischemia, but the intracellular weighted ²³Na signal slowly decreased. For the second part, significant variations of total sodium content, sodium distribution, and T₁ and T2∗ of ²³Na signal were observed between different vascular filling conditions. The measured sodium distribution correlates significantly with sodium T₁ and with the short and long T2∗ fractions. In contrast, significant changes in the proton T₂w signal were observed only in three muscles. Altogether, the mean T₂w signal intensity of all muscles as well as their mean T₂ did not vary significantly with the extracellular volume changes. In conclusion, at the expense of giving up spatial resolution, the proposed ²³Na spectroscopic method proved to be more sensitive than standard ¹H T₂ approach to monitor acute extracellular compartment changes within muscle tissue.

Low-power suppression of fast-motion spin 3/2 signals

Triple Quantum Filters (TQFs) are frequently used for the selection of bi-exponentially relaxing spin 3/2 nuclei (in particular ²³Na) in ordered environments, such as biological tissues. These methods provide an excellent selection of slow-motion spins, but their sensitivity is generally low, and coherence selection requirements may lead to long experiments when applied in vivo. Alternative methods, such as 2P DIM, have demonstrated that the sensitivities of the signals from bi-exponentially relaxing sodium can be significantly increased using strategies other than TQFs. A shortcoming of the 2P DIM method is its strong dependence on B₀ inhomogeneities. We describe here a method, which is sensitive to the slow-motion regime, while the signal from spins in the fast-motion regime is suppressed. This method is shown to be more effective than TQFs, requires minimal phase cycling for the suppression of the influence of rf inhomogeneity, and has less dependence on resonance offsets and B₀-inhomogeneity than 2P DIM.

Quantitative Sodium MR Imaging at 7 T: Initial Results and Comparison with Diffusion-weighted Imaging in Patients with Breast Tumors

Purpose To investigate the clinical feasibility of a quantitative sodium 23 ((23)Na) magnetic resonance (MR) imaging protocol developed for breast tumor assessment and to compare it with 7-T diffusion-weighted imaging (DWI). Materials and Methods Written informed consent in this institutional review board-approved study was obtained from eight healthy volunteers and 17 patients with 20 breast tumors (five benign, 15 malignant). To achieve the best image quality and reproducibility, the (23)Na sequence was optimized and tested on phantoms and healthy volunteers. For in vivo quantification of absolute tissue sodium concentration (TSC), an external phantom was used. Static magnetic field, or B0, and combined transmit and receive radiofrequency field, or B1, maps were acquired, and image quality, measurement reproducibility, and accuracy testing were performed. Bilateral (23)Na and DWI sequences were performed before contrast material-enhanced MR imaging in patients with breast tumors. TSC and apparent diffusion coefficient (ADC) were calculated and correlated for healthy glandular tissue and benign and malignant lesions. Results The (23)Na MR imaging protocol is feasible, with 1.5-mm in-plane resolution and 16-minute imaging time. Good image quality was achieved, with high reproducibility (mean TSC values ± standard deviation for the test, 36 mmol per kilogram of wet weight ± 2 [range, 34-37 mmol/kg]; for the retest, 37 mmol/kg ± 1 [range, 35-39 mmol/kg]; P = .610) and accuracy (r = 0.998, P < .001). TSC values in normal glandular and adipose breast tissue were 35 mmol/kg ± 3 and 18 mmol/kg ± 3, respectively. In malignant lesions (mean size, 31 mm ± 24; range, 6-92 mm), the TSC of 69 mmol/kg ± 10 was, on average, 49% higher than that in benign lesions (mean size, 14 mm ± 12; range, 6-35 mm), with a TSC of 47 mmol/kg ± 8 (P = .002). There were similar ADC differences between benign ([1.78 ± 0.23] × 10(-3) mm(2)/sec) and malignant ([1.03 ± 0.23] × 10(-3) mm(2)/sec) tumors (P = .002). ADC and TSC were inversely correlated (r = -0.881, P < .001). Conclusion Quantitative (23)Na MR imaging is clinically feasible, may provide good differentiation between malignant and benign breast lesions, and demonstrates an inverse correlation with ADC. (©) RSNA, 2016 Online supplemental material is available for this article.

Sodium MRI of the human heart at 7.0 T: preliminary results

The objective of this work was to examine the feasibility of three-dimensional (3D) and whole heart coverage (23)Na cardiac MRI at 7.0 T including single-cardiac-phase and cinematic (cine) regimes. A four-channel transceiver RF coil array tailored for (23)Na MRI of the heart at 7.0 T (f = 78.5 MHz) is proposed. An integrated bow-tie antenna building block is used for (1)H MR to support shimming, localization and planning in a clinical workflow. Signal absorption rate simulations and assessment of RF power deposition were performed to meet the RF safety requirements. (23) Na cardiac MR was conducted in an in vivo feasibility study. 3D gradient echo (GRE) imaging in conjunction with Cartesian phase encoding (total acquisition time T(AQ)  = 6 min 16 s) and whole heart coverage imaging employing a density-adapted 3D radial acquisition technique (T(AQ)  = 18 min 20 s) were used. For 3D GRE-based (23)Na MRI, acquisition of standard views of the heart using a nominal in-plane resolution of (5.0 × 5.0) mm(2) and a slice thickness of 15 mm were feasible. For whole heart coverage 3D density-adapted radial (23)Na acquisitions a nominal isotropic spatial resolution of 6 mm was accomplished. This improvement versus 3D conventional GRE acquisitions reduced partial volume effects along the slice direction and enabled retrospective image reconstruction of standard or arbitrary views of the heart. Sodium cine imaging capabilities were achieved with the proposed RF coil configuration in conjunction with 3D radial acquisitions and cardiac gating. Cardiac-gated reconstruction provided an enhancement in blood-myocardium contrast of 20% versus the same data reconstructed without cardiac gating. The proposed transceiver array enables (23)Na MR of the human heart at 7.0 T within clinical acceptable scan times. This capability is in positive alignment with the needs of explorations that are designed to examine the potential of (23)Na MRI for the assessment of cardiovascular and metabolic diseases.

Anatomically weighted second-order total variation reconstruction of 23Na MRI using prior information from 1H MRI

Sodium ((23)Na) MRI is a noninvasive tool to assess cell viability, which is linked to the total tissue sodium concentration (TSC). However, due to low in vivo concentrations, (23)Na MRI suffers from low signal-to-noise ratio (SNR) and limited spatial resolution. As a result, image quality is compromised by Gibbs ringing artifacts and partial volume effects. An iterative reconstruction algorithm that incorporates prior information from (1)H MRI is developed to reduce partial volume effects and to increase the SNR in non-proton MRI. Anatomically weighted second-order total variation (AnaWeTV) is proposed as a constraint for compressed sensing reconstruction of 3D projection reconstruction (3DPR) data. The method is evaluated in simulations and a MR measurement of a multiple sclerosis (MS) patient by comparing it to gridding and other reconstruction techniques. AnaWeTV increases resolution of known structures and reduces partial volume effects. In simulated MR brain data (nominal resolution Δx(3) = 3 × 3 × 3 mm(3)), the intensity error of four small MS lesions was reduced from (6.9 ± 3.8)% (gridding) to (2.8 ± 1.4)% (AnaWeTV with T2-weighted reference images). Compared to gridding, a substantial SNR increase of 130% was found in the white matter of the MS patient. The algorithm is robust against misalignment of the prior information on the order of the (23)Na image resolution. Features without prior information are still reconstructed with high contrast. AnaWeTV allows a more precise quantification of TSC in structures with prior knowledge. Thus, the AnaWeTV algorithm is in particular beneficial for the assessment of tissue structures that are visible in both (23)Na and (1)H MRI.

Short-T2 imaging for quantifying concentration of sodium (23 Na) of bi-exponential T2 relaxation

PURPOSE

This work intends to demonstrate a new method for quantifying concentration of sodium (²³ Na) of bi-exponential T₂ relaxation in patients on MRI scanners at 3.0 Tesla.

THEORY AND METHODS

Two single-quantum (SQ) sodium images acquired at very-short and short echo times (TE = 0.5 and 5.0 ms) are subtracted to produce an image of the short-T₂ component of the bi-exponential (or bound) sodium. An integrated calibration on the SQ and short-T₂ images quantifies both total and bound sodium concentrations. Numerical models were used to evaluate signal response of the proposed method to the short-T₂ components. MRI scans on agar phantoms and brain tumor patients were performed to assess accuracy and performance of the proposed method, in comparison with a conventional method of triple-quantum filtering.

RESULTS

A good linear relation (R²  = 0.98) was attained between the short-T₂ image intensity and concentration of bound sodium. A reduced total scan time of 22 min was achieved under the SAR restriction for human studies in quantifying both total and bound sodium concentrations.

CONCLUSION

The proposed method is feasible for quantifying bound sodium concentration in routine clinical settings at 3.0 Tesla. Magn Reson Med 74:162-174, 2015. © 2014 Wiley Periodicals, Inc.

Two-pulse biexponential-weighted 23Na imaging

A new method is proposed for acquiring 3D biexponential-weighted sodium images with two instead of three RF pulses to allow for shorter repetition time at high magnetic fields (B0≥7 T) and reduced SAR. The second pulse converts single- into triple-quantum coherences in regions containing sodium ions which are restricted in mobility. Since only single-quantum coherences can be detected, an image acquired after the second pulse is intrinsically single-quantum-filtered and can be used to generate a biexponential-weighted sodium image by a weighted subtraction with the spin-density-weighted image acquired between the pulses. The proposed sequence generates biexponential-weighted sodium images of in vivo human brain with 140% higher SNR than triple-quantum-filtered sodium images and 4% higher SNR than a biexponential-weighted sequence with three RF pulses at equal acquisition time and with 1/3 lower SAR. As SAR is reduced, accordingly repetition time can be spared to obtain even higher SNR-time efficiency. In comparison to a difference image generated from two images of a double-readout sequence, the proposed two-pulse sequence yields about 14% higher SNR. Our new two-pulse biexponential-weighted sequence allows for acquisition of full 3D data sets of the human brain in vivo with a nominal resolution of (5 mm)(3) in about 10 min.

Measurement techniques for magnetic resonance imaging of fast relaxing nuclei

In this review article, techniques for sodium ((23)Na) magnetic resonance imaging (MRI) are presented. These techniques can also be used to image other nuclei with short relaxation times (e.g., (39)K, (35)Cl, (17)O). Twisted projection imaging, density-adapted 3D projection reconstruction, and 3D cones are preferred because of uniform k-space sampling and ultra-short echo times. Sampling density weighted apodization can be applied if intrinsic filtering is desired. This approach leads to an increased signal-to-noise ratio compared to postfiltered acquisition in cases of short readout durations relative to T 2 (*) relaxation time. Different MR approaches for anisotropic resolution are presented, which are important for imaging of thin structures such as myocardium, cartilage, and skin. The third part of this review article describes different methods to put more weighting either on the intracellular or the extracellular sodium signal by means of contrast agents, relaxation-weighted imaging, or multiple-quantum filtering.