Dynamic Imaging of Glucose and Lactate Metabolism by 13C-MRS without Hyperpolarization

Metabolic reprogramming is one of the defining features of cancer and abnormal metabolism is associated with many other pathologies. Molecular imaging techniques capable of detecting such changes have become essential for cancer diagnosis, treatment planning, and surveillance. In particular, 18F-FDG (fluorodeoxyglucose) PET has emerged as an essential imaging modality for cancer because of its unique ability to detect a disturbed molecular pathway through measurements of glucose uptake. However, FDG-PET has limitations that restrict its usefulness in certain situations and the information gained is limited to glucose uptake only.13C magnetic resonance spectroscopy theoretically has certain advantages over FDG-PET, but its inherent low sensitivity has restricted its use mostly to single voxel measurements unless dissolution dynamic nuclear polarization (dDNP) is used to increase the signal, which brings additional complications for clinical use. We show here a new method of imaging glucose metabolism in vivo by MRI chemical shift imaging (CSI) experiments that relies on a simple, but robust and efficient, post-processing procedure by the higher dimensional analog of singular value decomposition, tensor decomposition. Using this procedure, we achieve an order of magnitude increase in signal to noise in both dDNP and non-hyperpolarized non-localized experiments without sacrificing accuracy. In CSI experiments an approximately 30-fold increase was observed, enough that the glucose to lactate conversion indicative of the Warburg effect can be imaged without hyper-polarization with a time resolution of 12s and an overall spatial resolution that compares favorably to 18F-FDG PET.

[1]  M. Rosay,et al.  Brute-Force Hyperpolarization for NMR and MRI. , 2015, Journal of the American Chemical Society.

[2]  Rasmus Bro,et al.  The N-way Toolbox for MATLAB , 2000 .

[3]  John Kurhanewicz,et al.  Analysis of cancer metabolism by imaging hyperpolarized nuclei: prospects for translation to clinical research. , 2011, Neoplasia.

[4]  Manuel Martín-Pastor,et al.  A new general-purpose fully automatic baseline-correction procedure for 1D and 2D NMR data. , 2006, Journal of magnetic resonance.

[5]  Qianjin Feng,et al.  Denoise diffusion-weighted images using higher-order singular value decomposition , 2017, NeuroImage.

[6]  Minh N. Do,et al.  Spatiotemporal denoising of MR spectroscopic imaging data by low-rank approximations , 2011, 2011 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[7]  R L Wahl,et al.  Pitfalls in oncologic diagnosis with FDG PET imaging: physiologic and benign variants. , 1999, Radiographics : a review publication of the Radiological Society of North America, Inc.

[8]  T. Jang,et al.  Hyperpolarized 13C‐lactate to 13C‐bicarbonate ratio as a biomarker for monitoring the acute response of anti‐vascular endothelial growth factor (anti‐VEGF) treatment , 2016, NMR in biomedicine.

[9]  Zhi-Pei Liang,et al.  Performance analysis of denoising with low-rank and sparsity constraints , 2013, 2013 IEEE 10th International Symposium on Biomedical Imaging.

[10]  Mark F. Lythgoe,et al.  In vivo imaging of glucose uptake and metabolism in tumors , 2012, Nature Medicine.

[11]  Robert D. Nowak,et al.  Wavelet-based Rician noise removal for magnetic resonance imaging , 1999, IEEE Trans. Image Process..

[12]  R. Haiges,et al.  Increasing sensitivity in determining chemical shifts in one dimensional Lorentzian NMR spectra. , 2013, The journal of physical chemistry. A.

[13]  H. Gudbjartsson,et al.  The rician distribution of noisy mri data , 1995, Magnetic resonance in medicine.

[14]  P. Eilers A perfect smoother. , 2003, Analytical chemistry.

[15]  Niamh M. Long,et al.  Causes and imaging features of false positives and false negatives on 18F-PET/CT in oncologic imaging , 2011, Insights into imaging.

[16]  M. Pomper,et al.  Molecular Imaging of Metastatic Potential , 2008, Journal of Nuclear Medicine.

[17]  Zhi-Pei Liang,et al.  Improved image reconstruction for subspace-based spectroscopic imaging using non-quadratic regularization , 2014, 2014 36th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[18]  Jacco D. van Beek,et al.  matNMR: A flexible toolbox for processing, analyzing and visualizing magnetic resonance data in Matlab® , 2007 .

[19]  Peter C M van Zijl,et al.  MRI detection of glycogen in vivo by using chemical exchange saturation transfer imaging (glycoCEST) , 2007, Proceedings of the National Academy of Sciences.

[20]  M. O'Doherty,et al.  Radiation dose rates from patients undergoing positron emission tomography: implications for technologists and waiting areas , 2000, European Journal of Nuclear Medicine.

