4003

Optimization of spin-lock preparation pulses for B1 and B0 insensitive cardiac T1ρ mapping
Haikun Qi1,2, Zhenfeng Lv1,2, Jian Xu3, and Peng Hu1,2
1School of Biomedical Engineering, ShanghaiTech University, Shanghai, China, 2Shanghai Clinical Research and Trial Center, Shanghai, China, 3UIH America, Inc., Houston, TX, United States

Synopsis

Keywords: Myocardium, Tissue Characterization

Cardiac T1ρ mapping is a promising technique for assessment of myocardial fibrosis without exogenous contrast agent. However, its wide application is hindered by the sensitivity of T1ρ preparation to B1 and B0 inhomogeneities. In this study, the state-of-the-art constant spin-lock methods including the composite and adiabatic excitation continuous-wave spin-lock methods were investigated using numerical simulations to assess their robustness to field inhomogeneities, and validated in phantoms and a preliminary subject. Two T1ρ preparation modules were found to generate superior T1ρ mapping quality in the presences of B0 and B1 inhomogeneities, indicating the potential of clinical application of cardiac T1ρ MRI.

Introduction

Cardiac T1ρ mapping is an emerging technique for assessment of focal and diffuse myocardial fibrosis without exogenous contrast agent 1-5. However, T1ρ CMR is challenging due to the sensitivity of T1ρ preparation (T1ρ prep) to B1 RF and B0 inhomogeneities. Various spin-lock techniques have been proposed to mitigate the influence of B1 and/or B0 field inhomogeneities, including the composite spin-lock (CSL) pulses 6,7 and the adiabatically excited continuous-wave spin-lock (ACW-SL) method 8. Especially, the ACW-SL has recently shown promising results for cardiac T1ρ mapping at 3T 9, while the key component, adiabatic excitation pulse, has not been systematically optimized. In this study, we optimized the pulse parameters for two commonly used adiabatic excitation pulses using numerical simulations and compared the optimized ACW-SL with the state-of-the-art composite spin-lock methods in phantoms and in vivo cardiac T1ρ mapping.

Methods

Four T1ρ prep modules were considered (Fig. 1): the composite spin-lock methods 7 with one (CSL-1) or two (CSL-2) 180° refocusing pulses for simultaneous B1 and B0 inhomogeneity compensation; the ACW-SL which consists of the adiabatic half-passage (AHP) tip-down pulse, the spin-lock pulse and the reverse AHP tip-up pulse (rAHP); the composite ACW-SL which has an additional refocusing pulse. The hyperbolic secant (HS) and tangent/hyperbolic tangent (tan/tanh) pulses were considered for the AHP. The amplitude and phase formulations of the two pulses are provided in Table 1. The pulse parameters were optimized in the following Bloch simulations.
Numerical simulations: Bloch simulations were performed for B0 inhomogeneity ranging from -200Hz to 200Hz and for B1 inhomogeneity ranging from -50% to 50% 10. For each combination of B0 and B1 inhomogeneities, the longitudinal magnetizations (Mz) were calculated for a range of spin-lock durations (TSL) from 0 to T1ρ. The spin-lock frequency was set to 350Hz 9. The myocardial T1ρ was to 50ms, while T2ρ was assumed to be 75ms 11. The residual sum of squares (RSS) between Mz and the mono-exponentially decayed signal Mz,d was calculated to assess the banding susceptibility: $$$RSS=\sum_n(M_{z,d}(n)-M_{z}(n))^{2}$$$. For quantification error (ΔE) analysis, the simulated longitudinal magnetizations were divided into four equally spaced segments and one sample was randomly selected in each segment for T1ρ fitting. The random T1ρ fitting process was repeated for 100 times to calculate the mean quantification error.
Optimization of the AHP pulse: The AHP pulse duration was set to 4ms for a tradeoff between satisfying the adiabatic condition and reducing relaxation during the AHP. Table 1 summarizes the optimization settings for the HS and tan/tanh pulses. For each pulse design considered, the RSS was calculated for 30 frequencies across ±150 Hz and 20 B1 field amplitudes between 70% and 130% of the maximum constraint. The brute-force search selected the set of design parameters that had the lowest RSS.
Imaging Experiments: The T1ρ prep methods of CSL-1, CSL-2 and ACW-SL with the optimized AHP pulse were implemented in a 3T United Imaging MR system, and compared in phantoms made of agarose and gadolinium contrast and in one healthy volunteer. The 2D T1ρ prep bSSFP 9 was performed with spin-lock frequency of 350Hz and TSL of 2, 16, 30 and 50ms under a single breath-hold of 10 heartbeats (two dummy cardiac cycles between TSL acquisitions). Volume shimming mode was adopted.

Results

Seeing from the RSS and ΔE maps in Fig. 2, the CSL-1 and ACW-SL with the optimized tan/tanh pulse are the two most robust methods, immune to a wide range of field inhomogeneities. The introduction of additional refocusing pulse to the ACW-SL degraded its performance. Figure 3 and 4 respectively shows the phantom and in vivo T1ρ mapping results and the B1 and B0 field maps. There are severe B0 inhomogeneities in the top and bottom phantoms, while the CSL-1 and tan/tanh-SL provided more reliable T1ρ-weighted images than the CSL-2, as evident by the T1ρ and RSS maps.

