4299

A Quadrature Detection Technique to Resolve J-Modulated [2-13C]Lactate in Hyperpolarized [2-13C] Pyruvate Imaging
Keshav Datta1,2 and Daniel Mark Spielman1,2

1Radiology, Stanford University, Stanford, CA, United States, 2Electrical Engineering, Stanford University, Stanford, CA, United States

Synopsis

For the imaging of hyperpolarized [2-13C]pyruvate, we propose a hybrid approach in which the closely spaced metabolic peaks of pyruvate, glutamate, citrate, acetoacetate, and acetyl-carnitine are measured using spectroscopic imaging, while the highly shifted [2-13C]Lactate peak is measured using selective excitation with a conventional fast imaging readout. We further demonstrate that a quadrature image reconstruction algorithm combining data from two fast imaging acquisitions, with the readout shifted by 1/2J, eliminates severe artifacts that normally arise when imaging a J-coupled system for which 1/J is comparable to the readout window.

Introduction

Hyperpolarized [1-13C]Pyruvate (Pyr) has shown considerable promise for measuring in-vivo glucose metabolism with the resulting [1-13C]Lactate (Lac) signal reflecting glycolysis and that from 13C-bicarbonate (Bic) reflecting Pyr conversion to Acetyl-coA. In contrast, hyperpolarized [2-13C]Pyr offers the potential to measure both lactate and TCA cycle intermediates. While the two most popular dynamic imaging methods for hyperpolarized [1-13C]Pyr are fast spectroscopic imaging1 and interleaved metabolite-selective imaging2, unfortunately, neither of these approaches is directly applicable for [2-13C]Pyr studies. The large chemical shift of [2-13C]Lac (almost 150 ppm relative to [2-13C]Pyr) results in a large chemical shift spatial misregistration error in slice-selective spectroscopic imaging acquisitions3 (Fig. 1). While this misregistration error can be eliminated by sequential selective imaging of [2-13C]Lac, the density of spectral peaks from pyruvate, glutamate, citrate, acetoacetate, and acetyl-carnitine are too close together to allow efficient selective excitation. Here we propose a hybrid approach for which crowded spectral peaks around 170-210 ppm are measured using spectroscopic imaging, while the highly shifted [2-13C]Lac peak is measured using selective excitation with a conventional fast imaging readout.

Methods

Interleaving a spectroscopic readout for the peaks around 170- 210 ppm with an imaging readout for selectively excited [2-13C]Lac is relatively straightforward. The major complication is that the [2-13C]Lac peak is a doublet with a 140 Hz J-coupling constant, which results in an unacceptable k-space weighting during readout. For example, as shown in Fig. 2 and Fig.3, there will be a cosine weighting primarily in the phase encode direction(ky) due to the C2-H J-coupling during a bidirectional echo planar imaging(EPI) acquisition. In general, this k-space weighting can be modeled as: I1(kx,ky) = I(kx,ky) cos(πJt) . Simlarly, a sin(πJt) weighted image can be generated using a readout shifted by (1/2J). The complex combination of these two “Quadrature” components results in: I’(kx,ky) = I(kx,ky) cos(πJt) + i.I(kx,ky) cos(πJ(t+1/2J)) = I(kx,ky) { cos(πJt) + i.sin(πJt) } = I(kx,ky) ei(πJt). Hence, combining two acquisitions with 1/2J shifted readouts (1/2J = 3.6 ms for [2-13C]Lac), an artifact-free, but spatially shifted, image can be reconstructed, with the spatial shift then easily removed during final image reconstruction. Simulations were performed in MATLAB (MathWorks Inc. MA, USA) by tracing an EPI k-space trajectory of a [2-13C]Lac Shepp-Logan phantom with parameters suitable for hyperpolarized pyruvate imaging in human brain (FOV=32cm, 32x32 imaging matrix, 10mmx10mmx20mm voxel size,) (Fig. 4). A Lac image, via quadrature combination of acquisitions from two echo times, TE1=10.6ms and TE2=14.2ms, was then computed to recover a shifted version of the original image. A correction for this spatial shift, computed as (Treadout*JCH/2), was then applied as the final step in image reconstruction.

Results

The cosine modulation of k-space due to J-coupling creates severe ghosting and blurring artifacts for metabolite-selective [2-13C]Lac imaging. These artifacts were effectively removed by quadrature combination of TE1=10.6 ms and TE2=14.2 ms acquisitions. The reconstruction from simulated single-shot EPI data sets collected at TE1=10.6ms and TE2=14.2ms of a [2-13C]Lactate Shepp-Logan phantom are presented in Fig. 4a & 4b respectively. Comparison of the original phantom (Fig. 4c) and the result of quadrature reconstruction (Fig. 4d) highlights the effectiveness of this new technique in recovering the original lactate image in the presence of J-modulation, even in the case of long readout times (Fig. 5).

