Data Format Standardization and DICOM Integration for Hyperpolarized 13C MRI

Hyperpolarized (HP) 13C MRI has shown promise as a valuable modality for in vivo measurements of metabolism and is currently in human trials at 15 research sites worldwide. With this growth, it is important to adopt standardized data storage practices as it will allow sites to meaningfully compare data. In this paper, we (1) describe data that we believe should be stored and (2) demonstrate pipelines and methods that utilize the Digital Imaging and Communications in Medicine (DICOM) standard. This includes proposing a set of minimum set of information that is specific to HP 13C MRI studies. We then show where the majority of these can be fit into existing DICOM attributes, primarily via the "Contrast/Bolus" module. We also demonstrate pipelines for utilizing DICOM for HP 13C MRI. DICOM is the most common standard for clinical medical image storage and provides the flexibility to accommodate the unique aspects of HP 13C MRI, including the HP agent information but also spectroscopic and metabolite dimensions. The pipelines shown include creating DICOM objects for studies on human and animal imaging systems with various pulse sequences. We also show a python-based method to efficiently modify DICOM objects to incorporate the unique HP 13C MRI information that is not captured by existing pipelines. Moreover, we propose best practices for HP 13C MRI data storage that will support future multi-site trials, research studies, and technical developments of this imaging technique.

Protocol for producing hyperpolarized 13C-bicarbonate for clinical MRI of extracellular pH in aggressive tumors

Tumor acidosis is one of the hallmarks indicating the initiation and progression of various cancers. Here, we present a protocol for preparing a hyperpolarized (HP) 13C-bicarbonate tissue pH MRI imaging contrast agent to detect aggressive tumors. We describe the steps for the formulation and polarization of a precursor molecule 13C-glycerol carbonate (13C-GLC), the post-dissolution reaction, and converting HP 13C-GLC to an injectable HP 13C-bicarbonate solution. We then detail procedures for MRI data acquisition to generate tumor pH maps for assessing tumor aggressiveness. For complete details on the use and execution of this protocol, please refer to Mu et al.1.

Current methods for hyperpolarized [1-13C]pyruvate MRI human studies

MRI with hyperpolarized (HP) 13C agents, also known as HP 13C MRI, can measure processes such as localized metabolism that is altered in numerous cancers, liver, heart, kidney diseases, and more. It has been translated into human studies during the past 10 years, with recent rapid growth in studies largely based on increasing availability of HP agent preparation methods suitable for use in humans. This paper aims to capture the current successful practices for HP MRI human studies with [1-13C]pyruvate-by far the most commonly used agent, which sits at a key metabolic junction in glycolysis. The paper is divided into four major topic areas: (1) HP 13C-pyruvate preparation; (2) MRI system setup and calibrations; (3) data acquisition and image reconstruction; and (4) data analysis and quantification. In each area, we identified the key components for a successful study, summarized both published studies and current practices, and discuss evidence gaps, strengths, and limitations. This paper is the output of the "HP 13C MRI Consensus Group" as well as the ISMRM Hyperpolarized Media MR and Hyperpolarized Methods and Equipment study groups. It further aims to provide a comprehensive reference for future consensus, building as the field continues to advance human studies with this metabolic imaging modality.

Data Format Standardization and DICOM Integration for Hyperpolarized 13C MRI

Hyperpolarized (HP) 13C MRI has shown promise as a valuable modality for in vivo measurements of metabolism and is currently in human trials at 15 research sites worldwide. With this growth it is important to adopt standardized data storage practices as it will allow sites to meaningfully compare data. In this paper we (1) describe data that we believe should be stored and (2) demonstrate pipelines and methods that utilize the Digital Imaging and Communications in Medicine (DICOM) standard. This includes proposing a set of minimum set of information that is specific to HP 13C MRI studies. We then show where the majority of these can be fit into existing DICOM Attributes, primarily via the "Contrast/Bolus" module. We also demonstrate pipelines for utilizing DICOM for HP 13C MRI. DICOM is the most common standard for clinical medical image storage and provides the flexibility to accommodate the unique aspects of HP 13C MRI, including the HP agent information but also spectroscopic and metabolite dimensions. The pipelines shown include creating DICOM objects for studies on human and animal imaging systems with various pulse sequences. We also show a python-based method to efficiently modify DICOM objects to incorporate the unique HP 13C MRI information that is not captured by existing pipelines. Moreover, we propose best practices for HP 13C MRI data storage that will support future multi-site trials, research studies and technical developments of this imaging technique.

Distinguishing metabolic signals of liver tumors from surrounding liver cells using hyperpolarized 13 C MRI and gadoxetate

PURPOSE

To use the hepatocyte-specific gadolinium-based contrast agent gadoxetate combined with hyperpolarized (HP) [1-13 C]pyruvate MRI to selectively suppress metabolic signals from normal hepatocytes while preserving the signals arising from tumors.

METHODS

Simulations were performed to determine the expected changes in HP 13 C MR signal in liver and tumor under the influence of gadoxetate. CC531 colon cancer cells were implanted into the livers of five Wag/Rij rats. Liver and tumor metabolism were imaged at 3 T using HP [1-13 C] pyruvate chemical shift imaging before and 15 min after injection of gadoxetate. Area under the curve for pyruvate and lactate were measured from voxels containing at least 75% of normal-appearing liver or tumor.

RESULTS

