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.

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.

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.

Spectroscopic imaging of D-2-hydroxyglutarate and other metabolites in pre-surgical patients with IDH-mutant lower-grade gliomas.

PURPOSE: Prognostically favorable IDH-mutant gliomas are known to produce oncometabolite D-2-hydroxyglutarate (2HG). In this study, we investigated metabolite-based features of patients with grade 2 and 3 glioma using 2HG-specific in vivo MR spectroscopy, to determine their relationship with image-guided tissue pathology and predictive role in progression-free survival (PFS). METHODS: Forty-five patients received pre-operative MRIs that included 3-D spectroscopy optimized for 2HG detection. Spectral data were reconstructed and quantified to compare metabolite levels according to molecular pathology (IDH1R132H, 1p/19q, and p53); glioma grade; histological subtype; and T2 lesion versus normal-appearing white matter (NAWM) ROIs. Levels of 2HG were correlated with other metabolites and pathological parameters (cellularity, MIB-1) from image-guided tissue samples using Pearson's correlation test. Metabolites predictive of PFS were evaluated with Cox proportional hazards models. RESULTS: Quantifiable levels of 2HG in 39/42 (93%) IDH+ and 1/3 (33%) IDH- patients indicated a 91.1% apparent detection accuracy. Myo-inositol/total choline (tCho) showed reduced values in astrocytic (1p/19q-wildtype), p53-mutant, and grade 3 (vs. 2) IDH-mutant gliomas (p < 0.05), all of which exhibited higher proportions of astrocytomas. Compared to NAWM, T2 lesions displayed elevated 2HG+ γ-aminobutyric acid (GABA)/total creatine (tCr) (p < 0.001); reduced glutamate/tCr (p < 0.001); increased myo-inositol/tCr (p < 0.001); and higher tCho/tCr (p < 0.001). Levels of 2HG at sampled tissue locations were significantly associated with tCho (R = 0.62; p = 0.002), total NAA (R = - 0.61; p = 0.002) and cellularity (R = 0.37; p = 0.04) but not MIB-1. Increasing levels of 2HG/tCr (p = 0.0007, HR 5.594) and thresholding (≥ 0.905, median value; p = 0.02) predicted adverse PFS. CONCLUSION: In vivo 2HG detection can reasonably be achieved on clinical scanners and increased levels may signal adverse PFS.

Hyperpolarized 13 C MRI data acquisition and analysis in prostate and brain at University of California, San Francisco.

Based on the expanding set of applications for hyperpolarized carbon-13 (HP-13 C) MRI, this work aims to communicate standardized methodology implemented at the University of California, San Francisco, as a primer for conducting reproducible metabolic imaging studies of the prostate and brain. Current state-of-the-art HP-13 C acquisition, data processing/reconstruction and kinetic modeling approaches utilized in patient studies are presented together with the rationale underpinning their usage. Organized around spectroscopic and imaging-based methods, this guide provides an extensible framework for handling a variety of HP-13 C applications, which derives from two examples with dynamic acquisitions: 3D echo-planar spectroscopic imaging of the human prostate and frequency-specific 2D multislice echo-planar imaging of the human brain. Details of sequence-specific parameters and processing techniques contained in these examples should enable investigators to effectively tailor studies around individual-use cases. Given the importance of clinical integration in improving the utility of HP exams, practical aspects of standardizing data formats for reconstruction, analysis and visualization are also addressed alongside open-source software packages that enhance institutional interoperability and validation of methodology. To facilitate the adoption and further development of this methodology, example datasets and analysis pipelines have been made available in the supporting information.

Characterization of serial hyperpolarized 13C metabolic imaging in patients with glioma.

