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

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.

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.

Quo Vadis Hyperpolarized 13C MRI?

Over the last two decades, hyperpolarized 13C MRI has gained significance in both preclinical and clinical studies, hereby relying on technologies like PHIP-SAH (ParaHydrogen-Induced Polarization-Side Arm Hydrogenation), SABRE (Signal Amplification by Reversible Exchange), and dDNP (dissolution Dynamic Nuclear Polarization), with dDNP being applied in humans. A clinical dDNP polarizer has enabled studies across 24 sites, despite challenges like high cost and slow polarization. Parahydrogen-based techniques like SABRE and PHIP offer faster, more cost-efficient alternatives but require molecule-specific optimization. The focus has been on imaging metabolism of hyperpolarized probes, which requires long T1, high polarization and rapid contrast generation. Efforts to establish novel probes, improve acquisition techniques and enhance data analysis methods including artificial intelligence are ongoing. Potential clinical value of hyperpolarized 13C MRI was demonstrated primarily for treatment response assessment in oncology, but also in cardiology, nephrology, hepatology and CNS characterization. In this review on biomedical hyperpolarized 13C MRI, we summarize important and recent advances in polarization techniques, probe development, acquisition and analysis methods as well as clinical trials. Starting from those we try to sketch a trajectory where the field of biomedical hyperpolarized 13C MRI might go.

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.

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.

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.

In vivo deuterium magnetic resonance imaging of xenografted tumors following systemic administration of deuterated water

In vivo deuterated water (2H2O) labeling leads to deuterium (2H) incorporation into biomolecules of proliferating cells and provides the basis for its use in cell kinetics research. We hypothesized that rapidly proliferating cancer cells would become preferentially labeled with 2H and, therefore, could be visualized by deuterium magnetic resonance imaging (dMRI) following a brief period of in vivo systemic 2H2O administration. We initiated systemic 2H2O administration in two xenograft mouse models harboring either human colorectal, HT-29, or pancreatic, MiaPaCa-2, tumors and 2H2O level of ~ 8% in total body water (TBW). Three schemas of 2H2O administration were tested: (1) starting at tumor seeding and continuing for 7 days of in vivo growth with imaging on day 7, (2) starting at tumor seeding and continuing for 14 days of in vivo growth with imaging on day 14, and (3) initiation of labeling following a week of in vivo tumor growth and continuing until imaging was performed on day 14. Deuterium chemical shift imaging of the tumor bearing limb and contralateral control was performed on either day 7 of 14 after tumor seeding, as described. After 14 days of in vivo tumor growth and 7 days of systemic labeling with 2H2O, a clear deuterium contrast was demonstrated between the xenografts and normal tissue. Labeling in the second week after tumor implantation afforded the highest contrast between neoplastic and healthy tissue in both models. Systemic labeling with 2H2O can be used to create imaging contrast between tumor and healthy issue, providing a non-radioactive method for in vivo cancer imaging.

Hyperpolarized δ-[1- 13C]gluconolactone imaging visualizes response to TERT or GABPB1 targeting therapy for glioblastoma

TERT promoter mutations are a hallmark of glioblastoma (GBM). Accordingly, TERT and GABPB1, a subunit of the upstream mutant TERT promoter transcription factor GABP, are being considered as promising therapeutic targets in GBM. We recently reported that the expression of TERT or GABP1 modulates flux via the pentose phosphate pathway (PPP). Here, we investigated whether 13C magnetic resonance spectroscopy (MRS) of hyperpolarized (HP) δ- [1-13C]gluconolactone can serve to image the reduction in PPP flux following TERT or GABPB1 silencing. We investigated two different human GBM cell lines stably expressing shRNAs targeting TERT or GABPB1, as well as doxycycline-inducible shTERT or shGABPB1cells. MRS studies were performed on live cells and in vivo tumors, and dynamic sets of 13C MR spectra were acquired following injection of HP δ-[1-13C]gluconolactone. HP 6-phosphogluconolactone (6PG), the product of δ-[1-13C]gluconolactone via the PPP, was significantly reduced in TERT or GABPB1-silenced cells or tumors compared to controls in all our models. Furthermore, a positive correlation between TERT expression and 6PG levels was observed. Our data indicate that HP δ-[1-13C]gluconolactone, an imaging tool with translational potential, could serve to monitor TERT expression and its silencing with therapies that target either TERT or GABPB1 in mutant TERT promoter GBM patients.