Numerical simulations predicted a 36% decrease in lactate-to-pyruvate (L/P) ratio in liver and 16% decrease in tumor. In vivo, baseline L/P ratio was 0.44 ± 0.25 in tumors versus 0.21 ± 0.08 in liver (p = 0.09). Following administration of gadoxetate, mean L/P ratio decreased by an average of 0.11 ± 0.06 (p < 0.01) in normal-appearing liver. In tumors, mean L/P ratio post-gadoxetate did not show a statistically significant change from baseline. Compared to baseline levels, the relative decrease in L/P ratio was significantly greater in liver than in tumors (-0.52 ± 0.16 vs. -0.19 ± 0.25, p < 0.05).

CONCLUSIONS

The intracellular hepatobiliary contrast agent showed a greater effect suppressing HP 13 C MRI metabolic signals (through T1 shortening) in normal-appearing liver when compared to tumors. The combined use of HP MRI with selective gadolinium contrast agents may allow more selective imaging in HP 13 C MRI.

Hyperpolarized 13 C metabolic imaging of the human abdomen with spatiotemporal denoising

PURPOSE

Improving the quality and maintaining the fidelity of large coverage abdominal hyperpolarized (HP) 13 C MRI studies with a patch based global-local higher-order singular value decomposition (GL-HOVSD) spatiotemporal denoising approach.

METHODS

Denoising performance was first evaluated using the simulated [1-13 C]pyruvate dynamics at different noise levels to determine optimal kglobal and klocal parameters. The GL-HOSVD spatiotemporal denoising method with the optimized parameters was then applied to two HP [1-13 C]pyruvate EPI abdominal human cohorts (n = 7 healthy volunteers and n = 8 pancreatic cancer patients).

RESULTS

The parameterization of kglobal  = 0.2 and klocal  = 0.9 denoises abdominal HP data while retaining image fidelity when evaluated by RMSE. The kPX (conversion rate of pyruvate-to-metabolite, X = lactate or alanine) difference was shown to be <20% with respect to ground-truth metabolic conversion rates when there is adequate SNR (SNRAUC  > 5) for downstream metabolites. In both human cohorts, there was a greater than nine-fold gain in peak [1-13 C]pyruvate, [1-13 C]lactate, and [1-13 C]alanine apparent SNRAUC . The improvement in metabolite SNR enabled a more robust quantification of kPL and kPA . After denoising, we observed a 2.1 ± 0.4 and 4.8 ± 2.5-fold increase in the number of voxels reliably fit across abdominal FOVs for kPL and kPA quantification maps.

CONCLUSION

Spatiotemporal denoising greatly improves visualization of low SNR metabolites particularly [1-13 C]alanine and quantification of [1-13 C]pyruvate metabolism in large FOV HP 13 C MRI studies of the human abdomen.

Probing human heart TCA cycle metabolism and response to glucose load using hyperpolarized [2-13 C]pyruvate MRS

INTRODUCTION

The healthy heart has remarkable metabolic flexibility that permits rapid switching between mitochondrial glucose oxidation and fatty acid oxidation to generate ATP. Loss of metabolic flexibility has been implicated in the genesis of contractile dysfunction seen in cardiomyopathy. Metabolic flexibility has been imaged in experimental models, using hyperpolarized (HP) [2-13 C]pyruvate MRI, which enables interrogation of metabolites that reflect tricarboxylic acid (TCA) cycle flux in cardiac myocytes. This study aimed to develop methods, demonstrate feasibility for [2-13 C]pyruvate MRI in the human heart for the first time, and assess cardiac metabolic flexibility.

METHODS

Good manufacturing practice [2-13 C]pyruvic acid was polarized in a 5 T polarizer for 2.5-3 h. Following dissolution, quality control parameters of HP pyruvate met all safety and sterility criteria for pharmacy release, prior to administration to study subjects. Three healthy subjects each received two HP injections and MR scans, first under fasting conditions, followed by oral glucose load. A 5 cm axial slab-selective spectroscopy approach was prescribed over the left ventricle and acquired at 3 s intervals on a 3 T clinical MRI scanner.

RESULTS

The study protocol, which included HP substrate injection, MR scanning, and oral glucose load, was performed safely without adverse events. Key downstream metabolites of [2-13 C]pyruvate metabolism in cardiac myocytes include the glycolytic derivative [2-13 C]lactate, TCA-associated metabolite [5-13 C]glutamate, and [1-13 C]acetylcarnitine, catalyzed by carnitine acetyltransferase (CAT). After glucose load, 13 C-labeling of lactate, glutamate, and acetylcarnitine from 13 C-pyruvate increased by an average of 39.3%, 29.5%, and 114% respectively in the three subjects, which could result from increases in lactate dehydrogenase, pyruvate dehydrogenase, and CAT enzyme activity as well as TCA cycle flux (glucose oxidation).

CONCLUSIONS

HP [2-13 C]pyruvate imaging is safe and permits noninvasive assessment of TCA cycle intermediates and the acetyl buffer, acetylcarnitine, which is not possible using HP [1-13 C]pyruvate. Cardiac metabolite measurement in the fasting/fed states provides information on cardiac metabolic flexibility and the acetylcarnitine pool.

Dynamic T 2 * relaxometry of hyperpolarized [1-13 C]pyruvate MRI in the human brain and kidneys

PURPOSE

This study aimed to quantifyT 2 * $$ {T}_2^{\ast } $$ for hyperpolarized [1-13 C]pyruvate and metabolites in the healthy human brain and renal cell carcinoma (RCC) patients at 3 T.

METHODS

DynamicT 2 * $$ {T}_2^{\ast } $$ values were measured with a metabolite-specific multi-echo spiral sequence. The dynamicT 2 * $$ {T}_2^{\ast } $$ of [1-13 C]pyruvate, [1-13 C]lactate, and 13 C-bicarbonate was estimated in regions of interest in the whole brain, sinus vein, gray matter, and white matter in healthy volunteers, as well as in kidney tumors and the contralateral healthy kidneys in a separate group of RCC patients.T 2 * $$ {T}_2^{\ast } $$ was fit using a mono-exponential function; and metabolism was quantified using pyruvate-to-lactate conversion rate maps and lactate-to-pyruvate ratio maps, which were compared with and without an estimatedT 2 * $$ {T}_2^{\ast } $$ correction.

