In order to use the transient signal of hyperpolarized tracers and their metabolites efficiently, dedicated imaging sequences are required. Here, we present a multi-echo bSSFP sequence with Dixon-based iterative reconstruction to obtain metabolite maps of hyperpolarized [1‑13C]pyruvate and the product of an enzymatic conversion [1-13C]lactate on a human 3T PET-MRI system in vitro and in vivo. When comparing to other methods (i.e. CSI and non-localized NMR spectra) we found that me-bSSFP provides good metabolite separation and reliable quantitative kinetic data more than 16 times faster than CSI (350 ms vs. 5.8 s), while consuming a similar amount of hyperpolarized magnetization.
[1-13C]pyruvate (HP-py) hyperpolarized with by DNP (Hypersense, Oxford Instruments) was dissolved in 100 mM PBS buffer and added either to 5 mL glass vials or 5 mm NMR tubes filled with aqueous and buffered solutions containing constant amounts of nicotinamide adenine dinucleotide (NADH) and increasing amounts of LDH according to the scheme in Table 1. For in vivo experiments, HP-py was injected approx. 24 s after dissolution via tail-vein into MAT-B-III bearing rats (female Fischer344, 1 mL, 90 mM HP-py, pH7.3)).
NMR 13C-spectroscopy was performed at a benchtop NMR spectrometer (Spinsolve, Magritek) with a pulse-acquire sequence. 13C-MRI with me-bSSFP and CSI (Table 2) was conducted on a clinical 3T PET-MRI system (Biograph, Siemens) equipped with a 1H-13C surface coil (Rapid Biomedical). A thermally polarized model solution (MS) containing [1-13C]lactate (TP-lac): 2.23 M, 0.15%(w/v) NaN3, was used for quality control and adjustments.
Zero-filling (2x) and Gaussian windowing was applied to CSI and me-bSSFP. Me-bSSFP images were reconstructed with a Dixon based iterative reconstruction4–6, and CSI data was Fourier transformed in time domain and peaks were integrated, to obtain frequency-selective metabolite maps at -186 Hz (HP-py) and +186 Hz (HP-lac). Signal-to-noise (SNR) ratios were calculated by dividing the mean signal by the standard deviation of selected regions of interest (ROI). NMR spectra were apodized by 5 Hz and quantified by integration of HP-py and HP-lac peaks. The py→lac conversion was quantified as area under the curve (AUC) ratio7 by dividing the AUC of the lactate signal by that of pyruvate.
In vitro assays with increasing LDH concentrations were performed with time-resolved 1D NMR (control) and time- and spatially-resolved me-bSSFP and CSI. We were able to obtain a two-fold better spatial resolution and a six-fold higher temporal resolution with me-bSSFP than with CSI (0.35s vs. 5.8s per image). Although me-bSSFP used a two-fold-higher flip angle, SNR after 50 s acquisition and later was better than with CSI (Figure 1).
All methods exhibited a linear response of the AUC ratio to LDH concentration (Figure 2). However, AUC ratios acquired with me-bSSFP did not show a significant difference to control measurements with NMR, while this was not the case for CSI. This observation is also apparent by eye in the dynamic data (Figure 2).
In vivo time-resolved axial metabolic maps of tumor-bearing rats were successfully acquired with CSI and me-bSSFP after the injection of HP-py at a clinical PET/MR scanner (Figure 3). As observed for in vitro experiments, me-bSSFP provided higher temporal and spatial resolution. HP-py and HP-lac signals are inhomogeneously distributed in the tumor, suggesting tumor heterogeneity. Spatial separation of tumor and healthy tissue was difficult for CSI because of the limited spatial resolution. The difference of the AUC ratio measured in vivo with me-bSSFP and CSI (0.34 and 0.17) replicates the discrepancy of the quantitative results from in vitro.
1. Leupold, J. et al. Fast chemical shift mapping with multiecho balanced SSFP. Magn. Reson. Mater. Phys. Biol. Med. 19, 267–273 (2006).
2. Leupold, J. et al. Fast multiecho balanced SSFP metabolite mapping of 1H and hyperpolarized 13C compounds. Magn. Reson. Mater. Phys. Biol. Med. 22, 251–256 (2009).
3. Abragam, A. & Goldman, M. Principles of dynamic nuclear polarisation. Rep. Prog. Phys. 41, 395 (1978).
4. Reeder, S. B. et al. Multicoil Dixon chemical species separation with an iterative least-squares estimation method. Magn. Reson. Med. 51, 35–45 (2004).
5. Reeder, S. B. et al. Iterative decomposition of water and fat with echo asymmetry and least-squares estimation (IDEAL): Application with fast spin-echo imaging. Magn. Reson. Med. 54, 636–644 (2005).
6. Yu, H. et al. Multiecho reconstruction for simultaneous water-fat decomposition and T2* estimation. J. Magn. Reson. Imaging 26, 1153–1161 (2007).
7. Daniels, C. J. et al. A comparison of quantitative methods for clinical imaging with hyperpolarized 13C-pyruvate. NMR Biomed. 29, 387–399 (2016).
Table 2: Sequence parameters for 13C NMR at 1 T and for 13C MRI at 3 T.