Functional non-invasive detection of glycolytic pancreatic ductal adenocarcinoma

doi:10.21203/rs.3.rs-1753218/v1

Download PDF

Research Article

Functional non-invasive detection of glycolytic pancreatic ductal adenocarcinoma

https://doi.org/10.21203/rs.3.rs-1753218/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Background: 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.

Methods: We used different patient-derived PDAC model systems (conventional and primary patient derived cells, patient derived xenografts (PDX) and patient FFPE 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-¹³C]pyruvate and [1-¹³C]lactate magnetic resonance spectroscopy (HP-MRS) in PDAC xenografts.

Results: 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 in 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 hyperpolarized ¹³C-magnetic resonance spectroscopy (HP-MRS), we were able to non-invasively identify highly glycolytic PDAC xenografts by detecting the last glycolytic enzymatic step and prominent intra-tumoral [1‑¹³C]pyruvate and [1-¹³C]lactate interconversion in vivo.

Conclusion: Our study adds functional metabolic phenotyping to transcriptome-based analysis and proposes a functional approach to identify highly glycolytic PDACs as candidates for anti-metabolic therapeutic avenues.

glycolysis

PDAC

lactate

molecular subtype

hyperpolarized magnetic resonance spectroscopy

Despite enormous research efforts, pancreatic ductal adenocarcinoma (PDAC) remains a fatal disease with marginal clinical advancement (1). Although genomic and transcriptional profiles of PDAC have been studied in great detail (2–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 or pancreatic progenitor with more epithelial differentiated tumor and quasi-mesenchymal (QM; term used hereafter) or squamous or basal-like (1, 5–7) with more mesenchymal differentiated tumor respectively. Especially QM PDACs demonstrate very aggressive phenotypes with shorter median survival and resistance to first-line chemotherapy with FOLFIRINOX (1). 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 (8). This allows cells not only to adapt but to thrive on particularly scarce conditions of hypoxia and nutrient limitations typically observed in PDAC (9). Glycolysis is the most prominent cancer associated metabolic pathway. Although high cancer dependency on glucose was described nearly 100 years ago by Otto Warburg (10), 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 (11). 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 (12). 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 primary cells, patient derived xenografts 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 non-invasive hyperpolarized ¹³C-magnetic resonance spectroscopy, we were able to detect the final glycolytic step in vivo, namely intratumoral conversion of [1-¹³C]pyruvate to [1-¹³C]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 [1-¹³C]lactate to [1-¹³C]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 a non-invasive 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 #11966025 and DMEM #A1443001, Thermo Fisher Scientific, Waltham, USA) adapted to final concentrations of 5 mM D-glucose (Thermo Fisher Scientific, Waltham, USA), 2 mM Lglutamine, 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)

For all metabolic analysis, PDC cell lines were cultivated in a 1:1 mixture of Keratinocyte-SF medium (#17005075, Thermo Fisher Scientific, Waltham, USA) and RPMI 1640 ( #11879020, Thermo Fisher Scientific, Waltham, USA) adapted to final concentrations of 5mM D-glucose, 4.5mM L-glutamine, 0.26mM sodium pyruvate, 6%v/vFBS and 1% v/v final mixture of penicillin/streptomycin (Thermo Fisher Scientific, Waltham, USA) and Antimycotic/Antibiotic (cat. nu. 15240-062, Thermo Fisher Scientific) if not stated otherwise.

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: MCT4, Atlas Antibodies (Sigma Aldrich, Cat#HPA021451, RRID:AB_1853663); HIF1a, BD Transduction laboratories #610959 (RRID:AB_398272); MCT1, Abcam, #ab85021(RRID: AB_10674945); KRT81, Santa Cruz, #sc-100929 (RRID: AB_2132772); panCytokeratin, Abcam #ab6401(RRID: AB_305450);

Hyperpolarized Magnetic Resonance Spectroscopy (HP-MRS)

Animal handling

All experiments were carried out in adherence to pertinent laws and regulations.

Detailed explanations of all experiments can be found in supplementary material and methods section.

Glycolytic metabolic gene transcripts are present in QM PDAC subtype

To evaluate whether glycolytic transcripts are omnipresent in the QM PDAC subtype (11), 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 primary patient derived cells (PDC, n = 11). Transcriptomes from bulk tissue of 204 PDAC samples from previously published resource was utilized (E-MTAB-1791). All of the samples were then subtyped to QM or classical group. For tumor subtype determination, negative matrix factorization (NMF) (13) was used, after median centering of the data. ConcensusClusterPlus (14) 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 (15)); PDAC xenograft (E-MTAB-4029 (16)); and bulk PDAC tissue (GSE16515 and GSE15471 (17, 18)) 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 was performed. By using our subtyping platform, in the PDAC patient sample cohort (204 samples) 88 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) (19). 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 Fig. 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 (20), glycolysis/glucose metabolism transcripts were also enriched in the QM patient PDAC samples, PDX and PDAC cell line datasets, 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 5mM glucose with addition of 2mM 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 to 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 what 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 (21). In publicly available TCGA RNA expression datasets, high expression of both lactate transporters MCT4 and MCT1 correlated with worse survival in PDAC patients associating 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 Fig. 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 (16), 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 (22). However, whether this effect is especially attributable to lactate producing high glycolytic QM PDAC cells is still not known. To investigate this, we cultivated PDAC cells were cultivated for 7 hours in i) “basal” DMEM or RPMI media without glucose or glutamine supplementation or in ii) “basal” media supplemented with lactate (basal + 10mM 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 Fig. 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 Fig. 2b).

