1473

R Dispersion imaging in human skeletal muscle at 3 Tesla
Fatemeh Adelnia1, Zhongliang Zu1,2, Feng Wang1,2, Kevin D Harkins3, and John C Gore1,2,3,4,5
1Vanderbilt University Institute of Imaging Science, Nashville, TN, United States, 2Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, United States, 3Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, United States, 4Department of Physics and Astronomy, Vanderbilt University, Nashville, TN, United States, 5Department of Molecular Physiology and Biophysics, Vanderbilt University, Nashville, TN, United States

Synopsis

R dispersion over a range of weak locking fields has the potential to reveal information on microvascular geometry and density such as microvascular spacing. This work presents in-vivo results supporting the application of R dispersion at low locking fields to the measurement of microvascular sizes and spacings in skeletal muscle. We present model fit parameters from measurements of R dispersion of human skeletal muscle that is close to ex-vivo data.

INTRODUCTION:

Measurements of the spin-lattice relaxation rate in the rotating frame (R=1/T), using spin-locking pulse sequences, are sensitive to local magnetic field fluctuations around the frequency of an applied locking pulse.1 R values vary as a function of the magnitude of the locking field frequency (FSL), and measurements of this variation (or dispersion) allow the derivation of intrinsic tissue properties. Recently, it has been shown that R dispersion at high static magnetic fields is dominated by chemical exchange (REx) and diffusion (RDiff) effects.2 Because the time scales of REx and RDiff are so different, these two processes are readily separated, as shown in Figure 1.3 Most in vivo studies of R dispersion have been performed at relatively high FSL, where they emphasize the sensitivity of R dispersion to chemical exchange.4,5 Here we provide in-vivo results supporting the application of R dispersion at lower FSL to quantify the signal dephasing caused by diffusion of tissue water molecules within field gradients produced by local magnetic field inhomogeneities.3,6 We present model fit parameters from R dispersion measurements of human skeletal muscle to estimate vascular spacing.

MATERIALS & METHODS:

Theory of R dispersion at low FSL: Diffusion within spatially varying magnetic field gradients induced by susceptibility variations produces a time-dependent modulation of spin frequencies that results in a net loss of transverse magnetization. These losses may be partially reduced by application of a locking field. For example, if the local field gradient varies sinusoidally with position x, it can be formulated as B0loc=(√2g/q)sin(qx), where g is the mean gradient strength and q is the spatial frequency of the gradient field. The contribution of this local field gradient to the relaxation rate from random diffusion is RDiff2g2D/((q2D)212), where D is the self-diffusion coefficient and ω1 is the FSL.3 In tissues, microscopic field gradients may be produced by, for example, the susceptibility differences between tissue and blood, so that q reflects the effective width of the spatial frequency spectrum of inhomogeneities corresponding to vasculature. Images acquired with different values of ω1 can be combined to specifically portray τc= 1/q2D, a correlation time which is a direct measure of the sizes and spacings of the capillaries, arterioles, and venules. Variations such as those produced by small micro-vessels of dimension d correspond to values of q≈π/d.
Numerical simulation: Monte-Carlo simulations were used to assess the effect of noise and residual unsuppressed fat signal on the bias and dispersion of measured T values in different conditions. The signals were generated using a monoexponential model with an additional term for the residual fat signal; ff=exp(-bDf). Different selections of TSL values were simulated along with the effects of added zero-mean Gaussian noise to vary SNR values.
MRI acquisition: 2D and 3D T images were acquired from the lower extremity skeletal muscle in healthy young subjects, positioned supine, feet first, using a 3T Ingenia-CX MRI scanner (Philips Healthcare). Spin-lock images were acquired with a TSE readout following a 90x-τ/2y-180y-τ/2−y-90x pulse preparation with SPIR fat suppression and PB volume shim.

RESULTS & DISCUSSION:

Fig.2 shows calculated T maps of the leg. There is a slight increase in T values at increasing spin locking fields as predicted. Fig.3 shows the quantitative R dispersion curves and the fitted parameters from selected ROIs in calf and thigh muscles. Both dispersions reveal similar spatial frequency (q) and magnetic gradient strength (g). Moreover, if D=2x10-5cm2/s, then the derived length scale characterizing the intrinsic gradients is d≈10μm which is close to the dimensions of microvasculature obtained from skeletal muscle in animals.7 The bias and dispersion of T were determined for a range of expected tissue properties and acquisition parameters for optimized scan times (see Fig.4a). The measured SNR data depicted in Fig.4b can be matched with the simulation results to estimate the bias and dispersion of the T value at different sets of TSL. Simulation results demonstrate that there is no substantial difference between different sets of TSL when the SNR is higher than 40. Furthermore, shorter TSL minimizes the adverse contributions of residual fat signal and B0-offset frequency. 3D-T images were collected with the optimized series of TSL and a smaller number of FSL in 3 slices from the mid thigh. The measured R dispersion within a clinically achievable scan time shows similar behavior to the 2D-R dispersion measured in the thigh and calf skeletal muscles.

CONCLUSION:

R dispersion imaging at low locking field amplitudes along with appropriate analyses may be applied to derive new types of parametric image information based on diffusion effects. The in-vivo results suggest we can use this approach to measure novel aspects of tissue microstructure including geometrical properties of the microvasculature. However, because of the complexity and contribution of both chemical exchange and diffusion effects to the R dispersion, a definite conclusion cannot be stated, and further investigation such as ex-vivo validation studies should be performed.

Acknowledgements

No acknowledgement found.

References

1. Davis DG, Perlman ME, London RE. Direct measurements of the dissociation-rate constant for inhibitor-enzyme complexes via the T1 rho and T2 (CPMG) methods. J Magn Reson B. Jul 1994;104(3):266-275.

2. Spear JT, Gore JC. New insights into rotating frame relaxation at high field. NMR Biomed. Sep 2016;29(9):1258-1273.

3. Spear JT, Zu Z, Gore JC. Dispersion of relaxation rates in the rotating frame under the action of spin-locking pulses and diffusion in inhomogeneous magnetic fields. Magn Reson Med. May 2014;71(5):1906-1911.

4. Cobb JG, Xie J, Gore JC. Contributions of chemical exchange to T1rho dispersion in a tissue model. Magn Reson Med. Dec 2011;66(6):1563-1571.

5. Wang F, Colvin DC, Wang S, et al. Spin-lock relaxation rate dispersion reveals spatiotemporal changes associated with tubulointerstitial fibrosis in murine kidney. Magn Reson Med. Oct 2020;84(4):2074-2087.

6. Zu Z, Janve V, Gore JC. Spin-lock imaging of intrinsic susceptibility gradients in tumors. Magn Reson Med. May 2020;83(5):1587-1595.

7. Krogh A. The number and distribution of capillaries in muscles with calculations of the oxygen pressure head necessary for supplying the tissue. J Physiol. May 20 1919;52(6):409-415.

Figures

Fig 1. a) Schematic diagram of an imaging voxel illustrating the chemical exchange of labile hydrogen, between water and hydroxyls, occurring in the presence of susceptibility gradients generated by the red blood cells containing paramagnetic deoxyhemoglobin. b) Simulated R dispersions for a spherical structure with a radius of 10 μm. The low-frequency dispersion corresponds to diffusion effects (blue), the high frequency corresponds to exchange (red), and the composite curve “double dispersion” added two rates linearly R = Rdipolar + RDiff + REx (adapted from Ref 3).

Fig 2. a) 2D T images collected with knee coil, TR = 3000 ms, TE = 10 ms, SENSE factor 2, a FOV of 200 x 200 mm, and a voxel size of 1 x1 x 5 mm with 8 TSL; 2, 12, 22, 32, 43, 52, 62, and 72ms resulted in a scan time of 2 min and 27 sec for each FSL. b) 2D T images collected with the flexible anterior coil, TR = 3000 ms, TE = 10 ms, SENSE factor 2, a FOV of 224 x 224 mm, and a voxel size of 1 x1 x 8 mm with 4 TSL; 2, 12, 28, 64 ms resulted in a scan time of 1 min and 30 sec for each FSL. Each map was calculated by fitting the signal decay to a mono-exponential function. ROIs placed in the areas that have approximately zero B0 off-set.

Fig 3. (a) quantitative R dispersion values from the selected ROIs (blue) for the thigh and (orange) for the calf muscle shown in Fig 2. The fitted curves calculated using R equation i.e. R = R2 + γ2g2D/((q2D)212), where ω1 is the FSLs and γ is the gyromagnetic ratio for proton (b) Fitted parameters extracted from R dispersion of skeletal muscles.

Fig 4. (a) Bias and dispersion (%) of T from Monte-Carlo simulations are shown for two different T values corresponding to the beginning and end of T dispersion: 25 and 35 ms. Columns denote different set of TSL: 1st, 2nd, and 3rd. Fat fraction (ff) varied in steps of 2% as shown on the x-axis. The SNR was incremented in different steps as follows: 10, 20, 40, 70, 100, 1000 as shown in the y-axis. (b) SNR map corresponding to slice# 3 of 26 yrs old subject (see Fig.5). (c) T signal decay simulated by solving the Bloch equation and adding a B0 off-set of 20 Hz with the T input value of 35ms.

Fig 5. 3D T images collected with the anterior coil, Slice N=3, TR = 1250 ms, shortest TE, SENSE factor 2, a FOV of 224 x 224 mm, and a voxel size of 1 x1 x 8 mm with TSL; 2, 12, 36 ms and FSL=0, 60, 120, 200, 300, 480 Hz resulted in a total scan time of 7 min and 56 sec. (a) Example of anatomy and R=1/T maps from slice #3. ROI placed in the area that has approximately B0 off-set less than 30Hz. (b) Averaged R values from selected ROI. The error bars show the standard deviation. The fitted curve calculated using R equation (c) Fitted parameters extracted from R dispersion depicted in (b).

Proc. Intl. Soc. Mag. Reson. Med. 29 (2021)
1473