RESULTS

TheT 2 * $$ {T}_2^{\ast } $$ of pyruvate was shown to vary during the acquisition, whereas theT 2 * $$ {T}_2^{\ast } $$ of lactate and bicarbonate were relatively constant through time and across the organs studied. TheT 2 * $$ {T}_2^{\ast } $$ of lactate was similar in gray matter (29.75 ± 1.04 ms), white matter (32.89 ± 0.9 ms), healthy kidney (34.61 ± 4.07 ms), and kidney tumor (33.01 ± 2.31 ms); and theT 2 * $$ {T}_2^{\ast } $$ of bicarbonate was different between whole-brain (108.17 ± 14.05 ms) and healthy kidney (58.45 ± 6.63 ms). TheT 2 * $$ {T}_2^{\ast } $$ of pyruvate had similar trends in both brain and RCC studies, reducing from 75.56 ± 2.23 ms to 22.24 ± 1.24 ms in the brain and reducing from 122.72 ± 9.86 ms to 57.38 ± 7.65 ms in the kidneys.

CONCLUSION

Multi-echo dynamic imaging can quantifyT 2 * $$ {T}_2^{\ast } $$ and metabolism in a single integrated acquisition. Clear differences were observed in theT 2 * $$ {T}_2^{\ast } $$ of metabolites and in their behavior throughout the timecourse.

Advanced Hyperpolarized 13C Metabolic Imaging Protocol for Patients with Gliomas: A Comprehensive Multimodal MRI Approach

This study aimed to implement a multimodal 1H/HP-13C imaging protocol to augment the serial monitoring of patients with glioma, while simultaneously pursuing methods for improving the robustness of HP-13C metabolic data. A total of 100 1H/HP [1-13C]-pyruvate MR examinations (104 HP-13C datasets) were acquired from 42 patients according to the comprehensive multimodal glioma imaging protocol. Serial data coverage, accuracy of frequency reference, and acquisition delay were evaluated using a mixed-effects model to account for multiple exams per patient. Serial atlas-based HP-13C MRI demonstrated consistency in volumetric coverage measured by inter-exam dice coefficients (0.977 ± 0.008, mean ± SD; four patients/11 exams). The atlas-derived prescription provided significantly improved data quality compared to manually prescribed acquisitions (n = 26/78; p = 0.04). The water-based method for referencing [1-13C]-pyruvate center frequency significantly reduced off-resonance excitation relative to the coil-embedded [13C]-urea phantom (4.1 ± 3.7 Hz vs. 9.9 ± 10.7 Hz; p = 0.0007). Significantly improved capture of tracer inflow was achieved with the 2-s versus 5-s HP-13C MRI acquisition delay (p = 0.007). This study demonstrated the implementation of a comprehensive multimodal 1H/HP-13C MR protocol emphasizing the monitoring of steady-state/dynamic metabolism in patients with glioma.

A user independent denoising method for x-nuclei MRI and MRS

PURPOSE

X-nuclei (also called non-proton MRI) MRI and spectroscopy are limited by the intrinsic low SNR as compared to conventional proton imaging. Clinical translation of x-nuclei examination warrants the need of a robust and versatile tool improving image quality for diagnostic use. In this work, we compare a novel denoising method with fewer inputs to the current state-of-the-art denoising method.

METHODS

Denoising approaches were compared on human acquisitions of sodium (23 Na) brain, deuterium (2 H) brain, carbon (13 C) heart and brain, and simulated dynamic hyperpolarized 13 C brain scans, with and without additional noise. The current state-of-the-art denoising method Global-local higher order singular value decomposition (GL-HOSVD) was compared to the few-input method tensor Marchenko-Pastur principal component analysis (tMPPCA). Noise-removal was quantified by residual distributions, and statistical analyses evaluated the differences in mean-square-error and Bland-Altman analysis to quantify agreement between original and denoised results of noise-added data.

RESULTS

GL-HOSVD and tMPPCA showed similar performance for the variety of x-nuclei data analyzed in this work, with tMPPCA removing ˜5% more noise on average over GL-HOSVD. The mean ratio between noise-added and denoising reproducibility coefficients of the Bland-Altman analysis when compared to the original are also similar for the two methods with 3.09 ± 1.03 and 2.83 ± 0.79 for GL-HOSVD and tMPPCA, respectively.

CONCLUSION

The strength of tMPPCA lies in the few-input approach, which generalizes well to different data sources. This makes the use of tMPPCA denoising a robust and versatile tool in x-nuclei imaging improvements and the preferred denoising method.

Investigating cerebral perfusion with high resolution hyperpolarized [1-13 C]pyruvate MRI

PURPOSE

To investigate high-resolution hyperpolarized (HP) 13 C pyruvate MRI for measuring cerebral perfusion in the human brain.

METHODS

HP [1-13 C]pyruvate MRI was acquired in five healthy volunteers with a multi-resolution EPI sequence with 7.5 × 7.5 mm2 resolution for pyruvate. Perfusion parameters were calculated from pyruvate MRI using block-circulant singular value decomposition and compared to relative cerebral blood flow calculated from arterial spin labeling (ASL). To examine regional perfusion patterns, correlations between pyruvate and ASL perfusion were performed for whole brain, gray matter, and white matter voxels.

RESULTS