To substantiate this finding, we cultivated PSN1 (QM), PaTu8988T (QM) and PaTu8988S (classical) cells in DMEM medium with 5mM glucose and 2mM glutamine without media change for 24-48-72-96 hours. 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 hours), as expected due to glucose consumption and lactate production and accumulation. Once the glucose was consumed from the medium (approx. after 72 hours in PaTu8988T/PSN1 cells), lactate concentration in the media decreased, indicating that in absence of other resources, PDAC cells start consuming self-produced lactate (supplementary Fig. 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 cell titer 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-¹³C]pyruvate and [1-¹³C] 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 (¹⁸F-FDG) positron emission tomography (PET) is clinically established method for tumor detection based on high glucose uptake into the cancer. However, ¹⁸F-FDG-PET detects only the very first step of glycolysis since ¹⁸F-FDG gets phosphorylated by the first glycolytic enzyme hexokinase or glucokinase and does not enter the further metabolic processing. Thus FDG-PET rather 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 hyperpolarized magnetic resonance spectroscopy (HP-MRS) with hyperpolarized (HP) [1-¹³C]pyruvate and [1-¹³C]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 Fig. 3a). Once the tumors reached a minimal size of 5x5 mm, metabolic spectroscopy was performed. HP-[1-¹³C]pyruvate was i.v. injected into the tail vein and intratumoral accumulation of HP-[1-¹³C]lactate was followed in real-time. Using MRS, significantly more HP-[1-¹³C]lactate was detected in PSN1 compared to HPAC tumors, supporting higher ¹³C label exchange between pyruvate and lactate specifically in PSN1 tumors (Fig. 5a and 5b).