Conclusion

In this work we show that by acquiring EPI images at two different echo times separated by 1/2J, J-coupling artifacts can be eliminated using a quadrature detection reconstruction technique. With this method, a hybrid imaging sequence that incorporates spectroscopic imaging readout for pyruvate, glutamate, citrate, acetoacetate, and acetyl-carnitine and spectrally selective fast EPI (or other fast imaging sequence such as spirals) for [2-13C]Lac can be developed for robust dynamic imaging of hyperpolarized [2-13C]Pyr and downstream products.

Acknowledgements

The Lucas Foundation, GE Healthcare, and NIH grants R01 CA17683603, R01 EB01901802, R01 MH110683, and P41 EB01589121.

References

1. Larson PE, Hu S, Lustig M, Kerr AB, Nelson SJ, Kurhanewicz J, Pauly JM, Vigneron DB. Fast dynamic 3D MR spectroscopic imaging with compressed sensing and multiband excitation pulses for hyperpolarized 13C studies. Magn Reson Med. 2011 Mar;65(3):610-9. doi: 10.1002/mrm.22650. Epub 2010 Oct 11. PubMed PMID: 20939089; PubMed Central PMCID: PMC3021589

2. Lau, A. Z., Chen, A. P., Ghugre, N. R., Ramanan, V. , Lam, W. W., Connelly, K. A., Wright, G. A. and Cunningham, C. H. (2010), Rapid multislice imaging of hyperpolarized 13C pyruvate and bicarbonate in the heart. Magn. Reson. Med., 64: 1323-1331. doi:10.1002/mrm.22525

3. Park JM, Josan S, Jang T, et al. Volumetric spiral chemical shift imaging of hyperpolarized [2-(13) c]pyruvate in a rat c6 glioma model. Magn Reson Med. 2015;75(3):973-84.

Figures

Figure 1: Chemical shift displacement artifacts. A) An averaged spectrum measured from a rat brain in-vivo using a 2D slice-selective FIDCSI after a bolus injection of 125-mM hyperpolarized [2-13C]pyruvate(post-DCA) at 3T. B) Zoomed region containing most of mitochondrialmetabolite peaks. C) Shifted slice excitation profiles, shown on sagittal 1H image of rat brain, for [2-13C]lactate, [2-13C]pyruvatehydrate, [5-13C]glutamate, and [2-13C]pyruvate due to the wide spectral distribution of the metabolites when an 8-mm slice(shaded-region) is prescribed at the resonance frequency of [5-13C]glutamate.Reconstructed metabolite maps of [2-13C]lactate (D), [2-13C]pyruvatehydrate (E), [5-13C]glutamate (F), [2-13C]pyruvate (G), overlaid on 1H MR images at the corresponding slice locations.

Figure 2: Simulations of the 800us Sinc radio frequency excitation pulse (blue) and the EPI gradient waveform (orange, blips for phase encode not shown) for typical hyperpolarized pyruvate imaging parameters- FOV=32cm, 32x32 imaging matrix, total readout time=10ms. The cosine temporal weighting function due to C2-H J-coupling (JCH=140Hz) in [2-13C]Lac is shown in red for two different acquisition times (a) TE1=10.6ms and (b) TE2=14.2ms

Figure 3: The 2D weighting functions in k-space (kx=readout & ky=phase encode) arising from the effect of J-coupling on the EPI readout for two acquisition times separated by 1/2J (a) TE1=10.6ms and (b) TE2=14.2ms

Figure 4: Simulation of reconstructions from the quadrature components (a) sine (im_sin) and (b) cosine (im_cos) for Treadout = 10ms; The result of complex combination (im_cos + i.im_sin) of the quadrature components is shown in (d) in comparison to the original phantom (c).

Figure 5: The effectiveness of the quadrature detection method is highlighted in the case of longer readout time (Treadout=80ms). The severe ghosting and blurring artifacts in the modulated images (a) TE1=50ms (im_sin) and (b) TE2=53.6ms (im_cos) are eliminated in the combined image (d), as can be seen by comparing it to the original (c). The resulting shift is easily corrected in the final reconstruction step.

Proc. Intl. Soc. Mag. Reson. Med. 27 (2019)
4299