Next Article in Journal
Prognosis of Single Early-Stage Hepatocellular Carcinoma (HCC) with CEUS Inconclusive Imaging (LI-RADS LR-3 and LR-4) Is No Better than Typical HCC (LR-5)
Next Article in Special Issue
In Vivo MRI-CEST Tumor pH Imaging Detects Resistance to Proton Pump Inhibitors in Human Prostate Cancer Murine Models
Previous Article in Journal
The Role of N6-Methyladenosine (m6A) Methylation Modifications in Hematological Malignancies
Previous Article in Special Issue
Impact of Inhibition of the Mitochondrial Pyruvate Carrier on the Tumor Extracellular pH as Measured by CEST-MRI
 
 
cancers-logo
Article Menu

Article Menu

Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Stephan Ursprung
1,2,
Ramona Woitek
1,2,
Mary A. McLean
1,2,
Andrew N. Priest
2,3,
Mireia Crispin-Ortuzar
1,
Cara R. Brodie
1,
Andrew B. Gill
1,2,
Marcel Gehrung
1,
Lucian Beer
1,2,
Antony C. P. Riddick
4,
Johanna Field-Rayner
1,2,
James T. Grist
2,
Surrin S. Deen
1,2,3,
Frank Riemer
1,2,
Joshua D. Kaggie
1,2,
Fulvio Zaccagna
1,2,
Joao A. G. Duarte
1,2,
Matthew J. Locke
1,2,
Amy Frary
1,2,
Tevita F. Aho
4,
James N. Armitage
4,
Ruth Casey
5,
Iosif A. Mendichovszky
1,2,3,
Sarah J. Welsh
1,6,7,
Tristan Barrett
1,2,
Martin J. Graves
3,
Tim Eisen
1,6,
Thomas J. Mitchell
1,4,7,8,
Anne Y. Warren
1,9,
Kevin M. Brindle
1,
Evis Sala
1,2,
Grant D. Stewart
1,4,7,† and
Ferdia A. Gallagher
1,2,*,†
add Show full author list remove Hide full author list
1
Cancer Research UK Cambridge Centre, University of Cambridge, Cambridge CB2 0QQ, UK
2
Department of Radiology, University of Cambridge, Cambridge CB2 0QQ, UK
3
Department of Radiology, Addenbrooke’s Hospital, Cambridge University Hospitals NHS Foundation Trust, Cambridge CB2 0QQ, UK
4
Department of Urology, Addenbrooke’s Hospital, Cambridge University Hospitals NHS Foundation Trust, Cambridge CB2 0QQ, UK
5
Department of Endocrinology, Addenbrooke’s Hospital, Cambridge University Hospitals NHS Foundation Trust, Cambridge CB2 0QQ, UK
6
Department of Oncology, Addenbrooke’s Hospital, Cambridge University Hospitals NHS Foundation Trust, Cambridge CB2 0QQ, UK
7
Department of Surgery, University of Cambridge, Cambridge CB2 0QQ, UK
8
Wellcome Sanger Institute, Hinxton CB10 1RQ, UK
9
Department of Pathology, Addenbrooke’s Hospital, Cambridge University Hospitals NHS Foundation Trust, Cambridge CB2 0QQ, UK
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Cancers 2022, 14(2), 335; https://doi.org/10.3390/cancers14020335
Submission received: 15 November 2021 / Revised: 4 January 2022 / Accepted: 6 January 2022 / Published: 11 January 2022
(This article belongs to the Special Issue Imaging Cancer Metabolism)

Abstract

:

Simple Summary

We evaluated renal cancer with varying aggressive appearances on histology, using an emerging form of non-invasive metabolic MRI. This imaging technique assesses the uptake and metabolism of a breakdown product of glucose (pyruvate) labelled with hyperpolarized carbon-13. We show that pyruvate metabolism is dependent on the aggressiveness of an individual tumor and we provide a mechanism for this finding from tissue analysis of molecules influencing pyruvate metabolism, suggesting a role for its membrane transporter.

Abstract

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.

Graphical Abstract

1. Introduction

Renal cell carcinoma (RCC) is the twelfth commonest solid cancer globally and the sixth in the United Kingdom [1]. Clear cell RCC (ccRCC) is the most frequent histological subtype found in 75–80% of patients with RCC [2]. In the localized disease stage, partial or radical nephrectomy is the treatment of choice [3]. However, with the increasing detection of incidental small renal masses by cross-sectional imaging, the risk of overtreatment is increasing, because neither conventional imaging nor tissue biopsies predict aggressive tumor behavior reliably [4,5]. Although tumor grade is a leading predictor of progression, recurrence, and survival, tumor biopsies may underestimate the grade owing to sampling error [6,7]. An imaging biomarker that could non-invasively grade ccRCC and exclude differential diagnoses would reduce overtreatment and could have significant clinical impact. This imaging biomarker may be derived from novel image analysis methods and machine learning approaches such as those in radiomics: the extraction and evaluation of high-dimensional quantitative data from images have shown promise in the grading of ccRCC [8,9]. Alternatively, new imaging techniques may provide a direct readout of tumor biology and in this study, we explored whether novel methods for probing tumor metabolism could be used to phenotype renal cancer more accurately.
Renal cell cancer is driven by dysregulation of several pathways involving cellular energy metabolism [10]. Dysregulation of hypoxia-inducible factor 1-alpha (HIF-1α) pathways resulting from VHL mutation or promotor hypermethylation is present in 91% of ccRCC leading to profound alterations in glucose metabolism [11]. Increasing tumor grade is associated with increased lactate concentrations and expression of proteins involved in glycolysis and lactic acid fermentation [12,13]. Metabolomic analyses require invasive sampling of vascularized tissue and are subject to sampling error. In contrast, non-invasive imaging enables repeated, non-destructive assessment of comparatively large volumes of tissue. For example, positron emission tomography (PET) using the glucose analog [18F]FDG, can indirectly probe these alterations. However, [18F]FDG-PET cannot distinguish metabolites, provides no direct information on glycolytic flux, and is limited in RCC by the background urinary signal [14].
Hyperpolarized carbon-13 magnetic resonance imaging (HP-13C-MRI) is a new clinical method to quantify tissue metabolism [15]. The technique can image the conversion of HP-[1-13C]pyruvate to [1-13C]lactate, catalyzed by lactate dehydrogenase (LDH), and is a promising pre-clinical and clinical imaging biomarker of aggressiveness and early treatment response [16]. HP-13C-MRI has been used to image the normal human brain and heart, as well as cancers of the prostate, breast, kidney, pancreas, and brain [17,18,19,20,21]. Monocarboxylate transporters (MCTs) are involved in cellular energy metabolism through the transport of pyruvate, lactate, and ketone bodies across the plasma membrane [22]. Findings from human breast and prostate cancer have identified MCT1 expression as a rate-limiting factor for 13C-lactate formation in HP-13C-MRI [20,23]. In this proof-of-principle study, we demonstrate that HP-13C-MRI provides a non-invasive measure of tumor aggressiveness in ccRCC from multiregional biopsies and show that the imaging measures of metabolism correlate with the expression of MCT1 for the first time in participants with RCC. We subsequently validated these changes in gene expression and patient outcome data from The Cancer Genome Atlas Kidney Renal Clear Cell Carcinoma (TCGA-KIRC) database [24].

2. Materials and Methods

2.1. Recruitment and Ethics

Study participants were recruited at Addenbrooke’s Hospital, Cambridge University Hospitals NHS Foundation Trust, Cambridge UK and provided written informed consent for this ethically approved study protocol: MISSION Substudy in Renal Cancer (Molecular Imaging and Spectroscopy with Stable Isotopes in Oncology and Neurology; Cambridge South Research Ethics Committee, 15/EE/0378) and DIAMOND (Evaluation of biomarkers in urological diseases; 03/018). The inclusion criteria were: age ≥18 years, clinical suspicion of T1b + RCC, candidate for surgery, and Eastern Cooperative Oncology Group (ECOG) performance status ≤1. Key exclusion criteria were metabolic disorders, allergy to gadolinium-based contrast agents, pregnancy, and unsuitability for MRI. Patients may have harbored localized or metastatic disease as long as surgery was clinically indicated, and they have not received any pre-surgical treatment.