High resolution 7.5 × 7.5 mm2 pyruvate images were used to obtain relative cerebral blood flow (rCBF) values that were significantly positively correlated with ASL rCBF values (r = 0.48, 0.20, 0.28 for whole brain, gray matter, and white matter voxels respectively). Whole brain voxels exhibited the highest correlation between pyruvate and ASL perfusion, and there were distinct regional patterns of relatively high ASL and low pyruvate normalized rCBF found across subjects.

CONCLUSION

Acquiring HP 13 C pyruvate metabolic images at higher resolution allows for finer spatial delineation of brain structures and can be used to obtain cerebral perfusion parameters. Pyruvate perfusion parameters were positively correlated to proton ASL perfusion values, indicating a relationship between the two perfusion measures. This HP 13 C study demonstrated that hyperpolarized pyruvate MRI can assess cerebral metabolism and perfusion within the same study.

Clinically Translatable Hyperpolarized 13C Bicarbonate pH Imaging Method for Use in Prostate Cancer

Solid tumors such as prostate cancer (PCa) commonly develop an acidic microenvironment with pH 6.5-7.2, owing to heterogeneous perfusion, high metabolic activity, and rapid cell proliferation. In preclinical prostate cancer models, disease progression is associated with a decrease in tumor extracellular pH, suggesting that pH imaging may reflect an imaging biomarker to detect aggressive and high-risk disease. Therefore, we developed a hyperpolarized carbon-13 MRI method to image the tumor extracellular pH (pHe) and prepared it for clinical translation for detection and risk stratification of PCa. This method relies on the rapid breakdown of hyperpolarized (HP) 1,2-glycerol carbonate (carbonyl-13C) via base-catalyzed hydrolysis to produce HP 13CO32-, which is neutralized and converted to HP H13CO3-. After injection, HP H13CO3- equilibrates with HP 13CO2 in vivo and enables the imaging of pHe. Using insights gleaned from mechanistic studies performed in the hyperpolarized state, we solved issues of polarization loss during preparation in a clinical polarizer system. We successfully customized a reaction apparatus suitable for clinical application, developed clinical standard operating procedures, and validated the radiofrequency pulse sequence and imaging data acquisition with a wide range of animal models. The results demonstrated that we can routinely produce a highly polarized and safe HP H13CO3- contrast agent suitable for human injection. Preclinical imaging studies validated the reliability and accuracy of measuring acidification in healthy kidney and prostate tumor tissue. These methods were used to support an Investigational New Drug application to the U.S. Food and Drug Administration. This methodology is now ready to be implemented in human trials, with the ultimate goal of improving the management of PCa.

Current Methods for Hyperpolarized [1-13C]pyruvate MRI Human Studies

MRI with hyperpolarized (HP) 13C agents, also known as HP 13C MRI, can measure processes such as localized metabolism that is altered in numerous cancers, liver, heart, kidney diseases, and more. It has been translated into human studies during the past 10 years, with recent rapid growth in studies largely based on increasing availability of hyperpolarized agent preparation methods suitable for use in humans. This paper aims to capture the current successful practices for HP MRI human studies with [1-13C]pyruvate - by far the most commonly used agent, which sits at a key metabolic junction in glycolysis. The paper is divided into four major topic areas: (1) HP 13C-pyruvate preparation, (2) MRI system setup and calibrations, (3) data acquisition and image reconstruction, and (4) data analysis and quantification. In each area, we identified the key components for a successful study, summarized both published studies and current practices, and discuss evidence gaps, strengths, and limitations. This paper is the output of the HP 13C MRI Consensus Group as well as the ISMRM Hyperpolarized Media MR and Hyperpolarized Methods & Equipment study groups. It further aims to provide a comprehensive reference for future consensus building as the field continues to advance human studies with this metabolic imaging modality.

Parallel detection of chemical reactions in a microfluidic platform using hyperpolarized nuclear magnetic resonance

The sensitivity of NMR may be enhanced by more than four orders of magnitude via dissolution dynamic nuclear polarization (dDNP), potentially allowing real-time, in situ analysis of chemical reactions. However, there has been no widespread use of the technique for this application and the major limitation has been the low experimental throughput caused by the time-consuming polarization build-up process at cryogenic temperatures and fast decay of the hyper-intense signal post dissolution. To overcome this limitation, we have developed a microfluidic device compatible with dDNP-MR spectroscopic imaging methods for detection of reactants and products in chemical reactions in which up to 8 reactions can be measured simultaneously using a single dDNP sample. Multiple MR spectroscopic data sets can be generated under the same exact conditions of hyperpolarized solute polarization, concentration, pH, and temperature. A proof-of-concept for the technology is demonstrated by identifying the reactants in the decarboxylation of pyruvate via hydrogen peroxide (e.g. 2-hydroperoxy-2-hydroxypropanoate, peroxymonocarbonate and CO2). dDNP-MR allows tracing of fast chemical reactions that would be barely detectable at thermal equilibrium by MR. We envisage that dDNP-MR spectroscopic imaging combined with microfluidics will provide a new high-throughput method for dDNP enhanced MR analysis of multiple components in chemical reactions and for non-destructive in situ metabolic analysis of hyperpolarized substrates in biological samples for laboratory and preclinical research.

Probing Human Heart TCA Cycle Metabolism and Response to Glucose Load using Hyperpolarized [2-13C]Pyruvate MR Spectroscopy

Introduction

