Functional noninvasive detection of glycolytic pancreatic ductal adenocarcinoma
Cancer & Metabolism volume 10, Article number: 24 (2022)
Pancreatic ductal adenocarcinoma (PDAC) lacks effective treatment options beyond chemotherapy. Although molecular subtypes such as classical and QM (quasi-mesenchymal)/basal-like with transcriptome-based distinct signatures have been identified, deduced therapeutic strategies and targets remain elusive. Gene expression data show enrichment of glycolytic genes in the more aggressive and therapy-resistant QM subtype. However, whether the glycolytic transcripts are translated into functional glycolysis that could further be explored for metabolic targeting in QM subtype is still not known.
We used different patient-derived PDAC model systems (conventional and primary patient-derived cells, patient-derived xenografts (PDX), and patient samples) and performed transcriptional and functional metabolic analysis. These included RNAseq and Illumina HT12 bead array, in vitro Seahorse metabolic flux assays and metabolic drug targeting, and in vivo hyperpolarized [1-13C]pyruvate and [1-13C]lactate magnetic resonance spectroscopy (HP-MRS) in PDAC xenografts.
We found that glycolytic metabolic dependencies are not unambiguously functionally exposed in all QM PDACs. Metabolic analysis demonstrated functional metabolic heterogeneity in patient-derived primary cells and less so in conventional cell lines independent of molecular subtype. Importantly, we observed that the glycolytic product lactate is actively imported into the PDAC cells and used in mitochondrial oxidation in both classical and QM PDAC cells, although more actively in the QM cell lines. By using HP-MRS, we were able to noninvasively identify highly glycolytic PDAC xenografts by detecting the last glycolytic enzymatic step and prominent intra-tumoral [1-13C]pyruvate and [1-13C]lactate interconversion in vivo.
Our study adds functional metabolic phenotyping to transcriptome-based analysis and proposes a functional approach to identify highly glycolytic PDACs as candidates for antimetabolic therapeutic avenues.
Despite enormous research efforts, pancreatic ductal adenocarcinoma (PDAC) remains a fatal disease with marginal clinical advancement . Although genomic and transcriptional profiles of PDAC have been studied in great detail [2,3,4], effective targeting strategies remain scarce. Sequencing efforts in large patient cohorts have identified distinct molecular PDAC subtypes in several independent studies with two dominant subgroups. Those are termed classical (term used hereafter) or pancreatic progenitor with more epithelial differentiated tumor and quasi-mesenchymal (QM; term used hereafter) or squamous or basal like [1, 5,6,7] with more mesenchymal tumor respectively. Especially, QM PDACs demonstrate very aggressive phenotypes with shorter median survival and resistance to first-line chemotherapy with FOLFIRINOX . Which cancer cell features contribute to the aggressive and therapy-resistant QM phenotype remains unknown.
Metabolic rewiring, i.e., an individual cell’s ability to use different metabolic pathways depending on alternating growth conditions including oxygen and nutrient availability, has been implicated as a major cause of therapy resistance in cancers and aggravates clinically successful targeting . This allows cells not only to adapt but also to thrive on particularly scarce conditions of hypoxia and nutrient limitations typically observed in PDAC . Glycolysis is the most prominent cancer-associated metabolic pathway. Although high cancer dependency on glucose was described nearly 100 years ago by Otto Warburg , glycolytic targeting is still not widely therapeutically exploited. Recently, expression of glycolytic metabolic transcripts has been associated with the resistant QM PDAC subtype in patients . Work in PDAC mouse models demonstrated that glycolysis is the major metabolic effector of oncogenic KRAS, the leading PDAC driver, and that co-targeting of RAS-RAF-MEK-MAPK cascade and glycolysis may be an effective approach in PDAC . However, functional evidence that glycolysis is indeed significantly operable in human QM PDACs is missing.
Here, we addressed this missing link and analyzed functional exposure of glycolysis in different clinically relevant PDAC samples ranging from long-term cultured PDAC cell lines to patient-derived xenografts and primary cells and bulk PDAC probes. We found considerable heterogeneity in the glycolytic behavior especially among patient-derived PDAC cells. However, individual representatives of the QM subtype were indeed functionally highly glycolytic what was preserved even in the in vivo xenograft setting. By using a noninvasive hyperpolarized 13C-magnetic resonance spectroscopy (HP-MRS), we were able to detect the final glycolytic step in vivo, namely intratumoral conversion of HP-[1-13C]pyruvate to HP-[1-13C]lactate. Importantly, QM PDAC cells actively consumed the final glycolytic product, lactate, in mitochondrial oxidative phosphorylation in vitro, what was further in vivo translated and detected as HP-[1-13C]lactate to HP-[1-13C]pyruvate conversion in QM PDAC xenografts. This suggested that glycolytic QM cells not only actively produce lactate but also metabolically use it. Our work opens a perspective for noninvasive detection of glycolytic PDACs and monitoring of individualized anti-glycolytic targeting approaches.
PDAC cell lines
All PDAC cell lines have been obtained from the ATCC and regularly externally authenticated (at least once a year). PDAC cell lines (PSN1 (RRID: CVCL_1644); Kp4 (RRID: CVCL_1338), PaTu8988T (CVCL_1847); MiaPaca2 (RRID: CVCL_0428), PaTu8988S (RRID: CVCL_1846), HPAC (RRID: CVCL_3517), HPAFII (RRID: CVCL_0313), and HupT4 (RRID: CVCL_1300) were grown in Dulbecco’s Modified Eagle Medium (1:1 mix of DMEM no. 11966025 and DMEM no. A1443001, Thermo Fisher Scientific, Waltham, USA) adapted to final concentrations of 5 mM D-glucose (Thermo Fisher Scientific, Waltham, USA), 2 mM L glutamine, 5% v/v fetal bovine serum (FBS, Thermo Fisher Scientific, Waltham, USA), and 1% v/v penicillin/streptomycin (P/S, Thermo Fisher Scientific, Waltham, USA) if not stated otherwise.
Patient-derived cells (PDCs)
From 11 PDX samples, we were able to isolate and cultivate cancer cells (PDCs) for further analysis. For all metabolic analysis, PDC cell lines were cultivated in a 1:1 mixture of Keratinocyte-SF medium (no. 17005075, Thermo Fisher Scientific, Waltham, USA) and RPMI 1640 (no. 11879020, Thermo Fisher Scientific, Waltham, USA) adapted to final concentrations of 5 mM D-glucose, 4.5 mM L-glutamine, 0.26 mM sodium pyruvate, 6%v/v FBS, and 1% v/v final mixture of penicillin/streptomycin (Thermo Fisher Scientific, Waltham, USA), and antimycotic/antibiotic (cat. no. 15240–062, Thermo Fisher Scientific) if not stated otherwise.
RNA isolation and gene expression analysis
Established/PDC cells were cultivated for 48 h in the respective media. At confluence of 70–80%, cells were placed on ice and washed twice with ice-cold PBS, mechanically scratched from the plate in 1 ml of ice-cold PBS, and centrifuged at 4 °C/400 g for 5 min. Pelleted cells were stored in − 80 °C until all cells were collected for RNA isolation. RNA was isolated using the Maxwell RSC simplyRNA Cells Kit (no. AS1390, Promega, Germany). Cell RNA isolation kit was used according to the manufacturer’s instructions. Total RNA was stored at − 80 °C until further processing and gene expression analysis. For PDX samples, RNA was isolated from fresh-frozen PDX tumor tissue using the PARIS (Ambion) isolation kit.
PDX samples preparation
Establishment of the PDX mouse model was performed using surgically resected PDAC tissues collected from patients.
Seahorse metabolic flux assays
All assays were performed following the manufacturer’s instructions (Agilent Technologies).
Immunohistochemistry and immunofluorescence
Immunohistochemistry was performed according to standard laboratory procedures on PFA fixed, FFPE tissue samples. Antibodies used in this study are as follows: MCT4, Atlas Antibodies (Sigma-Aldrich, Cat no. HPA021451, RRID:AB_1853663); HIF1a, BD Transduction laboratories no. 610959 (RRID: AB_398272); MCT1, Abcam, no. ab85021(RRID: AB_10674945); KRT81, Santa Cruz, no. sc-100929 (RRID: AB_2132772); and pancytokeratin, Abcam no. ab6401(RRID: AB_305450).
Hyperpolarized magnetic resonance spectroscopy (HP-MRS)
PSN1/HPAC cells were implanted subcutaneously (s.c.) into the back of male or female 6-week-old Crl:NIH-Foxn1rnu rats (Charles River). Pyruvate-lactate metabolism was measured with multi-frame slice spectroscopy (MRS, 15 mm slice thickness) using alternating metabolite-frequency-selective excitation (flip angle 30°, 250 Hz transmit bandwidth, 2 kHz receiver bandwidth, both metabolites separately excited and measured every 2 s) while injecting hyperpolarized HP-[1-13C]pyruvate or HP-[1-13C]lactate. Procedure optimization is described in detail in our previous study . For data analysis after HP-[1-13C]pyruvate injection, lactate and pyruvate spectral peak heights were summed over all time points and presented as the ratio of these areas under curves (AUClac/AUCpyr) . For data analysis after HP-[1-13C]lactate injection, due to low pyruvate signal, the signal intensities of lactate and pyruvate spectra were averaged over 10 time points near the maximum intensity and then fit with a constant offset plus a Lorentzian function with fixed 30 Hz full width at half maximum to determine the peak area (PA) , using the least-squares curve fit function in MatLab.
Glycolytic metabolic gene transcripts are present in QM PDAC subtype
To evaluate whether glycolytic transcripts are omnipresent in the QM PDAC subtype , we first performed gene expression analysis in multiple preclinical and clinical samples (Fig. 1a). RNA-seq or Illumina HT12 gene expression analysis was performed for conventional PDAC cell lines (n = 8), patient-derived xenografts (PDX, n = 34), and PDX-derived cancer cells (PDC, n = 11). Transcriptomes from bulk tissue of 204 PDAC samples from previously published resource were utilized (E-MTAB-1791). All samples were then subtyped to QM or classical group. For tumor subtype determination, nonnegative matrix factorization (NMF)  was used, after median centering of the data. ConsensusClusterPlus  was used to validate class assignment from NMF. In a first benchmarking step, we used publicly available transcriptionally subtyped PDAC cohorts (PDAC cell lines (GSE21654 ), PDAC xenograft (E-MTAB-4029 ), and bulk PDAC tissue (GSE16515 and GSE15471 [20, 21]) to verify the robustness of our classification pipeline. With this approach, we were able to reproduce more than 90% of the reported subtypes in our test datasets.
After establishing the pipeline, the same parameters were used for subtyping of patient PDAC samples, PDX cohort, and PDC samples to QM, and classical group and gene set enrichment analysis (GSEA) comparing the QM and classical groups were performed. By using our subtyping platform, in the PDAC patient sample cohort (204 samples), 88 are classified as classical and 116 as QM subtype. Samples clustering to the QM subtype presented significant enrichment of selected QM and squamous subtype assigner gene sets previously described [5, 6] (Fig. 1b and supplementary Table 1) supporting correct subtype assignment.
In the PDX cohort, 22 classical and 12 QM tumors were identified. Among PDCs, 6 classical (PDC44, 58, 59, 62, 70, 89) and 5 QM (PDC 34, 57, 69, 78, 80) were identified. The 8 PDAC cell lines used in this study were previously classified as QM (KP4, PSN1, MIAPaca2, PaTu8988T) and classical (PaTu8988S, HUPT4, HPAFII, HPAC) . We analyzed gene expression of vimentin (VIM) and E-cadherin (CDH1) as markers of mesenchymal and epithelial status respectively. As expected, general trend towards higher VIM expression in QM and CDH1 expression in classical PDAC cells was observed (supplementary Figure. 1a).
After classification, QM and classical groups were compared by GSEA for HALLMARK, REACTOME, and KEGG collections in all datasets. A full list of all enriched gene sets with respective normalized enrichment score (NES) and false discovery rate (FDR) values is given in supplementary Table 2. As expected for the mesenchymal phenotype, enrichment of the epithelial-to-mesenchymal transition (EMT) gene set was observed in the QM group in PDX, PDC, and patient PDAC samples (Fig. 1c) supporting correct assignment of the subtypes. In the QM samples, transcripts for glycolysis, hypoxia, and MYC-target genes were well preserved throughout different sample collections (Fig. 1c). The hypoxia gene set was enriched in QM bulk PDAC tissue, PDX, and PDC data sets, even though PDC cells were cultured under common laboratory normoxic conditions. Concordantly with the well-described correlation of hypoxia and glycolysis , glycolysis/glucose metabolism transcripts were also enriched in the QM patient PDAC samples, PDX and PDAC cell line datasets, and MYC target gene sets as well. Interestingly, in the QM PDCs, the glycolysis gene set was not enriched, possibly due to low sample numbers but also suggesting no unambiguous assignment of glycolytic genes to the QM subtype at least in PDCs. In summary, we observed strong transcriptional association of QM subtype with glycolysis in different preclinical and clinical samples.
Glycolytic pathway activity is high in individual QM PDAC cells
To investigate whether glycolysis is indeed functionally active in QM PDAC cells, we performed Seahorse metabolic flux assays and evaluated the extracellular acidification rate (ECAR) and oxygen consumption rate (OCR) as readouts of two major energy-supplying processes, glycolysis and oxidative phosphorylation respectively. ECAR and OCR levels were measured in cell lines and PDCs in media with physiological concentrations of 5 mM glucose with addition of 2 mM glutamine. Under these conditions, PSN1 and PDC69, both QM, presented the highest ECAR/OCR ratios and glycolytic energy phenotype among cell lines and PDCs, respectively (Fig. 2a and supplementary Table 3). However, energetic phenotype of other cell lines and PDCs was rather heterogeneous and independent of their transcriptional subtype, being QM or classical.
Notably, qPCR analysis revealed that the relative gene expression of the last glycolytic enzyme lactate dehydrogenase A (LDHA), lactate exporter MCT4 (SLC16A3), and importer MCT1 (SLC16A1) is high in PSN1 and PDC69 cells in comparison with other cells. The same was true for HIF1a, a central cellular regulator of hypoxia and glycolysis (Fig. 2b). Moreover, hierarchical clustering of transcriptome data showed generally higher expression of several glycolytic genes (e.g., HK1, HK2, ENO1, ENO2, PGK1) in QM cell lines and PDCs, especially in PSN1 and PDC69 (Fig. 2c). Taken together, active glycolysis was observed in some QM PDAC cells that correlated well with the high expression of glycolytic genes but was not unambiguously connected to QM subtype.
Lactate exporter MCT4 has previously been suggested to be a good marker of glycolytic PDACs . In publicly available TCGA RNA expression datasets, high expression of both lactate transporters MCT4 and MCT1 correlated with worse survival in PDAC patients associated with the resistant QM subtype (Fig. 3a). Immunohistochemical analysis of MCT1 and MCT4 in FFPE samples of 30 human PDACs suggested that both MCT4 and MCT1 were expressed on cancer and stromal cells with MCT4 being more prominently expressed on cancer cells, while MCT1 was often prominently expressed in the surrounding stroma as well (Fig. 3b). In our cohorts, we observed that MCT4 (SLC16A3) gene expression levels were higher than MCT1 (SLC16A1) in both bulk PDAC and PDX tissue samples (Fig. 3c and supplementary Figure. 2a), suggesting a lead role of MCT4 as lactate transporter in tissue context. Furthermore, multiplex immunofluorescence for pancytokeratin (PanCK), MCT4 and KRT81, an established QM marker , in 6 PDAC FFPE specimens, showed that the proportion of MCT4-positive cells was higher among KRT81 positive (30–50%) than KRT81-negative cancer cells (< 20%) (Fig. 3d). Our data support the use of MCT4 as a surrogate marker of QM PDACs with activated glycolysis.
PDAC cells actively use lactate as oxidative fuel
Intrigued by high expression of lactate transporters detected in some of the QM PDAC cells, we aimed to investigate lactate metabolism in PDAC. It is now well accepted that lactate is not only the end waste product of glycolysis but is also actively used in metabolic processes in cancer as well. Lactate conversion to pyruvate and subsequent oxidation in the mitochondria has been suggested in murine PDAC . However, whether this effect is especially attributable to lactate-producing high-glycolytic QM PDAC cells is still not known. To investigate this, PDAC cells were cultivated for 7 h in (i) “basal” DMEM or RPMI media without glucose or glutamine supplementation or in (ii) “basal” media supplemented with lactate (basal + 10 mM L-lactate). Consequently, Seahorse metabolic flux measurement was performed, and OCR values measured in media with and without lactate were compared. In basal media, classical PDAC cell lines (HPAC, HPAFII, and HupT4) presented generally higher basal OCR levels than the QM cell lines, indicating that oxidative phosphorylation (OXPHOS) is well supplied by alternative fuels other than glucose or glutamine in the classical cell lines (supplementary Figure 2b). Interestingly, lactate supplementation to the medium led to OCR boost in all cells of both QM, and classical subtypes with however more pronounced OCR increase in the QM PDAC cell lines (Fig. 4a). In PDCs, lactate treatment led to an OCR increase in all cells, without pronounced subtype-specific effect (Fig. 4a and supplementary Figure 2b).
To substantiate this finding, we cultivated PSN1 (QM), PaTu8988T (QM), and PaTu8988S (classical) cells in DMEM medium with 5 mM glucose and 2 mM glutamine without media change for 24–48 to 72–96 h. Glucose and lactate concentrations in the media were measured at given time points. With time, glucose concentration in the media decreased, and lactate increased (0–72 h), as expected due to glucose consumption and lactate production and accumulation. Once the glucose was consumed from the medium (approx. after 72 h in PaTu8988T/PSN1 cells), lactate concentration in the media decreased, indicating that in the absence of other resources, PDAC cells start consuming self-produced lactate (supplementary Figure 2c).
We also challenged the detected changes in glycolysis and lactate metabolism with the inhibitor of lactate dehydrogenase GNE-140 and followed the concentration-dependent inhibition of metabolic activity in cells via CellTiter-Glo assay (Fig. 4b). GNE-140 treatment indeed induced a decrease in cell viability especially in the QM cell lines, being most effective in PSN1, MIAPaca2, and PaTu8988T cells. PDCs were in general less sensitive to GNE-140, and the observed inhibitory effects were, as expected from Seahorse lactate supplementations assays, not subtype dependent.
In conclusion, PDAC cells, regardless of subtype, not only actively produce and excrete glycolytically produced lactate but also actively use it potentially as an oxidative fuel. This phenomenon is more pronounced in QM than in classical PDAC cell lines and is exposed to specific metabolic targeting with lactate dehydrogenase inhibitors.
Hyperpolarized magnetic resonance spectroscopy of [1-13C]pyruvate and [1-13C] lactate identifies highly glycolytic tumors
Pharmacological inhibition suggested efficacy of GNE-140 in glycolytic cells arguing for the need of unequivocal identification of highly glycolytic PDACs for successful metabolic targeting. However, detection of dominant metabolic pathways driving tumor phenotypes remains a highly challenging task. Fluorodeoxyglucose (18F-FDG) positron emission tomography (PET) is clinically established method for tumor detection based on high glucose uptake into the cancer. However, 18F-FDG-PET detects only the very first step of glycolysis since 18F-FDG gets phosphorylated by the first glycolytic enzyme hexokinase or glucokinase and does not enter the further metabolic processing. Thus, FDG-PET detects the glucose trapping in the cell, rather than the real glycolytic activity of the tumor. We sought to evaluate the last step of glycolysis, pyruvate to lactate conversion, and explored HP-MRS with [1-13C]pyruvate and [1-13C]lactate in PDAC in an in vivo approach. For this purpose, rats were subcutaneously implanted with glycolytic QM PSN1 or classical HPAC PDAC cell lines. Consistent with the respective molecular subtype, PSN1 tumors presented an undifferentiated mesenchymal histology, while HPAC tumors showed a more differentiated epithelial tumor (supplementary Figure 3a). Once the tumors reached a minimal size of 5 × 5 mm, metabolic spectroscopy was performed. HP-[1-13C]pyruvate was i.v. injected into the tail vein, and intra-tumoral accumulation of HP-[1-13C]lactate was followed in real time. Using MRS, significantly more HP-[1-13C]lactate was detected in PSN1 compared to HPAC tumors, supporting higher 13C label exchange between pyruvate and lactate specifically in PSN1 tumors (Fig. 5 a and b, supplementary Figure 3b).
To evaluate whether lactate is imported and used by tumors in vivo as observed in vitro in Seahorse experiments, we also performed the reverse experiment and intravenously injected HP-[1-13C]lactate in PSN1 and HPAC tumor rats in vivo. High HP-[1-13C]lactate uptake and intratumoral HP-[1-13C]pyruvate were clearly detected in PSN1 compared to very low HP-[1-13C]pyruvate signal in HPAC xenografts (Fig. 5c, supplementary Figure 3b). Accordingly, significantly higher peak area (PA) PApyr/PAlac ratios were measured for PSN1 than HPAC tumors (Fig. 5d).
Taken together, highly glycolytic PSN1 xenografts could readily be discriminated based on high HP-[1-13C]pyruvate to HP-[1-13C]lactate interconversion observed in HP-MRS.
Lactate dehydrogenase (LDH) enzymatic activity measured ex vivo after the spectroscopy experiment in snap frozen tissues was also higher in PSN1 compared to HPAC tumors (Fig. 5e) consistent with the in vivo finding. We further confirmed the highly glycolytic nature of rat PSN1 tumors by immunohistochemical analysis of glycolytic markers HIF1A and MCT4 in rat PSN1 and HPAC xenografts used in MRS experiments. MCT4 showed the typical membrane-associated expression in cancer cells in both xenografts, with only slightly stronger staining intensity in PSN1 than in HPAC tumors (Fig. 5f) but not prominently different. Intriguingly, staining for HIF1A, a major glycolytic regulator in cancers, was found exclusively in the PSN1 tumors with typical nuclear expression pattern in the cancer cells (Fig. 5f). Next, we also analyzed HIF1A and MCT4 expression in murine xenografts of human PDAC cell lines (supplementary Figure 4). Indeed, MCT4 staining intensity was also only lightly stronger in the murine QM xenografts compared to xenografts of classical cell lines. Specific nuclear HIF1A expression was limited to QM tumors only (PSN1, KP4, MIAPaCa2, PaTu8988T), and not detected in classical tumors (HPAFII, PaTu8988S, HUPT4, HPAC) (supplementary Fig. 4).
The challenge in PDAC is its enormous therapy resistance due to the evolution of aggressive cancer cells driven by oncogenic KRAS and loss of key tumor suppressors in a complex adapting microenvironment with various signaling effectors and biophysical and hypoxic restraints. Despite considerable genetic homogeneity with regard to oncogenic KRAS as lead driver, many studies support the existence of several molecular PDAC subtypes, including classical/progenitor and QM/squamous/basal-like and hybrid states with more or less pronounced subtype-specific transcriptional programs [3,4,5]. Though indisputably present, functional aspects and phenotypic cues of the defined transcriptional subtypes are less described. One key feature of PDAC is the metabolic rewiring that may lead to phenotypic features not entirely captured by transcriptomic signatures. PDACs identified as transcriptionally glycolytic show amplification of KRAS and MYC genes and are associated with a worse prognosis both in resectable and metastatic setting . In patients, targeting of KRAS-MEK-MAPK pathway in a monotherapy approach only is not successful potentially due to activation of escape routes such as PI3K-AKT. However, blocking glycolysis with 2-deoxyglucose in combination with MAPK inhibitor is at least in mice effective and leads to apoptosis induction and reduction in tumor volume, suggesting high potential of this co-targeting . Functional identification of patients with highly glycolytic PDACs can lead clinical decision-making and introduction of anti-glycolytic drugs in the clinic. In this work, we pursued two aims: (i) to evaluate whether the presence of glycolytic transcripts is indeed translated into operable glycolysis in PDAC QM subtype and (ii) offer a noninvasive imaging-based approach for detection of highly glycolytic tumors. We focused this analysis on patient-derived model systems including PDX and PDCs to value the molecular and metabolic heterogeneity in primary PDAC model systems. Furthermore, all our metabolic assays are performed under supplementation with physiological levels of glucose (5 mM), thus omitting the metabolic artifacts that can be caused by the usage of typical high glucose media.
Gene expression analysis in four different model systems (cell lines, PDC, PDX, and bulk tissue samples) indeed identified glucose metabolism/glycolysis/hypoxia/MYC targets as dominating metabolic transcripts of the QM subtype. This is in line with the previously observed “glycolytic” subtype in mesenchymal PDAC cell lines  and the recently reported “glycolytic” transcriptional PDAC subtypes in patients . However, glycolysis was not unambiguously functionally dominant in all cells of QM subtype, being cell lines or primary. In functional assays, we observed notable heterogeneity in metabolic behavior especially in patient-derived cells. We found active functional glycolysis in single representatives of the QM subtype, such as PSN1 and PDC69 cells. In Seahorse assays, these cells demonstrated high ECAR to OCR ratios, suggesting that cell intrinsic energy metabolism relies rather on glycolysis than on OXPHOS. It should however be noted that Seahorse assays evaluate ECAR and OCR values in in vitro conditions and are highly dependent on cell culture features such as current cellular density, growth pattern, cell cycle, and current mitochondrial number  and should be interpreted only as indication of the cellular energetic status. Glycolytic energetic status of PSN1 and PDC69 correlated well with high gene expression of the lactate producer and transporters LDHA and MCT1/4, respectively, supporting the translation of transcripts in active glucose metabolism. Interestingly, HIF1A, a major transcriptional regulator of glycolysis and cellular response to hypoxia , was also well expressed in the identified glycolytic cells here grown in typical in vitro normoxic conditions, supporting intrinsic gene expression programs well preserved in QM cells. Furthermore, by using multiplex immunofluorescence approach, we observed MCT4 expression in KRT81-positive cells in human PDAC samples, further suggesting correlation of QM subtype and glycolytic phenotype. In line with our observations, MCT4 has already been suggested as marker of glycolytic PDACs with poor prognosis .
The heterogeneity observed in our results suggests that rigid transcriptome-based classification of PDAC subtypes may not be sufficient as the basis for clinical decisions regarding metabolic targeting approaches. Rather, individual PDACs may often present a continuum of different metabolic states that are more or less phenotypically presented depending on various cell-autonomous and non-cell-autonomous cues. Hybrid PDAC subtypes with transcriptomic signatures in between the classical and QM/basal-like states have been highlighted recently [4, 11]. Similar to our study, a correlation of molecular cues and functional oxidative phosphorylation was very recently reported for PDAC cells . The authors emphasize on metabolic heterogeneity and flexibility and shifts from OXPHOS or glycolysis when necessary, supporting the existence of plastic metabolic states dependent on the environmental challenges. It is reasonable to assume that among PDAC cells, a whole spectrum of weakly to highly glycolytic QM PDAC cells exists. The exclusive dependency on the one or the other metabolic pathway is an unlikely scenario. However, individual tumors with high activity of specific metabolic pathway may exist, and their identification will be the key to successful targeting. We show here that both PSN1 and HPAC xenografted tumors import HP-[1-13C]pyruvate, however show different conversion rates to HP-[1-13C]lactate. Glycolytic PSN1 tumors were readily detectable by HP-MRS due to higher 13C-label exchange among pyruvate and lactate, indicating high activity of the last glycolytic enzyme LDH and high intratumoral pyruvate to lactate conversion. Similarly, in breast cancer patients, high HP-[1-13C]pyruvate to HP-[1-13C]lactate conversion rates identified strongly glycolytic aggressive triple negative breast cancer with high HIF1a and MCT1 tissue expression  and high-grade lesions in prostate cancer with increased MCT4 expression . This approach is already being used in personalized therapy monitoring in prostate and brain cancer [30, 31]. Additionally, alanine-to-lactate signal ratio upon injection of HP-[1-13C]pyruvate has also been shown to distinguish well between preneoplastic lesions and PDAC in mouse models [32, 33]. In our study, we focused on potentials of imaging pyruvate-to-lactate interconversions for detection of a subpopulation of highly glycolytic PDACs.
Furthermore, we also confirmed in vivo that HP-[1-13C]lactate enters the PSN1 and HPAC tumors, with however lower uptake in classical HPAC xenografts. As a result, we clearly detected HP-[1-13C]pyruvate in PSN1 tumors and very low levels in HPAC tumors, what suggests potential subsequent use of pyruvate in the TCA cycle and high metabolic flexibility of PSN1 tumors for pyruvate and lactate as oxidative fuels. Lactate has recently been considered as one of the important actors in tumor metabolism . Tumors use the advantage of lactate being the second most abundant metabolite in the systemic circulation and readily feed the TCA cycle with pyruvate generated from lactate in a reverse lactate dehydrogenase reaction [25, 35, 36]. Indeed, in Seahorse experiments, we also observed OXPHOS activation with lactate in PDAC cells, especially in the QM cell lines. However, it should be noted that Seahorse experiments were performed under deprivation of main OXPHOS fuels, glucose and glutamine. This may potentially lead to overemphasis of lactate usage in OXPHOS under these in vitro conditions.
We hypothesize that the hypoxic microenvironment of the tumor favors the epithelial-to-mesenchymal transformation (EMT) of the cancer cells and appearance of the glycolytic QM tumors. These tumors potentially adapt their oxidative metabolism to fuels which are then locally produced, either by themselves or by neighboring cancer, stromal, or immune cells. In some cases, this will likely be through oxidation of lactate, although lactate oxidation by itself may not be a sufficient marker for the classical to QM transition.
Our HP-MRS experiments provide evidence for the concept that PDACs with high reliance on glycolysis are potentially detectable via HP-[1-13C]pyruvate/lactate MRS imaging (MRSI) in clinical practice. Thus, identification of highly glycolytic, aggressive PDACs by HP-[1-13C]pyruvate and HP-[1-13C]lactate MRSI may be used to guide and monitor tumor treatment with anti-glycolytic therapies. The limitation for clinical translation of this method is the low pyruvate SNR measured after the HP-lactate injection, leading to an uncertainty on the performed quantifications. Therefore, future improved methods of higher sensitivity of the HP-[1-13C]lactate experiment such as employment of cryogenically cooled receiver coils would be of high interest and may even allow to quantify metabolic fluxes within the different subtypes of PDAC.
In contrast to biopsy-based tumor characterization, metabolic imaging allows dynamic evaluation of the whole tumor limiting sampling bias and addressing tumor heterogeneity . Although likely not all QM tumors are potentially extremely glycolytic, noninvasive detection of highly glycolytic PDACs by HP-[1-13C]pyruvate/lactate MRS is one of the first methods for successful individual metabolic approaches.
Availability of data and materials
All data generated in this paper are available from the corresponding author upon a reasonable request.
Pancreatic ductal adenocarcinoma
Hyperpolarized magnetic resonance spectroscopy (imaging)
Formalin-fixed paraffin embedded
Patient derived cells
Monocarboxylate transporter 4
Monocarboxylate transporter 1
Normalized enrichment score
False discovery rate
Oxygen consumption rate
Extracellular acidification rate
The Cancer Genome Atlas
Dulbecco’s Modified Eagle Medium
Roswell Park Memorial Institute
Positron emission tomography
Aung KL, Fischer SE, Denroche RE, Jang GH, Dodd A, Creighton S, et al. Genomics-driven precision medicine for advanced pancreatic cancer: early results from the COMPASS trial. Clin Cancer Res. 2018;24(6):1344–54.
Waddell N, Pajic M, Patch AM, Chang DK, Kassahn KS, Bailey P, et al. Whole genomes redefine the mutational landscape of pancreatic cancer. Nature. 2015;518(7540):495–501.
Moffitt RA, Marayati R, Flate EL, Volmar KE, Loeza SG, Hoadley KA, et al. Virtual microdissection identifies distinct tumor- and stroma-specific subtypes of pancreatic ductal adenocarcinoma. Nat Genet. 2015;47(10):1168–78.
Chan-Seng-Yue M, Kim JC, Wilson GW, Ng K, Figueroa EF, O’Kane GM, et al. Transcription phenotypes of pancreatic cancer are driven by genomic events during tumor evolution. Nat Genet. 2020;52(2):231–40.
Collisson EA, Sadanandam A, Olson P, Gibb WJ, Truitt M, Gu S, et al. Subtypes of pancreatic ductal adenocarcinoma and their differing responses to therapy. Nat Med. 2011;17(4):500–3.
Bailey P, Chang DK, Nones K, Johns AL, Patch AM, Gingras MC, et al. Genomic analyses identify molecular subtypes of pancreatic cancer. Nature. 2016;531(7592):47–52.
Cancer Genome Atlas Research Network. Electronic address aadhe, Cancer Genome Atlas Research N. Integrated genomic characterization of pancreatic ductal adenocarcinoma. Cancer Cell. 2017;32(2):185–203 e13.
DeBerardinis RJ, Chandel NS. Fundamentals of cancer metabolism. Sci Adv. 2016;2(5): e1600200.
Biancur DE, Kimmelman AC. The plasticity of pancreatic cancer metabolism in tumor progression and therapeutic resistance. Biochim Biophys Acta Rev Cancer. 2018;1870(1):67–75.
Koppenol WH, Bounds PL, Dang CV. Otto Warburg’s contributions to current concepts of cancer metabolism. Nat Rev Cancer. 2011;11(5):325–37.
Karasinska JM, Topham JT, Kalloger SE, Jang GH, Denroche RE, Culibrk L, et al. Altered gene expression along the glycolysis-cholesterol synthesis axis is associated with outcome in pancreatic cancer. Clin Cancer Res. 2020;26(1):135–46.
Yan L, Tu B, Yao J, Gong J, Carugo A, Bristow CA, et al. Targeting glucose metabolism sensitizes pancreatic cancer to MEK inhibition. Cancer Res. 2021;81(15):4054–65.
Topping GJ, Heid I, Trajkovic-Arsic M, Kritzner L, Grashei M, Hundshammer C, et al. Hyperpolarized (13)C spectroscopy with simple slice-and-frequency-selective excitation. Biomedicines. 2021;9(2):121. https://www.mdpi.com/2227-9059/9/2/121.
Hill DK, Orton MR, Mariotti E, Boult JK, Panek R, Jafar M, et al. Model free approach to kinetic analysis of real-time hyperpolarized 13C magnetic resonance spectroscopy data. PLoS One. 2013;8(9): e71996.
Hyacinthe JN, Buscemi L, Le TP, Lepore M, Hirt L, Mishkovsky M. Evaluating the potential of hyperpolarised [1-(13)C] L-lactate as a neuroprotectant metabolic biosensor for stroke. Sci Rep. 2020;10(1):5507.
Gaujoux R, Seoighe C. A flexible R package for nonnegative matrix factorization. BMC Bioinformatics. 2010;11:367.
Wilkerson MD, Hayes DN. ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinformatics. 2010;26(12):1572–3.
Maupin KA, Sinha A, Eugster E, Miller J, Ross J, Paulino V, et al. Glycogene expression alterations associated with pancreatic cancer epithelial-mesenchymal transition in complementary model systems. PLoS One. 2010;5(9): e13002.
Noll EM, Eisen C, Stenzinger A, Espinet E, Muckenhuber A, Klein C, et al. CYP3A5 mediates basal and acquired therapy resistance in different subtypes of pancreatic ductal adenocarcinoma. Nat Med. 2016;22(3):278–87.
Pei H, Li L, Fridley BL, Jenkins GD, Kalari KR, Lingle W, et al. FKBP51 affects cancer cell response to chemotherapy by negatively regulating Akt. Cancer Cell. 2009;16(3):259–66.
Badea L, Herlea V, Dima SO, Dumitrascu T, Popescu I. Combined gene expression analysis of whole-tissue and microdissected pancreatic ductal adenocarcinoma identifies genes specifically overexpressed in tumor epithelia. Hepatogastroenterology. 2008;55(88):2016–27.
Daemen A, Peterson D, Sahu N, McCord R, Du X, Liu B, et al. Metabolite profiling stratifies pancreatic ductal adenocarcinomas into subtypes with distinct sensitivities to metabolic inhibitors. Proc Natl Acad Sci U S A. 2015;112(32):E4410–7.
Semenza GL. Defining the role of hypoxia-inducible factor 1 in cancer biology and therapeutics. Oncogene. 2010;29(5):625–34.
Baek G, Tse YF, Hu Z, Cox D, Buboltz N, McCue P, et al. MCT4 defines a glycolytic subtype of pancreatic cancer with poor prognosis and unique metabolic dependencies. Cell Rep. 2014;9(6):2233–49.
Hui S, Ghergurovich JM, Morscher RJ, Jang C, Teng X, Lu W, et al. Glucose feeds the TCA cycle via circulating lactate. Nature. 2017;551(7678):115–8.
Little AC, Kovalenko I, Goo LE, Hong HS, Kerk SA, Yates JA, et al. High-content fluorescence imaging with the metabolic flux assay reveals insights into mitochondrial properties and functions. Commun Biol. 2020;3(1):271.
Masoud R, Reyes-Castellanos G, Lac S, Garcia J, Dou S, Shintu L, et al. Targeting mitochondrial complex I overcomes chemoresistance in high OXPHOS pancreatic cancer. Cell Rep Med. 2020;1(8): 100143.
Gallagher FA, Woitek R, McLean MA, Gill AB, Manzano Garcia R, Provenzano E, et al. Imaging breast cancer using hyperpolarized carbon-13 MRI. Proc Natl Acad Sci U S A. 2020;117(4):2092–8.
Sushentsev N, McLean MA, Warren AY, Benjamin AJV, Brodie C, Frary A, et al. Hyperpolarised (13)C-MRI identifies the emergence of a glycolytic cell population within intermediate-risk human prostate cancer. Nat Commun. 2022;13(1):466.
Park I, Larson PEZ, Gordon JW, Carvajal L, Chen HY, Bok R, et al. Development of methods and feasibility of using hyperpolarized carbon-13 imaging data for evaluating brain metabolism in patient studies. Magn Reson Med. 2018;80(3):864–73.
Aggarwal R, Vigneron DB, Kurhanewicz J. Hyperpolarized 1-[(13)C]-pyruvate magnetic resonance imaging detects an early metabolic response to androgen ablation therapy in prostate cancer. Eur Urol. 2017;72(6):1028–9.
Dutta P, Pando SC, Mascaro M, Riquelme E, Zoltan M, Zacharias NM, et al. Early detection of pancreatic intraepithelial neoplasias (PanINs) in transgenic mouse model by hyperpolarized (13)C metabolic magnetic resonance spectroscopy. Int J Mol Sci. 2020;21(10):3722. https://pubmed.ncbi.nlm.nih.gov/32466260/.
Serrao EM, Kettunen MI, Rodrigues TB, Dzien P, Wright AJ, Gopinathan A, et al. MRI with hyperpolarised [1-13C]pyruvate detects advanced pancreatic preneoplasia prior to invasive disease in a mouse model. Gut. 2016;65(3):465–75.
Brooks GA. The science and translation of lactate shuttle theory. Cell Metab. 2018;27(4):757–85.
Faubert B, Li KY, Cai L, Hensley CT, Kim J, Zacharias LG, et al. Lactate metabolism in human lung tumors. Cell. 2017;171(2):358-71 e9.
Draoui N, Feron O. Lactate shuttles at a glance: from physiological paradigms to anti-cancer treatments. Dis Model Mech. 2011;4(6):727–32.
Hayashi A, Yavas A, McIntyre CA, Ho YJ, Erakky A, Wong W, et al. Genetic and clinical correlates of entosis in pancreatic ductal adenocarcinoma. Mod Pathol. 2020;33(9):1822–31.
Open Access funding enabled and organized by Projekt DEAL. This work was supported by the grant of Wilhelm-Sander Stiftung (grant number: 2019.008.1) to M. T. A. and J. S. and by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) through 405344257 (SI 1549/3–2), SI1549/4–1, and CRC824 “Imaging for Selection, Monitoring and Individualization of Cancer Therapies” (SFB824; projects C4, C6, Z2, and A7); work in the lab of J. T. S. is further supported by the German Cancer Consortium (DKTK) and by the German Cancer Aid (no. 70112505/PIPAC, no. 70113834/PREDICT-PACA); and J. T. S. and W. W. are supported by the German Federal Ministry of Education and Research (BMBF; 01KD2206A/SATURN3). P. F Y. C is supported by the DFG (CH 2320/2–3).
Ethics approval and consent to participate
All relevant ethic approvals concerning use of patient material or mouse studies are appropriately mentioned in the supplementary methods section. Establishment of the PDX mouse model was performed using surgically resected PDAC tissues collected from patients at the Ruhr-University Bochum Comprehensive Cancer Center. Informed and written consent was obtained from all patients. The study was approved by the ethics committee of the Ruhr-University Bochum (permission nos. 3534–9, 3841–10, 16–5792). Patient tumor tissues were xenografted in both flanks of nude mice and expanded, isolated, and reimplanted for at least three generations. All animal experiments were performed according to the guidelines of the local Animal Use and Care Committees at the Ruhr-University Bochum (8.87–50.10.32.09.018, 84–02.04.2012.A328, and 81–02.04.2017.A423). For immunohistochemistry and multiplex immunofluorescence on human PDAC samples, a cohort of 31 patient samples from Radboud University Medical Center in Nijmegen, the Netherlands, was used. Given the retrospective nature of this study and the anonymized handling of data, informed consent was waived by the institutional review board (protocol CMO2018-4420). For HP-MRS animal experiments, approval of the animal protection and welfare review board was received prior to study initiation (ROB-55.2–2532.Vet_02-18–91). All experiments were carried out in adherence to pertinent laws and regulations.
Consent for publication
J. T. S. receives honoraria as consultant or for continuing medical education presentations from AstraZeneca, Bayer, Bristol-Myers Squibb, Immunocore, MSD Sharp Dohme, Novartis, Roche/Genentech, and Servier. His institution receives research funding from Abalos Therapeutics, Bristol-Myers Squibb, Celgene, Eisbach Bio, and Roche/Genentech; he holds ownership and serves on the Board of Directors of Pharma15, all outside the submitted work. DR received consultant and lecture fees from AstraZeneca, Merck Serono, Takeda, Pfizer, Novartis, Boehringer Ingelheim, Sanofi-Aventis, and BMS. DR received consultant and lecture fees from AstraZeneca, Merck Serono, Takeda, Pfizer, Novartis, Boehringer Ingelheim, Sanofi-Aventis, and BMS. DR is a founder and consultants of PearlRiver Bio GmbH and shareholder of Centessa Pharmaceuticals plc. DV reports the following: Gilead, Pfizer (advisory board), Bristol Myers Squibb (speaker’s honoraria), and AbbVie (travel support and congress registration fees). WW reports the following: advisory boards and speaker for Roche, MSD, BMS, AstraZeneca, Pfizer, Merck, Lilly, Boehringer, Novartis, Takeda, Bayer, Amgen, Astellas, Eisai, Illumina, Siemens, Agilent, ADC, GSK, and Molecular Health. Research funding is from Roche, MSD, BMS, and AstraZeneca. RV reports “pro bono” lectures and provides expert opinions at German court of law. The other authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Heid, I., Münch, C., Karakaya, S. et al. Functional noninvasive detection of glycolytic pancreatic ductal adenocarcinoma. Cancer Metab 10, 24 (2022). https://doi.org/10.1186/s40170-022-00298-5