BACKGROUND: Hyperpolarized carbon-13 (HP-13C) MRI is a non-invasive imaging technique for probing brain metabolism, which may improve clinical cancer surveillance. This work aimed to characterize the consistency of serial HP-13C imaging in patients undergoing treatment for brain tumors and determine whether there is evidence of aberrant metabolism in the tumor lesion compared to normal-appearing tissue. METHODS: Serial dynamic HP [1-13C]pyruvate MRI was performed on 3 healthy volunteers (6 total examinations) and 5 patients (21 total examinations) with diffuse infiltrating glioma during their course of treatment, using a frequency-selective echo-planar imaging (EPI) sequence. HP-13C imaging at routine clinical timepoints overlapped treatment, including radiotherapy (RT), temozolomide (TMZ) chemotherapy, and anti-angiogenic/investigational agents. Apparent rate constants for [1-13C]pyruvate conversion to [1-13C]lactate (kPL) and [13C]bicarbonate (kPB) were simultaneously quantified based on an inputless kinetic model within normal-appearing white matter (NAWM) and anatomic lesions defined from 1H MRI. The inter/intra-subject consistency of kPL-NAWM and kPB-NAWM was measured in terms of the coefficient of variation (CV). RESULTS: When excluding scans following anti-angiogenic therapy, patient values of kPL-NAWM and kPB-NAWM were 0.020 s-1 ± 23.8% and 0.0058 s-1 ± 27.7% (mean ± CV) across 17 HP-13C MRIs, with intra-patient serial kPL-NAWM/kPB-NAWM CVs ranging 6.8-16.6%/10.6-40.7%. In 4/5 patients, these values (0.018 s-1 ± 13.4% and 0.0058 s-1 ± 24.4%; n = 13) were more similar to those from healthy volunteers (0.018 s-1 ± 5.0% and 0.0043 s-1 ± 12.6%; n = 6) (mean ± CV). The anti-angiogenic agent bevacizumab was associated with global elevations in apparent rate constants, with maximum kPL-NAWM in 2 patients reaching 0.047 ± 0.001 and 0.047 ± 0.003 s-1 (±model error). In 3 patients with progressive disease, anatomic lesions showed elevated kPL relative to kPL-NAWM of 0.024 ± 0.001 s-1 (±model error) in the absence of gadolinium enhancement, and 0.032 ± 0.008, 0.040 ± 0.003 and 0.041 ± 0.009 s-1 with gadolinium enhancement. The lesion kPB in patients was reduced to unquantifiable values compared to kPB-NAWM. CONCLUSION: Serial measures of HP [1-13C]pyruvate metabolism displayed consistency in the NAWM of healthy volunteers and patients. Both kPL and kPB were globally elevated following bevacizumab treatment, while progressive disease demonstrated elevated kPL in gadolinium-enhancing and non-enhancing lesions. Larger prospective studies with homogeneous patient populations are planned to evaluate metabolic changes following treatment.

Tensor image enhancement and optimal multichannel receiver combination analyses for human hyperpolarized 13 C MRSI.

PURPOSE: With the initiation of human hyperpolarized 13 C (HP-13 C) trials at multiple sites and the development of improved acquisition methods, there is an imminent need to maximally extract diagnostic information to facilitate clinical interpretation. This study aims to improve human HP-13 C MR spectroscopic imaging through means of Tensor Rank truncation-Image enhancement (TRI) and optimal receiver combination (ORC). METHODS: A data-driven processing framework for dynamic HP 13 C MR spectroscopic imaging (MRSI) was developed. Using patient data sets acquired with both multichannel arrays and single-element receivers from the brain, abdomen, and pelvis, we examined the theory and application of TRI, as well as 2 ORC techniques: whitened singular value decomposition (WSVD) and first-point phasing. Optimal conditions for TRI were derived based on bias-variance trade-off. RESULTS: TRI and ORC techniques together provided a 63-fold mean apparent signal-to-noise ratio (aSNR) gain for receiver arrays and a 31-fold gain for single-element configurations, which particularly improved quantification of the lower-SNR-[13 C]bicarbonate and [1-13 C]alanine signals that were otherwise not detectable in many cases. Substantial SNR enhancements were observed for data sets that were acquired even with suboptimal experimental conditions, including delayed (114 s) injection (8× aSNR gain solely by TRI), or from challenging anatomy or geometry, as in the case of a pediatric patient with brainstem tumor (597× using combined TRI and WSVD). Improved correlation between elevated pyruvate-to-lactate conversion, biopsy-confirmed cancer, and mp-MRI lesions demonstrated that TRI recovered quantitative diagnostic information. CONCLUSION: Overall, this combined approach was effective across imaging targets and receiver configurations and could greatly benefit ongoing and future HP 13 C MRI research through major aSNR improvements.

Coil combination methods for multi-channel hyperpolarized 13C imaging data from human studies.