The normal heart has remarkable metabolic flexibility that permits rapid switching between mitochondrial glucose oxidation and fatty acid (FA) oxidation to generate ATP. Loss of metabolic flexibility has been implicated in the genesis of contractile dysfunction seen in cardiomyopathy. Metabolic flexibility has been imaged in experimental models, using hyperpolarized (HP) [2-13C]pyruvate MRI, which enables interrogation of metabolites that reflect tricarboxylic acid (TCA) cycle flux in cardiac myocytes. This study aimed to develop methods, demonstrate feasibility for [2-13C]pyruvate MRI in the human heart for the first time, and assess cardiac metabolic flexibility.

Methods

Good Manufacturing Practice [2-13C]pyruvic acid was polarized in a 5T polarizer for 2.5-3 hours. Following dissolution, QC parameters of HP pyruvate met all safety and sterility criteria for pharmacy release, prior to administration to study subjects. Three healthy subjects each received two HP injections and MR scans, first under fasting conditions, followed by oral glucose load. A 5cm axial slab-selective spectroscopy approach was prescribed over the left ventricle and acquired at 3s intervals on a 3T clinical MRI scanner.

Results

The study protocol which included HP substrate injection, MR scanning and oral glucose load, was performed safely without adverse events. Key downstream metabolites of [2-13C]pyruvate metabolism in cardiac myocytes include the glycolytic derivative [2-13C]lactate, TCA-associated metabolite [5-13C]glutamate, and [1-13C]acetylcarnitine, catalyzed by carnitine acetyltransferase (CAT). After glucose load, 13C-labeling of lactate, glutamate, and acetylcarnitine from 13C-pyruvate increased by 39.3%, 29.5%, and 114%, respectively in the three subjects, that could result from increases in lactate dehydrogenase (LDH), pyruvate dehydrogenase (PDH), and CAT enzyme activity as well as TCA cycle flux (glucose oxidation).

Conclusions

HP [2-13C]pyruvate imaging is safe and permits non-invasive assessment of TCA cycle intermediates and the acetyl buffer, acetylcarnitine, which is not possible using HP [1-13C]pyruvate. Cardiac metabolite measurement in the fasting/fed states provides information on cardiac metabolic flexibility and the acetylcarnitine pool.

Hyperpolarized [2-13C]pyruvate MR molecular imaging with whole brain coverage

Hyperpolarized (HP) 13C Magnetic Resonance Imaging (MRI) was applied for the first time to image and quantify the uptake and metabolism of [2-13C]pyruvate in the human brain to provide new metabolic information on cerebral energy metabolism. HP [2-13C]pyruvate was injected intravenously and imaged in 5 healthy human volunteer exams with whole brain coverage in a 1-minute acquisition using a specialized spectral-spatial multi-slice echoplanar imaging (EPI) pulse sequence to acquire 13C-labeled volumetric and dynamic images of [2-13C]pyruvate and downstream metabolites [5-13C]glutamate and [2-13C]lactate. Metabolic ratios and apparent conversion rates of pyruvate-to-lactate (kPL) and pyruvate-to-glutamate (kPG) were quantified to investigate simultaneously glycolytic and oxidative metabolism in a single injection.

Multi-parametric hyperpolarized 13C/1H imaging reveals Warburg-related metabolic dysfunction and associated regional heterogeneity in high-grade human gliomas

BACKGROUND

Dynamic hyperpolarized (HP)-13C MRI has enabled real-time, non-invasive assessment of Warburg-related metabolic dysregulation in glioma using a [1-13C]pyruvate tracer that undergoes conversion to [1-13C]lactate and [13C]bicarbonate. Using a multi-parametric 1H/HP-13C imaging approach, we investigated dynamic and steady-state metabolism, together with physiological parameters, in high-grade gliomas to characterize active tumor.

METHODS

Multi-parametric 1H/HP-13C MRI data were acquired from fifteen patients with progressive/treatment-naïve glioblastoma [prog/TN GBM, IDH-wildtype (n = 11)], progressive astrocytoma, IDH-mutant, grade 4 (G4AIDH+, n = 2) and GBM manifesting treatment effects (n = 2). Voxel-wise regional analysis of the cohort with prog/TN GBM assessed imaging heterogeneity across contrast-enhancing/non-enhancing lesions (CEL/NEL) and normal-appearing white matter (NAWM) using a mixed effects model. To enable cross-nucleus parameter association, normalized perfusion, diffusion, and dynamic/steady-state (HP-13C/spectroscopic) metabolic data were collectively examined at the 13C resolution. Prog/TN GBM were similarly compared against progressive G4AIDH+ and treatment effects.

RESULTS

Regional analysis of Prog/TN GBM metabolism revealed statistically significant heterogeneity in 1H choline-to-N-acetylaspartate index (CNI)max, [1-13C]lactate, modified [1-13C]lactate-to-[1-13C]pyruvate ratio (CELval > NELval > NAWMval); [1-13C]lactate-to-[13C]bicarbonate ratio (CELval > NELval/NAWMval); and 1H-lactate (CELval/NELval > NAWMundetected). Significant associations were found between normalized perfusion (cerebral blood volume, nCBV; peak height, nPH) and levels of [1-13C]pyruvate and [1-13C]lactate, as well as between CNImax and levels of [1-13C]pyruvate, [1-13C]lactate and modified ratio. GBM, by comparison to G4AIDH+, displayed lower perfusion %-recovery and modeled rate constants for [1-13C]pyruvate-to-[1-13C]lactate conversion (kPL), and higher 1H-lactate and [1-13C]pyruvate levels, while having higher nCBV, %-recovery, kPL, [1-13C]pyruvate-to-[1-13C]lactate and modified ratios relative to treatment effects.

CONCLUSIONS

GBM consistently displayed aberrant, Warburg-related metabolism and regional heterogeneity detectable by novel HP-13C/1H imaging techniques.