Tensor denoising of multidimensional MRI data

PURPOSE

To develop a denoising strategy leveraging redundancy in high-dimensional data.

THEORY AND METHODS

The SNR fundamentally limits the information accessible by MRI. This limitation has been addressed by a host of denoising techniques, recently including the so-called MPPCA: principal component analysis of the signal followed by automated rank estimation, exploiting the Marchenko-Pastur distribution of noise singular values. Operating on matrices comprised of data patches, this popular approach objectively identifies noise components and, ideally, allows noise to be removed without introducing artifacts such as image blurring, or nonlocal averaging. The MPPCA rank estimation, however, relies on a large number of noise singular values relative to the number of signal components to avoid such ill effects. This condition is unlikely to be met when data patches and therefore matrices are small, for example due to spatially varying noise. Here, we introduce tensor MPPCA (tMPPCA) for the purpose of denoising multidimensional data, such as from multicontrast acquisitions. Rather than combining dimensions in matrices, tMPPCA uses each dimension of the multidimensional data's inherent tensor-structure to better characterize noise, and to recursively estimate signal components.

RESULTS

Relative to matrix-based MPPCA, tMPPCA requires no additional assumptions, and comparing the two in a numerical phantom and a multi-TE diffusion MRI data set, tMPPCA dramatically improves denoising performance. This is particularly true for small data patches, suggesting that tMPPCA can be especially beneficial in such cases.

CONCLUSIONS

The MPPCA denoising technique can be extended to high-dimensional data with improved performance for smaller patch sizes.

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.

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.

Imaging biomarkers of TERT or GABPB1 silencing in TERT-positive glioblastoma

BACKGROUND

TERT promoter mutations are observed in 80% of wild-type IDH glioblastoma (GBM). Moreover, the upstream TERT transcription factor GABPB1 was recently identified as a cancer-specific therapeutic target for tumors harboring a TERT promoter mutation. In that context, noninvasive imaging biomarkers are needed for the detection of TERT modulation.

METHODS

Multiple GBM models were investigated as cells and in vivo tumors and the impact of TERT silencing, either directly or by targeting GABPB1, was determined using 1H and hyperpolarized 13C magnetic resonance spectroscopy (MRS). Changes in associated metabolic enzymes were also investigated.

RESULTS

1H-MRS revealed that lactate and glutathione (GSH) were the most significantly altered metabolites when either TERT or GABPB1 was silenced, and lactate and GSH levels were correlated with cellular TERT expression. Consistent with the drop in lactate, 13C-MRS showed that hyperpolarized [1-13C]lactate production from [1-13C]pyruvate was also reduced when TERT was silenced. Mechanistically, the reduction in GSH was associated with a reduction in pentose phosphate pathway flux, reduced activity of glucose-6-phosphate dehydrogenase, and reduced NADPH. The drop in lactate and hyperpolarized lactate were associated with reductions in glycolytic flux, NADH, and expression/activity of GLUT1, monocarboxylate transporters, and lactate dehydrogenase A.

CONCLUSIONS

Our study indicates that MRS-detectable GSH, lactate, and lactate production could serve as metabolic biomarkers of response to emerging TERT-targeted therapies for GBM with activating TERT promoter mutations. Importantly these biomarkers are readily translatable to the clinic, and thus could ultimately improve GBM patient management.

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.

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.

Post-Acquisition Hyperpolarized 29Silicon Magnetic Resonance Image Processing for Visualization of Colorectal Lesions Using a User-Friendly Graphical Interface