Effective coil combination methods for human hyperpolarized 13C spectroscopy multi-channel data had been relatively unexplored. This study implemented and tested several coil combination methods, including (1) the sum-of-squares (SOS), (2) singular value decomposition (SVD), (3) Roemer method by using reference peak area as a sensitivity map (RefPeak), and (4) Roemer method by using ESPIRiT-derived sensitivity map (ESPIRiT). These methods were evaluated by numerical simulation, thermal phantom experiments, and human cancer patient studies. Overall, the SVD, RefPeak, and ESPIRiT methods demonstrated better accuracy and robustness than the SOS method. Extracting complex pyruvate signal provides an easy and excellent approximation of the coil sensitivity map while maintaining valuable phase information of the coil-combined data.

Comparison between 8- and 32-channel phased-array receive coils for in vivo hyperpolarized 13 C imaging of the human brain.

PURPOSE: To compare the performance of an 8-channel surface coil/clamshell transmitter and 32-channel head array coil/birdcage transmitter for hyperpolarized 13 C brain metabolic imaging. METHODS: To determine the field homogeneity of the radiofrequency transmitters, B1 + mapping was performed on an ethylene glycol head phantom and evaluated by means of the double angle method. Using a 3D echo-planar imaging sequence, coil sensitivity and noise-only phantom data were acquired with the 8- and 32-channel receiver arrays, and compared against data from the birdcage in transceiver mode. Multislice frequency-specific 13 C dynamic echo-planar imaging was performed on a patient with a brain tumor for each hardware configuration following injection of hyperpolarized [1-13 C]pyruvate. Signal-to-noise ratio (SNR) was evaluated from pre-whitened phantom and temporally summed patient data after coil combination based on optimal weights. RESULTS: The birdcage transmitter produced more uniform B1 + compared with the clamshell: 0.07 versus 0.12 (fractional error). Phantom experiments conducted with matched lateral housing separation demonstrated 8- versus 32-channel mean transceiver-normalized SNR performance: 0.91 versus 0.97 at the head center; 6.67 versus 2.08 on the sides; 0.66 versus 2.73 at the anterior; and 0.67 versus 3.17 on the posterior aspect. While the 8-channel receiver array showed SNR benefits along lateral aspects, the 32-channel array exhibited greater coverage and a more uniform coil-combined profile. Temporally summed, parameter-normalized patient data showed SNRmean,slice ratios (8-channel/32-channel) ranging 0.5-2.00 from apical to central brain. White matter lactate-to-pyruvate ratios were conserved across hardware: 0.45 ± 0.12 (8-channel) versus 0.43 ± 0.14 (32-channel). CONCLUSION: The 8- and 32-channel hardware configurations each have advantages in particular brain anatomy.

Measuring Tumor Metabolism in Pediatric Diffuse Intrinsic Pontine Glioma Using Hyperpolarized Carbon-13 MR Metabolic Imaging.

Objective: The purpose of this study was to demonstrate the feasibility of using hyperpolarized carbon-13 (13C) metabolic imaging with [1-13C]-labeled pyruvate for evaluating real-time in vivo metabolism of orthotopic diffuse intrinsic pontine glioma (DIPG) xenografts. Materials and Methods: 3D 13C magnetic resonance spectroscopic imaging (MRSI) data were acquired on a 3T scanner from 8 rats that had been implanted with human-derived DIPG cells in the brainstem and 5 healthy controls, following injection of 2.5 mL (100 mM) hyperpolarized [1-13C]-pyruvate. Results: Anatomical images from DIPG-bearing rats characteristically exhibited T2-hyperintensity throughout the cerebellum and pons that was not accompanied by contrast enhancement. Evaluation of real-time in vivo13C spectroscopic data revealed ratios of lactate-to-pyruvate (p < 0.002), lactate-to-total carbon (p < 0.002), and normalized lactate (p < 0.002) that were significantly higher in T2 lesions harboring tumor relative to corresponding values of healthy normal brain. Elevated levels of lactate in lesions demonstrated a distinct metabolic profile that was associated with infiltrative, viable tumor recapitulating the histopathology of pediatric DIPG. Conclusions: Results from this study characterized pyruvate and lactate metabolism in orthotopic DIPG xenografts and suggest that hyperpolarized 13C MRSI may serve as a noninvasive imaging technique for in vivo monitoring of biochemical processes in patients with DIPG.