A referenceless Nyquist ghost correction workflow for echo planar imaging of hyperpolarized [1-13 C]pyruvate and [1-13 C]lactate

MR-Nucleomics: The study of pathological cellular processes with multinuclear magnetic resonance spectroscopy and imaging in vivo

Clinical medicine has experienced a rapid development in recent decades, during which therapies targeting specific cellular signaling pathways, or specific cell surface receptors, have been increasingly adopted. While these developments in clinical medicine call for improved precision in diagnosis and treatment monitoring, modern medical imaging methods are restricted mainly to anatomical imaging, lagging behind the requirements of precision medicine. Although positron emission tomography and single photon emission computed tomography have been used clinically for studies of metabolism, their applications have been limited by the exposure risk to ionizing radiation, the subsequent limitation in repeated and longitudinal studies, and the incapability in assessing downstream metabolism. Magnetic resonance spectroscopy (MRS) or spectroscopic imaging (MRSI) are, in theory, capable of assessing molecular activities in vivo, although they are often limited by sensitivity. Here, we review some recent developments in MRS and MRSI of multiple nuclei that have potential as molecular imaging tools in the clinic.

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.

Hyperpolarized Metabolic MRI-Acquisition, Reconstruction, and Analysis Methods

Hyperpolarized metabolic MRI with 13C-labeled agents has emerged as a powerful technique for in vivo assessments of real-time metabolism that can be used across scales of cells, tissue slices, animal models, and human subjects. Hyperpolarized contrast agents have unique properties compared to conventional MRI scanning and MRI contrast agents that require specialized imaging methods. Hyperpolarized contrast agents have a limited amount of available signal, irreversible decay back to thermal equilibrium, bolus injection and perfusion kinetics, cellular uptake and metabolic conversion kinetics, and frequency shifts between metabolites. This article describes state-of-the-art methods for hyperpolarized metabolic MRI, summarizing data acquisition, reconstruction, and analysis methods in order to guide the design and execution of studies.

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-¹³ 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-¹³ 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-¹³ 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

•

Serial HP ¹³C MRI was evaluated for data consistency and abnormal metabolism.

•

Metabolism of [1-¹³C]pyruvate to lactate and bicarbonate was kinetically modeled.

•

Conversion rates within NAWM were consistent in healthy volunteer and patient scans

•

Progressed tumor lesions showed higher relative conversion rates to [1-¹³C]lactate.

•

Globally elevated rate constants were observed with anti-angiogenic treatment.

Background

Hyperpolarized carbon-13 (HP-¹³C) 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-¹³C 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-¹³C]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-¹³C imaging at routine clinical timepoints overlapped treatment, including radiotherapy (RT), temozolomide (TMZ) chemotherapy, and anti-angiogenic/investigational agents. Apparent rate constants for [1-¹³C]pyruvate conversion to [1-¹³C]lactate (k_PL) and [¹³C]bicarbonate (k_PB) were simultaneously quantified based on an inputless kinetic model within normal-appearing white matter (NAWM) and anatomic lesions defined from ¹H MRI. The inter/intra-subject consistency of k_PL-NAWM and k_PB-NAWM was measured in terms of the coefficient of variation (CV).

Results

When excluding scans following anti-angiogenic therapy, patient values of k_PL-NAWM and k_PB-NAWM were 0.020 s⁻¹ ± 23.8% and 0.0058 s⁻¹ ± 27.7% (mean ± CV) across 17 HP-¹³C MRIs, with intra-patient serial k_PL-NAWM/k_PB-NAWM CVs ranging 6.8–16.6%/10.6–40.7%. In 4/5 patients, these values (0.018 s⁻¹ ± 13.4% and 0.0058 s⁻¹ ± 24.4%; n = 13) were more similar to those from healthy volunteers (0.018 s⁻¹ ± 5.0% and 0.0043 s⁻¹ ± 12.6%; n = 6) (mean ± CV). The anti-angiogenic agent bevacizumab was associated with global elevations in apparent rate constants, with maximum k_PL-NAWM in 2 patients reaching 0.047 ± 0.001 and 0.047 ± 0.003 s⁻¹ (±model error). In 3 patients with progressive disease, anatomic lesions showed elevated k_PL relative to k_PL-NAWM of 0.024 ± 0.001 s⁻¹ (±model error) in the absence of gadolinium enhancement, and 0.032 ± 0.008, 0.040 ± 0.003 and 0.041 ± 0.009 s⁻¹ with gadolinium enhancement. The lesion k_PB in patients was reduced to unquantifiable values compared to k_PB-NAWM.

Conclusion

Serial measures of HP [1-¹³C]pyruvate metabolism displayed consistency in the NAWM of healthy volunteers and patients. Both k_PL and k_PB were globally elevated following bevacizumab treatment, while progressive disease demonstrated elevated k_PL in gadolinium-enhancing and non-enhancing lesions. Larger prospective studies with homogeneous patient populations are planned to evaluate metabolic changes following treatment.

Hyperpolarized 13C MRI: A novel approach for probing cerebral metabolism in health and neurological disease