Model-constrained reconstruction accelerated with Fourier-based undersampling for hyperpolarized [1-13 C] pyruvate imaging

PURPOSE

Model-constrained reconstruction with Fourier-based undersampling (MoReFUn) is introduced to accelerate the acquisition of dynamic MRI using hyperpolarized [1-13 C]-pyruvate.

METHODS

The MoReFUn method resolves spatial aliasing using constraints introduced by a pharmacokinetic model that describes the signal evolution of both pyruvate and lactate. Acceleration was evaluated on three single-channel data sets: a numerical digital phantom that is used to validate the accuracy of reconstruction and model parameter restoration under various SNR and undersampling ratios, prospectively and retrospectively sampled data of an in vitro dynamic multispectral phantom, and retrospectively undersampled imaging data from a prostate cancer patient to test the fidelity of reconstructed metabolite time series.

RESULTS

All three data sets showed successful reconstruction using MoReFUn. In simulation and retrospective phantom data, the restored time series of pyruvate and lactate maintained the image details, and the mean square residual error of the accelerated reconstruction increased only slightly (< 10%) at a reduction factor up to 8. In prostate data, the quantitative estimation of the conversion-rate constant of pyruvate to lactate was achieved with high accuracy of less than 10% error at a reduction factor of 2 compared with the conversion rate derived from unaccelerated data.

CONCLUSION

The MoReFUn technique can be used as an effective and reliable imaging acceleration method for metabolic imaging using hyperpolarized [1-13 C]-pyruvate.

Monitoring response to a clinically relevant IDH inhibitor in glioma-Hyperpolarized 13C magnetic resonance spectroscopy approaches

Background

Mutant isocitrate dehydrogenase (IDHmut) catalyzes 2-hydroxyglutarate (2HG) production and is considered a therapeutic target for IDHmut tumors. However, response is mostly associated with inhibition of tumor growth. Response assessment via anatomic imaging is therefore challenging. Our goal was to directly detect IDHmut inhibition using a new hyperpolarized (HP) 13C magnetic resonance spectroscopy-based approach to noninvasively assess α-ketoglutarate (αKG) metabolism to 2HG and glutamate.

Methods

We studied IDHmut-expressing normal human astrocyte (NHAIDH1mut) cells and rats with BT257 tumors, and assessed response to the IDHmut inhibitor BAY-1436032 (n ≥ 4). We developed a new 13C Echo Planar Spectroscopic Imaging sequence with an optimized RF pulse to monitor the fate of HP [1-13C]αKG and [5-12C,1-13C]αKG with a 2.5 × 2.5 × 8 mm3 spatial resolution.

Results

Cell studies confirmed that BAY-1436032-treatment leads to a drop in HP 2HG and an increase in HP glutamate detectable with both HP substrates. Data using HP [5-12C,1-13C]αKG also demonstrated that its conversion to 2HG is detectable without the proximal 1.1% natural abundance [5-13C]αKG signal. In vivo studies showed that glutamate is produced in normal brains but no 2HG is detectable. In tumor-bearing rats, we detected the production of both 2HG and glutamate, and BAY-1436032-treatment led to a drop in 2HG and an increase in glutamate. Using HP [5-12C,1-13C]αKG we detected metabolism with an signal-to-noise ratio of 23 for 2HG and 17 for glutamate.

Conclusions

Our findings point to the clinical potential of HP αKG, which recently received FDA investigational new drug approval for research, for noninvasive localized imaging of IDHmut status.

Whole-Abdomen Metabolic Imaging of Healthy Volunteers Using Hyperpolarized [1-13 C]pyruvate MRI

BACKGROUND

Hyperpolarized 13 C MRI quantitatively measures enzyme-catalyzed metabolism in cancer and metabolic diseases. Whole-abdomen imaging will permit dynamic metabolic imaging of several abdominal organs simultaneously in healthy and diseased subjects.

PURPOSE

Image hyperpolarized [1-13 C]pyruvate and products in the abdomens of healthy volunteers, overcoming challenges of motion, magnetic field variations, and spatial coverage. Compare hyperpolarized [1-13 C]pyruvate metabolism across abdominal organs of healthy volunteers.

STUDY TYPE

Prospective technical development.

SUBJECTS

A total of 13 healthy volunteers (8 male), 21-64 years (median 36).

FIELD STRENGTH/SEQUENCE

A 3 T. Proton: T1 -weighted spoiled gradient echo, T2 -weighted single-shot fast spin echo, multiecho fat/water imaging. Carbon-13: echo-planar spectroscopic imaging, metabolite-specific echo-planar imaging.

ASSESSMENT

Transmit magnetic field was measured. Variations in main magnetic field (ΔB0 ) determined using multiecho proton acquisitions were compared to carbon-13 acquisitions. Changes in ΔB0 were measured after localized shimming. Improvements in metabolite signal-to-noise ratio were calculated. Whole-organ regions of interests were drawn over the liver, spleen, pancreas, and kidneys by a single investigator. Metabolite signals, time-to-peak, decay times, and mean first-order rate constants for pyruvate-to-lactate (kPL ) and alanine (kPA ) conversion were measured in each organ.

STATISTICAL TESTS

Linear regression, one-sample Kolmogorov-Smirnov tests, paired t-tests, one-way ANOVA, Tukey's multiple comparisons tests. P ≤ 0.05 considered statistically significant.

RESULTS

Proton ΔB0 maps correlated with carbon-13 ΔB0 maps (slope = 0.93, y-intercept = -2.88, R2  = 0.73). Localized shimming resulted in mean frequency offset within ±25 Hz for all organs. Metabolite SNR significantly increased after denoising. Mean kPL and kPA were highest in liver, followed by pancreas, spleen, and kidneys (all comparisons with liver were significant).