Medical imaging devices often use automated processing that creates and displays a self-normalized image. When improperly executed, normalization can misrepresent information or result in an inaccurate analysis. In the case of diagnostic imaging, a false positive in the absence of disease, or a negative finding when disease is present, can produce a detrimental experience for the patient and diminish their health prospects and prognosis. In many clinical settings, a medical technical specialist is trained to operate an imaging device without sufficient background information or understanding of the fundamental theory and processes involved in image creation and signal processing. Here, we describe a user-friendly image processing algorithm that mitigates user bias and allows for true signal to be distinguished from background. For proof-of-principle, we used antibody-targeted molecular imaging of colorectal cancer (CRC) in a mouse model, expressing human MUC1 at tumor sites. Lesion detection was performed using targeted magnetic resonance imaging (MRI) of hyperpolarized silicon particles. Resulting images containing high background and artifacts were then subjected to individualized image post-processing and comparative analysis. Post-acquisition image processing allowed for co-registration of the targeted silicon signal with the anatomical proton magnetic resonance (MR) image. This new methodology allows users to calibrate a set of images, acquired with MRI, and reliably locate CRC tumors in the lower gastrointestinal tract of living mice. The method is expected to be generally useful for distinguishing true signal from background for other cancer types, improving the reliability of diagnostic MRI.

Hyperpolarized 13C-Pyruvate Metabolism as a Surrogate for Tumor Grade and Poor Outcome in Renal Cell Carcinoma-A Proof of Principle Study

Differentiating aggressive clear cell renal cell carcinoma (ccRCC) from indolent lesions is challenging using conventional imaging. This work prospectively compared the metabolic imaging phenotype of renal tumors using carbon-13 MRI following injection of hyperpolarized [1-13C]pyruvate (HP-13C-MRI) and validated these findings with histopathology. Nine patients with treatment-naïve renal tumors (6 ccRCCs, 1 liposarcoma, 1 pheochromocytoma, 1 oncocytoma) underwent pre-operative HP-13C-MRI and conventional proton (1H) MRI. Multi-regional tissue samples were collected using patient-specific 3D-printed tumor molds for spatial registration between imaging and molecular analysis. The apparent exchange rate constant (kPL) between 13C-pyruvate and 13C-lactate was calculated. Immunohistochemistry for the pyruvate transporter (MCT1) from 44 multi-regional samples, as well as associations between MCT1 expression and outcome in the TCGA-KIRC dataset, were investigated. Increasing kPL in ccRCC was correlated with increasing overall tumor grade (ρ = 0.92, p = 0.009) and MCT1 expression (r = 0.89, p = 0.016), with similar results acquired from the multi-regional analysis. Conventional 1H-MRI parameters did not discriminate tumor grades. The correlation between MCT1 and ccRCC grade was confirmed within a TCGA dataset (p < 0.001), where MCT1 expression was a predictor of overall and disease-free survival. In conclusion, metabolic imaging using HP-13C-MRI differentiates tumor aggressiveness in ccRCC and correlates with the expression of MCT1, a predictor of survival. HP-13C-MRI may non-invasively characterize metabolic phenotypes within renal cancer.

Hypoxia Imaging As a Guide for Hypoxia-Modulated and Hypoxia-Activated Therapy

Significance: Oxygen imaging techniques, which can probe the spatiotemporal heterogeneity of tumor oxygenation, could be of significant clinical utility in radiation treatment planning and in evaluating the effectiveness of hypoxia-activated prodrugs. To fulfill these goals, oxygen imaging techniques should be noninvasive, quantitative, and capable of serial imaging, as well as having sufficient temporal resolution to detect the dynamics of tumor oxygenation to distinguish regions of chronic and acute hypoxia. Recent Advances: No current technique meets all these requirements, although all have strengths in certain areas. The current status of positron emission tomography (PET)-based hypoxia imaging, oxygen-enhanced magnetic resonance imaging (MRI), ¹⁹F MRI, and electron paramagnetic resonance (EPR) oximetry are reviewed along with their strengths and weaknesses for planning hypoxia-guided, intensity-modulated radiation therapy and detecting treatment response for hypoxia-targeted prodrugs. Critical Issues: Spatial and temporal resolution emerges as a major concern for these areas along with specificity and quantitative response. Although multiple oxygen imaging techniques have reached the investigative stage, clinical trials to test the therapeutic effectiveness of hypoxia imaging have been limited. Future Directions: Imaging elements of the redox environment besides oxygen by EPR and hyperpolarized MRI may have a significant impact on our understanding of the basic biology of the reactive oxygen species response and may extend treatment possibilities.