To evaluate whether lactate is used by tumors in vivo as observed in vitro in seahorse experiments, we performed the reverse experiment and intravenously injected HP-[1-¹³C]lactate in PSN1 and HPAC tumor rats in vivo. Intratumoral HP-[1-¹³C]pyruvate was detected in PSN1 but not HPAC xenografts (Fig. 5c). 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-¹³C]pyruvate to HP-[1-¹³C]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 xenografts by immunohistochemical analysis of glycolytic markers HIF1A and MCT4. 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). We next analyzed HIF1A and MCT4 expression in murine xenografts of human PDAC cell lines (supplementary Fig. 3b). Indeed, MCT4 staining intensity was also only lightly stronger in the QM xenografts. However, 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. 3b).

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, QM/squamous/basal-like and hybrid states with more or less pronounced subtype specific transcriptional programs (3–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 the worse prognosis both in resectable and metastatic setting (11). 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 (12). 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 presence of glycolytic transcripts is indeed translated into operable glycolysis in PDAC QM subtype and ii) offer a non-invasive 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 (5mM), thus omitting the metabolic artefacts 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 (19) and the recently reported “glycolytic” transcriptional PDAC subtypes in patients (11). 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 (23) and should be interpreted only as indication of the cellular energetic status. Glycolytic energetic status of PSN1 and PD69 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 (20), 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 (21).

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 (24). 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 mesenchymal and 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 glycolytic PSN1 tumors were readily detectable with HP-MRS due to higher ¹³C-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-¹³C]pyruvate to HP-[1-¹³C]lactate conversion rates identified strongly glycolytic aggressive triple negative breast cancer with high HIF1a and MCT1 tissue expression (25). This approach is already being used in personalized therapy monitoring in prostate and brain cancer (26, 27).

We were also able to confirm in vivo the active import and conversion of HP-[1-¹³C]lactate into HP-[1-¹³C]pyruvate in PSN1 QM-type but not in HPAC classical-type xenografts. Lactate has recently been considered as one of the important actors in tumor metabolism (28). 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 (22, 29). Indeed, we also observed OxPhos activation with lactate in PDAC cells, especially in the QM cell lines.

We speculate 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.

Our HP-MRS experiments provide evidence for the concept that PDACs with high reliance on glycolysis are potentially detectable via HP-[1-¹³C]pyruvate/lactate spectroscopy in clinical practice. Thus, identification of highly glycolytic, aggressive PDACs by HP-[1-¹³C]pyruvate and HP-[1-¹³C]lactate spectroscopy may be used to guide and monitor tumor treatment with anti-glycolytic therapies. In contrast to biopsy-based tumor characterization, metabolic imaging allows dynamic evaluation of the whole tumor limiting sampling bias and addressing tumor heterogeneity (30). Although likely not all QM tumors are potentially extremely glycolytic, non-invasive detection of highly glycolytic PDACs by HP-[1-¹³C]pyruvate/lactate MRS is one of the first methods for successful individual metabolic approaches.

PDAC-Pancreatic Ductal Adenocarcinoma

QM-quasi mesenchymal

HP-MRS- hyperpolarized magnetic resonance spectroscopy

PDX-patient derived xenografts

FFPE-formalin fixed paraffin embded

PDC-patient derived cells

MCT4-monocarboxylate transporter 4

MCT1-monocarboxylate transporter 4

LDH-lactate dehydrogenase

HIF-hypoxia inducible factor

NES-normalized enrichment score

FDR-false discovery rate

OCR- Oxygen Consumption Rate

ECAR- Extracellular Acidification Rate

TCGA-The Cancer Genome Atlas

DMEM-Dulbeco’s Modified Eagle Medium

RPMI-Roswell Park Memorial Institute

OxPhos-Oxidative Phosporylation

¹⁸F-FDG-18F-fluorodeoxyglucose

PET-positron emission tomography

Ethics approvals: 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 no. 3534-9, 3841-10, 16-5792). Patient tumor tissues were xenografted in both flanks of nude mice and expanded, isolated and re-implanted 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: Not applicable.

Availability of data and material: All data generated in this paper are available from the corresponding author upon a reasonable request.

Competing interests: JTS receives honoraria as consultant or for continuing medical education presentations from AstraZeneca, Bayer, Immunocore, Novartis, Roche/Genentech, Servier. His institution receives research funding from Bristol-Myers Squibb, Celgene, 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 Astra-Zeneca, Merck-Serono, Takeda, Pfizer, Novartis, Boehringer Ingelheim, Sanofi-Aventis and BMS. DR received consultant and lecture fees from Astra-Zeneca, 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: Gilead, Pfizer (Advisory Board), Pfizer, Bristol Myers Squibb (Speaker’s honoraria), Abbvie (Travel support and Congress registration fees). WW reports: 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 from Roche, MSD, BMS and AstraZeneca. RV reports “pro-bono” lectures and provides expert opinions at German court of law.