Discussion

The state-of-the-art constant spin-lock methods were investigated using numerical simulations, and validated in phantoms and a preliminary subject. The CSL-1 and optimized tan/tanh-SL methods are able to generate superior T1ρ mapping quality in the presences of B0 and B1 inhomogeneities. The optimized B1 amplitude for the tan/tanh pulse is the same to the spin-lock pulse, which echoes the previous finding 8 that the ACW-SL is robust to both B1 and B0 inhomogeneities when the B1 amplitude of the AHP and the rAHP pulse is the same to that of the spin-lock pulse.

Acknowledgements

No acknowledgement found.

References

1. Bustin A, Toupin S, Sridi S, Yerly J, Bernus O, Labrousse L, Quesson B, Rogier J, Haissaguerre M, van Heeswijk R, Jais P, Cochet H, Stuber M. Endogenous assessment of myocardial injury with single-shot model-based non-rigid motion-corrected T1 rho mapping. J Cardiovasc Magn Reson 2021;23(1):119.

2. Han Y, Liimatainen T, Gorman RC, Witschey WR. Assessing Myocardial Disease Using T1rho MRI. Curr Cardiovasc Imaging Rep 2014;7(2):9248.

3. Qi H, Bustin A, Kuestner T, Hajhosseiny R, Cruz G, Kunze K, Neji R, Botnar RM, Prieto C. Respiratory motion-compensated high-resolution 3D whole-heart T1rho mapping. J Cardiovasc Magn Reson 2020;22(1):12.

4. Thompson EW, Kamesh Iyer S, Solomon MP, Li Z, Zhang Q, Piechnik S, Werys K, Swago S, Moon BF, Rodgers ZB, Hall A, Kumar R, Reza N, Kim J, Jamil A, Desjardins B, Litt H, Owens A, Witschey WRT, Han Y. Endogenous T1rho cardiovascular magnetic resonance in hypertrophic cardiomyopathy. J Cardiovasc Magn Reson 2021;23(1):120.

5. Witschey WR, Zsido GA, Koomalsingh K, Kondo N, Minakawa M, Shuto T, McGarvey JR, Levack MM, Contijoch F, Pilla JJ, Gorman JH, 3rd, Gorman RC. In vivo chronic myocardial infarction characterization by spin locked cardiovascular magnetic resonance. J Cardiovasc Magn Reson 2012;14:37.

6. Gram M, Seethaler M, Gensler D, Oberberger J, Jakob PM, Nordbeck P. Balanced spin-lock preparation for B1 -insensitive and B0 -insensitive quantification of the rotating frame relaxation time T1rho. Magn Reson Med 2021;85(5):2771-2780.

7. Witschey WRT, Borthakur A, Elliott MA, Mellon E, Niyogi S, Wallman DJ, Wang CY, Reddy R. Artifacts in T1p-weighted imaging: Compensation for B-1 and B-0 field imperfections. J Magn Reson 2007;186(1):75-85.

8. Chen W. Artifacts correction for T1rho imaging with constant amplitude spin-lock. J Magn Reson 2017;274:13-23.

9. Qi H, Lv Z, Hu J, Xu J, Botnar R, Prieto C, Hu P. Accelerated 3D free-breathing high-resolution myocardial T1rho mapping at 3 Tesla. Magn Reson Med 2022;88(6):2520-2531.

10. Kellman P, Herzka DA, Hansen MS. Adiabatic Inversion Pulses for Myocardial T1 Mapping. Magnetic Resonance in Medicine 2014;71(4):1428-1434.

11. Wheaton AJ, Borthakur A, Kneeland JB, Regatte RR, Akella SV, Reddy R. In vivo quantification of T1rho using a multislice spin-lock pulse sequence. Magn Reson Med 2004;52(6):1453-1458.

Figures

Fig. 1: A, B: the composite spin-lock methods with one (CSL-1) or two (CSL-2) refocusing pulses. C: the adiabatically excited continuous-wave spin-lock (ACW-SL) module. D: the composite ACW-SL with one additional refocusing pulse.

Table 1: Amplitude and frequency formulations of the HS and tan/tanh pulses, as well as the settings of design parameters used in the optimization process.

Fig. 2: The RSS and ΔE maps for the five spin-lock methods. The optimized pulse parameters were used in the formulations of HS and tan/tanh given in Table 1. HS: B1,max, 14.1uT; ∆𝑓, 700Hz; 𝜇, 2.7. tan/tanh: B1,max, 8.2uT; fmax, 3000Hz; 𝜉, 20; tan(k), 38.

Fig. 3: Phantom results. A: The B0 and B1 field maps. B: The T1ρ-weighted images, T1ρ and RSS maps for the three spin-lock methods of CSL-2, CSL-1 and tan/tanh-SL.

Fig. 4: Results of a healthy subject with heart rate of ~80bpm. A: The B0 and B1 field maps. B: The T1ρ-weighted images, T1ρ and RSS maps for the three spin-lock methods of CSL-2, CSL-1 and tan/tanh-SL.

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