3044

Simulation of Hyperpolarized Perfusion MRI with a Segmented Snapshot Acquisition
Keith A Michel1, Christopher M Walker1, Yunyun Chen2, Jorge Delacerda1, Stephen Lai2, and James A Bankson1

1Department of Imaging Physics, UT MD Anderson Cancer Center, Houston, TX, United States, 2Department of Head & Neck Surgery, UT MD Anderson Cancer Center, Houston, TX, United States

Synopsis

Metabolically inert hyperpolarized tracers can be used for perfusion measurements in vivo. In this work we present a simulation study to compare the performance for measurement of tissue perfusion of a conventional imaging approach to a segmented snapshot acquisition where each segment ends with a 90 degree RF pulse.

Introduction

The hyperpolarization of metabolically inert molecules provides unique possibilities for quantitative perfusion imaging1. The magnitude of the hyperpolarized (HP) signal from such molecules provides a direct signal source that is virtually background-free and therefore not dependent upon the changes in the proton relaxation rates of tissues that underlie the quantification of physiological parameters in other methods of perfusion MRI. Furthermore, hyperpolarized MRI allows for the selection of imaging sequence parameters that can precisely control the observable signal of the perfusion agent in space and time. These acquisition parameters must be chosen very carefully since the HP signal is non-renewable and is depleted by RF excitation. A kinetic model for the evolution of the HP signal is therefore needed to evaluate and optimize imaging strategies for accurate and efficient quantification of tissue perfusion with HP agents in vivo2.

One unique advantage of HP MRI is the ability to expunge longitudinal magnetization within a slice with a 90 degree RF pulse. As Figure 1 shows, longitudinal magnetization may be replenished by delivery of fresh HP spins in blood, and one technique to exploit this potential is a segmented acquisition wherein a 90 degree pulse is applied at the end of each segment of repetitions. This technique may allow for improved measurements of perfusion parameters since the uncertainty in the amount of residual HP spins within the extravascular space immediately following a perfect 90 degree pulse is significantly reduced. Here we employ a kinetic model for perfusion imaging with HP agents to compare the performance of a conventional snapshot acquisition to the segmented snapshot acquisition described above.

Methods

Tissue perfusion was modeled using a two-compartment kinetic model with reversible exchange:

$$X \left( t \right) = e^{ A \left( t - \Delta \right) } v_e X_e \left( \Delta \right) + k_{ve} \int_{\Delta}^{t} e^{ A \left( t - \tau \right)} X_b \left( \tau \right) d \tau + v_b X_b \left( t \right)$$

$$A= - \left( \frac{k_{ve}}{v_e} + R_X \right)$$

$$M_{xy} \left( t + \delta \right) = M_z \left( t \right) \sin \left( \alpha \right) \propto X \left( t \right) \sin \left( \alpha \right)$$

$$M_{z} \left( t + \delta \right) = M_z \left( t \right) \cos \left( \alpha \right) \propto X \left( t \right) \cos \left( \alpha \right)$$

where X is the total amount of HP agent; Xe and Xb are the amounts of HP agent, and ve and vb are volume fractions for the extravascular and blood compartments, respectively; kve is the rate of spin exchange; RX is the rate of HP spin loss to factors other than reverse exchange and RF excitation; Mxy and Mz are the transverse and longitudinal magnetizations, respectively; α is the RF tip angle; Δ is the time at which the previous RF excitation ends and δ is the RF pulse duration. This model is similar to the well-established extended Tofts model and models proposed previously for HP pyruvate2, but without chemical exchange of HP spins.

The accuracy of estimates for kve, vb and ve was assessed for both acquisition strategies depicted in Figure 1 with a range of repetition times and flip angles by generating simulated HP time curves, adding a fixed level of noise and fitting the perfusion model to the noisy curves with these three perfusion parameters as fit parameters using a trust-region reflective method. Ground truth values of 0.02 s-1, 0.09 and 0.2 were utilized for kve, vb and ve, respectively.

Within the kinetic model presented above, the VIF is assumed to be unperturbed by RF excitation, however when 90 degree pulses are applied to thick slices HP spins in major vessels may be lost, severely limiting the imaging window. We investigated this effect in an orthotopic mouse model of thyroid cancer by acquiring single-shot EPI images of HP 13C urea for a 1 cm slice covering the tumor and major vessels in the neck with a TR of 2 s and a flip angle of 90 degrees.

Results

The simulation results show improvement in the accuracy of estimates for all three perfusion parameters obtained at low flip angles and repetition times with the segmented acquisition compared to the conventional approach. The in vivo images show a duration of HP urea signal that is brief but sufficiently long to support the 30 s duration of our simulations.

Discussion & Conclusion

The improvement in perfusion parameter accuracy for the segmented acquisition at low flip angles could preserve valuable HP signal for longer scan times or multislice imaging. Overall, this segmented approach to perfusion imaging holds promise to provide superior perfusion information compared to conventional acquisitions.

Acknowledgements

This work was supported in part by the by the National Institutes of Health (P30-CA016672) and the Cancer Prevention & Research Institute of Texas (RP140021-P5).

References

1. Ardenkjær-Larsen JH, et al. Increase in signal-to-noise ratio of > 10,000 times in liquid-state NMR. Proc Natl Acad Sci USA. 100(18), 10158-10163. doi:10.1073/pnas.1733835100 (2003).

2. Bankson, J. A. et al. Kinetic Modeling and Constrained Reconstruction of Hyperpolarized [1-13C]-Pyruvate Offers Improved Metabolic Imaging of Tumors. Cancer Res. 75, 4708-4717. doi:10.1158/0008-5472.CAN-15-0171 (2015).

Figures

Figure 1: Simulated HP perfusion agent time curves generated with a gamma-variate vascular input function for a conventional acquisition with TR of 1 s, flip angle of 20 degrees, kve of 0.02 s-1, vb of 0.09 and ve of 0.2 (Left), and for a segmented acquisition of three 10 s segments each ending in a 90 degree pulse, otherwise using the same acquisition and physiological parameters (Right). Times at which 90 degree pulses occur are indicated with dotted lines.

Figure 2: Absolute percent error as a function of TR and flip angle for physiological parameters fit to simulated HP agent time curves with noise for a conventional acquisition (Top: A-C) and for a segmented acquisition consisting of three 10 s segments each ending in a 90 degree pulse (Middle: D-F). Relative change in absolute percent error for physiological parameters fit to simulated noisy data for the conventional acquisition relative to the otherwise equivalent segmented approach (Bottom: G-I).

Figure 3: Averaged HP urea image overlaid on anatomical T2-weighted image of a mouse with orthotopic thyroid tumor (Left) and HP signal time curves for single voxels in the tumor and a blood vessel on the contralateral side of the animal’s neck. HP Urea images were obtained using a multi-repetition single-shot EPI sequence with a TR of 2 s and flip angle of 90 degrees.

Proc. Intl. Soc. Mag. Reson. Med. 25 (2017)
3044