Acquisition and quantification pipeline for in vivo hyperpolarized 13 C MR spectroscopy

PURPOSE

The goal of this study was to combine a specialized acquisition method with a new quantification pipeline to accurately and efficiently probe the metabolism of hyperpolarized ¹³ C-labeled compounds in vivo. In this study, we tested our approach on [2-¹³ C]pyruvate and [1-¹³ C]α-ketoglutarate data in rat orthotopic brain tumor models at 3T.

METHODS

We used a multiband metabolite-specific radiofrequency (RF) excitation in combination with a variable flip angle scheme to minimize substrate polarization loss and measure fast metabolic processes. We then applied spectral-temporal denoising using singular value decomposition to enhance spectral quality. This was combined with LCModel-based automatic ¹³ C spectral fitting and flip angle correction to separate overlapping signals and rapidly quantify the different metabolites.

RESULTS

Denoising improved the metabolite signal-to-noise ratio (SNR) by approximately 5. It also improved the accuracy of metabolite quantification as evidenced by a significant reduction of the Cramer Rao lower bounds. Furthermore, the use of the automated and user-independent LCModel-based quantification approach could be performed rapidly, with the kinetic quantification of eight metabolite peaks in a 12-spectrum array achieved in less than 1 minute.

CONCLUSION

The specialized acquisition method combined with denoising and a new quantification pipeline using LCModel for the first time for hyperpolarized ¹³ C data enhanced our ability to monitor the metabolism of [2-¹³ C]pyruvate and [1-¹³ C]α-ketoglutarate in rat orthotopic brain tumor models in vivo. This approach could be broadly applicable to other hyperpolarized agents both preclinically and in the clinical setting.

Current human brain applications and challenges of dynamic hyperpolarized carbon-13 labeled pyruvate MR metabolic imaging

The ability of hyperpolarized carbon-13 MR metabolic imaging to acquire dynamic metabolic information in real time is crucial to gain mechanistic insights into metabolic pathways, which are complementary to anatomic and other functional imaging methods. This review presents the advantages of this emerging functional imaging technology, describes considerations in clinical translations, and summarizes current human brain applications. Despite rapid development in methodologies, significant technological and physiological related challenges continue to impede broader clinical translation.

Specialized computational methods for denoising, B1 correction, and kinetic modeling in hyperpolarized 13 C MR EPSI studies of liver tumors

PURPOSE

To develop a novel post-processing pipeline for hyperpolarized (HP) ¹³ C MRSI that integrates tensor denoising and B 1 + correction to measure pyruvate-to-lactate conversion rates (k_PL ) in patients with liver tumors.

METHODS

Seven HP ¹³ C MR scans of progressing liver tumors were acquired using a custom ¹³ C surface transmit/receive coil and the echo-planar spectroscopic imaging (EPSI) data analysis included B₀ correction, tensor rank truncation, and zero- and first-order phase corrections to recover metabolite signals that would otherwise be obscured by spectral noise as well as a correction for inhomogeneous transmit ( B 1 + ) using a B 1 + map aligned to the coil position for each patient scan. Processed HP data and corrected flip angles were analyzed with an inputless two-site exchange model to calculate k_PL .

RESULTS

Denoising averages SNR increases of pyruvate, lactate, and alanine were 37.4-, 34.0-, and 20.1-fold, respectively, with lactate and alanine dynamics most noticeably recovered and better defined. In agreement with Monte Carlo simulations, over-flipped regions underestimated k_PL and under-flipped regions overestimated k_PL . B 1 + correction addressed this issue.

CONCLUSION

The new HP ¹³ C EPSI post-processing pipeline integrated tensor denoising and B 1 + correction to measure k_PL in patients with liver tumors. These technical developments not only recovered metabolite signals in voxels that did not receive the prescribed flip angle, but also increased the extent and accuracy of k_PL estimations throughout the tumor and adjacent regions including normal-appearing tissue and additional lesions.