2.2. Hyperpolarized [1-13C]Pyruvate MRI Acquisition

The preparation of the 13C-pyruvate has been described previously [20] and are summarized in the Supplementary Methods S1. Hyperpolarized 13C images were acquired on a clinical 3 T MR system (MR750, GE Healthcare, Waukesha, WI, USA) using a clamshell coil for 13C-transmission. A 13C-tuned 8-channel array coil (Rapid Biomedical, Rimpar, Germany) was centered over the tumor. An Iterative Decomposition with Shifted Echo times and Least Squares Estimation (IDEAL) spiral chemical shift imaging (CSI) sequence [25] was employed to acquire five axial slices of 3 cm thickness with a 5 mm gap. The center of one slice was placed through the center of the tumor and as much of the tumor and contralateral normal kidney as possible were covered (Supplementary Figure S1). Radiofrequency pulses with a nominal flip angle (α) of 15° were applied to acquire a 34 cm × 34 cm field-of-view (FOV) with a matrix size of 40 × 40 voxels and a temporal resolution of 4 s for 20 time points (repetition time 0.5 s), starting 12 s after the injection of HP-pyruvate. Images had a true resolution of 17 mm × 17 mm and were reconstructed to 128 × 128 voxels. The patients were instructed to hold their breath for approximately the first 20 s and to breathe gently thereafter.

2.3. 13C-MRI Data Analysis

The IDEAL spiral data were reconstructed in MATLAB. Complex imaging data acquired over the eight channels of the abdominal array coil were combined using a singular value decomposition approach similar to Chen et al. [26]. For metabolite signals, the complex data were summed over all time points prior to coil combination to minimize noise propagation.
Polygonal regions of interest (ROIs) of the tumor and normal kidney were drawn manually on co-registered T1-weighted images and propagated to metabolic maps using OsiriX 10.0 (Pixemo SARL, Bernex, Switzerland). For individual biopsies, cylindrical ROIs with a diameter and height of 3 cm were centered on the biopsy. Voxel-wise intensity values were exported from OsiriX using the JavaScript Object Notation (JSON) format. Estimation of the mean and standard deviation (S.D.) of the noise was performed in an extracorporeal region of the image where spiral artifacts were absent. The signal-to-noise ratio (SNRmetabolite) was calculated from maps of pyruvate, lactate, and pyruvate+lactate (total carbon) within each ROI was defined as follows (see Formula (1)):
S N R m e t a b o l i t e = mean ( SI R O I ) mean ( SI n o i s e )   2   S . D .   ( SI n o i s e )
where mean (SIROI) is the mean signal intensity in the ROI, and the mean and S.D. of SInoise are computed over the ROI containing background only. The factor of √2 accounts for the narrowed Rayleigh distribution of magnitude noise, with an approximate adjustment for the use of multiple receivers.
The apparent reaction rate constant for the exchange of the HP-13C label between pyruvate and lactate (kPL) was calculated using a two-site exchange model in a frequency-domain approach and linear least-squares fitting. Back-conversion of lactate-to-pyruvate (kPL) and spin-lattice relaxation effects (T1) were combined as an effective relaxation term T1eff [27].
Masks containing only voxels with a summed carbon signal-to-noise ratio (SNR) of ≥5 were created to avoid voxels with a poor fit.

2.4. Proton (1H) MRI

The multiparametric 1H-MRI was acquired immediately after the HP-13C-MRI, following repositioning of the patient into a 32-channel cardiac array coil (GE Healthcare, Waukesha, WI, USA) on the same 3 T MRI scanner. The protocol included diffusion weighted, blood-oxygenation level dependent, and dynamic contrast-enhanced MRI in addition to T1- and T2-weighted sequences. The sequence parameters in Supplementary Table S1 were used for all except one participant where the presence of tumor thrombus in the intrahepatic inferior vena cava required a focused pre-surgical assessment (Supplementary Table S2). The sequences and image post-processing are described in the Supplementary Methods S1. All imaging data were analyzed blinded to the clinical outcome.

2.5. D-Printed Patient Specific Tumor Molds

Tumor molds were created as described previously [28]. In brief: the tumor, adjacent normal appearing kidney, and the perirenal fat including Gerota’s fascia were segmented manually in OsiriX (Pixmeo SARL, Switzerland). The segmentation was performed on the T1w LavaFlex sequence because this provided the highest spatial resolution. The segmentation was re-oriented using in-house MATLAB code and the mold generated using in-house code implemented in Python, interfacing with Meshlab and OpenSCAD. The 3D model was sliced in PrusaSlicer and a Prusa i3 MK3S was used to print the mold in polylactic acid (PLA).

2.6. Histology and Immunohistochemistry

Following nephrectomy, the specimen was sliced using a patient-specific 3D-printed tumor mold to allow co-registration of tissue samples and imaging [28]. The location of biopsies was recorded on maps rendered to represent the cross-section of the tissue at the location of the tissue slice and then transferred to the corresponding location in the T1w LavaFlex reference sequence (Supplementary Tables S1 and S2). Three to nine representative tumor samples (6 mm punches of fresh tissue) were obtained depending on tumor size, as well as samples from the ipsilateral normal appearing renal tissue. These were formalin-fixed and paraffin embedded. A subspecialized pathologist determined tumor grade on hematoxylin and eosin (H&E) stained slides according to the WHO/ISUP scale [6], assigning a grade to each tissue sample. Tissue samples with <75% viable tumor tissue were excluded.
Immunohistochemical (IHC) staining for MCT1 and MCT4 was performed using Leica’s Polymer Refine Detection System (DS9800) in combination with their Bond automated system (Leica Biosystems Newcastle Ltd., Newcastle upon Tyne, UK).
Sections were cut to 4 μm thickness and incubated with Tris EDTA for 20 min for antigen retrieval. Endogenous peroxidase activity was quenched using 3–4% (v/v) hydrogen peroxide. The primary antibodies were HPA003324 for MCT1 and HPA021451 for MCT4 (both Atlas Antibodies, Bromma, Sweden) both at a dilution of 1:500. The sections were then incubated in Anti-rabbit Poly-HRP-IgG polymer (<25 μg/mL) containing 10% (v/v) animal serum in Tris-buffered saline /0.09% ProClin 950 (Sigma-Aldrich, Gillingham, UK). The complex was visualized using 66 mM 3,3′-Diaminobenzidine tetrahydrochloride hydrate in a stabilizer solution and ≤0.1% (v/v) Hydrogen Peroxide. Leica DAB Enhancer was added to enhance staining. Cell nuclei were counterstained with <0.1% hematoxylin.
HALO v2.2.1870.15 (Indica Labs, Albuquerque, NM, USA) and the area quantification v2.1.11 module were used for automated analysis of scanned sections. Areas of weak, moderate, and strong staining were summed and divided by the total tissue area to obtain the percentage of positive tissue for MCT1 and MCT4 expression. Optical densities for weak, moderate, and strong stains used for the automated quantitative analysis of scanned sections were 0.2164, 0.3274 and 0.4938, respectively. The median percentage of positive tissue of all biopsies from one patient was taken to compensate for spatially heterogeneous marker expression within tumors. Similarly, cell density was quantified using HALO by dividing the number of cells counted in each biopsy by the tissue area. The median for all biopsies from one patient was used for patient-level comparisons.

2.7. TCGA-KIRC Data