Cerebral metabolism is tightly regulated and fundamental for healthy neurological function. There is increasing evidence that alterations in this metabolism may be a precursor and early biomarker of later stage disease processes. Proton magnetic resonance spectroscopy (1H-MRS) is a powerful tool to non-invasively assess tissue metabolites and has many applications for studying the normal and diseased brain. However, the technique has limitations including low spatial and temporal resolution, difficulties in discriminating overlapping peaks, and challenges in assessing metabolic flux rather than steady-state concentrations. Hyperpolarized carbon-13 magnetic resonance imaging is an emerging clinical technique that may overcome some of these spatial and temporal limitations, providing novel insights into neurometabolism in both health and in pathological processes such as glioma, stroke and multiple sclerosis. This review will explore the growing body of pre-clinical data that demonstrates a potential role for the technique in assessing metabolism in the central nervous system. There are now a number of clinical studies being undertaken in this area and this review will present the emerging clinical data as well as the potential future applications of hyperpolarized 13C magnetic resonance imaging in the brain, in both clinical and pre-clinical studies.

Fast Imaging for Hyperpolarized MR Metabolic Imaging

MRI with hyperpolarized carbon-13 agents has created a new type of noninvasive, in vivo metabolic imaging that can be applied in cell, animal, and human studies. The use of ¹³ C-labeled agents, primarily [1-¹³ C]pyruvate, enables monitoring of key metabolic pathways with the ability to image substrate and products based on their chemical shift. Over 10 sites worldwide are now performing human studies with this new approach for studies of cancer, heart disease, liver disease, and kidney disease. Hyperpolarized metabolic imaging studies must be performed within several minutes following creation of the hyperpolarized agent due to irreversible decay of the net magnetization back to equilibrium, so fast imaging methods are critical. The imaging methods must include multiple metabolites, separated based on their chemical shift, which are also undergoing rapid metabolic conversion (via label exchange), further exacerbating the challenges of fast imaging. This review describes the state-of-the-art in fast imaging methods for hyperpolarized metabolic imaging. This includes the approach and tradeoffs between three major categories of fast imaging methods-fast spectroscopic imaging, model-based strategies, and metabolite specific imaging-as well additional options of parallel imaging, compressed sensing, tailored RF flip angles, refocused imaging methods, and calibration methods that can improve the scan coverage, speed, signal-to-noise ratio (SNR), resolution, and/or robustness of these studies. To date, these approaches have produced extremely promising initial human imaging results. Improvements to fast hyperpolarized metabolic imaging methods will provide better coverage, SNR, resolution, and reproducibility for future human imaging studies. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY STAGE: 1.

Partial Fourier reconstruction for improved resolution in 3D hyperpolarized 13 C EPI

PURPOSE

Asymmetric in-plane k-space sampling of EPI can reduce the minimum achievable TE in hyperpolarized 13 C with spectral-spatial radio frequency pulses, thereby reducing T 2 * weighting and signal-losses. Partial Fourier image reconstruction exploits the approximate Hermitian symmetry of k-space data and can be applied to asymmetric data sets to synthesize unmeasured data. Here we tested whether the application of partial Fourier image reconstruction would improve spatial resolution from hyperpolarized [1- 13 C ]pyruvate scans in the human brain.

METHODS

Fifteen healthy control subjects were imaged using a volumetric dual-echo echo-planar imaging sequence with spectral-spatial radio frequency excitation. Images were reconstructed by zero-filling as well as with the partial Fourier reconstruction algorithm projection-on-convex-sets. Resulting images were quantitatively evaluated with a no-reference image quality assessment.

RESULTS

The no-reference image sharpness metric agreed with perceived improvements in image resolution and contrast. The [1- 13 C ]lactate images benefitted most, followed by the [1- 13 C ]pyruvate images. The 13 C -bicarbonate images were improved by the smallest degree, likely owing to relatively lower SNR.

CONCLUSIONS

Partial Fourier imaging and reconstruction were shown to improve the sharpness and contrast of human HP 13 C brain data and is a viable method for enhancing resolution.

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

Effective coil combination methods for human hyperpolarized ¹³C 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.

Translation of Carbon-13 EPI for hyperpolarized MR molecular imaging of prostate and brain cancer patients

PURPOSE

To develop and translate a metabolite-specific imaging sequence using a symmetric echo planar readout for clinical hyperpolarized (HP) Carbon-13 (¹³ C) applications.

METHODS

Initial data were acquired from patients with prostate cancer (N = 3) and high-grade brain tumors (N = 3) on a 3T scanner. Samples of [1-¹³ C]pyruvate were polarized for at least 2 h using a 5T SPINlab system operating at 0.8 K. Following injection of the HP substrate, pyruvate, lactate, and bicarbonate (for brain studies) were sequentially excited with a singleband spectral-spatial RF pulse and signal was rapidly encoded with a single-shot echo planar readout on a slice-by-slice basis. Data were acquired dynamically with a temporal resolution of 2 s for prostate studies and 3 s for brain studies.

RESULTS

High pyruvate signal was seen throughout the prostate and brain, with conversion to lactate being shown across studies, whereas bicarbonate production was also detected in the brain. No Nyquist ghost artifacts or obvious geometric distortion from the echo planar readout were observed. The average error in center frequency was 1.2 ± 17.0 and 4.5 ± 1.4 Hz for prostate and brain studies, respectively, below the threshold for spatial shift because of bulk off-resonance.

CONCLUSION

This study demonstrated the feasibility of symmetric EPI to acquire HP ¹³ C metabolite maps in a clinical setting. As an advance over prior single-slice dynamic or single time point volumetric spectroscopic imaging approaches, this metabolite-specific EPI acquisition provided robust whole-organ coverage for brain and prostate studies while retaining high SNR, spatial resolution, and dynamic temporal resolution.