Denoising of hyperpolarized 13 C MR images of the human brain using patch-based higher-order singular value decomposition

PURPOSE

To improve hyperpolarized ¹³ C (HP-¹³ C) MRI by image denoising with a new approach, patch-based higher-order singular value decomposition (HOSVD).

METHODS

The benefit of using a patch-based HOSVD method to denoise dynamic HP-¹³ C MR imaging data was investigated. Image quality and the accuracy of quantitative analyses following denoising were evaluated first using simulated data of [1-¹³ C]pyruvate and its metabolic product, [1-¹³ C]lactate, and compared the results to a global HOSVD method. The patch-based HOSVD method was then applied to healthy volunteer HP [1-¹³ C]pyruvate EPI studies. Voxel-wise kinetic modeling was performed on both non-denoised and denoised data to compare the number of voxels quantifiable based on SNR criteria and fitting error.

RESULTS

Simulation results demonstrated an 8-fold increase in the calculated SNR of [1-¹³ C]pyruvate and [1-¹³ C]lactate with the patch-based HOSVD denoising. The voxel-wise quantification of k_PL (pyruvate-to-lactate conversion rate) showed a 9-fold decrease in standard errors for the fitted k_PL after denoising. The patch-based denoising performed superior to the global denoising in recovering k_PL information. In volunteer data sets, [1-¹³ C]lactate and [¹³ C]bicarbonate signals became distinguishable from noise across captured time points with over a 5-fold apparent SNR gain. This resulted in >3-fold increase in the number of voxels quantifiable for mapping k_PB (pyruvate-to-bicarbonate conversion rate) and whole brain coverage for mapping k_PL .

CONCLUSIONS

Sensitivity enhancement provided by this denoising significantly improved quantification of metabolite dynamics and could benefit future studies by improving image quality, enabling higher spatial resolution, and facilitating the extraction of metabolic information for clinical research.

Schulte R, Laustsen C. Comprehensive Literature Review of Hyperpolarized Carbon-13 MRI: The Road to Clinical Application

This review provides a comprehensive assessment of the development of hyperpolarized (HP) carbon-13 metabolic MRI from the early days to the present with a focus on clinical applications. The status and upcoming challenges of translating HP carbon-13 into clinical application are reviewed, along with the complexity, technical advancements, and future directions. The road to clinical application is discussed regarding clinical needs and technological advancements, highlighting the most recent successes of metabolic imaging with hyperpolarized carbon-13 MRI. Given the current state of hyperpolarized carbon-13 MRI, the conclusion of this review is that the workflow for hyperpolarized carbon-13 MRI is the limiting factor.

PCA denoising and Wiener deconvolution of 31 P 3D CSI data to enhance effective SNR and improve point spread function

Purpose

This study evaluates the performance of 2 processing methods, that is, principal component analysis‐based denoising and Wiener deconvolution, to enhance the quality of phosphorus 3D chemical shift imaging data.

Methods

Principal component analysis‐based denoising increases the SNR while maintaining spectral information. Wiener deconvolution reduces the FWHM of the voxel point spread function, which is increased by Hamming filtering or Hamming‐weighted acquisition. The proposed methods are evaluated using simulated and in vivo 3D phosphorus chemical shift imaging data by 1) visual inspection of the spatial signal distribution; 2) SNR calculation of the PCr peak; and 3) fitting of metabolite basis functions.

Results

With the optimal order of processing steps, we show that the effective SNR of in vivo phosphorus 3D chemical shift imaging data can be increased. In simulations, we show we can preserve phosphorus‐containing metabolite peaks that had an SNR < 1 before denoising. Furthermore, using Wiener deconvolution, we were able to reduce the FWHM of the voxel point spread function with only partially reintroducing Gibb‐ringing artifacts while maintaining the SNR. After data processing, fitting of the phosphorus‐containing metabolite signals improved.

Conclusion

In this study, we have shown that principal component analysis‐based denoising in combination with regularized Wiener deconvolution allows increasing the effective spectral SNR of in vivo phosphorus 3D chemical shift imaging data, with reduction of the FWHM of the voxel point spread function. Processing increased the effective SNR by at least threefold compared to Hamming weighted acquired data and minimized voxel bleeding. With these methods, fitting of metabolite amplitudes became more robust with decreased fitting residuals.

