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.