DATA CONCLUSION

Whole-abdomen coverage with hyperpolarized carbon-13 MRI was feasible despite technical challenges. Multiecho gradient echo 1 H acquisitions accurately predicted chemical shifts observed using carbon-13 spectroscopy. Carbon-13 acquisitions benefited from local shimming. Metabolite energetics in the abdomen compiled for healthy volunteers can be used to design larger clinical trials in patients with metabolic diseases.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 1.

Improving multiparametric MR-transrectal ultrasound guided fusion prostate biopsies with hyperpolarized 13 C pyruvate metabolic imaging: A technical development study

PURPOSE

To develop techniques and establish a workflow using hyperpolarized carbon-13 (13 C) MRI and the pyruvate-to-lactate conversion rate (kPL ) biomarker to guide MR-transrectal ultrasound fusion prostate biopsies.

METHODS

The integrated multiparametric MRI (mpMRI) exam consisted of a 1-min hyperpolarized 13 C-pyruvate EPI acquisition added to a conventional prostate mpMRI exam. Maps of kPL values were calculated, uploaded to a picture archiving and communication system and targeting platform, and displayed as color overlays on T2 -weighted anatomic images. Abdominal radiologists identified 13 C research biopsy targets based on the general recommendation of focal lesions with kPL  >0.02(s-1 ), and created a targeting report for each study. Urologists conducted transrectal ultrasound-guided MR fusion biopsies, including the standard 1 H-mpMRI targets as well as 12-14 core systematic biopsies informed by the research 13 C-kPL targets. All biopsy results were included in the final pathology report and calculated toward clinical risk.

RESULTS

This study demonstrated the safety and technical feasibility of integrating hyperpolarized 13 C metabolic targeting into routine 1 H-mpMRI and transrectal ultrasound fusion biopsy workflows, evaluated via 5 men (median age 71 years, prostate-specific antigen 8.4 ng/mL, Cancer of the Prostate Risk Assessment score 2) on active surveillance undergoing integrated scan and subsequent biopsies. No adverse event was reported. Median turnaround time was less than 3 days from scan to 13 C-kPL targeting, and scan-to-biopsy time was 2 weeks. Median number of 13 C targets was 1 (range: 1-2) per patient, measuring 1.0 cm (range: 0.6-1.9) in diameter, with a median kPL of 0.0319 s-1 (range: 0.0198-0.0410).

CONCLUSIONS

This proof-of-concept work demonstrated the safety and feasibility of integrating hyperpolarized 13 C MR biomarkers to the standard mpMRI workflow to guide MR-transrectal ultrasound fusion biopsies.

Multiparametric Magnetic Resonance Imaging and Metabolic Characterization of Patient-Derived Xenograft Models of Clear Cell Renal Cell Carcinoma

Patient-derived xenografts (PDX) are high-fidelity cancer models typically credentialled by genomics, transcriptomics and proteomics. Characterization of metabolic reprogramming, a hallmark of cancer, is less frequent. Dysregulated metabolism is a key feature of clear cell renal cell carcinoma (ccRCC) and authentic preclinical models are needed to evaluate novel imaging and therapeutic approaches targeting metabolism. We characterized 5 PDX from high-grade or metastatic ccRCC by multiparametric magnetic resonance imaging (MRI) and steady state metabolic profiling and flux analysis. Similar to MRI of clinical ccRCC, T2-weighted images of orthotopic tumors of most PDX were homogeneous. The increased hyperintense (cystic) areas observed in one PDX mimicked the cystic phenotype typical of some RCC. The negligible hypointense (necrotic) areas of PDX grown under the highly vascularized renal capsule are beneficial for preclinical studies. Mean apparent diffusion coefficient (ADC) values were equivalent to those of ccRCC in human patients. Hyperpolarized (HP) [1-13C]pyruvate MRI of PDX showed high glycolytic activity typical of high-grade primary and metastatic ccRCC with considerable intra- and inter-tumoral variability, as has been observed in clinical HP MRI of ccRCC. Comparison of steady state metabolite concentrations and metabolic flux in [U-13C]glucose-labeled tumors highlighted the distinctive phenotypes of two PDX with elevated levels of numerous metabolites and increased fractional enrichment of lactate and/or glutamate, capturing the metabolic heterogeneity of glycolysis and the TCA cycle in clinical ccRCC. Culturing PDX cells and reimplanting to generate xenografts (XEN), or passaging PDX in vivo, altered some imaging and metabolic characteristics while transcription remained like that of the original PDX. These findings show that PDX are realistic models of ccRCC for imaging and metabolic studies but that the plasticity of metabolism must be considered when manipulating PDX for preclinical studies.

Kinetic analysis of multi-resolution hyperpolarized 13 C human brain MRI to study cerebral metabolism

PURPOSE

To investigate multi-resolution hyperpolarized (HP) 13 C pyruvate MRI for measuring kinetic conversion rates in the human brain.

METHODS

HP [1-13 C]pyruvate MRI was acquired in 6 subjects with a multi-resolution EPI sequence at 7.5 × 7.5 mm2 resolution for pyruvate and 15 × 15 mm2 resolution for lactate and bicarbonate. With the same lactate data, 2 quantitative maps of pyruvate-to-lactate conversion (kPL ) maps were generated: 1 using 7.5 × 7.5 mm2 resolution pyruvate data and the other using synthetic 15 × 15 mm2 resolution pyruvate data to simulate a standard constant resolution acquisition. To examine local kPL values, 4 voxels were manually selected in each study representing brain tissue near arteries, brain tissue near veins, white matter, and gray matter.