[21]  Prateek Jain,et al.  Tensor vs. Matrix Methods: Robust Tensor Decomposition under Block Sparse Perturbations , 2015, AISTATS.

[22]  D. van Ormondt,et al.  SVD-based quantification of magnetic resonance signals , 1992 .

[23]  Joos Vandewalle,et al.  On the Best Rank-1 and Rank-(R1 , R2, ... , RN) Approximation of Higher-Order Tensors , 2000, SIAM J. Matrix Anal. Appl..

[24]  A. Bax,et al.  Improved linear prediction of damped NMR signals using modified forward-backward linear prediction , 1992 .

[25]  C. Eckart,et al.  The approximation of one matrix by another of lower rank , 1936 .

[26]  C Chauvin,et al.  Wavelets and related time‐frequency techniques in magnetic resonance spectroscopy , 2001, NMR in biomedicine.

[27]  Pernille R. Jensen,et al.  Magnetic resonance imaging of pH in vivo using hyperpolarized 13C-labelled bicarbonate , 2008, Nature.

[28]  Zhiqiang Weng,et al.  Communication An efficient algorithm for automatic phase correction of NMR spectra based on entropy minimization , 2002 .

[29]  James A. Cadzow,et al.  Signal enhancement-a composite property mapping algorithm , 1988, IEEE Trans. Acoust. Speech Signal Process..

[30]  A S Stern,et al.  Does the maximum entropy method improve sensitivity? , 1990, Proceedings of the National Academy of Sciences of the United States of America.

[31]  J. Bomanji,et al.  18F-FDG PET/CT Imaging In Oncology , 2011, Annals of Saudi medicine.

[32]  J. V. van Beek,et al.  matNMR: a flexible toolbox for processing, analyzing and visualizing magnetic resonance data in Matlab. , 2007, Journal of magnetic resonance.

[33]  Y. Yen,et al.  Dynamic metabolic imaging of hyperpolarized [2‐13C]pyruvate using spiral chemical shift imaging with alternating spectral band excitation , 2014, Magnetic resonance in medicine.

[34]  M. Mattingly,et al.  Transport and imaging of brute-force (13)C hyperpolarization. , 2015, Journal of magnetic resonance.

[35]  N. Hay,et al.  Reprogramming glucose metabolism in cancer: can it be exploited for cancer therapy? , 2016, Nature Reviews Cancer.

[36]  L M Hamberg,et al.  The dose uptake ratio as an index of glucose metabolism: useful parameter or oversimplification? , 1994, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[37]  I. Daubechies Orthonormal bases of compactly supported wavelets , 1988 .

[38]  Bruce D Cheson,et al.  Progress and Promise of FDG-PET Imaging for Cancer Patient Management and Oncologic Drug Development , 2005, Clinical Cancer Research.

[39]  H. Chen,et al.  A Signal-Enhancement Algorithm for the Quantification of NMR Data in the Time Domain , 1994 .

[40]  Tinsu Pan,et al.  Frequent Diagnostic Errors in Cardiac PET/CT Due to Misregistration of CT Attenuation and Emission PET Images: A Definitive Analysis of Causes, Consequences, and Corrections , 2007, Journal of Nuclear Medicine.

[41]  Yan Liu,et al.  Improved Low-Rank Filtering of Magnetic Resonance Spectroscopic Imaging Data Corrupted by Noise and $B_0$ Field Inhomogeneity , 2016, IEEE Transactions on Biomedical Engineering.

[42]  VandewalleJoos,et al.  On the Best Rank-1 and Rank-(R1,R2,. . .,RN) Approximation of Higher-Order Tensors , 2000 .

[43]  J. Ardenkjær-Larsen,et al.  Increase in signal-to-noise ratio of > 10,000 times in liquid-state NMR , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[44]  Ray Freeman,et al.  Supercycles for broadband heteronuclear decoupling , 1982 .

[45]  Sabine Van Huffel,et al.  Parameter Estimation with Prior Knowledge of Known Signal Poles for the Quantification of NMR Spectroscopy Data in the Time Domain , 1996 .

[46]  K. Brindle,et al.  17: Magnetic resonance imaging of tumour metabolism , 2014 .

[47]  Tamara G. Kolda,et al.  Tensor Decompositions and Applications , 2009, SIAM Rev..

[48]  R Gruetter,et al.  Automatic, localized in Vivo adjustment of all first‐and second‐order shim coils , 1993, Magnetic resonance in medicine.

[49]  Bryan A. Clifford,et al.  High‐resolution dynamic 31P‐MRSI using a low‐rank tensor model , 2017, Magnetic resonance in medicine.