Funding: This work was supported by the grant of Wilhelm-Sander Stiftung (grant number: 2019.008.1) to M.T-A and J.S and the German Research Foundation (DFG) within the SFB-Initiative 824 (collaborative research center), “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), by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) - 405344257 (SI 1549/3-2) and SI1549/4-1 and by the German Cancer Aid (#70112505/PIPAC, #70113834/PREDICT-PACA); P.FY.C is supported by the DFG (CH 2320/2-3).

Author contributions: Heid I- In vivo study design, data collection and analysis, drafting and critical revision of the manuscript; Münch C, Karakaya S- data collection and analysis, critical revision of the manuscript; Lueong SS, Winkelkotte A, Liffers ST, Godfrey L, Cheung PFY-data collection and analysis, critical revision of the manuscript; Savvatakis K-material preparation; Topping G-method optimization and critical revision of the manuscript; Englert F, Kritzner L, Grashei M-method optimization, data collection; Tannapfel A, Viebahn R, Wolters H, Uhl W, Vangala D, Smeets EMM, Aarntzen EHJG-material availability; SH-material availability and critical revision of the manuscript; DR, WW-critical revision of the manuscript; JDH-material and data availability; FS, RB- resources, critical revision of the manuscript; MTA- study design, data collection and analysis, resources, drafting and critical revision of the manuscript; JT-study design, resources, drafting and critical revision of the manuscript;

All authors confirm access to the data and agree to manuscript submission.

Acknowldgements: Not applicable.

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.
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.
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.
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.
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.

No competing interests reported.

Download PDF

Editorial decision: Major revision
16 Aug, 2022
Reviews received at journal
07 Jul, 2022
Reviewers agreed at journal
05 Jul, 2022
Reviewers invited by journal
01 Jul, 2022
Editor assigned by journal
15 Jun, 2022
Submission checks completed at journal
15 Jun, 2022
First submitted to journal
13 Jun, 2022

You are reading this latest preprint version

Functional non-invasive detection of glycolytic pancreatic ductal adenocarcinoma

Status:

Version 1

Abstract

Figures

Background

Methods

PDAC cell lines

Patient Derived Cells (PDCs)

PDX samples preparation

Seahorse metabolic flux assays

Immunohistochemistry and immunofluorescence

Hyperpolarized Magnetic Resonance Spectroscopy (HP-MRS)

Animal handling

Results

Glycolytic metabolic gene transcripts are present in QM PDAC subtype

Glycolytic pathway activity is high in individual QM PDAC cells

PDAC cells actively use lactate as oxidative fuel

Hyperpolarized magnetic resonance spectroscopy of [1-¹³C]pyruvate and [1-¹³C] lactate identifies highly glycolytic tumors

Discussion

Abbreviations

Declarations

References

Additional Declarations

Supplementary Files

Status:

Version 1

Functional non-invasive detection of glycolytic pancreatic ductal adenocarcinoma

Status:

Version 1

Abstract

Figures

Background

Methods

PDAC cell lines

Patient Derived Cells (PDCs)

PDX samples preparation

Seahorse metabolic flux assays

Immunohistochemistry and immunofluorescence

Hyperpolarized Magnetic Resonance Spectroscopy (HP-MRS)

Animal handling

Results

Glycolytic metabolic gene transcripts are present in QM PDAC subtype

Glycolytic pathway activity is high in individual QM PDAC cells

PDAC cells actively use lactate as oxidative fuel

Hyperpolarized magnetic resonance spectroscopy of [1-13C]pyruvate and [1-13C] lactate identifies highly glycolytic tumors

Discussion

Abbreviations

Declarations

References

Additional Declarations

Supplementary Files

Status:

Version 1

Hyperpolarized magnetic resonance spectroscopy of [1-¹³C]pyruvate and [1-¹³C] lactate identifies highly glycolytic tumors