RESULTS

High resolution 7.5 × 7.5 mm2 pyruvate images increased the spatial delineation of brain structures and decreased partial volume effects compared to coarser resolution 15 × 15 mm2 pyruvate images. Voxels near arteries, veins and in white matter exhibited higher calculated kPL for multi-resolution images.

CONCLUSION

Acquiring HP 13 C pyruvate metabolic data with a multi-resolution approach minimized partial volume effects from vascular pyruvate signals while maintaining the SNR of downstream metabolites. Higher resolution pyruvate images for kinetic fitting resulted in increased kinetic rate values, particularly around the superior sagittal sinus and cerebral arteries, by reducing extracellular pyruvate signal contributions from adjacent blood vessels. This HP 13 C study showed that acquiring pyruvate with finer resolution improved the quantification of kinetic rates throughout the human brain.

Development of specialized magnetic resonance acquisition techniques for human hyperpolarized [13 C,15 N2 ]urea + [1-13 C]pyruvate simultaneous perfusion and metabolic imaging

PURPOSE

This study aimed to develop and demonstrate the in vivo feasibility of a 3D stack-of-spiral balanced steady-state free precession(3D-bSSFP) urea sequence, interleaved with a metabolite-specific gradient echo (GRE) sequence for pyruvate and metabolic products, for improving the SNR and spatial resolution of the first hyperpolarized 13 C-MRI human study with injection of co-hyperpolarized [1-13 C]pyruvate and [13 C,15 N2 ]urea.

METHODS

A metabolite-specific bSSFP urea imaging sequence was designed using a urea-specific excitation pulse, optimized TR, and 3D stack-of-spiral readouts. Simulations and phantom studies were performed to validate the spectral response of the sequence. The image quality of urea data acquired by the 3D-bSSFP sequence and the 2D-GRE sequence was evaluated with 2 identical injections of co-hyperpolarized [1-13 C]pyruvate and [13 C,15 N2 ]urea formula in a rat. Subsequently, the feasibility of the acquisition strategy was validated in a prostate cancer patient.

RESULTS

Simulations and phantom studies demonstrated that 3D-bSSFP sequence achieved urea-only excitation, while minimally perturbing other metabolites (<1%). An animal study demonstrated that compared to GRE, bSSFP sequence provided an ∼2.5-fold improvement in SNR without perturbing urea or pyruvate kinetics, and bSSFP approach with a shorter spiral readout reduced blurring artifacts caused by J-coupling of [13 C,15 N2 ]urea. The human study demonstrated the in vivo feasibility and data quality of the acquisition strategy.

CONCLUSION

The 3D-bSSFP urea sequence with a stack-of-spiral acquisition demonstrated significantly increased SNR and image quality for [13 C,15 N2 ]urea in co-hyperpolarized [1-13 C]pyruvate and [13 C,15 N2 ]urea imaging studies. This work lays the foundation for future human studies to achieve high-quality and high-SNR metabolism and perfusion images.

Assessment of higher-order singular value decomposition denoising methods on dynamic hyperpolarized [1-13C]pyruvate MRI data from patients with glioma

BACKGROUND

Real-time metabolic conversion of intravenously-injected hyperpolarized [1-13C]pyruvate to [1-13C]lactate and [13C]bicarbonate in the brain can be measured using dynamic hyperpolarized carbon-13 (HP-13C) MRI. However, voxel-wise evaluation of metabolism in patients with glioma is challenged by the limited signal-to-noise ratio (SNR) of downstream 13C metabolites, especially within lesions. The purpose of this study was to evaluate the ability of higher-order singular value decomposition (HOSVD) denoising methods to enhance dynamic HP [1-13C]pyruvate MRI data acquired from patients with glioma.

METHODS

Dynamic HP-13C MRI were acquired from 14 patients with glioma. The effects of two HOSVD denoising techniques, tensor rank truncation-image enhancement (TRI) and global-local HOSVD (GL-HOSVD), on the SNR and kinetic modeling were analyzed in [1-13C]lactate data with simulated noise that matched the levels of [13C]bicarbonate signals. Both methods were then evaluated in patient data based on their ability to improve [1-13C]pyruvate, [1-13C]lactate and [13C]bicarbonate SNR. The effects of denoising on voxel-wise kinetic modeling of kPL and kPB was also evaluated. The number of voxels with reliable kinetic modeling of pyruvate-to-lactate (kPL) and pyruvate-to-bicarbonate (kPB) conversion rates within regions of interest (ROIs) before and after denoising was then compared.

RESULTS

Both denoising methods improved metabolite SNR and regional signal coverage. In patient data, the average increase in peak dynamic metabolite SNR was 2-fold using TRI and 4-5 folds using GL-HOSVD denoising compared to acquired data. Denoising reduced kPL modeling errors from a native average of 23% to 16% (TRI) and 15% (GL-HOSVD); and kPB error from 42% to 34% (TRI) and 37% (GL-HOSVD) (values were averaged voxelwise over all datasets). In contrast-enhancing lesions, the average number of voxels demonstrating within-tolerance kPL modeling error relative to the total voxels increased from 48% in the original data to 84% (TRI) and 90% (GL-HOSVD), while the number of voxels showing within-tolerance kPB modeling error increased from 0% to 15% (TRI) and 8% (GL-HOSVD).

CONCLUSION

Post-processing denoising methods significantly improved the SNR of dynamic HP-13C imaging data, resulting in a greater number of voxels satisfying minimum SNR criteria and maximum kinetic modeling errors in tumor lesions. This enhancement can aid in the voxel-wise analysis of HP-13C data and thereby improve monitoring of metabolic changes in patients with glioma following treatment.