[50]  C. Kuntner Kinetic modeling in pre-clinical positron emission tomography. , 2014, Zeitschrift fur medizinische Physik.

[51]  Shingo Matsumoto,et al.  Pyruvate sensitizes pancreatic tumors to hypoxia-activated prodrug TH-302 , 2015, Cancer & metabolism.

[52]  S. Confort-Gouny,et al.  Optimization of residual water signal removal by HLSVD on simulated short echo time proton MR spectra of the human brain. , 2001, Journal of magnetic resonance.

[53]  Zhipeng Yang,et al.  Global denoising for 3D MRI , 2016, Biomedical engineering online.

[54]  Gil Navon,et al.  Molecular imaging of tumors and metastases using chemical exchange saturation transfer (CEST) MRI , 2013, Scientific Reports.

[55]  Angel Moreno,et al.  Alternative 1‐13C glucose infusion protocols for clinical 13C MRS examinations of the brain , 2001, Magnetic resonance in medicine.

[56]  D. Mankoff,et al.  Hypoxia imaging-directed radiation treatment planning , 2006, European Journal of Nuclear Medicine and Molecular Imaging.

[57]  H. Degani,et al.  Glycolysis as a metabolic marker in orthotopic breast cancer, monitored by in vivo (13)C MRS. , 2002, American journal of physiology. Endocrinology and metabolism.

[58]  Minh N. Do,et al.  Denoising MR Spectroscopic Imaging Data With Low-Rank Approximations , 2013, IEEE Transactions on Biomedical Engineering.

[59]  Kieran Clarke,et al.  Hyperpolarized Magnetic Resonance: A Novel Technique for the In Vivo Assessment of Cardiovascular Disease , 2011, Circulation.

[60]  I. Johnstone,et al.  Ideal spatial adaptation by wavelet shrinkage , 1994 .

[61]  Jie Zheng,et al.  Energy metabolism of cancer: Glycolysis versus oxidative phosphorylation (Review). , 2012, Oncology letters.

[62]  Zheng-Hua Wu,et al.  Improved image reconstruction for subspace-based spectroscopic imaging using non-quadratic regularization. , 2014, Conference proceedings : ... Annual International Conference of the IEEE Engineering in Medicine and Biology Society. IEEE Engineering in Medicine and Biology Society. Annual Conference.

[63]  E. Henry,et al.  [8] Singular value decomposition: Application to analysis of experimental data , 1992 .

[64]  E. Küstermann,et al.  Simultaneous in vivo monitoring of hepatic glucose and glucose‐6‐phosphate by 13C‐NMR spectroscopy , 2000, Magnetic resonance in medicine.

[65]  Shantanu H. Joshi,et al.  Denoising of MR spectroscopy signals using total variation and iterative Gauss-Seidel gradient updates , 2015, 2015 IEEE 12th International Symposium on Biomedical Imaging (ISBI).

[66]  Jean-Yves Tourneret,et al.  Hybrid sparse regularization for magnetic resonance spectroscopy , 2013, 2013 35th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC).

[67]  Markus Schwaiger,et al.  In Vivo Imaging of Tumor Metabolism and Acidosis by Combining PET and MRI-CEST pH Imaging. , 2016, Cancer research.

[68]  Justin P. Haldar,et al.  Denoising of MR spectroscopic imaging data with spatial-spectral regularization , 2010, 2010 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[69]  Xinyuan Zhang,et al.  Denoising of 3D magnetic resonance images by using higher-order singular value decomposition , 2015, Medical Image Anal..

[70]  S. Huffel,et al.  MR spectroscopy quantitation: a review of time‐domain methods , 2001, NMR in biomedicine.

[71]  Haim Azhari,et al.  Super-resolution in PET imaging , 2006, IEEE Transactions on Medical Imaging.

[72]  J Hennig,et al.  RARE imaging: A fast imaging method for clinical MR , 1986, Magnetic resonance in medicine.

[73]  W. Dietrich,et al.  Fast and precise automatic baseline correction of one- and two-dimensional nmr spectra , 1991 .

[74]  Fan Lam,et al.  A subspace approach to high‐resolution spectroscopic imaging , 2014, Magnetic resonance in medicine.

[75]  R. Hamanaka,et al.  Targeting glucose metabolism for cancer therapy , 2012, The Journal of experimental medicine.

[76]  Sabine Van Huffel,et al.  Subspace-based parameter estimation of exponentially damped sinusoids using prior knowledge of frequency and phase , 1997, Signal Process..

[77]  Ray Freeman,et al.  Broadband Decoupling in High-Resolution Nuclear Magnetic Resonance Spectroscopy , 1984 .

[78]  N. Freedman,et al.  Physiologic and Benign Variants , 2014 .