Pilot Study of Hyperpolarized 13C Metabolic Imaging in Pediatric Patients with Diffuse Intrinsic Pontine Glioma and Other CNS Cancers

BACKGROUND AND PURPOSE

Pediatric CNS tumors commonly present challenges for radiographic interpretation on conventional MR imaging. This study sought to investigate the safety and tolerability of hyperpolarized carbon-13 (HP-¹³C) metabolic imaging in pediatric patients with brain tumors.

MATERIALS AND METHODS

Pediatric patients 3 to 18 years of age who were previously diagnosed with a brain tumor and could undergo MR imaging without sedation were eligible to enroll in this safety study of HP [1-¹³C]pyruvate. Participants received a one-time injection of HP [1-¹³C]pyruvate and were imaged using dynamic HP-¹³C MR imaging. We assessed 2 dose levels: 0.34 mL/kg and the highest tolerated adult dose of 0.43 mL/kg. Participants were monitored throughout imaging and for 60 minutes postinjection, including pre- and postinjection electrocardiograms and vital sign measurements.

RESULTS

Between February 2017 and July 2019, ten participants (9 males; median age, 14 years; range, 10-17 years) were enrolled, of whom 6 completed injection of HP [1-¹³C]pyruvate and dynamic HP-¹³C MR imaging. Four participants failed to undergo HP-¹³C MR imaging due to technical failures related to generating HP [1-¹³C]pyruvate or MR imaging operability. HP [1-¹³C]pyruvate was well-tolerated in all participants who completed the study, with no dose-limiting toxicities or adverse events observed at either 0.34 (n = 3) or 0.43 (n = 3) mL/kg. HP [1-¹³C]pyruvate demonstrated characteristic conversion to [1-¹³C]lactate and [¹³C]bicarbonate in the brain. Due to poor accrual, the study was closed after only 3 participants were enrolled at the highest dose level.

CONCLUSIONS

Dynamic HP-¹³C MR imaging was safely performed in 6 pediatric patients with CNS tumors and demonstrated HP [1-¹³C]pyruvate brain metabolism.

Comparison of hyperpolarized 13 C and non-hyperpolarized deuterium MRI approaches for imaging cerebral glucose metabolism at 4.7 T

PURPOSE

The purpose of this study was to directly compare two isotopic metabolic imaging approaches, hyperpolarized (HP) ¹³ C MRI and deuterium metabolic imaging (DMI), for imaging specific closely related segments of cerebral glucose metabolism at 4.7 T.

METHODS

Comparative HP-¹³ C and DMI neuroimaging experiments were conducted consecutively in normal rats during the same scanning session. Localized conversions of [1-¹³ C]pyruvate and [6,6-² H₂ ]glucose to their respective downstream metabolic products were measured by spectroscopic imaging, using an identical 2D-CSI sequence with parameters optimized for the respective experiments. To facilitate direct comparison, a pair of substantially equivalent 2.5-cm double-tuned X/¹ H RF surface coils was developed. For improved results, multidimensional low-rank reconstruction was applied to denoise the raw DMI data.

RESULTS

Localized conversion of HP [1-¹³ C]pyruvate to [1-¹³ C]lactate, and [6,6-² H₂ ]glucose to [3,3-² H₂ ]lactate and Glx-d (glutamate and glutamine), was detected in rat brain by spectroscopic imaging at 4.7 T. The SNR and spatial resolution of HP-¹³ C MRI was superior to DMI but limited to a short time window, whereas the lengthy DMI acquisition yielded maps of not only lactate, but also Glx production, albeit with relatively poor spectral discrimination between metabolites at this field strength. Across the individual rats, there was an apparent inverse correlation between cerebral production of HP [1-¹³ C]lactate and Glx-d, along with a trend toward increased [3,3-² H₂ ]lactate.

CONCLUSION

The HP-¹³ C MRI and DMI methods are both feasible at 4.7 T and have significant potential for metabolic imaging of specific segments of glucose metabolism.