Findings in the prospective dataset were validated using normalized RNAseq by Expectation Maximization (RSEM) data from the TCGA-KIRC dataset of ccRCC (accessed through Broad Institute, https://gdac.broadinstitute.org/, accessed on 24 January 2021). Gene expression data was z-transformed and the Kruskal–Wallis test used to test if the expression was dependent on tumor grade, complemented by pairwise Wilcoxon rank-sum tests with correction for multiple comparisons (Benjamini–Hochberg [29]). The proportion of variance in the expression of genes explained by another gene was measured by linear regression. Associations of gene expression and overall or disease-free survival (OS and DFS) were explored with Kaplan–Meier analysis at a pre-specified cut-off at the 85th percentile. Univariate and multivariate Cox regression models including clinical variables were fit.

2.8. Statistical Analysis

Statistical analyses were conducted using R (version 4.0.0, R foundation for statistical computing, Vienna, Austria). Mean and standard deviation were used as summary statistics for (approximatively) normally distributed continuous variables. Median and range were preferred with skewed data. Comparisons between unpaired samples used the Wilcoxon rank-sum test, paired samples the Wilcoxon Signed-Rank test and correlations the Spearman correlation for skewed data. Comparisons between approximately normally distributed data samples used Student’s t-test and Pearson’s correlation. Linear mixed models and two-way ordinal analysis of variance (ANOVA) were used to account for dependencies of biopsy-level data. The dplyr (version 1.0.2), ggplot2 (3.3.3), GGally (2.1.0) survival (3.2–7), survminer (0.4.8) packages were used for data processing, figure design and survival analysis.

3. Results

3.1. Participants

Nine participants with suspected renal masses were recruited as part of this physiological study (Supplementary Figure S2). Participant characteristics are reported in Table 1. All patients were imaged with HP-13C-MRI and multiparametric 1H-MRI shortly before tissue sampling (median of 12 days; interquartile range: 10–21). Diagnoses of the malignant tumors were confirmed on post-nephrectomy histology. Six patients were diagnosed with clear cell renal cell carcinoma (ccRCC); a dedifferentiated renal liposarcoma and a pheochromocytoma were imaged as they were thought to be RCCs based on pre-operative imaging. One participant, recruited for comparison purposes, harbored a benign oncocytoma which was diagnosed with multi-core needle biopsy. Research tissue was available for all malignant tumors except the pheochromocytoma where all samples were used to establish the clinical diagnosis.

3.2. Hyperpolarized 13C-MRI

Time-summed spectra of the slice containing the largest cross-section of the tumor and overlays of the [1-13C]pyruvate and [1-13C]lactate signals summed over time on T1w images for all nine patients are displayed in Figure 1. Quality control parameters of the 13C-pyruvate preparation are summarized in Supplementary Table S3. A median signal-to-noise ratio (SNR) for 13C-pyruvate of 26.7 (16.8–61.3) was measured in ccRCC (n = 6), similar to the contralateral normal kidney with a median SNR of 30.1 (14.8–64.7; n = 8). There was no significant difference in the median 13C-lactate SNR between ccRCC and the normal kidney: 5.7 (1.9–9.6) and 3.4 (2.0–9.3; p = 0.51), respectively. The dedifferentiated liposarcoma, pheochromocytoma, and oncocytoma showed a median pyruvate SNR of 31.9, 34.0, and 12.3, and a lactate SNR of 10.5, 6.0, and 1.0, respectively.
There was no significant difference in the median lactate-to-pyruvate ratio (Lac/Pyr) between ccRCC (0.13; 0.10–0.42) and normal kidney (0.14; 0.12–0.22; p = 0.40). The Lac/Pyr was 0.35 for the liposarcoma, 0.17 for the intrarenal pheochromocytoma, and 0.14 for the oncocytoma. The median kPL was 0.0065 (0.0024–0.0151) in ccRCC and 0.0043 (0.0028–0.0076) in the normal kidney. The liposarcoma was the most metabolically active tumor (kPL: 0.0152; pheochromocytoma: 0.0086). The oncocytoma presented with the lowest kPL (0.0022). The pyruvate and lactate SNR were not correlated, however the Lac/Pyr and kPL correlated strongly (r = 0.88, p = 0.03, Figure 2).
An increasing median kPL, in contrast to Lac/Pyr, correlated with a higher WHO/ISUP tumor grade (ρ = 0.92, p = 0.009; Figure 3a; Supplementary Figure S3a). The median kPL for individual biopsies (n = 44) was strongly patient-dependent (p < 0.001). A higher kPL in individual biopsies remained significantly predictive of higher tumor grade after accounting for dependence (p = 0.03, Figure 3b). Similarly, the Lac/Pyr was significantly increased in higher-grade tumor samples (p = 0.005, Supplementary Figure S3b).
Supplementary Figure S4 shows an example of 13C- and 1H-MRI obtained from one participant. Comparison of 1H-MRI parameters with WHO/ISUP tumor grade showed that perfusion fraction, a measure of capillary volume, trended towards a negative association with grade, but failed to reach statistical significance (Supplementary Figure S3d). No other 1H-MRI parameter correlated with tumor grade. The Lac/Pyr negatively correlated with the apparent diffusion coefficient D0 (p = 0.003, r = −0.98; Figure 2b, Supplementary Figure S5). The perfusion fraction showed a trend for correlation with the median kPL (p = 0.30, r = −0.58; Figure 2c).
Tumor volume showed no association with WHO/ISUP grade (Supplementary Figure S3c) and 13C-MRI.

3.3. Immunohistochemistry

MCT1 expression in ccRCC was positively correlated with the kPL (r = 0.89, p = 0.016; Figure 4), Lac/Pyr (r = 0.94, p = 0.005; Supplementary Figure S6) and grade. The kPL showed borderline significance in predicting the variance in tumor MCT1 expression at the biopsy level (p = 0.052; Supplementary Figure S8) and was independent of the cell density in tumor biopsies (p = 0.95, Supplementary Figure S7). MCT4 expression was not correlated with metabolic activity (Supplementary Figure S6).
The median MCT1 for all samples in each patient showed a trend towards decreased expression in ccRCC compared to normal renal tissue which was not significant (4.0 ± 6.0% vs. 6.2 ± 3.4%, p = 0.17, n = 6). Median MCT4 showed a higher positive pixel count in ccRCC compared to normal renal tissue (37 ± 10.5% vs. 10.2 ± 6.0%, mean ± S.D.; p = 0.009, n = 6).
The dedifferentiated liposarcoma showed a lower expression of MCT4 compared to ccRCC (16% vs. 37%) and an MCT1 expression comparable to WHO/ISUP grade 4 ccRCC (5.4% vs. 6.2%).
The concentration of serum LDH measured immediately before the HP-13C-MRI was independent of the tumor aggressiveness (Supplementary Figure S9). Measurements were not available for two participants owing to hemolyzed samples.

3.4. TCGA-KIRC RNA Expression

Analysis of RNA expression from 526 patients with ccRCC in the TCGA showed that MCT1 and MCT4 expression was grade dependent (p < 0.001 and 0.008 respectively) with overexpression in higher grade tumors (Figure 4 and Supplementary Figure S10). LDHA expression was independent of tumor grade (p = 0.06, Supplementary Figure S11). While the expression of all three genes was highly correlated (all p < 0.001), R2adj was low (MCT1/MCT4: 0.03, MCT1/LDHA: 0.16, MCT4/LDHA: 0.22; Supplementary Figure S12).
150 patients reached the overall survival endpoint while the surviving patients were censored after a median of 49 months (inter quartile range (IQR): 27–73 months). One-hundred and fifty patients reached the overall survival endpoint while the surviving patients were censored after a median of 49 months (inter quartile range (IQR): 27–73 months). A total of 127 patients progressed while the non-progressive patients were followed-up for a median of 45 months (IQR: 20–64 months). Univariate Cox and Kaplan–Meier analysis demonstrated that MCT1 and MCT4 expression was predictive of OS, while MCT1 and LDHA were predictive of progression-free survival (PFS; Figure 4, Supplementary Figures S10 and S11). In a multivariate Cox regression model, the presence of metastasis, age at diagnosis, and MCT1 expression, were independent predictors of OS (Table 2). In contrast, only metastasis, nodal stage and grade were predicted DFS (Supplementary Table S4).

4. Discussion

This proof of principle study has demonstrated the positive correlation between the apparent kinetic rate constant for the conversion of pyruvate to lactate (kPL) measured using HP-13C-MRI, the WHO/ISUP tumor grade and MCT1 expression in patients with ccRCC.
Previous work in RCC cell lines has shown increased 13C-lactate efflux in a more aggressive cell, which was attributed to increased MCT4 expression as the predominate mediator of lactate export and accompanied by a reduced MCT1 expression [30]. In vivo experiments showed an association between increasing HP pyruvate-to-lactate flux and higher LDHA expression [30,31,32]. Proteomic and metabolomic analyses of human RCCs confirm the relationship between lactate and tumor grade [12,13], and demonstrate the wide range of potential molecular drivers of HP-13C-lactate signal, raising the question of what may be the predominant factor in human ccRCC in vivo. A recent clinical study of HP-13C-pyruvate MRI in six patients with ccRCC (two additional patients had chromophobe RCCs), showed an increasing overall Lac/Pyr was associated with a higher tumor grade [21]. Here we have analyzed 44 multiregional biopsies from six ccRCC patients comparing biopsy grade with the co-registered metabolic imaging using a 3D tumor mold and compared immunohistochemical measures of pyruvate metabolism with imaging to explain the underlying molecular mechanism. The present study also assessed the complementary value of HP-13C-MRI and 1H-MRI and the value of the absolute quantification of the metabolic activity using kPL.
Multiparametric 1H-MRI was used to characterize the tumors morphologically and physiologically. The perfusion fraction was most strongly associated with tumor grade, but this association was weaker than for kPL, and did not reach statistical significance. The independence of perfusion fraction and kPL suggests that the metabolic changes detected between tumor grades are likely to reflect alterations in cellular metabolism rather than a perfusion effect. Although RCC volume was not related to metabolic activity in this study, larger and more hypoxic tumors have been shown to be more metabolically active in breast cancer [20].
The expression of MCT1 on IHC, the membrane transporter facilitating pyruvate influx [33], correlated with kPL and Lac/Pyr. The latter association has also been demonstrated in breast and prostate cancer [20,23]. The overexpression of MCT1 in grade 4 ccRCC in the TCGA-KIRC database strengthens this finding and, as MCT1 also predicts OS, it underlines its potential prognostic significance which is in agreement with existing literature [34,35]. Other biological factors may also control the formation of the labelled lactate, such as the expression of the enzyme LDH, the concentration of cofactor NADH/NAD+, as well as endogenous lactate levels [36].
The results from 1H-MRI supports the more glycolytic phenotype of aggressive tumors: the decreased perfusion fraction may indicate that aggressive tumors outgrow their blood supply and rely more on glycolysis. However, while the negative correlation of the perfusion fraction with tumor grade has been described previously, there is debate in the literature with other studies finding no association [37,38].
ccRCC and the contralateral normal kidney showed similar levels of HP-13C-pyruvate metabolism, which is supported by the similar MCT1 expression in RCC and normal kidney. MCT1 expression has been shown to be an important determinant of pyruvate-to-lactate conversion in preclinical models of pancreatic and breast cancer [39]. In this study, MCT4 was overexpressed in the tumor compared to normal renal tissue, in keeping with existing data [40]. MCT4 exports lactate into the extracellular space, which is central to maintaining a high glycolytic activity and promoting tumor cell and invasion [41,42]. Owing to the high metabolic activity of the normal kidney, HP-13C-MRI is unlikely to be used for detection of renal tumors determine tumor aggressiveness, as these tumors will usually have been detected using conventional imaging techniques as part of routine workflow; major clinical unmet needs include improved non-invasive tumor and treatment stratification, particularly if a biopsy cannot be obtained easily. Therefore, imaging tumor metabolism using HP-13C-MRI, will be complementary to the structural and functional information provided by imaging modalities such as CT and 1H-MRI.
Alterations in tumor metabolism, and cellular energy metabolism in particular, are a hallmark of cancer and have been exploited as biomarkers of treatment response [43,44]. In vivo studies of HP-13C-pyruvate MRI in many cancer models, including ccRCC, have shown treatment-induced metabolic changes as early as 24 h after the start of treatment [45]. Recent evidence from a clinical study has shown that early treatment response to neoadjuvant therapy can be assessed with HP-13C-pyruvate MRI after 7–11 days, which is associated with pathological complete response in patients with breast cancer [36]. Detecting successful treatment response earlier will play an increasingly important role in ccRCC with the advent of an increasing number of treatment options in advanced disease.
Future clinical research will assess the clinical utility of HP-13C-pyruvate MRI for these key questions of tumor grading and treatment response detection in ccRCC. For example, metabolic imaging could be used to prioritize small but metabolically active renal lesions for surgery rather than surveillance.
Although the number of participants in this study is relatively small, the use of multiregional biopsies has enhanced the power to detect significance between the imaging metrics and the histological and immunohistochemical measurements. The multiregional analyzes employing a 3D-printed tumor mold for accurate co-registration of imaging and tissue samples, confirming the association between kPL, Lac/Pyr, and tumor grade, as well as kPL and MCT1 expression. Two of the patients whose pre-surgical diagnosis suggested the presence of RCC were diagnosed with an intrarenal pheochromocytoma and a dedifferentiated liposarcoma, presenting the first opportunity to investigate the metabolism of these tumors with HP-13C-MRI. RCC exhibits marked intratumoral heterogeneity [46] and therefore tissue samples and molecular biomarkers may be subject to sampling error. The non-invasive identification of this intratumoral metabolic heterogeneity represents one of the strengths of this emerging technique and could provide important new measures to assist tumor stratification in the future.

5. Conclusions

The apparent kinetic rate constant for the conversion of pyruvate-to-lactate derived from HP-13C-MRI (kPL) differentiated overall ccRCC tumor grade in each tumor, as well as the individual tumor grade acquired from multiregional biopsies in each tumor, in this proof of principle study. This kPL was correlated with the expression of the MCT1 pyruvate transporter, which is an independent predictor of survival in patients with ccRCC. Therefore, metabolic imaging may have important prognostic implications for stratifying ccRCCs, could be used to identify tumor metabolic heterogeneity, and could be used to monitor therapy response in the future.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cancers14020335/s1, Methods S1: Hyperpolarized [1-13C]pyruvate sample preparation, proton MRI technique, and proton image processing, Figure S1: Positioning of the slices for the acquisition of the HP-13C-pyruvate MRI. The blue slice indicates the largest cross-section of the tumor which is also shown in Figure 1, Figures S4 and S5, Figure S2: Patient flow diagram, Figure S3: Correlation between WHO/ISUP tumor grade and imaging parameters, Figure S4: Representative example images of a patient with grade IV ccRCC (Patient ccRCC 4), Figure S5: Diffusivity (D0) maps corresponding to the metabolic images presented in Figure 1 of the manuscript, Figure S6: Immunohistochemistry for the monocarboxylate transporters MCT1 and MCT4 from biopsy samples, Figure S7: Correlation between kPL and histological cell density, Figure S8: Correlation between kPL and MCT1 expression at the single biopsy level, Figure S9: Association between serum LDH levels and tumor aggressiveness, Figure S10: Representative micrographs of a clear cell renal cell carcinoma with high and low expression of MCT4, Figure S11: Association between LDHA expression, survival and tumor aggressiveness in the TCGA-KIRC dataset, Figure S12: Correlation between gene expression parameters in the TCGA-KIRC dataset. Gene expressions were z-score transformed, Table S1: Proton MR parameters for all patients except ccRCC 1, Table S2: Proton MR parameters for ccRCC 1, Table S3: 13C pyruvate preparation QC parameters, Table S4: Multivariate logistic regression for progression-free survival in the TCGA KIRC cohort.

Author Contributions

Conceptualization: S.U., M.A.M., T.E., E.S., K.M.B., G.D.S., F.A.G.; methodology: S.U., R.W., M.A.M., A.N.P., M.C.-O., C.R.B., A.B.G., M.G., J.T.G., S.S.D., F.R., J.D.K., F.Z., J.A.G.D., M.J.L., R.C., I.A.M., S.J.W., M.J.G., T.E., T.J.M., A.Y.W., K.M.B., E.S., G.D.S., F.A.G.; software S.U., M.C.-O., M.A.M., A.N.P., C.R.B.; formal analysis, S.U., R.W., M.A.M., A.N.P., C.R.B., A.B.G., A.Y.W., F.A.G.; investigation, S.U., R.W., L.B., A.C.P.R., J.F.-R., A.F., T.F.A., J.N.A., R.C., I.A.M., T.B., T.J.M., A.Y.W., G.D.S., F.A.G.; resources, F.A.G.; data curation, S.U., M.A.M.; writing—original draft preparation, S.U.; writing—review and editing, all authors; visualization, S.U.; supervision, E.S., G.D.S., F.A.G.; project administration, M.J.L., A.F., J.F.-R.; funding acquisition, F.A.G.; Edited and approved the manuscript: all co-authors. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Cancer Research UK (CRUK; C19212/A16628, C19212/A29082, C19212/A27150), Wellcome Trust (095962), The Mark Foundation for Cancer Research and the CRUK Cambridge Centre (C9685/A25177), Addenbrooke’s Charitable Trust, the Cambridge National Institute for Health Research Biomedical Research Centre (NIHR BRC), Cambridge University Hospitals NHS Foundation Trust and the Austrian Science Fund (J4025-B26). The Urological Malignancies Programme is part of the CRUK Cambridge Centre, funded by a CRUK Major Centre Award, and supported by the NIHR Cambridge BRC. Infrastructure for the Cambridge Urological Biorepository received funding from the NIHR Cambridge BRC and CRUK Cambridge Centre. The Human Research Tissue Bank is supported by the NIHR Cambridge BRC. AYW is supported by the NIHR Cambridge Biomedical Research Centre. JTG was funded by a Medical Research Council award. This project received research support from GE Healthcare.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Research Ethics Committee of Cambridge South (UK) 15/EE/0378 on 14/12/2015.

Informed Consent Statement

Written informed consent was obtained from all subjects involved in the study. Written informed consent has been obtained from the patients to publish this paper.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to ethical restrictions.

Acknowledgments

The results published here are in part based upon data generated by the TCGA Research Network: https://www.cancer.gov/tcga, accessed on 24 January 2021. We are grateful to Dominique-Laurent Couturier, senior biostatistician at the Cancer Research UK Cambridge Centre, for his statistical advice. We are grateful to Vincent Gnanapragasam, chief investigator of the DIAMOND study (REC 03/018), for access to tissue samples. We thank the tissue doners and research groups involved in the generation of these data for their contribution. S.J.W. is funded by a grant from the M.R.C. We acknowledge the help of Gaspar Delso (GE Healthcare) and Dattesh Shanbhag (GE Global Research) for the use and ongoing support of their magnetic resonance imaging motion-correction programming code. The views expressed are those of the authors and not necessarily those of the NIHR or the Department of Health and Social Care. Open Access Funding by the Austrian Science Fund (FWF).

Conflicts of Interest

The authors declare ongoing research collaboration with GE Healthcare. G.D.S. has received educational grants from Pfizer, AstraZeneca and Intuitive Surgical; consultancy fees from Pfizer, Merck, EUSA Pharma and CMR Surgical; Travel expenses from Pfizer and Speaker fees from Pfizer. T.E. has received research support from AstraZeneca, Bayer and Pfizer; he was employed by AstraZeneca until March 2020 and is now employed by Roche; he has stock in AstraZeneca and Roche. S.J.W. has received travel funding from Ipsen. E.S. is a co-founder and shareholder of Lucida Medical, receives consulting fees from Amazon and speaker’s fees from G.S.K. and Siemens, none of these relate to the work presented here. F.A.G. has grants from G.S.K, research support from GE Healthcare and has consulted for AstraZeneca. J.T.G. holds a grant from Bristol Myers Squibb, receives research support from GE Healthcare, PulseTeq, and Polarean, and receives consultancy fees from the Journal of Magnetic Resonance Imaging. None of these relate to the work presented here.

References

  1. Ferlay, J.; Ervik, M.; Lam, F.; Colombet, M.; Mery, L.; Piñeros, M.; Znaor, A.; Soerjomataram, I.; Bray, F. Global Cancer Observatory: Cancer Today. Available online: https://gco.iarc.fr/today/home (accessed on 24 January 2021).
  2. Cheville, J.C.; Lohse, C.M.; Zincke, H.; Weaver, A.L.; Blute, M.L. Comparisons of Outcome and Prognostic Features Among Histologic Subtypes of Renal Cell Carcinoma. Am. J. Surg. Pathol. 2003, 27, 612–624. [Google Scholar] [CrossRef] [PubMed]
  3. Ljungberg, B.; Albiges, L.; Abu-Ghanem, Y.; Bensalah, K.; Dabestani, S.; Montes, S.F.P.; Giles, R.H.; Hofmann, F.; Hora, M.; Kuczyk, M.A.; et al. European Association of Urology Guidelines on Renal Cell Carcinoma: The 2019 Update. Eur. Urol. 2019, 75, 799–810. [Google Scholar] [CrossRef] [PubMed]
  4. Gumundsson, E.; Hellborg, H.; Lundstam, S.; Erikson, S.; Ljungberg, B. Metastatic potential in renal cell carcinomas <7 cm: Swedish kidney cancer quality register data. Eur. Urol. 2011, 60, 975–982. [Google Scholar] [CrossRef]
  5. Marconi, L.; Dabestani, S.; Lam, T.B.; Hofmann, F.; Stewart, F.; Norrie, J.; Bex, A.; Bensalah, K.; Canfield, S.E.; Hora, M.; et al. Systematic review and meta-analysis of diagnostic accuracy of percutaneous renal tumour biopsy. Eur. Urol. 2016, 69, 660–673. [Google Scholar] [CrossRef]
  6. Delahunt, B.; Cheville, J.C.; Martignoni, G.; Humphrey, P.A.; Magi-Galluzzi, C.; McKenney, J.; Egevad, L.; Algaba, F.; Moch, H.; Grignon, D.J.; et al. The International Society of Urological Pathology (ISUP) grading system for renal cell carcinoma and other prognostic parameters. Am. J. Surg. Pathol. 2013, 37, 1490–1504. [Google Scholar] [CrossRef] [Green Version]
  7. Dagher, J.; Delahunt, B.; Rioux-Leclercq, N.; Egevad, L.; Srigley, J.R.; Coughlin, G.; Dunglinson, N.; Gianduzzo, T.; Kua, B.; Malone, G.; et al. Clear cell renal cell carcinoma: Validation of World Health Organization/International Society of Urological Pathology grading. Histopathology 2017, 71, 918–925. [Google Scholar] [CrossRef]
  8. Stanzione, A.; Ricciardi, C.; Cuocolo, R.; Romeo, V.; Petrone, J.; Sarnataro, M.; Mainenti, P.P.; Improta, G.; De Rosa, F.; Insabato, L.; et al. MRI Radiomics for the Prediction of Fuhrman Grade in Clear Cell Renal Cell Carcinoma: A Machine Learning Exploratory Study. J. Digit. Imaging 2020, 33, 879–887. [Google Scholar] [CrossRef]
  9. Ursprung, S.; Beer, L.; Bruining, A.; Woitek, R.; Stewart, G.D.; Gallagher, F.A.; Sala, E. Radiomics of computed tomography and magnetic resonance imaging in renal cell carcinoma—a systematic review and meta-analysis. Eur. Radiol. 2020, 30, 1–9. [Google Scholar] [CrossRef] [Green Version]
  10. Linehan, W.M.; Srinivasan, R.; Schmidt, L.S. The genetic basis of kidney cancer: A metabolic disease. Nat. Rev. Urol. 2010, 7, 277–285. [Google Scholar] [CrossRef] [PubMed]
  11. Nickerson, M.L.; Jaeger, E.; Shi, Y.; Durocher, J.A.; Mahurkar, S.; Zaridze, D.; Matveev, V.; Janout, V.; Kollarova, H.; Bencko, V.; et al. Improved Identification of von Hippel-Lindau Gene Alterations in Clear Cell Renal Tumors. Clin. Cancer Res. 2008, 14, 4726–4734. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Wettersten, H.I.; Hakimi, A.A.; Morin, D.; Bianchi, C.; Johnstone, M.E.; Donohoe, D.R.; Trott, J.F.; Abu Aboud, O.; Stirdivant, S.; Neri, B.; et al. Grade-dependent metabolic reprogramming in kidney cancer revealed by combined proteomics and metabolomics analysis. Cancer Res. 2015, 75, 2541–2552. [Google Scholar] [CrossRef] [Green Version]
  13. Jing, L.; Guigonis, J.M.; Borchiellini, D.; Durand, M.; Pourcher, T.; Ambrosetti, D. LC-MS based metabolomic profiling for renal cell carcinoma histologic subtypes. Sci. Rep. 2019, 9, 1–10. [Google Scholar] [CrossRef]
  14. Witney, T.H.; Kettunen, M.I.; Day, S.E.; Hu, D.; Neves, A.A.; Gallagher, F.A.; Fulton, S.M.; Brindle, K.M. A comparison between radiolabeled fluorodeoxyglucose uptake and hyperpolarized 13C-labeled pyruvate utilization as methods for detecting tumor response to treatment. Neoplasia 2009, 11, 574–582. [Google Scholar] [CrossRef] [Green Version]
  15. Zaccagna, F.; Grist, J.T.; Deen, S.S.; Woitek, R.; Lechermann, L.M.; McLean, M.A.; Basu, B.; Gallagher, F.A. Hyperpolarized carbon-13 magnetic resonance spectroscopic imaging: A clinical tool for studying tumour metabolism. Br. J. Radiol. 2018, 91, 20170688. [Google Scholar] [CrossRef]
  16. Woitek, R.; McLean, M.A.; Gill, A.B.; Grist, J.T.; Provenzano, E.; Patterson, A.J.; Ursprung, S.; Torheim, T.; Zaccagna, F.; Locke, M.; et al. Hyperpolarized 13C MRI of Tumor Metabolism Demonstrates Early Metabolic Response to Neoadjuvant Chemotherapy in Breast Cancer. Radiol. Imaging Cancer 2020, 2, e200017. [Google Scholar] [CrossRef] [PubMed]
  17. Albers, M.J.; Bok, R.; Chen, A.P.; Cunningham, C.H.; Zierhut, M.L.; Zhang, V.Y.; Kohler, S.J.; Tropp, J.; Hurd, R.E.; Yen, Y.F.; et al. Hyperpolarized 13C lactate, pyruvate, and alanine: Noninvasive biomarkers for prostate cancer detection and grading. Cancer Res. 2008, 68, 8607–8615. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Grist, J.T.; McLean, M.A.; Riemer, F.; Schulte, R.F.; Deen, S.S.; Zaccagna, F.; Woitek, R.; Daniels, C.J.; Kaggie, J.D.; Matyz, T.; et al. Quantifying normal human brain metabolism using hyperpolarized [1–13C]pyruvate and magnetic resonance imaging. Neuroimage 2019, 189, 171–179. [Google Scholar] [CrossRef] [PubMed]
  19. Cunningham, C.H.; Lau, J.Y.C.; Chen, A.P.; Geraghty, B.J.; Perks, W.J.; Roifman, I.; Wright, G.A.; Connelly, K.A. Hyperpolarized 13 C Metabolic MRI of the Human Heart. Circ. Res. 2016, 119, 1177–1182. [Google Scholar] [CrossRef] [PubMed]
  20. Gallagher, F.A.; Woitek, R.; McLean, M.A.; Gill, A.B.; Garcia, R.M.; Provenzano, E.; Riemer, F.; Kaggie, J.; Chhabra, A.; Ursprung, S.; et al. Imaging breast cancer using hyperpolarized carbon-13 MRI. Proc. Natl. Acad. Sci. USA 2020, 117, 2092–2098. [Google Scholar] [CrossRef] [Green Version]
  21. Tang, S.; Meng, M.V.; Slater, J.B.; Gordon, J.W.; Vigneron, D.B.; Stohr, B.A.; Larson, P.E.Z.; Wang, Z.J. Metabolic imaging with hyperpolarized 13C pyruvate magnetic resonance imaging in patients with renal tumors—Initial experience. Cancer 2021, 127, 2693–2704. [Google Scholar] [CrossRef] [PubMed]
  22. Felmlee, M.A.; Jones, R.S.; Rodriguez-Cruz, V.; Follman, K.E.; Morris, M.E. Monocarboxylate transporters (SLC16): Function, regulation, and role in health and disease. Pharmacol. Rev. 2020, 72, 466–485. [Google Scholar] [CrossRef] [PubMed]
  23. Granlund, K.L.; Tee, S.S.; Vargas, H.A.; Lyashchenko, S.K.; Reznik, E.; Fine, S.; Laudone, V.; Eastham, J.A.; Touijer, K.A.; Reuter, V.E.; et al. Hyperpolarized MRI of Human Prostate Cancer Reveals Increased Lactate with Tumor Grade Driven by Monocarboxylate Transporter 1. Cell Metab. 2020, 31, 105–114. [Google Scholar] [CrossRef] [PubMed]
  24. Creighton, C.J.; Morgan, M.; Gunaratne, P.H.; Wheeler, D.A.; Gibbs, R.A.; Robertson, G.; Chu, A.; Beroukhim, R.; Cibulskis, K.; Signoretti, S.; et al. Comprehensive molecular characterization of clear cell renal cell carcinoma. Nature 2013, 499, 43–49. [Google Scholar] [CrossRef] [Green Version]
  25. Wiesinger, F.; Weidl, E.; Menzel, M.I.; Janich, M.A.; Khegai, O.; Glaser, S.J.; Haase, A.; Schwaiger, M.; Schulte, R.F. IDEAL spiral CSI for dynamic metabolic MR imaging of hyperpolarized [1-13C]pyruvate. Magn. Reson. Med. 2012, 68, 8–16. [Google Scholar] [CrossRef] [PubMed]
  26. Chen, H.Y.; Autry, A.W.; Brender, J.R.; Kishimoto, S.; Krishna, M.C.; Vareth, M.; Bok, R.A.; Reed, G.D.; Carvajal, L.; Gordon, J.W.; et al. Tensor image enhancement and optimal multichannel receiver combination analyses for human hyperpolarized 13C MRSI. Magn. Reson. Med. 2020, 84, 3351–3365. [Google Scholar] [CrossRef]
  27. Khegai, O.; Schulte, R.F.; Janich, M.A.; Menzel, M.I.; Farrell, E.; Otto, A.M.; Ardenkjaer-Larsen, J.H.; Glaser, S.J.; Haase, A.; Schwaiger, M.; et al. Apparent rate constant mapping using hyperpolarized [1-13C]pyruvate. NMR Biomed. 2014, 27, 1256–1265. [Google Scholar] [CrossRef] [PubMed]
  28. Crispin-Ortuzar, M.; Gehrung, M.; Ursprung, S.; Gill, A.B.; Warren, A.Y.; Beer, L.; Gallagher, F.A.; Mitchell, T.J.; Mendichovszky, I.A.; Priest, A.N.; et al. Three-Dimensional Printed Molds for Image-Guided Surgical Biopsies: An Open Source Computational Platform. JCO Clin. Cancer Inform. 2020, 4, 736–748. [Google Scholar] [CrossRef]
  29. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B 1995, 57, 289–300. [Google Scholar] [CrossRef]
  30. Keshari, K.R.; Sriram, R.; Koelsch, B.L.; Van Criekinge, M.; Wilson, D.M.; Kurhanewicz, J.; Wang, Z.J. Hyperpolarized 13C-pyruvate magnetic resonance reveals rapid lactate export in metastatic renal cell carcinomas. Cancer Res. 2013, 73, 529–538. [Google Scholar] [CrossRef] [Green Version]
  31. Sriram, R.; Van Criekinge, M.; Santos, J.D.; Keshari, K.R.; Wilson, D.M.; Peehl, D.; Kurhanewicz, J.; Wang, Z.J. Non-Invasive Differentiation of Benign Renal Tumors from Clear Cell Renal Cell Carcinomas Using Clinically Translatable Hyperpolarized C Pyruvate Magnetic Resonance. Tomography 2016, 2, 35–42. [Google Scholar] [CrossRef]
  32. Sriram, R.; Gordon, J.; Baligand, C.; Ahamed, F.; Santos, J.D.; Qin, H.; Bok, R.A.; Vigneron, D.B.; Kurhanewicz, J.; Larson, P.E.Z.; et al. Non-Invasive Assessment of Lactate Production and Compartmentalization in Renal Cell Carcinomas Using Hyperpolarized 13 C Pyruvate MRI. Cancers 2018, 10, 313. [Google Scholar] [CrossRef] [Green Version]
  33. Parks, S.K.; Chiche, J.; Pouysségur, J. Disrupting proton dynamics and energy metabolism for cancer therapy. Nat. Rev. Cancer 2013, 13, 611–623. [Google Scholar] [CrossRef]
  34. Kim, Y.S.; Choi, J.W.; Lee, J.H.; Kim, Y.S. Expression of lactate/H+ symporters MCT1 and MCT4 and their chaperone CD147 predicts tumor progression in clear cell renal cell carcinoma: Immunohistochemical and the Cancer Genome Atlas data analyses. Hum. Pathol. 2015, 46, 104–112. [Google Scholar] [CrossRef] [PubMed]
  35. Ambrosetti, D.; Dufies, M.; Dadone, B.; Durand, M.; Borchiellini, D.; Amiel, J.; Pouyssegur, J.; Rioux-Leclercq, N.; Pages, G.; Burel-Vandenbos, F.; et al. The two glycolytic markers GLUT1 and MCT1 correlate with tumor grade and survival in clear-cell renal cell carcinoma. PLoS ONE 2018, 13, e0193477. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Woitek, R.; McLean, M.A.; Ursprung, S.; Rueda, O.M.; Manzano Garcia, R.; Locke, M.J.; Beer, L.; Baxter, G.; Rundo, L.; Provenzano, E.; et al. Hyperpolarized Carbon-13 MRI for Early Response Assessment of Neoadjuvant Chemotherapy in Breast Cancer Patients. Cancer Res. 2021, 81, 6004–6017. [Google Scholar] [CrossRef]
  37. Yu, X.; Qin, Y.; Liu, Y.; Li, C.; Xie, Y.; Sun, Y. Value of intravoxel incoherent motion diffusion-weighted imaging in differentiating the pathological grade of clear cell renal cell carcinoma. Int. J. Clin. Exp. Med. 2018, 11, 10923–10928. [Google Scholar]
  38. Ye, J.; Xu, Q.; Wang, S.A.; Zheng, J.; Dou, W.Q. Quantitative Evaluation of Intravoxel Incoherent Motion and Diffusion Kurtosis Imaging in Assessment of Pathological Grade of Clear Cell Renal Cell Carcinoma. Acad. Radiol. 2020, 27, e176–e182. [Google Scholar] [CrossRef]
  39. Rao, Y.; Gammon, S.; Zacharias, N.M.; Liu, T.; Salzillo, T.; Xi, Y.; Wang, J.; Bhattacharya, P.; Piwnica-Worms, D. Hyperpolarized [1-13C]pyruvate-to-[1-13C]lactate conversion is rate-limited by monocarboxylate transporter-1 in the plasma membrane. Proc. Natl. Acad. Sci. USA 2020, 117, 22378–22389. [Google Scholar] [CrossRef]
  40. Fisel, P.; Kruck, S.; Winter, S.; Bedke, J.; Hennenlotter, J.; Nies, A.T.; Scharpf, M.; Fend, F.; Stenzl, A.; Schwab, M.; et al. DNA methylation of the SLC16A3 promoter regulates expression of the human lactate transporter MCT4 in renal Cancer with consequences for clinical outcome. Clin. Cancer Res. 2013, 19, 5170–5181. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Dimmer, K.S.; Friedrich, B.; Lang, F.; Deitmer, J.W.; Broer, S. The low-affinity monocarboxylate transporter MCT4 is adapted to the export of lactate in highly glycolytic cells. Biochem. J. 2000, 350, 219–227. [Google Scholar] [CrossRef]
  42. Gerlinger, M.; Santos, C.R.; Spencer-Dene, B.; Martinez, P.; Endesfelder, D.; Burrell, R.A.; Vetter, M.; Jiang, M.; Saunders, R.E.; Kelly, G.; et al. Genome-wide RNA interference analysis of renal carcinoma survival regulators identifies MCT4 as a Warburg effect metabolic target. J. Pathol. 2012, 227, 146–156. [Google Scholar] [CrossRef] [Green Version]
  43. Hanahan, D.; Weinberg, R.A. Hallmarks of cancer: The next generation. Cell 2011, 144, 646–674. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Johnson, S.; Kumar, A.; Matasar, M.; Schöder, H.; Rademaker, J. Imaging for Staging and Response Assessment in Lymphoma. Radiology 2015, 276, 323–338. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Dong, Y.; Eskandari, R.; Ray, C.; Granlund, K.L.; Dos Santos-Cunha, L.; Miloushev, V.Z.; Tee, S.S.; Jeong, S.; Aras, O.; Chen, Y.B.; et al. Hyperpolarized MRI visualizes Warburg effects and predicts treatment response to mTOR inhibitors in patient-derived CCRCC xenograft models. Cancer Res. 2019, 79, 242–250. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Gerlinger, M.; Rowan, A.J.; Horswell, S.; Larkin, J.; Endesfelder, D.; Gronroos, E.; Martinez, P.; Matthews, N.; Stewart, A.; Tarpey, P.; et al. Intratumor Heterogeneity and Branched Evolution Revealed by Multiregion Sequencing. N. Engl. J. Med. 2012, 366, 883–892. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Summary of all tumors included in this study. Axial T1w images as a reference. Non-localized 13C spectra of the axial slice containing the largest tumor cross-section summed over all timepoints. 13C-Pyruvate and 13C-lactate signal summed over all time points superimposed on an axial T1w image of the largest tumor cross-section. The border of the tumor is outlined in blue. A: Aorta, Lac: Lactate, Pyr: Pyruvate, PyrH: Pyruvate hydrate. A PASS score (Thompson, The American Journal of Surgical Pathology: May 2002) of ≥4 is associated with a potential for clinically aggressive behavior in pheochromocytoma. 1H-MRI.
Figure 1. Summary of all tumors included in this study. Axial T1w images as a reference. Non-localized 13C spectra of the axial slice containing the largest tumor cross-section summed over all timepoints. 13C-Pyruvate and 13C-lactate signal summed over all time points superimposed on an axial T1w image of the largest tumor cross-section. The border of the tumor is outlined in blue. A: Aorta, Lac: Lactate, Pyr: Pyruvate, PyrH: Pyruvate hydrate. A PASS score (Thompson, The American Journal of Surgical Pathology: May 2002) of ≥4 is associated with a potential for clinically aggressive behavior in pheochromocytoma. 1H-MRI.
Cancers 14 00335 g001
Figure 2. Relationships between hyperpolarized 13C-pyruvate MRI parameters as well as between HP-13C-MRI and proton MRI parameters. (a) The median lactate-to-pyruvate ratio was strongly correlated to the median kPL. The volume of the tumors, represented by the area of the points on the plot, ranged from 53 to 1350 cm3 (53–908 cm3 for ccRCC). There was no correlation between ccRCC volume and lactate-to pyruvate ratio and kPL. (b) The median lactate-to-pyruvate ratio was negatively correlated with the median diffusivity. Correlation coefficients calculated for ccRCC only. (c) Association of the median perfusion fraction from IVIM-type DWI with the median kPL from 13C-MRI.
Figure 2. Relationships between hyperpolarized 13C-pyruvate MRI parameters as well as between HP-13C-MRI and proton MRI parameters. (a) The median lactate-to-pyruvate ratio was strongly correlated to the median kPL. The volume of the tumors, represented by the area of the points on the plot, ranged from 53 to 1350 cm3 (53–908 cm3 for ccRCC). There was no correlation between ccRCC volume and lactate-to pyruvate ratio and kPL. (b) The median lactate-to-pyruvate ratio was negatively correlated with the median diffusivity. Correlation coefficients calculated for ccRCC only. (c) Association of the median perfusion fraction from IVIM-type DWI with the median kPL from 13C-MRI.
Cancers 14 00335 g002
Figure 3. Comparison of the ISUP/WHO Grade with median kPL. (a) Increasing tumor grade was associated with increased metabolic activity in ccRCC. The liposarcoma and pheochromocytoma showed metabolic activity comparable to grade IV ccRCC. The benign renal oncocytoma showed the lowest metabolic activity. (b) An increase in kPL values was associated with an increasing tumor grade determined for the individual tumor biopsy. LS: Liposarcoma, Onc: Oncocytoma, Pheo: Pheochromocytoma, RCC: clear cell renal cell carcinoma.
Figure 3. Comparison of the ISUP/WHO Grade with median kPL. (a) Increasing tumor grade was associated with increased metabolic activity in ccRCC. The liposarcoma and pheochromocytoma showed metabolic activity comparable to grade IV ccRCC. The benign renal oncocytoma showed the lowest metabolic activity. (b) An increase in kPL values was associated with an increasing tumor grade determined for the individual tumor biopsy. LS: Liposarcoma, Onc: Oncocytoma, Pheo: Pheochromocytoma, RCC: clear cell renal cell carcinoma.
Cancers 14 00335 g003
Figure 4. Representative micrographs of a clear cell renal cell carcinoma with high and low expression of MCT1. (a) The left panels show a hematoxylin and eosin stain, and the right panels the corresponding immunohistochemical stain for MCT1 on an adjacent tissue slice. (b) Correlation between MCT1 IHC expression and kPL in the patient dataset imaged as part of this study (MISSION, dots represent ccRCC, the triangle the liposarcomas) (c) Box plot comparing z-transformed MCT1 expression as a function of histological tumor grade from a Cancer Genome Atlas dataset (TCGA-KIRC, yellow). (d) Kaplan–Meier plots for the association of MCT1 expression with overall and (e) progression/recurrence-free survival using an expression cut-off at the 85th percentile in the TCGA-IRC dataset. MCT1 expression was significantly associated with overall and progression/recurrence-free survival.
Figure 4. Representative micrographs of a clear cell renal cell carcinoma with high and low expression of MCT1. (a) The left panels show a hematoxylin and eosin stain, and the right panels the corresponding immunohistochemical stain for MCT1 on an adjacent tissue slice. (b) Correlation between MCT1 IHC expression and kPL in the patient dataset imaged as part of this study (MISSION, dots represent ccRCC, the triangle the liposarcomas) (c) Box plot comparing z-transformed MCT1 expression as a function of histological tumor grade from a Cancer Genome Atlas dataset (TCGA-KIRC, yellow). (d) Kaplan–Meier plots for the association of MCT1 expression with overall and (e) progression/recurrence-free survival using an expression cut-off at the 85th percentile in the TCGA-IRC dataset. MCT1 expression was significantly associated with overall and progression/recurrence-free survival.
Cancers 14 00335 g004
Table 1. Patient characteristics.
Table 1. Patient characteristics.
Patient CharacteristicsDistribution
Patients (male/female)9 (8/1)
Patient Age (median ± IQR) (years)59.5 ± 8.7
Histology6 clear cell renal cell carcinoma
1 Pheochromocytoma
1 Dedifferentiated Liposarcoma
1 Renal Oncocytoma
Tumor Stage at Surgery (RCC only)2 pT1b pNx cM0
3 pT3a pN0 cM0
1 pT3b pNx cM1
WHO/ISUP Tumor Grade at Surgery (RCC only)1 Grade 2
2 Grade 3
3 Grade 4
Location of Metastasis1 Lung
Patient Weight (median ± IQR) (kg)90.1 ± 13.5
Time between imaging and surgery (median ± IQR) (days)12 ± 11
Laterality5 left/4 right
Plasma glucose (median ± IQR) (mmol/l)5.0 ± 0.3
Patients (male/female)9 (8/1)
S.D.: standard deviation, IQR: Inter-quartile range.
Table 2. Multivariate logistic regression for overall survival in the TCGA-IRC dataset.
Table 2. Multivariate logistic regression for overall survival in the TCGA-IRC dataset.
Covariatesp ValueHR95% Confidence Interval
LowerUpper
Age0.0021.0341.0121.058
Female Sex0.170.6960.4161.162
LN0.301.5790.6653.748
Metastasis<0.0013.1791.8345.510
Size0.361.1910.8181.734
Grade
Grade 10.990.0010
Grade 20.150.5910.2871.218
Grade 30.500.7990.4141.540
MCT10.0101.3091.0651.609
MCT40.551.0730.8501.355
LDHA0.0820.7970.6171.029
HR: Hazard ratio, LDHA: Lactate dehydrogenase A, LN: Lymph node, MCT1/4: Monocarboxylate transporter ¼. Concordance index = 0.783.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ursprung, S.; Woitek, R.; McLean, M.A.; Priest, A.N.; Crispin-Ortuzar, M.; Brodie, C.R.; Gill, A.B.; Gehrung, M.; Beer, L.; Riddick, A.C.P.; et al. Hyperpolarized 13C-Pyruvate Metabolism as a Surrogate for Tumor Grade and Poor Outcome in Renal Cell Carcinoma—A Proof of Principle Study. Cancers 2022, 14, 335. https://doi.org/10.3390/cancers14020335

AMA Style

Ursprung S, Woitek R, McLean MA, Priest AN, Crispin-Ortuzar M, Brodie CR, Gill AB, Gehrung M, Beer L, Riddick ACP, et al. Hyperpolarized 13C-Pyruvate Metabolism as a Surrogate for Tumor Grade and Poor Outcome in Renal Cell Carcinoma—A Proof of Principle Study. Cancers. 2022; 14(2):335. https://doi.org/10.3390/cancers14020335

Chicago/Turabian Style

Ursprung, Stephan, Ramona Woitek, Mary A. McLean, Andrew N. Priest, Mireia Crispin-Ortuzar, Cara R. Brodie, Andrew B. Gill, Marcel Gehrung, Lucian Beer, Antony C. P. Riddick, and et al. 2022. "Hyperpolarized 13C-Pyruvate Metabolism as a Surrogate for Tumor Grade and Poor Outcome in Renal Cell Carcinoma—A Proof of Principle Study" Cancers 14, no. 2: 335. https://doi.org/10.3390/cancers14020335

APA Style

Ursprung, S., Woitek, R., McLean, M. A., Priest, A. N., Crispin-Ortuzar, M., Brodie, C. R., Gill, A. B., Gehrung, M., Beer, L., Riddick, A. C. P., Field-Rayner, J., Grist, J. T., Deen, S. S., Riemer, F., Kaggie, J. D., Zaccagna, F., Duarte, J. A. G., Locke, M. J., Frary, A., ... Gallagher, F. A. (2022). Hyperpolarized 13C-Pyruvate Metabolism as a Surrogate for Tumor Grade and Poor Outcome in Renal Cell Carcinoma—A Proof of Principle Study. Cancers, 14(2), 335. https://doi.org/10.3390/cancers14020335

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop