Main
Previously, we reported that 37.2% of primary PDAC tumours exhibit a loss of PHGDH expression6. These PHGDH-deficient PDAC cells lack the capacity for de novo serine synthesis and therefore rely on exogenous serine (exSer) for their survival and proliferation. In the primary tumour setting, these exSer-dependent PDAC cells use serine that is released from sensory nerves, and this constitutes a distinct and therapeutically targetable metabolic vulnerability6. However, it remains unclear how exSer-dependent PDAC cells obtain serine in the metastatic setting, and particularly in liver metastases, where the tumour microenvironment (TME) differs markedly from that of the primary tumour7. This lack of knowledge is due in part to a limited availability of clinical specimens, because patients with liver metastases are ineligible for surgical resection1.
Hepatocyte PHGDH is increased near PDAC
To examine the serine synthesis pathway (SSP) in PDAC liver metastases, both in cancer cells and in the surrounding metastatic microenvironment, we analysed the expression of PHGDH in human PDAC liver-metastasis samples collected through a rapid-autopsy programme8. We found that 10 out of 32 cases (31.3%) were PHGDH-negative in tumour cells (Fig. 1a,b), a frequency comparable to that observed in primary tumours6. Notably, hepatocytes adjacent to metastatic lesions, particularly those with weak or negative expression of PHGDH in tumour cells, exhibited strong PHGDH expression, whereas hepatocytes in non-metastatic liver tissue showed weak or negative staining (Fig. 1a,b). Dual immunofluorescence staining for PHGDH and albumin (a hepatocyte marker) confirmed this reciprocal pattern, with strong PHGDH expression in hepatocytes adjacent to PHGDH-negative PDAC metastases but little or no PHGDH expression around PHGDH-strong metastases (Fig. 1c). To validate this finding in vivo, we established liver metastases in nude mice using human PDAC cell lines with differential PHGDH expression: MiaPaCa2 (PHGDH-high) and PaTu8902 (PHGDH-null)6. Consistent with the human data, dual immunofluorescence staining revealed minimal expression of PHGDH in hepatocytes surrounding MiaPaCa2 metastases but a marked upregulation of PHGDH in those adjacent to PaTu8902 metastases (Extended Data Fig. 1a). Given the central role of hepatocytes in maintaining metabolic homeostasis, we hypothesized that hepatocytes upregulate serine synthesis to metabolically support exSer-dependent PDAC cells in the liver metastatic niche.
a, Representative haematoxylin and eosin (H&E) and PHGDH immunohistochemistry (IHC) of normal human livers and human PDAC liver metastases (32 cases each) stratified by tumour-cell PHGDH staining (strong, weak or negative; 10/32, 12/32 and 10/32 cases, respectively). Boxes indicate magnified regions; T, tumour; H, hepatocytes. Scale bars as indicated. b, Quantification of hepatocyte PHGDH intensity in a. Hepatocyte PHGDH staining (strong versus negative or weak) was compared between normal livers and each metastasis group stratified by tumour-cell PHGDH status. Numbers in bars indicate the number of cases. P values by two-sided Fisher’s exact test. c, Representative immunofluorescence of albumin (green), PHGDH (red) and DAPI (blue) in human PDAC liver metastases from a. Dashed lines indicate T–H boundaries. Scale bar, 20 μm. d, Scheme for hepatocyte-specific deletion of Phgdh in Cas9 mice. e, Tumour area in liver sections from Cas9 mice that were fed a −SG diet, quantified by MetFinder (n = 6, 7, 8 and 8 mice, from left to right). Data are mean ± s.e.m. P values by one-way ANOVA with Holm–Šidák’s multiple-comparisons test. f, Representative MetFinder images from e. g, Representative immunofluorescence of CK19 (green), PHGDH (red) and DAPI (blue) in sections from e. Dashed lines indicate T–H boundaries. Scale bars, 50 μm.
Source data
PHGDH loss in PDAC cells induces hepatocyte PHGDH
To determine the causal role of PDAC PHGDH loss in driving hepatocyte PHGDH upregulation, we knocked out Phgdh in HY19636 mouse PDAC cells. HY19636 cells were derived from a KPC pancreatic tumour (LSL-KrasG12D/+; Trp53flox/+; p48-Cre)9, a model that retains intact PHGDH expression and an active SSP (ref. 10). After confirming efficient Phgdh deletion and extracellular serine dependence in HY19636 cells transduced with Phgdh-targeting single-guide RNAs (sgRNAs; sgPhgdh), but not in tdTomato-targeting control cells (sgTOM) (Extended Data Fig. 1b,c), we intrasplenically injected these cells into C57BL/6 mice to establish liver metastases. Hepatocytes adjacent to sgPhgdh tumours showed strong PHGDH staining, whereas those near sgTOM tumours or in normal livers showed weak or negative staining (Extended Data Fig. 1d–f). Quantitative PCR (qPCR) analysis of hepatocytes isolated from metastasis-bearing livers further supported this upregulation (Extended Data Fig. 1g). These findings suggest that metastatic PDAC cells, particularly exSer-dependent ones, induce the expression of PHGDH in hepatocytes.
Hepatocyte PHGDH supports metastasis
To assess the functional contribution of hepatocyte-derived PHGDH to PDAC liver metastasis, we generated mice with hepatocyte-specific knockout of Phgdh, by using adeno-associated virus (AAV)-based in vivo CRISPR editing in Rosa-LSL-Cas9-EGFP mice (Cas9 mice)11,12 (Extended Data Fig. 2a). We used the AAV2/8-U6-sgRNA-TBG-Cre (AAV-sgRNA) vector, which drives hepatocyte-specific expression of Cre recombinase under the thyroid-hormone-binding globulin (TBG) promoter12. Efficient deletion of Phgdh was confirmed in EGFP-positive hepatocytes isolated from mice injected with AAV-sgPhgdh (Extended Data Fig. 2b). Notably, plasma biochemistry showed no evidence of liver dysfunction in these mice under a serine- and glycine-deficient (–SG) diet (Extended Data Fig. 2c), ruling out impaired liver function as a confounding factor in subsequent metastasis assays. To evaluate the effects of hepatocyte Phgdh loss on metastatic colonization, we intrasplenically injected sgTOM or sgPhgdh cells and monitored metastatic outgrowth in the liver (Fig. 1d). Under a –SG diet, the liver metastatic burden of sgPhgdh cells was significantly reduced in hepatocyte-specific Phgdh-KO mice of both sexes (Fig. 1e,f and Extended Data Fig. 2d,e). This suppression was partially rescued by a complete amino acid (AA) diet (Extended Data Fig. 2f), indicating that dietary serine can compensate for the loss of hepatocyte-derived serine. Supporting the role of hepatocyte PHGDH in this context, microcolonies of sgPhgdh PDAC cells in hepatocyte-specific Phgdh-KO mice were predominantly located near hepatocytes that retained residual PHGDH expression, probably reflecting incomplete knockout (Fig. 1g and Extended Data Fig. 2g,h). These findings show that hepatocyte-derived serine, through PHGDH activity, is required to support the colonization and outgrowth of exSer-dependent PDAC cells in the liver under dietary SG restriction. Targeting serine biosynthesis in hepatocytes might therefore offer a therapeutic strategy to limit liver metastasis in exSer-dependent PDAC cells.
Hepatocytes rescue exSer-dependent PDAC cells
To test the hypothesis that hepatocytes support PDAC cell growth by supplying de-novo-synthesized serine, we co-cultured primary mouse hepatocytes with PDAC cells and quantified PDAC cell growth and proliferation (Fig. 2a–d and Extended Data Fig. 3a,b). Under −SG conditions, the growth of exSer-dependent PDAC cells was markedly impaired, but it was restored by co-culture with hepatocytes (Fig. 2b,d and Extended Data Fig. 3b). This hepatocyte-dependent growth increase was not observed in PHGDH-proficient PDAC cells or in SG-containing medium (Fig. 2b,d and Extended Data Fig. 3b), indicating that hepatocyte-derived serine specifically promotes the proliferation of exSer-dependent PDAC cells in conditions in which extracellular serine is limited. These in vitro findings were mirrored in vivo: sgPhgdh PDAC cells exhibited increased proliferation at the tumour–liver interface, compared with the tumour core, suggesting that there is a spatial relationship between tumour growth and hepatocyte proximity (Extended Data Fig. 3c,d).
a, Schematic of transwell co-culture of PDAC cells with primary mouse hepatocytes under +SG or −SG conditions. b, Growth and viability of PaTu8902 (PHGDH-null) and HY19636 PDAC cells, quantified by crystal violet staining. c, Schematic of direct co-culture of mouse PDAC cells with primary mouse hepatocytes. d, Percentage of EdU+ PDAC cells in c by flow cytometry. e, Schematic of serine synthesis from 13C6-glucose by the SSP. Red circles denote 13C. f, Design for 13C6-glucose tracing in PDAC cells cultured alone or with hepatocytes under −SG conditions. GC–MS, gas chromatography–mass spectrometry. g, Intracellular serine isotopologue distribution in HY19636 PDAC cells from f. h, qPCR analysis of SSP gene expression in hepatocytes cultured alone or co-cultured with sgTOM or sgPhgdh HY19636 cells. i, Schematic of the CM assay. CM was collected from cell-free dishes (control CM) or from dishes containing sgTOM or sgPhgdh HY19636 PDAC cells (sgTOM CM and sgPhgdh CM) after 24 h in SG-free DMEM (T0 CM). Primary mouse hepatocytes were pre-starved of SG and glucose for 24 h, then incubated with T0 CM supplemented with 13C6-glucose (25 mM) for 24 h. Hepatocytes and media were collected at 24 h (T24 CM) for the GC–MS measurements shown in j–l. j, Intracellular distribution of serine isotopologues in hepatocytes from i. k, Intracellular serine concentration in hepatocytes from i. l, Secreted serine in media from i after subtraction of baseline serine in T0 CM. Data are mean ± s.d. (b) or mean ± s.e.m. (d,g,j–l). n (biological replicates): b, 6; d, 4; g, 3; h, 2; j–l, 3. Statistical significance by one-way ANOVA with Tukey’s multiple-comparisons test (b,d,k,l). Exact P values are shown in the relevant panels.
Source data
To determine whether activation of the SSP leads to functional serine production and transfer to cancer cells, we performed co-culture assays with 13C6-glucose (Fig. 2e,f). Labelled M+3 serine, a product of the SSP, was detected in sgPhgdh PDAC cells only when co-cultured with hepatocytes, and not in monoculture (Fig. 2g and Extended Data Fig. 3e,f), suggesting that hepatocyte-derived serine is synthesized de novo, released into the medium and taken up by PDAC cells.
Consistently, qPCR and RNA sequencing (RNA-seq) analyses revealed an upregulation of SSP genes in hepatocytes that were co-cultured with sgPhgdh HY19636 cells (Fig. 2h), along with a trend towards enrichment of signatures of glycine–serine–threonine metabolism (Extended Data Fig. 3g). These transcriptional changes recapitulate what was observed in patients with liver metastases and in mouse models (Fig. 1 and Extended Data Fig. 1). Notably, −SG conditions alone did not clearly induce the expression of SSP genes in monocultured hepatocytes, whereas co-culture with exSer-dependent PDAC cells did (Fig. 2h), suggesting that PDAC-derived signals have a role in hepatocyte metabolic reprogramming.
PDAC-derived CM drives hepatocyte serine supply
To directly assess whether hepatocytes engage in de novo serine synthesis and actively secrete serine in response to tumour-derived cues, we used conditioned medium (CM) collected from either sgTOM or sgPhgdh HY19636 cells after 24 h of culture in Dulbecco’s modified Eagle’s medium (DMEM) lacking SG (T0 CM). Primary mouse hepatocytes were pre-starved of SG and glucose for 24 h, then switched to T0 CM supplemented with 13C6-glucose (Fig. 2i). Hepatocytes cultured in CM derived from sgPhgdh cells exhibited significantly higher fractional enrichment of M+3-labelled serine, as well as increased total intracellular levels of serine, compared with those that were cultured in CM from sgTOM cells or control medium (Fig. 2j,k). Notably, the concentration of serine in the culture medium was highest in hepatocytes grown in CM from sgPhgdh cells, after subtracting baseline serine concentrations (T0 CM) (Fig. 2l), reaching a concentration range that was sufficient to increase the growth of sgPhgdh cells (Extended Data Fig. 3h). These findings suggest that paracrine factors secreted by exSer-dependent PDAC cells promote the biosynthesis and secretion of serine in hepatocytes, consistent with a tumour–host metabolic cross-talk that supports the growth of serine-auxotrophic cancer cells.
CXCL5–CXCR2 drives serine cross-talk
We next sought to identify soluble factors secreted by exSer-dependent PDAC cells that induce PHGDH expression in hepatocytes. RNA-seq analysis identified Cxcl5 as a candidate gene, which was upregulated (1) in sgPhgdh HY19636 cells compared with sgTOM HY19636 cells under co-culture with hepatocytes, and (2) in sgPhgdh HY19636 cells grown in co-culture with hepatocytes compared with those grown in monoculture (Fig. 3a and Supplementary Table 1). This upregulation of CXCL5 was further confirmed at the protein level in co-culture systems (Extended Data Fig. 4a) and in tissue interstitial fluids (TIFs) from mouse liver metastases (Fig. 3b). Notably, the increased expression of CXCL5 persisted for at least 14 days in vitro in sgPhgdh HY19636 cells after isolation from liver metastases (Extended Data Fig. 4b,c). CXCL5 concentrations were also increased in human liver-metastasis samples, with the highest levels in tissues bearing PHGDH-weak or PHGDH-negative metastases, lower levels in PHGDH-strong metastases and undetectable levels in normal liver tissue (Fig. 3c).
a, RNA-seq workflow. The candidate gene Cxcl5 is shown in bold. FDR, false discovery rate; KO, knockout. b, CXCL5 concentrations in TIFs from mouse liver metastases, quantified by ELISA (n = 3 mice per group). c, CXCL5 concentrations in human liver homogenates, quantified by ELISA (n as indicated). d, Growth of sgPhgdh (SKO) or sgPhgdh/sgCxcl5 (DKO) HY19636 cells in the transwell assay ± rCXCL5, quantified by crystal violet staining (n = 4 biological replicates). e, Intracellular serine levels in HY19636 PDAC cells, measured using a fluorometric serine assay. Data are shown relative to sgTOM cells without rCXCL5 treatment (n = 3 biological replicates). f,g, Whole liver weight (f) and tumour area (g) in liver sections from mice (fed a −SG diet) after intrasplenic injection of SKO or DKO cells (n = 9 (left) and n = 10 (right) mice in f; n = 5 (left) and n = 7 (right) mice in g). Representative MetFinder images are shown on the right in g. Data in f,g are from independent experiments. h, Representative immunofluorescence of albumin (green), PHGDH (red) and DAPI (blue) in SKO or DKO liver metastases. Dashed lines indicate tumour (T)–hepatocyte (H) boundaries. Scale bars, 50 μm. i, PHGDH intensity in hepatocytes adjacent to SKO or DKO metastases (n = 4 mice per group). A.u., arbitrary units. j, Scheme for hepatocyte-specific knockout of Cxcr2 in Cas9 mice. k, Tumour area in liver sections from j (mice fed a −SG diet), quantified by MetFinder (n = 5, 5, 9 and 9 mice, from left to right). l, Representative MetFinder images from j. Data are mean ± s.e.m. (b,c,f,g,i,k) or mean ± s.d. (d,e). Statistical significance by one-way ANOVA with Tukey’s multiple-comparisons test (b–e,k) or by two-sided unpaired Welch’s t-test (f,g,i). Each symbol represents one independent sample. Exact P values are shown in the relevant panels.
Source data
We previously showed that exSer-dependent PDAC cells in primary tumours rely on nerve growth factor (NGF)-driven tumour innervation and neuronal serine supply6. Consistently, sgPhgdh HY19636 cells showed increased production of NGF under –SG conditions, but this induction was blunted by co-culture with hepatocytes (Extended Data Fig. 4d). Moreover, mouse liver metastases showed only a non-significant trend towards increased NGF in TIFs and no increase in tumour innervation (Extended Data Fig. 4e,f). Conversely, sgPhgdh tumours in the pancreas did not upregulate Cxcl5 after dietary serine restriction (Extended Data Fig. 4g). Together, these data suggest that exSer-dependent PDAC cells engage distinct tissue-dependent paracrine programs for serine acquisition, with the NGF–nerve axis predominating at the primary site and CXCL5 emerging as a candidate mediator of tumour–hepatocyte cross-talk in liver metastases.
CXCL5 enables hepatocyte-mediated serine supply
To evaluate the role of tumour-derived CXCL5 in hepatocyte-mediated serine supply, we deleted Cxcl5 in both sgTOM and sgPhgdh HY19636 cells (Extended Data Fig. 4h) and assessed their growth in co-culture. Ablation of Cxcl5 did not affect PDAC proliferation in monoculture, but abolished the hepatocyte-mediated rescue of sgPhgdh cell growth under −SG conditions (Fig. 3d). This effect was accompanied by reduced expression of PHGDH in hepatocytes (Extended Data Fig. 4i) and decreased intracellular serine levels in co-cultured PDAC cells (Fig. 3e). All of these effects were fully rescued by recombinant CXCL5 (rCXCL5). Pharmacological inhibition of CXCR2, the cognate receptor for CXCL5, suppressed the CXCL5-induced upregulation of Phgdh in hepatocytes to the level observed in monoculture (Extended Data Fig. 4j). Consistent with these findings, CXCR2 inhibition also downregulated SSP genes beyond Phgdh (Extended Data Fig. 4k).
CXCL5–CXCR2 supports PHGDH-null PDAC
To test this interaction in vivo, Phgdh single-knockout (sgPhgdh; SKO) or Phgdh/Cxcl5 double-knockout (sgPhgdh/sgCxcl5; DKO) HY19636 cells were injected intrasplenically into C57BL/6 mice maintained on a −SG diet. Mice that were injected with DKO cells showed a significantly lower liver metastatic burden than did those injected with SKO cells, as measured by whole liver weight and tumour area in the liver sections (Fig. 3f,g). Notably, PHGDH expression was reduced in hepatocytes adjacent to DKO metastases (Fig. 3h,i), and the proliferation of tumour cells also decreased (Extended Data Fig. 4l,m). Given the role of the CXCL5–CXCR2 axis in immune-cell recruitment13, we generated hepatocyte-specific Cxcr2-knockout mice using AAV-mediated in vivo CRISPR editing, to test hepatocyte CXCR2 signalling directly (Fig. 3j and Extended Data Figs. 2a and 4n). When sgTOM or sgPhgdh HY19636 cells were injected into these mice, ablation of Cxcr2 in hepatocytes significantly reduced the liver metastatic burden of sgPhgdh—but not sgTOM—PDAC cells (Fig. 3k,l). These findings mirror the effects observed in hepatocyte-specific Phgdh-knockout mice (Fig. 1), and collectively support a model in which CXCL5–CXCR2 signalling in hepatocytes drives de novo serine synthesis to support exSer-dependent PDAC outgrowth during liver metastasis.
CXCR2–AKT derepresses PHGDH through FOXO3A
To delineate the signalling axis that is activated in hepatocytes after co-culture with sgPhgdh HY19636 cells under CXCR2 signalling, we performed over-representation analysis (ORA) in RNA-seq data from primary mouse hepatocytes. Specifically, we compared hepatocytes cultured with or without sgPhgdh HY19636 cells and, under co-culture conditions, with or without the CXCR2 inhibitor SB225002. In both comparisons, the PI3K–AKT signalling pathway was significantly over-represented (Fig. 4a,b). Consistently, co-culture with PDAC cells increased AKT phosphorylation and PHGDH expression in hepatocytes (Extended Data Fig. 5a). Pharmacological inhibition of either CXCR2 or PI3K–AKT blocked this upregulation of PHGDH (Extended Data Fig. 5a) and suppressed de novo serine synthesis from 13C6-glucose (Figs. 2e and 4c,d and Extended Data Fig. 5b). Consequently, the growth-promoting effect of hepatocytes on exSer-dependent PDAC cells under −SG conditions was lost (Extended Data Fig. 5c,d). These results identify the PI3K–AKT pathway as a crucial node that mediates PHGDH upregulation in hepatocytes downstream of the CXCL5–CXCR2 axis.
a,b, Top KEGG pathways identified by ORA of RNA-seq data from primary mouse hepatocytes cultured alone or with sgPhgdh HY19636 cells (a), or co-cultured with sgPhgdh HY19636 cells with or without CXCR2i (SB225002, 100 nM) (b). Gene sets with adjusted P < 0.05 were considered significantly enriched. c, Experimental design for 13C6-glucose tracing in PDAC cells co-cultured with hepatocytes with or without the indicated inhibitors for 24 h. d, Intracellular serine fractional labelling measured by GC–MS. Stacked bars indicate mean ± s.d. from n = 3 biological replicates. e, Proposed model for PHGDH regulation mediated by the CXCL5–CXCR2–PI3K–AKT–FOXO3A axis in hepatocytes. f, Experimental design for AML12 hepatocytes co-cultured with sgPhgdh HY19636 cells treated with dimethyl sulfoxide (DMSO), CXCR2i (SB225002, 100 nM) or PI3Ki (BKM120, 50 nM), followed by immunofluorescence and chromatin immunoprecipitation (ChIP). g, Representative FOXO3A immunofluorescence in AML12 hepatocytes from f. White boxes indicate magnified regions. Scale bars, 20 μm. Right, FOXO3A nuclear localization quantified as the correlation coefficient between DAPI and FOXO3A signals (n = 3 biological replicates). Data are mean ± s.d. Statistical significance by one-way ANOVA with Dunnett’s multiple-comparisons test. h, Left, TGTTT-core FOXO consensus-like motifs in the Phgdh promoter (P1 and P2). P3, distal gene-body region (negative control). The TGTTT motif is shown in red for clarity in the bottom sequences. Right, FOXO3A ChIP–qPCR at the Phgdh locus in AML12 hepatocytes, shown as fold change (FOXO3A/IgG) normalized to 2% input. The dotted line indicates IgG control, set to 1 (n = 3 biological replicates). Data are mean ± s.d. Statistical significance by one-way ANOVA with Holm–Šidák’s multiple-comparisons test versus untreated co-culture. Exact P values are shown in the relevant panels.
Source data
Among the candidate transcription factors that might regulate Phgdh expression, we focused on FOXO3A for two reasons: (1) FOXO3 is inactivated by PI3K–AKT signalling through phosphorylation-dependent cytoplasmic sequestration and subsequent proteasomal degradation14,15; and (2) FOXO3 can function as both a repressor and an activator of transcription, depending on the target genes16,17,18. Consistently, gene set enrichment analysis (GSEA) revealed a significant negative enrichment of the FOXO signalling pathway in hepatocytes co-cultured with sgPhgdh PDAC cells, relative to monoculture, whereas CXCR2 inhibition reversed this pattern to positive enrichment (Extended Data Fig. 5e,f). Similarly, GSEA using transcription factor–target-gene sets showed a negative enrichment of FOXO3 target genes in co-cultured hepatocytes, with a trend towards positive enrichment after inhibition of CXCR2 (Extended Data Fig. 5g), further supporting FOXO3A as a plausible downstream effector of PI3K–AKT signalling in this context.
On the basis of these findings, we hypothesized that nuclear FOXO3 represses Phgdh expression in hepatocytes, and that PI3K–AKT activation downstream of the CXCL5–CXCR2 axis alleviates this repression by promoting FOXO3 phosphorylation and cytoplasmic retention14,15,19 (Fig. 4e). Supporting this model, pharmacological inhibition of CXCR2 or PI3K in hepatocytes co-cultured with sgPhgdh PDAC cells enhanced FOXO3 nuclear localization (Fig. 4f,g) and increased FOXO3 occupancy at TGTTT-core FOXO consensus-like motifs in the Phgdh promoter (P1 and P2), but not at a distal gene-body region that was used as a negative control (P3)18 (Fig. 4h). Of note, this regulatory relationship was also evident in human liver-metastasis samples: hepatocytes with strong PHGDH expression exhibited predominantly non-nuclear FOXO3, whereas hepatocytes with weak or absent PHGDH expression exhibited predominantly nuclear FOXO3 (Extended Data Fig. 5h). Consistent with this, inhibiting CXCR2 in mice bearing sgPhgdh PDAC liver metastases promoted FOXO3 nuclear localization in hepatocytes (Extended Data Fig. 5i). Together, these findings establish FOXO3 as a transcriptional repressor of PHGDH, whose activity is relieved by PI3K–AKT signalling downstream of the CXCL5–CXCR2 axis.
FOXO3A is phosphorylated by AKT at Thr32, Ser253 and Ser315, which promotes its cytoplasmic retention19. To investigate the functional importance of these phosphorylation sites, we used a non-phosphorylatable FOXO3A mutant (T32A/S253A/S315A; FOXO3A(AAA)) by replacing these three residues with alanine19. Endogenous FOXO3A was deleted in a mouse hepatocyte cell line (AML12) and replaced with either wild-type (WT) or mutant (AAA) constructs (Extended Data Fig. 6a). After inhibition of PI3K, FOXO3A(WT) translocated to the nucleus, whereas the non-phosphorylatable FOXO3A(AAA) mutant remained constitutively nuclear (Extended Data Fig. 6b,c). Consistent with its nuclear localization, FOXO3A(AAA) expression significantly reduced the activity of the Phgdh promoter, mimicking the effect of PI3K inhibition (Extended Data Fig. 6d,e). In co-culture with exSer-dependent PDAC cells, FOXO3A(AAA)-expressing hepatocytes showed enhanced FOXO3A binding to the Phgdh promoter (Extended Data Fig. 6f) and impaired serine production, as evidenced by reduced incorporation of 13C6 into co-cultured PDAC cells (Extended Data Fig. 6g,h). Moreover, both Cxcr2 deletion and FOXO3A(AAA) expression in AML12 cells similarly suppressed their capacity to support the growth of exSer-dependent PDAC cells under −SG conditions, supporting the functional relevance of FOXO3A-mediated PHGDH regulation downstream of CXCR2–PI3K–AKT signalling (Extended Data Fig. 6i–k). Together, these findings identify FOXO3A as a transcriptional repressor of Phgdh that is inactivated by CXCL5–CXCR2–PI3K–AKT signalling in hepatocytes, representing an important downstream component of the hepatocyte serine-supplying program that supports exSer-dependent PDAC cells during liver metastasis.
CXCR2 inhibition and –SG suppress PHGDH-null PDAC
Our findings thus far show that genetic ablation of Cxcl5 in PDAC cells, or targeted deletion of Phgdh or Cxcr2 in hepatocytes (Figs. 1 and 3), effectively suppresses the growth of exSer-dependent PDAC under dietary serine restriction. Several CXCR2 (or CXCR1 and CXCR2; CXCR1/2) inhibitors have entered clinical testing for PDAC, including SX-682 (NCT04477343). Furthermore, the safety and feasibility of a −SG diet have been demonstrated in a phase I clinical trial20. To evaluate the translational potential of our findings, we next tested the therapeutic efficacy of combining a pharmacological CXCR2 inhibitor (SB225002) with dietary serine restriction in preclinical models of exSer-dependent PDAC.
To this end, mice were intrasplenically injected with sgTOM or sgPhgdh HY19636 cells and subsequently maintained on either an AA control diet or a –SG diet, with concurrent administration of vehicle or a CXCR2 inhibitor (CXCRi; SB225002) (Fig. 5a). Notably, the combination of CXCR2 inhibition and a −SG diet did not affect liver weights in mice that underwent sham surgery (Extended Data Fig. 7a,b), suggesting the potential safety and translational relevance of this strategy for future clinical applications. CXCR2 inhibition significantly reduced the liver metastatic burden of sgPhgdh PDAC cells in mice that were fed a –SG diet, as measured by liver weight and metastatic area (Fig. 5b–d). This effect was not observed in mice injected with sgTOM cells or in those maintained on a complete AA control diet (Fig. 5b–d), suggesting that the therapeutic effect of CXCR2 inhibition is mediated mainly by the disruption of hepatocyte-derived serine supply rather than by immune modulation. Consistently, PHGDH expression was robustly upregulated in hepatocytes adjacent to sgPhgdh metastases, and this effect was diminished after CXCR2 inhibition (Extended Data Fig. 7c,d), further supporting a model in which CXCL5–CXCR2 signalling induces hepatocyte-derived serine synthesis to support the metastatic outgrowth of exSer-dependent PDAC when dietary serine is restricted. To assess therapeutic benefits beyond tumour burden, we repeated the intrasplenic injection experiment using a lower number of injected cells to monitor survival. Mice bearing sgPhgdh metastases and treated with SB225002 survived significantly longer than did untreated controls (Fig. 5e).
a, Experimental scheme for b–d. b, Whole liver weight (n = 5, 5, 5, 5, 7, 7, 8 and 8 mice, from left to right). Blue dotted line indicates the mean liver weight in sham-operated mice (958.2 mg; n = 17 mice; Extended Data Fig. 7a). c, Tumour area in liver sections, quantified by MetFinder (n = 5, 5, 5, 5, 6, 6, 6 and 6 mice, from left to right). d, Representative MetFinder images from c. e, Survival of mice bearing sgTOM or sgPhgdh HY19636 liver metastases under a −SG diet, with or without treatment with a CXCR2 inhibitor (n = 12 mice per group). Survival curves were compared by two-sided log-rank test. f, Experimental scheme for g–i. OT, orthotopic injection; PTR, primary tumour resection. g, Resected primary tumour weight (n = 4, 4, 4, 5, 5 and 5 mice, from left to right). h, Percentage of tumour area in liver sections, quantified by MetFinder (n = 4, 4, 4, 5, 5 and 5 mice, from left to right). i, Representative MetFinder images from h. j, Model of liver-metastasis-promoting cross-talk between exSer-dependent PDAC cells and hepatocytes. exSer-dependent PDAC cells activate CXCL5–CXCR2–PI3K–AKT signalling in hepatocytes, leading to FOXO3A phosphorylation, nuclear exclusion, PHGDH derepression and hepatocyte serine production. Right, the combination of a −SG diet and CXCR2 inhibition suppresses this axis, restores FOXO3A nuclear localization, reduces hepatocyte PHGDH expression and serine production and inhibits the proliferation of metastatic PDAC cells. Data are mean ± s.e.m. (b,c,g,h). Statistical significance by one-way ANOVA with Tukey’s multiple-comparisons test (b,c,g,h). Each symbol represents an individual mouse. Exact P values are shown in the relevant panels.
Source data
Although the intrasplenic injection model is well suited to the study of later stages of metastasis, including tumour-cell engraftment and expansion in the liver, it bypasses early steps of the metastatic cascade. To address this limitation and better recapitulate clinical progression after primary tumour resection, we used a spontaneous metastasis model21. sgTOM or sgPhgdh HY19636 cells were orthotopically implanted into the pancreas. Three weeks later, primary tumours were surgically resected, and mice were euthanized at week 7 (Fig. 5f). To assess the stage-dependent therapeutic effects of CXCR2 inhibition, SB225002 was administered either 2 days after (CXCR2i (during)) or 2 weeks after (CXCR2i (after)) tumour implantation. CXCR2 inhibition had no effect on primary tumour growth regardless of Phgdh status (Fig. 5g), but significantly reduced the liver metastasis of sgPhgdh PDAC cells when administered at either time point (Fig. 5h,i). These data show that the CXCL5–CXCR2 axis has an essential role in sustaining the liver metastatic outgrowth of exSer-dependent PDAC under serine restriction, but that it has little effect on primary tumour growth.
Discussion
How exSer-dependent PDAC cells secure serine within liver metastases, and whether this involves organ-specific stromal metabolic responses, has remained unclear. Here we show that loss of PHGDH and dependence on exSer reprogram AA metabolism in hepatocytes through a CXCL5–CXCR2–PI3K–AKT–FOXO3A–PHGDH axis, thereby supporting the outgrowth of PDAC in the liver metastatic niche (Fig. 5j).
We observed a loss of PHGDH in 31.3% of liver metastases—a level comparable with that reported in primary PDAC tumours (37.2%)6—which suggests that SSP deficiency does not prevent liver metastatic colonization (Fig. 1a,b). These SSP-deficient cells might compensate by acquiring extracellular serine, and might even gain a selective advantage from this dependency. As reported previously in brain tumours, the uptake of extracellular serine can spare glucose-derived carbon by reducing cells’ reliance on de novo serine synthesis, and can redirect it towards other biosynthetic processes22; a similar mechanism might support PDAC growth if glucose is constrained in a similar manner in liver metastases, as has been described in primary PDAC23.
Our study also indicates that serine acquisition is organ-dependent in this PDAC subtype. Whereas primary tumours rely on NGF-driven tumour innervation and neuronal serine supply6, liver metastases show no clear increase in innervation but instead engage a CXCL5-associated, hepatocyte-dependent serine supply program (Fig. 3b,c and Extended Data Fig. 4a–f). Consistently, dietary serine restriction did not induce the expression of Cxcl5 in pancreatic primary tumours (Extended Data Fig. 4g). Notably, hepatocyte co-culture attenuated NGF induction while enhancing CXCL5 expression in Phgdh-knockout PDAC cells (Fig. 3a and Extended Data Fig. 4d), further supporting the idea that distinct microenvironments favour different paracrine programs for serine acquisition. These findings also raise the possibility of bidirectional PDAC–hepatocyte communication, although the hepatocyte-derived factor(s) responsible remain unidentified. Together, these observations suggest that exSer-dependent PDAC cells use different organ-specific stromal circuits to satisfy the same metabolic need. Thus, CXCL5–CXCR2 signalling should not be viewed as a universal mechanism that supports all PHGDH-deficient lesions, but rather as a pathway directed at the liver metastatic niche, with particular therapeutic relevance in peri-operative or post-resection settings, in which suppression of liver metastatic recurrence is a major goal.
Because the CXCL5–CXCR2 axis is also implicated in immune-cell recruitment13, whether immune mechanisms contribute to the anti-metastatic effects observed here remains an important question. However, the selective reduction of Phgdh-deficient PDAC metastases in hepatocyte-specific Cxcr2-knockout mice, together with the loss of CXCR2 inhibition efficacy under a complete AA control diet, supports hepatocyte CXCR2-dependent metabolic responses as a major contributor in this model (Figs. 3j–l and 5b–d). Future studies will be needed to define how this metabolic axis interacts with immune responses.
Our data support combining dietary SG restriction with CXCR2 inhibition as a liver-metastasis-directed strategy in this PDAC subtype. Although this combination showed no clear efficacy against the pancreatic primary tumour within the time frame of our study, it suppressed liver metastatic burden, prolonged survival and reduced liver metastatic recurrence after primary tumour resection (Fig. 5g–i). Because CXCR1/2 inhibitors are in clinical testing for PDAC, including SX-682 (NCT04477343), and because an –SG diet has shown safety and feasibility in a phase I trial20, this approach might be practical as a time-limited peri-operative intervention for patients with PHGDH-deficient, exSer-dependent PDAC, given the physiological roles of serine in normal tissues24,25. PHGDH immunostaining of resected specimens, or of biopsy samples obtained pre-operatively by endoscopic ultrasound-guided fine-needle aspiration (EUS-FNA), could be used to help to identify patients who are most likely to benefit. Clinical implementation will require careful optimization of treatment duration, nutritional support and biomarker-guided monitoring.
Direct inhibition of PHGDH represents a key alternative strategy for targeting serine metabolism, and pharmacological PHGDH inhibition has been shown to cooperate with dietary SG restriction to suppress tumour growth in vivo26. However, PHGDH inhibitors that are currently available act systemically rather than selectively in hepatocytes, and might therefore suppress compensatory de novo serine biosynthesis in normal tissues under extracellular serine limitation. By contrast, our rationale for CXCR2 inhibition is to block the tumour-induced activation of a hepatocyte-derived serine supply while potentially preserving this compensatory pathway in normal tissues, thereby exploiting the differential vulnerability of exSer-dependent PDAC cells. Nevertheless, pharmacological inhibition of CXCR2 is not hepatocyte-specific and is likely to exert pleiotropic effects throughout the metastatic microenvironment, including on immune and stromal cells. Accordingly, the anti-metastatic effects observed with systemic CXCR2 inhibition in this study cannot be attributed solely to the blockade of CXCR2 signalling in hepatocytes. Thus, whether CXCR2 inhibition provides a therapeutic advantage over systemic PHGDH inhibition remains unresolved and will require direct comparisons of their efficacy and tolerability under dietary SG restriction.
Several limitations should be noted. First, because the human analyses used rapid-autopsy specimens, which reflect end-stage disease, tumour serine demand and hepatocyte serine supply might be accentuated by reduced food intake and cachexia-associated declines in circulating serine in the terminal phase of disease. Thus, the liver metastatic niche during early colonization warrants separate investigation, as does the conservation of PHGDH status between matched primary tumours and liver metastases. Second, the mechanisms that drive CXCL5 induction in PDAC cells and serine efflux from hepatocytes remain unresolved, although both represent processes that might be targetable. Third, the extent to which hepatocyte-derived serine supports tumour growth in human disease, and whether this dependency is even greater under standard dietary conditions than it is in mouse models, remains unclear. If the latter is the case, dietary serine restriction could prove dispensable in patients, and CXCR2 inhibition alone might suffice in selected settings. Although our absolute quantification experiments support the feasibility of a hepatocyte-derived serine supply under the tested culture conditions (Fig. 2i–l and Extended Data Fig. 3h), they do not directly compare hepatocyte serine export with tumour-cell uptake demand in vivo. Addressing these questions will require future metabolic flux analyses, including in vivo 13C-glucose tracing.
In summary, we identify a CXCL5–CXCR2–PI3K–AKT–FOXO3A–PHGDH axis through which PHGDH-deficient, exSer-dependent PDAC rewires hepatocyte AA metabolism in the liver metastatic niche to support metastatic outgrowth. These findings support a model in which tumour-intrinsic metabolic states can elicit organ-specific stromal programs, providing a basis for tailoring therapy according to both tumour metabolism and anatomical site.
Methods
Human samples
Snap-frozen and formalin-fixed paraffin-embedded (FFPE) liver tissue sections from de-identified patients with PDAC and individuals without PDAC were obtained from the Rapid Autopsy Program at the University of Nebraska Medical Center (UNMC), in accordance with institutional review board (IRB) approval (IRB 091-01). Written informed consent for the research use of these specimens was obtained from all participants in accordance with the approved protocol. To ensure specimen quality, organs were collected within 3 h post-mortem, after which tissues were either flash-frozen in liquid nitrogen or immediately fixed in formalin. For histological analyses, FFPE tissues were sectioned at a thickness of 4 μm and mounted onto charged glass slides. Tissue homogenates were prepared from snap-frozen samples using Lysis Buffer (R&D Systems, 895347) and a bead-based tissue lyser (TissueLyser LT; QIAGEN), according to the manufacturers’ protocols.
Immunohistochemistry and immunofluorescence (manual staining)
FFPE sections of mouse and human liver tissue (5 µm) were baked at 65 °C for 1 h, cooled to room temperature, deparaffinized in xylene and rehydrated through a graded ethanol series. Antigen retrieval was performed by heating sections in 0.05 M sodium citrate buffer (pH 6.0) for 20 min using a microwave oven. Endogenous peroxidase activity was quenched by incubation in 3% hydrogen peroxide in methanol for 10 min. Sections were blocked for 1 h at room temperature with 2.5% horse serum (Vector Laboratories, MP-7800-15) for immunohistochemistry (IHC) or with 10% normal goat serum (NGS; Vector Laboratories, S-1000-20) for immunofluorescence. Sections were incubated overnight at 4 °C in a humidified chamber with primary antibodies: anti-PHGDH (Sigma, HPA021241, 1:100), anti-albumin (Abcam, ab106582, 1:200), anti-CK19 (DSHB, TROMA-III, 1:100), anti-Ki67 (Abcam, ab15580, 1:400), and anti-FOXO3A (Cell Signaling Technology, 2497, 1:200).
For chromogenic IHC, signal detection was performed using the VECTASTAIN Elite ABC-HRP kit (Vector Laboratories, PK-6100) according to the manufacturer’s instructions, followed by development with the DAB Substrate Kit (Vector Laboratories, SK-4100). Sections were counterstained with haematoxylin (Vector Laboratories, H-3401), dehydrated through graded ethanol, cleared in xylene and mounted with Permount mounting medium (Thermo Fisher Scientific, SP15-100). Bright-field images were acquired using an Aperio digital slide scanner (Leica Biosystems) using a 20× objective lens. For immunofluorescence, sections were washed and incubated with fluorescently conjugated secondary antibodies (1:1,000, Supplementary Table 2) at room temperature for 1 h, followed by nuclear counterstaining with DAPI. Some samples that were included in analyses of liver weight or whole-section tumour burden were not included in immunostaining quantification when extensive necrosis, tissue disruption, poor section quality or staining artefacts prevented reliable assessment of marker-positive cells or staining intensity.
For quantification of PHGDH IHC intensity, chromogenic PHGDH IHC images were analysed using ImageScope (Leica Biosystems), as previously described6. Tumour and hepatocyte regions were annotated manually, and staining intensity was quantified using the Positive Pixel Count v.9 algorithm. Pixels were classified as negative (0–100), weak (100–175) or strong (175–200) on the basis of intensity thresholds. For quantification of PHGDH immunofluorescence intensity, images were analysed using ZEN lite (ZEISS, v.3.8.99) and Fiji (ImageJ 1.54g; Java 1.8.0_345) to measure fluorescence intensities. Colocalization analysis is described in ‘Colocalization analysis’ below. For mouse samples, five to eight randomly selected tissue sections per group were analysed for IHC. For human samples, all available patient specimens were included in the IHC analysis (n = 32 for normal liver samples and n = 32 for liver-metastasis samples). For immunofluorescence-based correlation analyses, three randomly selected individuals per group were analysed, with three fields per individual.
PGP9.5 immunohistochemistry (automated staining)
For PGP9.5 staining, 5-μm-thick sections were prepared and immunostained on a Leica BondRX automated stainer, according to the manufacturer’s instructions. In brief, tissues underwent deparaffinization online, followed by epitope retrieval for 20 min at 100 °C with Leica Biosystems ER1 (pH 6) solution (Leica, AR9961) and endogenous peroxidase activity blocking with H2O2. Sections were then incubated with primary antibodies against PGP9.5 (Abcam, ab108986) at a 1:350 dilution for 15 min at ambient temperature. Primary antibodies were detected with the anti-rabbit HRP-conjugated polymer and 3,3′-diaminobenzidine (DAB) substrate that are provided in the Leica BOND Polymer Refine Detection System (DS9800). Finally, samples were counterstained with haematoxylin, dehydrated and coverslipped with Permount. Slides were scanned at 40× on a Hamamatsu Nanozoomer (2.0-HT) whole-slide scanner using NDP.scan (v.3.4.2). For quantification of PGP9.5 IHC, 4–13 randomly selected fields per mouse at 40× magnification were analysed, depending on the size of the metastatic lesions, to determine the PGP9.5-positive area. Regions containing necrosis or overt non-specific staining were excluded from the analysis. Quantification of DAB-positive areas was performed using Fiji (ImageJ). Colour deconvolution (H-DAB) was applied to separate the DAB signal, followed by background subtraction and fixed thresholding applied uniformly across all images. The percentage of DAB-positive area relative to the total image area was calculated automatically. For each mouse, values are presented as the mean of the analysed fields. Statistical analysis was performed by one-way ANOVA followed by Holm–Šídák’s post-hoc multiple-comparison tests.
Colocalization analysis
Fluorescence imaging was performed on a ZEISS spinning-disk confocal microscope. DAPI was excited at 353 nm and emission was collected at 465 nm; albumin was excited at 488 nm and emission was collected at 509 nm; and FOXO3A was excited at 587 nm and emission was collected at 610 nm. Acquisition settings were kept constant across experimental conditions to enable direct comparison of signal intensity and colocalization. Colocalization between FOXO3A and nuclear DAPI staining was quantified using the Coloc 2 plug-in in Fiji (ImageJ, v.2.1.0). To minimize diffuse background signal, manual background subtraction was applied uniformly across images by subtracting an intensity value of 112 from the DAPI channel and 107 from the FOXO3A channel in human samples. In mouse samples, 73 was subtracted from the DAPI channel and 80 was subtracted from the FOXO3A channel across both images to reduce background and improve statistical colocalization. Regions of interest (ROIs) were drawn over nuclei with a strong DAPI signal to ensure analysis was restricted to areas with clear nuclear staining. Within each ROI, Coloc 2 was run to calculate colocalization statistics, including Pearson’s correlation coefficient and Spearman’s rank correlation. Default Coloc 2 settings were used unless otherwise noted. For mouse liver-metastasis samples treated with the CXCR2 inhibitor, DAPI staining occasionally appeared more diffuse than it did in vehicle-treated controls, which can reduce the dynamic range of pixel-intensity-based colocalization metrics (particularly Pearson’s correlation). Pearson’s and Spearman’s coefficients were therefore interpreted in conjunction with a qualitative assessment of FOXO3A signal distribution. In these samples, the FOXO3A signal in hepatocytes showed a round, punctate pattern consistent with nuclear localization.
Quantification of metastatic burden in the liver
Metastatic burden in the liver was quantified using H&E-stained sections. For each mouse, H&E-stained sections from multiple liver lobes were scanned at 20× magnification using an Aperio digital slide scanner (Leica Biosystems). Whole-slide H&E images were analysed using MetFinder, a deep-learning-based tool for automated quantification of metastatic tumour burden in mouse organs (https://metfinder.org/)27. MetFinder automatically annotated each whole-slide image into normal liver and metastatic tumour regions. The resulting heat maps were carefully inspected to verify appropriate annotation, and metastatic burden was calculated as the percentage of metastatic tumour area relative to the total liver section area.
Cell culture
The human pancreatic cancer cell lines PaTu8902 and MiaPaCa2 and the mouse hepatocyte cell line AML12 were obtained from the American Type Culture Collection (ATCC) or DSMZ. PaTu8902 and MiaPaCa2 cells were cultured in DMEM (Corning) supplemented with 10% fetal bovine serum (FBS; Atlanta Biologicals, S11550H). AML12 cells were maintained in DMEM/F12 (Corning) supplemented with 10% FBS, 10 µg ml−1 insulin, 5.5 µg ml−1 transferrin, 5 ng ml−1 selenium and 40 ng ml−1 dexamethasone (Gibco). All media were supplemented with 1% penicillin–streptomycin (Gibco). The primary mouse PDAC cell line HY19636 was established from pancreatic tumours derived from KPC genetically engineered mouse models (LSL-KrasG12D/+; Trp53flox/+; p48-Cre), as described previously9. HY19636 cells were genetically modified to silence or overexpress genes of interest using CRISPR–Cas9-based approaches. All mouse PDAC cell lines were cultured in DMEM (Corning) supplemented with 10% FBS and 1% penicillin–streptomycin (Gibco). Primary mouse hepatocytes were isolated from C57BL/6 mice as described elsewhere in the Methods. Isolated primary mouse hepatocytes were plated on 0.01% collagen-coated culture plates and maintained for 6–10 days in Williams’ E medium (Corning) supplemented with 1% GlutaMAX (Gibco) and 1% penicillin–streptomycin (Gibco). All cells were cultured at 37 °C in a humidified incubator with 5% CO2. Cultures were routinely confirmed to be mycoplasma free using the MycoAlert Detection Kit (Lonza) or PCR-based assays. Cell lines were periodically authenticated by short tandem repeat profiling, and low-passage cultures were maintained in a centralized laboratory cell bank. For experiments using media containing different concentrations of serine and glycine, basal DMEM containing 2.78 mM (0.5 g l−1) glucose, 4 mM glutamine and 1 mM sodium pyruvate was supplemented as indicated with serine and glycine at 400 µM (high SG), 150 µM (low SG), or 0 µM (−SG), and with 10% dialysed FBS, unless otherwise noted.
Reagents
For pharmacological inhibition of CXCR2 and PI3K, SB225002 (Selleck, S7651) and BKM120 (buparlisib; Selleck, S2247) were used, respectively. SB225002 was prepared as a 20 mg ml−1 stock solution in DMSO, and BKM120 was prepared in DMSO. For in vitro experiments, SB225002 and BKM120 were used at final concentrations of 100 nM and 50 nM, respectively, unless otherwise indicated. Vehicle control samples received an equivalent volume of DMSO. Recombinant mouse CXCL5 (LIX) (rCXCL5; R&D Systems, 433-MC-025) and a CXCL5-neutralizing antibody (R&D Systems, MAB433-100) were used at final concentrations of 500 ng ml−1 and 2,000 ng ml−1, respectively.
Isolation of mouse hepatocytes
Primary mouse hepatocytes were isolated using a two-step collagenase digestion method as described previously28. Mice were anaesthetized, and a midline abdominal incision was made. After placement of a medium-sized ligation clip (Teleflex, 002200) on the inferior vena cava (IVC) immediately below the diaphragm and above the liver, the IVC was cannulated with a 20-gauge infusion catheter. The liver was first perfused with 20 ml pre-warmed perfusion buffer consisting of Hank’s balanced salt solution (HBSS) without Ca2+, Mg2+ or phenol red, supplemented with 10 mM HEPES and 0.5 mM EGTA (pH 7.4), at 37 °C and a flow rate of 7 ml min−1. The portal vein was severed to allow drainage of blood and perfusion buffer. Subsequently, the liver was perfused with 30 ml of digestion buffer composed of HBSS containing Ca2+, Mg2+ and phenol red, supplemented with 10 mM HEPES and 25 µg ml−1 Liberase TM (Sigma-Aldrich, 5401127001), at a flow rate of 5 ml min−1 for 6 min. After digestion, the liver was carefully excised and transferred to a 10-cm culture dish containing HBSS without Ca2+ and Mg2+. The liver capsule was gently torn, and hepatocytes were released into the buffer by gentle agitation. The resulting cell suspension was collected using a 25-ml serological pipette, filtered through a 70-µm cell strainer and centrifuged at 40g for 1 min at 4 °C without deceleration. The supernatant was aspirated, and the cell pellet was resuspended in HBSS. This low-speed centrifugation step was repeated several times to remove debris and non-parenchymal cells. The final pellet was resuspended in 10 ml Williams’ E medium (Corning) and mixed with an equal volume of Percoll solution, followed by centrifugation at 200g for 10 min at 4 °C. The resulting hepatocyte-enriched pellet was used for primary culture or downstream cell-sorting experiments.
Plasmids and generation of stable cell lines
To generate control (sgTOM), Phgdh-knockout or Cxcl5-knockout PDAC cell lines, the lentiviral CRISPR–Cas9 vector lentiCRISPRv2 (a gift from F. Zhang; Addgene plasmid 52961) was used. Single-guide RNAs (sgRNAs) targeting the indicated genes were cloned into the lentiCRISPRv2 backbone according to the standard protocol. The sgRNA target sequences were as follows:
sgTOM: 5′-GGCCACGAGTTCGAGATCGA-3′
sgPhgdh_1: 5′-TGAGCCCGGAATACGAGCAG-3′
sgPhgdh_2: 5′-AGGTGCTCCCTACCAAGCCG-3′
sgCxcl5_1: 5′-ATGGCGAGATGGAACCGCTG-3′
sgCxcl5_2: 5′-TTCCTCAGTCATAGCCGCAA-3′.
After lentiviral transduction, single-cell-derived clones were generated using a MoFlo XDP cell sorter (Beckman Coulter) and subsequently screened for efficient gene knockout. To minimize potential off-target effects and clonal variability, multiple independently derived clones with validated knockout were pooled and used for downstream experiments. For Phgdh or Cxcl5 knockout, clones generated using two independent sgRNAs were pooled before analysis. Rescue experiments were performed using the doxycycline-inducible lentiviral expression vector pCW57.1 carrying either wild-type or catalytically inactive human PHGDH (Addgene plasmids 154916 and 154903), which were gifts from R. Possemato and M. Pacold. Lentiviral transduction was performed according to standard protocols, and transduced cells were selected with puromycin (2 µg ml−1). For the generation of GFP-expressing sgTOM and mCherry-expressing sgPhgdh PDAC cell lines, MIGR1 (Addgene, 27490) and pMSCV-IRES-mCherry FP (Addgene, 52114; a gift from W. Pear) were used. AML12 cells with stable knockout of Cxcr2 or Foxo3a were generated using sgRNAs cloned into the pSpCas9(BB)-2A-GFP (PX458) vector (Addgene plasmid 48138). The sgRNA target sequences were as follows: sgCxcr2, 5′-TACGCAGTACGACCCTCAAA-3′; sgFoxo3a-sg1, 5′-GGACTGTCGTCTGCCGACTC-3′; and sgFoxo3a-sg2, 5′-TCTCGATGGCGCGGGTGATC-3′. For Foxo3a knockout, clones generated using the two independent sgRNAs were pooled before analysis. For ectopic expression of human FOXO3A, AML12 cells were transfected with pcDNA3-Flag-FOXO3A/FKHRL1 (wild type; Addgene plasmid 10708) or pcDNA3-Flag-FOXO3A/FKHRL1-AAA, a non-phosphorylatable mutant (Addgene plasmid 10709).
For transfection of retroviral plasmids, HEK293FT cells were co-transfected with the retroviral vector and the packaging plasmids pHit60 and VSVG at a 0.5:0.25:0.25 ratio. For transfection of lentiviral vectors, HEK293FT cells were co-transfected with the lentiviral vector and the packaging plasmids psPAX2 (Addgene plasmid 12260) and pMD2.G (Addgene plasmid 12259) at a 0.5:0.25:0.25 ratio. Lipofectamine 3000 (Thermo Fisher Scientific) was used as the transfection reagent according to the manufacturer’s instructions. Viral supernatants were collected 48 h and 72 h after transfection, filtered through a 0.45-µm filter and used for infection in the presence of polybrene (8 µg ml−1; EMD Millipore).
Mouse experiments
Intrasplenic and orthotopic injections of PDAC cells were performed as described previously29. In brief, mice were anaesthetized by intraperitoneal administration of ketamine and xylazine. After anaesthesia, a small incision was made in the top-left region of the abdomen, and either the pancreas or the spleen was gently exteriorized. For intrasplenic injection, PDAC cells were suspended in 100 µl PBS and loaded into an insulin syringe (28-gauge needle; BD, 329461) pre-filled with 100 µl PBS. The externalized spleen was divided using ligating clips (Teleflex, 002200), and cells were injected into the distal hemispleen. After injection, the splenic vein was ligated at the splenic hilum using ligating clips (Teleflex, 001200), and the injected hemispleen was surgically removed. For orthotopic pancreatic implantation, PDAC cells were resuspended in 20 µl of a 1:1 mixture of Matrigel (Corning, 356231) and HBSS and injected into the pancreatic tail using insulin syringes fitted with a 29-gauge needle (BD, 324702). Unless otherwise specified, 5 × 105 cells (for end-point assays) or 2.5 × 105 cells (for survival analysis) were used for intrasplenic injections, whereas 5 × 104 cells were injected for orthotopic pancreatic implantation.
For the orthotopic implantation followed by primary tumour resection (PTR) model (Fig. 5f–i), primary tumours were surgically resected 3 weeks after orthotopic implantation by distal pancreatectomy with splenectomy, essentially as described previously21. In brief, under anaesthesia, a sharp midline laparotomy was performed, with an optional left subcostal ‘hockey-stick’ extension as needed. The peritoneum was entered with electrocautery, and the pancreatic tail tumour was elevated without direct manipulation of the tumour mass by handling non-tumour pancreatic tissue (no-touch technique). Adhesions to surrounding structures were divided sharply or with electrocautery. Normal pancreas was identified and separated from gastric tissue, and two medium-sized ligating clips (Horizon, 002200) were applied proximal to the gross tumour edge to include draining vasculature. The pancreas was transected between the clips, and dissection was continued along the connective tissue plane between the tumour or spleen and the stomach until the superior splenic pole vasculature was encountered. Small-sized clips (Horizon 001200) were applied to the superior attachment between spleen and stomach, and the specimen was removed. The surgical field was confirmed to be haemostatic before closing the peritoneum and skin. Mice were euthanized at week 7 (4 weeks after PTR), and livers were collected for downstream analyses.
When the procedure was complete, the peritoneum was closed with 3-0 VICRYL violet sutures (Ethicon, J311H), and the skin incision was closed using the BD AutoClip wound closure system (BD). Mice were euthanized at the experimental end point, and livers or pancreatic tumours were collected for downstream analyses. For survival studies, mice were monitored daily and euthanized after reaching predefined humane end points in accordance with the Institutional Animal Care and Use Committee (IACUC) of New York University (NYU) Grossman School of Medicine.
For mouse experiments, C57BL/6J and NCr nude mice aged 8–10 weeks were used. Female mice were used in all experiments unless otherwise specified. Mice were fed either an AA control diet (ENVIGO, TD.01084) or a −SG diet (ENVIGO, TD.180296) for 15 days before tumour-cell implantation and were maintained on the same diet until the experimental end point. For pharmacological inhibition of CXCR2, mice were treated with the CXCR2 antagonist SB225002 (Selleck, S7651) or vehicle control by oral gavage at a dose of 1 mg kg−1. The vehicle consisted of 2% DMSO, 30% polyethylene glycol (PEG) and 5% Tween 20. Drug administration began 48 h after tumour-cell implantation and continued throughout the study period until euthanasia, with dosing performed every other day.
The experiments in this study were performed in compliance with the NYU IACUC, under protocols IA16-00507 and IA16-01331. All mouse experiments were performed in a specific-pathogen-free conventional animal facility. Mice were maintained in pre-sterilized, disposable irradiated cages supplied with irradiated chow and acidified drinking water. Microisolator cages were housed on ventilated rack systems. Sample sizes were determined empirically on the basis of preliminary experiments, and no formal statistical power calculations were performed. Investigators were not blinded to group allocation, because knowledge of treatment groups was required for proper execution of the experimental procedures. Tumour weight, an objective end-point measurement, was assessed only at the experimental end point after euthanasia and tumour collection. In accordance with the approved IACUC protocol, mice bearing palpable tumours were regularly examined, and tumour dimensions were assessed by palpation. Mice were euthanized if a tumour exceeded 2 cm in any dimension. No pancreatic tumours reached this limit. For the liver-metastasis models, the IACUC protocol did not specify a maximum liver size or weight; mice were instead monitored for general condition and body-weight loss and were euthanized if their body weight decreased by more than 20% from baseline. This body-weight limit was not exceeded in any of the liver-metastasis experiments.
In vivo hepatocyte-specific gene knockout using AAV-based CRISPR
In vivo hepatocyte-specific knockout of Phgdh or Cxcr2 was performed as described previously12 (Extended Data Fig. 2a). AAV vectors expressing sgRNAs were generated using the AAV2/8-U6-sgRNA-TBG-iCre-WPRE plasmid as the backbone. Plasmid construction and sgRNA cloning were performed by GenScript. The sgRNA target sequences were as follows:
sgScr (control): 5′-TACACGTCGCTAGGTTGCCC-3′
sgPhgdh: 5′-AGGTGCTCCCTACCAAGCCG-3′
sgCxcr2: 5′- TACGCAGTACGACCCTCAAA-3′.
Plasmids were submitted to the Viral Vector Core at the University of Iowa for AAV production. Rosa26-LSL-Cas9-EGFP mice on a C57BL/6J background were obtained from the Jackson Laboratory (strain 026175). Male and female homozygous Rosa26-LSL-Cas9-EGFP mice aged 6–8 weeks were used. Male mice were used only in the experiment shown in Extended Data Fig. 2e; female mice were used in all other experiments. Mice were injected retro-orbitally with AAV particles carrying either control sgRNA (sgScr) or sgRNAs targeting Phgdh or Cxcr2 at a dose of 5 × 1011 genome copies per mouse. One week after AAV administration, Cas9-expressing mice were injected intrasplenically with PDAC cells for downstream analyses.
For validation of knockout efficiency, hepatocytes were isolated from mice seven days after AAV injection, and EGFP-positive hepatocytes were sorted using a MoFlo cell sorter (Beckman Coulter). Sorted cells were subsequently subjected to western blot analysis.
RNA extraction
Total RNA was extracted from cells cultured in vitro or isolated by cell sorting using the PureLink RNA Mini Kit (Thermo Fisher Scientific, 12183025), according to the manufacturer’s instructions. For transwell co-culture experiments, total RNA was extracted from hepatocytes and cancer cells plated in the bottom chamber of a six-well Boyden co-culture plate (Greiner Bio-One, 657610) using the same kit.
Immunoblotting
Whole-cell protein lysates were prepared from hepatocytes or cancer cells cultured under monoculture or co-culture conditions using radioimmunoprecipitation assay (RIPA) buffer (Sigma-Aldrich, 20-188) supplemented with protease and phosphatase inhibitor cocktails (Roche). Protein lysates were separated on 4–20% gradient gels (Bio-Rad, 4561096) and transferred to polyvinylidene fluoride (PVDF) membranes (Merck Millipore, IPVH00010) using Tris–glycine transfer buffer containing 10% methanol. Membranes were blocked with 3% bovine serum albumin (BSA; Sigma-Aldrich, A2058) for at least 1 h at room temperature and incubated with primary antibodies overnight at 4 °C at the following dilutions: PHGDH (1:1,000; Sigma-Aldrich, HPA021241), AKT (1:1,000; Cell Signaling Technology (CST), 4691), phospho-AKT (1:1,000; CST, 9271), CXCL5 (LIX) (mouse) (1:300; R&D Systems, MAB433-100), CXCR2 (1:1,000; Invitrogen, PA1-31217), FOXO3A (1:1,000; CST, FOXO3A (75D8) rabbit monoclonal antibody, 2497S), phospho-FOXO3A (Ser315) (1:500; Proteintech, 28755-1-AP) and phospho-FOXO3A (Ser253) (1:500; Invitrogen, PA5-36816) and β-actin (1:1,000; Sigma-Aldrich, A5441). Membranes were washed at least three times (15 min each) with Tris-buffered saline containing 0.1% Tween 20 (TBST) and incubated with horseradish peroxidase (HRP)-conjugated secondary antibodies (anti-mouse IgG (1:5,000; CST, 7076S) or anti-rabbit IgG (1:5,000; CST, 7074S)) for 1 h at room temperature. Signals were detected using an enhanced chemiluminescence (ECL) detection system (Bio-Rad, 1705061), and images were acquired using a ChemiDoc imaging system (Image Lab Touch Software, BIO-RAD, v.2.3.0.07). When reliable signals could not be obtained after membrane stripping and reprobing, equal amounts of protein lysate were loaded onto separate gels and processed in parallel. These instances are indicated in the legends to Supplementary Fig. 1, which contains the uncropped gel and blot source images.
qPCR
Total RNA was reverse-transcribed using SuperScript VILO IV Reverse Transcriptase (Thermo Fisher Scientific) with oligo(dT) primers to generate complementary DNA (cDNA). qPCR was performed using SYBR Green Supermix (Bio-Rad) using the CFX Manager software (Bio-Rad, v.3.1.1517.0823). Relative mRNA expression levels were calculated using the ΔCt method and normalized to Actb mRNA expression. The sequences of qPCR primers were as follows: Phgdh forward, 5′-GACGTGAACTTGGTGAACGC-3′; Phgdh reverse, 5′-AGACAGCTCCGTTGAGCATC-3′; Psat1 forward, 5′-GCCTGAGACAGCGAACCAATG-3′; Psat1 reverse, 5′-CATGGTGCTAAGGCGACAGC-3′; Psph forward, 5′-AACTGGTTCTCCCGTCATCG-3′; Psph reverse, 5′-CTCTTAAAAGCGCCGAACCG-3′; Actb forward, 5′-CGATATCGCTGCGCTGGTC-3′; Actb reverse, 5′-CCACGATGGAGGGGAATACAG-3′; Cxcl5 forward, 5′-TAAAAGGGGTGCAGTGGGTT-3′; Cxcl5 reverse, 5′-GAGCACCAGCTCGGGATATG-3′.
All primers were specific for mouse genes.
Bulk RNA-seq
Bulk RNA-seq library preparation, FASTQ generation, read alignment to the reference genome using HISAT2 (v.2.2.1), read mapping and differential gene-expression analysis using DESeq2 (v.1.44.0) were performed by Novogene. Differentially regulated genes (log2-transformed fold change ≤ −2 and P < 0.05) were subjected to GSEA for pathway analysis. Enrichment scores and corresponding P values are reported in the figures.
GSEA
Differential expression results were converted into ranked gene lists for GSEA by selecting all genes with available log2(fold change) values and assigning each a ranking metric equal to the log2(fold change) from the DE analysis. The ranked list was saved as a two-column, tab-delimited text file with the first column containing gene symbols and the second containing the ranking metric, with no header line and unique ranking values.
The KEGG pathway ‘FoxO signalling pathway - Mus musculus (house mouse)’ (mmu04068) gene set was obtained from the KEGG database (https://www.genome.jp/dbget-bin/www_bget?path:mmu04068). The list of Mus musculus gene symbols was extracted from the ‘Genes’ section of the pathway entry and formatted into Gene Matrix Transposed (GMT) format, consisting of a single line with the pathway ID in the first column, a brief description in the second column and the gene symbols separated by tab characters in subsequent columns.
GSEA was performed in R (v.4.4.2) using the clusterProfiler package (v.4.14.6). Ranked gene lists were read into R and converted into named numeric vectors sorted in decreasing order of the ranking metric. The mmu04068 GMT file was read using read.gmt() and supplied as the TERM2GENE argument to the GSEA() function. All available genes were tested (pvalueCutoff = 1) to enable subsequent filtering. Enrichment plots were generated with the gseaplot2() function from the enrichplot package (v.1.26.6).
ORA
Raw gene-expression data were processed and differential expression analysis was performed using DESeq2. Genes were annotated with Entrez Gene IDs using the bitr function from the clusterProfiler R package and the org.Mm.eg.db mouse genome annotation database (v.3.20.0). Genes with Padj < 0.05 and absolute log2(fold change) > 1 were considered significantly differentially expressed.
To identify enriched biological pathways, ORA was performed using the enrichKEGG function from clusterProfiler, specifying the mouse organism (mmu). Resulting enriched pathways were filtered for significance at an FDR threshold of 5% ( Padj < 0.05). Pathway annotations were converted to human-readable gene symbols using setReadable. Visualization of the top enriched KEGG pathways was performed with the dotplot function.
Cell sorting
For isolation of hepatocytes from mouse livers with or without liver metastases, crude hepatocyte suspensions obtained as described elsewhere in the Methods were further purified by fluorescence-activated cell sorting using a MoFlo cell sorter (Beckman Coulter). Dead cells were excluded by DAPI staining. In in vivo hepatocyte-specific knockout experiments, EGFP-positive hepatocytes that are indicative of successful Cas9-mediated gene editing were isolated by cell sorting for downstream analyses. For the isolation of GFP- or mCherry-expressing PDAC cells from mouse liver metastases, metastatic liver nodules were excised, mechanically minced with scissors and enzymatically digested in DMEM containing 1 mg ml−1 collagenase IV (Gibco), 100 µg ml−1 DNase I (Roche), 1% FBS, 10 mM HEPES and 2% antibiotic–antimycotic (Thermo Fisher Scientific) for 40 min at 37 °C in the dark, with gentle agitation every 10 min. Digested tissues were washed twice with DMEM containing 10% FBS and filtered through a 40-µm nylon mesh strainer (Corning). Cell suspensions were treated with ACK lysis buffer (Thermo Fisher Scientific) for 10 min at 4 °C in the dark to remove red blood cells, followed by sorting of GFP- or mCherry-positive PDAC cells using a MoFlo cell sorter (Beckman Coulter). In all sorting experiments, dead cells were excluded by DAPI staining. Sorted cells were used for protein extraction, RNA extraction or in vitro culture for the collection of CM.
Flow cytometry
EdU incorporation in PDAC cells grown in direct co-culture with mouse hepatocytes was assessed using the Click-iT Plus EdU Alexa Fluor 647 Flow Cytometry Assay Kit (Thermo Fisher Scientific, C10634), according to the manufacturer’s instructions. In brief, mCherry-labelled Phgdh-knockout HY19636 cells were mixed with mouse hepatocytes at a 1:1 ratio (1,500 cells each per well) and seeded into 12-well plates. After overnight incubation in complete DMEM, cells were washed with PBS and cultured in basal DMEM, with or without SG supplementation, for 48 h. Cells were then treated with 10 μM EdU for 60 min, collected by trypsinization and fixed. EdU incorporation was detected following the manufacturer’s protocol. Flow cytometric analysis was performed using a BD LSR II UV flow cytometer, and EdU uptake was quantified specifically in mCherry-positive PDAC cells.
Transwell co-culture experiments
For transwell co-culture experiments, primary mouse hepatocytes or AML12 cells were seeded onto collagen-coated plates at a density of 12,000 cells per well for 6-well plates or 3,000 cells per well for 24-well plates. PDAC cells (15,000 cells per well) were seeded separately onto 1.0-µm pore-size transwell inserts (Greiner Bio-One, 657610). In selected experiments, the configuration was reversed, with PDAC cells seeded in the bottom wells and hepatocytes seeded onto the transwell inserts.
After overnight attachment, cells were gently washed with PBS, and the culture medium was replaced with basal DMEM containing 2.78 mM glucose, 4 mM glutamine and 1 mM sodium pyruvate, supplemented with 10% dialysed FBS and either 400 µM (high), 150 µM (low) or 0 µM (null) serine and glycine. Co-cultures were maintained for the indicated time periods and subsequently underwent protein extraction, RNA extraction or intracellular metabolite analysis, as described below. For glucose tracing experiments, co-cultures were processed as described in ‘Metabolomics’ below.
Metabolomics
For 13C6-glucose tracing experiments, PDAC cells and hepatocytes were plated either as monocultures or in transwell co-culture configurations, as described previously30,31,32. After two washes with PBS, cells were cultured in basal DMEM without glucose for 12 h and subsequently incubated in basal DMEM supplemented with 25 mM 13C6-glucose (Cambridge Isotope Laboratories, CLM-1396) for an additional 24 h.
To remove residual medium and contaminants, cells were rinsed with 0.9% NaCl prepared in high-performance-liquid-chromatography-grade water. Metabolites were extracted using 80% methanol containing 1 µg norvaline as an internal standard. AA standards (Cambridge Isotope Laboratories, MSK-A2-1.2) were included for metabolite identification and quantification. Cell extracts were vortexed for 15 min at 4 °C and centrifuged at maximum speed for 10 min to separate insoluble material. Polar metabolites were collected from the aqueous phase, and 300 µl of each sample was transferred into polypropylene vials (Agilent Technologies, 5190-2243) and dried under vacuum using a SpeedVac concentrator (Thermo Fisher Scientific Savant, SPD111V). Dried metabolites were derivatized by incubation with 20 µl methoxyamine hydrochloride (20 mg ml−1 in pyridine; freshly prepared) for 60 min at 37 °C, followed by incubation with 20 µl MTBSTFA containing 1% tert-butyldimethylchlorosilane (t-BDMSC) for 30 min at 37 °C. Samples were analysed with a 7890B gas chromatograph coupled to a 5977B mass spectrometer (Agilent Technologies) equipped with a DB-35ms Ultra Inert column (Agilent Technologies, 122-3832UI) using Agilent MassHunter GC/MS Acquisition (v.B.07.05.2479). The GC–MS parameters, quantification and correction for natural isotope abundances were performed as described previously31,33,34.
Growth assays
PDAC cells were seeded either as monocultures or in transwell co-culture with hepatocytes in 24-well plates as described in the previous section. Cells were cultured for 5 days or for the indicated time periods in medium containing high or low SG concentrations. At the end of the culture period, PDAC cells were washed with PBS, fixed with ice-cold methanol (pre-chilled at −20 °C) and stained with crystal violet solution (0.5% w/v in methanol). After incubation for 20 min at room temperature with gentle rocking, plates were washed thoroughly under running tap water and air-dried for 2 h or overnight. Crystal violet staining was then solubilized with 300 µl methanol per well and transferred to 96-well plates, and absorbance was measured at 570 nm using a SpectraMax M5 plate reader (Agilent) with SoftMax Pro (Molecular Devices, v.7.0).
Chromatin immunoprecipitation
Chromatin immunoprecipitation (ChIP) was performed using the SimpleChIP Enzymatic Chromatin IP Kit (Magnetic Beads) (CST, 9003) according to the manufacturer’s instructions, with minor modifications. AML12 hepatocytes were cultured alone or co-cultured with PDAC cells (Phgdh-KO HY19636 cells) using transwell assays (see ‘Transwell co-culture experiments’ for details) and treated as indicated. Experimental conditions included AML12 monoculture, AML12–PDAC co-culture and co-culture in the presence of CXCR2 or PI3K inhibitors. After treatment, cells were cross-linked, nuclei were isolated and chromatin was digested as per the kit protocol. After resuspension, the samples were sonicated in Diagenode Bioruptor for 15 cycles with 60 s on and 60 s off at 4 °C. For each immunoprecipitation, 5–10 μg of digested chromatin was incubated overnight at 4 °C with rotation with an anti-FOXO3 antibody (1:200, CST, 2497). Normal IgG supplied with the kit was used as a negative control. A 2% aliquot of chromatin was reserved as input control. Immunocomplexes were captured and washed sequentially with low- and high-salt buffers according to the manufacturer’s instructions. Chromatin was eluted, cross-links were reversed and DNA was purified using the spin columns provided with the kit. Purified DNA was analysed by qPCR using primers targeting putative FOXO-binding regions (TGTTT-core FOXO consensus-like motifs) in the Phgdh promoter (P1 and P2) and a distal gene-body region (P3; negative control)18 (Fig. 4h). Primer sequences: P1, Fw CTGGGCCAGAGAAGGGAAAG, Rv GCCTGTGCTGTTACCTCCAT; P2, Fw TGAGATTTAATTCCCTCGTGGAG, Rv GTGAGCTTTAACACGCACGAT; and P3, Fw ACAGTAAGGCGCTCAGTCAC, Rv TGGCTGGATTCAGTAACGGC). ChIP–qPCR signals were normalized to the 2% input sample and enrichment was expressed as the FOXO3 immunoprecipitation relative to IgG control (FOXO3/IgG).
Dual-reporter analysis
PHGDH promoter activity was assessed using the Secrete-Pair Gaussia luciferase/secreted alkaline phosphatase (GLuc/SEAP) dual luminescence assay system (GeneCopoeia, LF031) together with a GLuc reporter construct containing the mouse Phgdh promoter (GeneCopoeia, MPRM60654). Mouse hepatocyte AML12 cells, in which endogenous Foxo3a was knocked out and reconstituted with human wild-type FOXO3A or a FOXO3A(AAA) mutant, were transfected with the GLuc/SEAP reporter vectors using Lipofectamine 3000 transfection reagent (Thermo Fisher Scientific) according to the manufacturer’s instructions. Forty-eight hours after transfection, cells were co-cultured with PDAC cells for an additional 48 h in 12-well plates under conditions described elsewhere. CM (200 µl per well) was collected and transferred to 96-well plates for luminescence measurements. GLuc and SEAP activities were measured using a SpectraMax M5 plate reader (Agilent) with SoftMax Pro (Molecular Devices, v.7.0) following the manufacturer’s instructions. Each experimental condition was analysed using three biological replicates per group. GLuc luminescence values were normalized to the corresponding SEAP signals to control for transfection efficiency, and results were expressed as fold change relative to the wild-type group. The experiment was independently repeated twice.
D-serine and L-serine assay
For Fig. 3e, concentrations of D-serine and L-serine were measured using a DL-Serine Assay kit (Abcam, 241027) according to the manufacturer’s instructions. In brief, PDAC cells cultured under monoculture or co-culture conditions were collected and centrifuged, followed by pretreatment with the sample clean-up mix provided in the kit. Samples were then deproteinized and clarified by centrifugation through 3-kDa molecular weight cut-off filters (Amicon Ultra-0.5 Centrifugal Filter; MilliporeSigma, UFC5003). Processed samples, together with d-serine standards supplied with the kit, were subjected to fluorescence-based quantification using a SpectraMax M5 plate reader (Agilent) with SoftMax Pro (Molecular Devices, v.7.0) (excitation/emission = 535/587 nm). Concentrations of d-serine and l-serine were calculated according to the manufacturer’s protocol.
Isolation of TIF from mouse livers
Isolation of TIF was adapted from previously established protocols35,36,37 with minor modifications. Freshly isolated mouse livers were cut into small pieces (1–3 mm3) in ice-cold PBS supplemented with a protease inhibitor cocktail. Tissue fragments were carefully transferred into 15-ml tubes using pre-chilled glass pipettes, washed briefly with PBS and resuspended in an equal volume of PBS relative to tissue volume. Samples were incubated at 37 °C in a CO2 incubator for 1 h with occasional gentle tapping, followed by centrifugation at 1,000g for 8 min at 4 °C. The supernatant was transferred to a pre-chilled microcentrifuge tube using a pre-chilled glass pipette and centrifuged at 2,000g for 15 min at 4 °C. The resulting supernatant was further clarified by centrifugation at 20,000g for 30 min at 4 °C. The clarified supernatant was concentrated using a vacuum concentrator for 1 h to obtain TIF at a final protein concentration of approximately 2–5 mg ml−1 per liver from C57BL/6 mice. TIF samples were isolated from three to five mice per group. Sample purity was assessed using a lactate dehydrogenase (LDH) assay kit (Abcam, ab102526) according to the manufacturer’s instructions. Samples exhibiting less than 10% cytotoxicity were snap-frozen, stored at −80 °C and used for subsequent analyses.
Enzyme-linked immunosorbent assay
The concentrations of CXCL5 and NGF in CM from PDAC cells, TIF derived from mouse tissues and tissue homogenates of human tissues were measured using ELISA kits for human CXCL5 (R&D Systems, DY254), mouse CXCL5 (R&D Systems, MX000) and mouse NGF (LSBio, LS-F5156), according to the manufacturers’ instructions. In brief, CM from mouse PDAC cells cultured with or without hepatocyte co-culture was collected 24 h after medium replacement and centrifuged at low speed to remove cellular debris. The resulting clarified supernatants were then subjected to enzyme-linked immunosorbent assay (ELISA). The preparation of TIF from mouse liver tissue and tissue homogenates from frozen human liver tissue is described elsewhere in the Methods.
For CXCL5 ELISA using homogenates of snap-frozen liver-metastasis tissues from autopsy samples (Fig. 3c), PHGDH status was assigned for each ELISA sample based on PHGDH immunostaining of corresponding metastatic lesions in available sections, and samples were grouped as PHGDH-strong versus PHGDH-weak/negative on the basis of PHGDH staining intensity in metastatic tumour cells. Because the ELISA samples were prepared from frozen tissue pieces, the set of samples used for ELISA did not necessarily match one-to-one with the FFPE specimens used for quantitative PHGDH evaluation in both tumour cells and hepatocytes (Fig. 1a,b). In some cases, although PHGDH staining was performed to categorize the frozen-tissue ELISA samples, the corresponding FFPE sections lacked sufficient hepatocytes to score hepatocyte PHGDH, and those cases were therefore excluded from the analyses in Fig. 1a,b.
Statistical analysis and reproducibility
Statistical analyses were performed using GraphPad Prism (v.9.0 and v.10.5.0). For comparisons between two groups, two-tailed unpaired Welch’s t-tests were used. For experiments involving more than two groups, one-way ANOVA followed by Tukey’s multiple-comparisons or Holm–Šídák’s post-hoc multiple-comparison tests was applied. P < 0.05 was considered statistically significant. For stacked isotopologue distribution plots, statistical significance was assessed for each isotopologue by two-way ANOVA followed by Tukey’s multiple-comparisons test. Exact P values are provided in the corresponding source data files. Data are presented as mean ± s.d. or as individual data points, as indicated in the figure legends. For bar graphs, each symbol represents one biologically independent sample unless otherwise stated. Immunoblot experiments were independently repeated at least twice with similar results. For in vivo studies, each intrasplenic or intrapancreatic injection cohort represents one independent experiment from multiple repeats; the number of biologically independent mice per group (n) is indicated in the corresponding figure legends, and each symbol represents one mouse.
Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
Data availability
RNA-seq data have been deposited in the Gene Expression Omnibus (GEO) data repository with accession numbers GSE324743 and GSE324997. Mass-spectrometry data have been deposited in the MassIVE repository under accession number MSV000102778. All other data supporting the findings of this study are available in the Article and its Supplementary Information. Source data are provided with this paper.
Code availability
No custom software, original algorithms or bespoke computational methods were developed for this study. RNA-seq analyses were performed using established, publicly available software and Bioconductor packages, as described in the Methods.
References
Inzunza, J. & Del Valle, A. C. Deciphering the liver’s role in pancreatic cancer metastasis: pathways and therapeutic approaches. npj Precis. Oncol. 9, 395 (2025).
Article CAS PubMed PubMed Central Google Scholar
Link, J. M. et al. Ongoing replication stress tolerance and clonal T cell responses distinguish liver and lung recurrence and outcomes in pancreatic cancer. Nat. Cancer 6, 123–144 (2025).
Article CAS PubMed PubMed Central Google Scholar
Ganguly, K. & Kimmelman, A. C. Reprogramming of tissue metabolism during cancer metastasis. Trends Cancer 9, 461–471 (2023).
Article CAS PubMed PubMed Central Google Scholar
Krieg, S., Fernandes, S. I., Kolliopoulos, C., Liu, M. & Fendt, S. M. Metabolic signaling in cancer metastasis. Cancer Discov. 14, 934–952 (2024).
Article CAS PubMed PubMed Central Google Scholar
Demicco, M., Liu, X. Z., Leithner, K. & Fendt, S. M. Metabolic heterogeneity in cancer. Nat. Metab. 6, 18–38 (2024).
Article PubMed Google Scholar
Banh, R. S. et al. Neurons release serine to support mRNA translation in pancreatic cancer. Cell 183, 1202–1218 (2020).
Article CAS PubMed PubMed Central Google Scholar
Zhang, S. et al. Single cell transcriptomic analyses implicate an immunosuppressive tumor microenvironment in pancreatic cancer liver metastasis. Nat. Commun. 14, 5123 (2023).
Article ADS CAS PubMed PubMed Central Google Scholar
Law, H. C. et al. The proteomic landscape of pancreatic ductal adenocarcinoma liver metastases identifies molecular subtypes and associations with clinical response. Clin. Cancer Res. 26, 1065–1076 (2020).
Article CAS PubMed Google Scholar
Bardeesy, N. et al. Both p16Ink4a and the p19Arf-p53 pathway constrain progression of pancreatic adenocarcinoma in the mouse. Proc. Natl Acad. Sci. USA 103, 5947–5952 (2006).
Article ADS CAS PubMed PubMed Central Google Scholar
Maddocks, O. D. K. et al. Modulating the therapeutic response of tumours to dietary serine and glycine starvation. Nature 544, 372–376 (2017).
Article ADS CAS PubMed Google Scholar
Platt, R. J. et al. CRISPR–Cas9 knockin mice for genome editing and cancer modeling. Cell 159, 440–455 (2014).
Article CAS PubMed PubMed Central Google Scholar
Katsuda, T. et al. Rapid in vivo multiplexed editing (RIME) of the adult mouse liver. Hepatology 78, 486–502 (2023).
Article PubMed Google Scholar
Steele, C. W. et al. CXCR2 inhibition profoundly suppresses metastases and augments immunotherapy in pancreatic ductal adenocarcinoma. Cancer Cell 29, 832–845 (2016).
Article CAS PubMed PubMed Central Google Scholar
Huang, H. & Tindall, D. J. Regulation of FOXO protein stability via ubiquitination and proteasome degradation. Biochim. Biophys. Acta 1813, 1961–1964 (2011).
Article CAS PubMed PubMed Central Google Scholar
Santo, E. E. et al. FOXO3a is a major target of inactivation by PI3K/AKT signaling in aggressive neuroblastoma. Cancer Res. 73, 2189–2198 (2013).
Article CAS PubMed Google Scholar
Ferber, E. C. et al. FOXO3a regulates reactive oxygen metabolism by inhibiting mitochondrial gene expression. Cell Death Differ. 19, 968–979 (2012).
Article CAS PubMed Google Scholar
Litvak, V. et al. A FOXO3–IRF7 gene regulatory circuit limits inflammatory sequelae of antiviral responses. Nature 490, 421–425 (2012).
Article ADS CAS PubMed PubMed Central Google Scholar
Tao, R., Xiong, X., DePinho, R. A., Deng, C. X. & Dong, X. C. FoxO3 transcription factor and Sirt6 deacetylase regulate low density lipoprotein (LDL)-cholesterol homeostasis via control of the proprotein convertase subtilisin/kexin type 9 (Pcsk9) gene expression. J. Biol. Chem. 288, 29252–29259 (2013).
Article CAS PubMed PubMed Central Google Scholar
Brunet, A. et al. Akt promotes cell survival by phosphorylating and inhibiting a Forkhead transcription factor. Cell 96, 857–868 (1999).
Article ADS CAS PubMed Google Scholar
Tong, H. et al. Dual impacts of serine/glycine-free diet in enhancing antitumor immunity and promoting evasion via PD-L1 lactylation. Cell Metab. 36, 2493–2510 (2024).
Article CAS PubMed Google Scholar
Sorrentino, A. et al. Surgical resection of murine PDAC alters hepatic metastases and immune microenvironment. Pancreas 54, e831–e841 (2025).
Article CAS Google Scholar
Scott, A. J. et al. Rewiring of cortical glucose metabolism fuels human brain cancer growth. Nature 646, 413–422 (2025).
Article ADS CAS PubMed PubMed Central Google Scholar
Kamphorst, J. J. et al. Human pancreatic cancer tumors are nutrient poor and tumor cells actively scavenge extracellular protein. Cancer Res. 75, 544–553 (2015).
Article CAS PubMed PubMed Central Google Scholar
Gheller, B. J. et al. Extracellular serine and glycine are required for mouse and human skeletal muscle stem and progenitor cell function. Mol. Metab. 43, 101106 (2021).
Article CAS PubMed Google Scholar
Lim, E. W. et al. Serine and glycine physiology reversibly modulate retinal and peripheral nerve function. Cell Metab. 36, 2315–2328 (2024).
Article CAS PubMed Google Scholar
Tajan, M. et al. Serine synthesis pathway inhibition cooperates with dietary serine and glycine limitation for cancer therapy. Nat. Commun. 12, 366 (2021).
Article ADS CAS PubMed PubMed Central Google Scholar
Karz, A. et al. MetFinder: a tool for automated quantitation of metastatic burden in histological sections from preclinical models. Pigment Cell Melanoma Res. 38, e13195 (2025).
Article CAS PubMed Google Scholar
Charni-Natan, M. & Goldstein, I. Protocol for primary mouse hepatocyte isolation. STAR Protoc. 1, 100086 (2020).
Article PubMed PubMed Central Google Scholar
Yamamoto, K. et al. Autophagy promotes immune evasion of pancreatic cancer by degrading MHC-I. Nature 581, 100–105 (2020).
Article ADS CAS PubMed PubMed Central Google Scholar
Sousa, C. M. et al. Pancreatic stellate cells support tumour metabolism through autophagic alanine secretion. Nature 536, 479–483 (2016).
Article ADS CAS PubMed PubMed Central Google Scholar
Parker, S. J. et al. Selective alanine transporter utilization creates a targetable metabolic niche in pancreatic cancer. Cancer Discov. 10, 1018–1037 (2020).
Article CAS PubMed PubMed Central Google Scholar
Encarnación-Rosado, J. et al. Targeting pancreatic cancer metabolic dependencies through glutamine antagonism. Nat. Cancer 5, 85–99 (2024).
Article PubMed Google Scholar
Hollinshead, K. E. R. et al. Respiratory supercomplexes promote mitochondrial efficiency and growth in severely hypoxic pancreatic cancer. Cell Rep. 33, 108231 (2020).
Article CAS PubMed PubMed Central Google Scholar
Parker, S. J. et al. Spontaneous hydrolysis and spurious metabolic properties of α-ketoglutarate esters. Nat. Commun. 12, 4905 (2021).
Article ADS CAS PubMed PubMed Central Google Scholar
Celis, J. E. et al. Proteomic characterization of the interstitial fluid perfusing the breast tumor microenvironment. Mol. Cell. Proteomics 3, 327–344 (2004).
Article CAS PubMed Google Scholar
Sun, W. et al. Characterization of the liver tissue interstitial fluid (TIF) proteome indicates potential for application in liver disease biomarker discovery. J. Proteome Res. 9, 1020–1031 (2010).
Article ADS CAS PubMed Google Scholar
Sun, W., Jiang, Y. & He, F. in Serum/Plasma Proteomics: Methods in Molecular Biology Vol. 728 (eds Simpson, R. J. & Greening, D. W.) 247–257 (2011).
Download references
Acknowledgements
We thank H. Ying for providing the HY19636 cells; the NYU Langone Health Experimental Pathology Laboratory and Flow Cytometry Core, each supported in part by a Cancer Center Support Grant (P30CA016087) at the Laura and Isaac Perlmutter Cancer Center; Novogene for RNA-seq and subsequent bioinformatic analysis; S. Kurosaki and S. Kawamura for assistance with isolating and culturing primary hepatocytes; and V. Costiniti and R. S. Banh for assistance with metabolomics experiments.
Funding
This work was supported by National Cancer Institute (NCI) grants P01CA117969, R35CA232124, P30CA016087-38 and 1R01CA251726-01A1 (A.C.K.); the Lustgarten Foundation (A.C.K.); SU2C (A.C.K.); the Yasuda Medical Foundation (K.Y.); the Tokyo Society of Medical Sciences (K.Y.); the MSD Life Science Foundation (K.Y.); the Astellas Foundation for Research on Metabolic Disorders (K.Y.); the Uehara Memorial Foundation (K.Y.); the Princess Takamatsu Cancer Research Fund (K.Y.); the Cell Science Research Foundation (K.Y.); the Chugai Foundation for Innovative Drug Discovery Science (C-FINDs) (K.Y.); the Takeda Science Foundation (K.Y.); the Foundation for Promotion of Cancer Research (K.Y.); JSPS KAKENHI grants 21K20828, 22H02898, 23K06647, 24K02302, 25K02642 and 25K02515 (K.Y.); NIH grant K00CA245822-06 (D.E.B.); the Howard Hughes Medical Institute Gilliam Fellowship for Advanced Study (J.E.-R.); JST Moonshot Research and Development Program grant JPMJMS2214-10 (M.F.); Kobayashi Foundation for Cancer Research (K.Y.) and the Naito Foundation (K.Y.). The UNMC Rapid Autopsy Program for Pancreas was supported by the SPORE in Pancreatic Cancer (P50CA127297), the Pancreatic Cancer Detection Consortium (U01CA210240), a NCI Cancer Center Support Grant (P30CA36727), a NCI Research Specialist Award (R50CA211462), U01CA284086 and U54CA274329.
Ethics declarations
Competing interests
A.C.K. is an inventor on patents pertaining to KRAS-regulated metabolic pathways; redox control pathways in pancreatic cancer; targeting GOT1 as a therapeutic approach; alanine import inhibition; and the autophagic control of iron metabolism. A.C.K. is on the scientific advisory board of Cornerstone Pharmaceuticals (formerly Rafael Pharmaceuticals). E.Z.W. completed this work during his postdoctoral fellowship at NYU Grossman School of Medicine and is currently an employee of Merck Sharp & Dohme, a subsidiary of Merck; the opinions and perspectives expressed herein are those of the author (E.Z.W.) and do not represent the opinions or perspectives of his current employer. The remaining authors declare no competing interests.
Peer review
Peer review information
Nature thanks the anonymous reviewers for their contribution to the peer review of this work.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Fig. 1 Loss of PHGDH in PDAC cells induces the upregulation of PHGDH in hepatocytes.
a, Representative immunofluorescence of albumin (green), PHGDH (red) and DAPI (blue) in liver metastases generated in nude mice by intrasplenic injection of human PDAC cell lines that are PHGDH-high (MiaPaCa2) or PHGDH-null (PaTu8902). Dashed lines indicate tumour (T)–hepatocyte (H) boundaries. Scale bar, 20 μm. Similar results were observed in five biologically independent mice per group. b, Immunoblot validation of CRISPR-mediated Phgdh knockout in mouse HY19636 PDAC cells (sgTOM, Phgdh WT; sgPhgdh, Phgdh KO) cultured in serine and glycine-replete (+SG) or serine- and glycine-depleted (−SG) medium. β-actin was used as a loading control. c, Proliferation of sgTOM or sgPhgdh HY19636 cells cultured in +SG or −SG medium, quantified by crystal violet staining and shown as fold change relative to day 0. The exact P value is shown for sgPhgdh +SG versus sgPhgdh −SG at day 5 (two-sided unpaired Student’s t-test). Data are mean ± s.d. (n = 3 biological replicates). d, Representative H&E and PHGDH IHC of normal mouse liver and livers bearing metastases derived from HY19636 sgTOM or sgPhgdh cells (n = 5, 6 and 6 mice, respectively). Boxes indicate tumour (T, red) and hepatocyte (H, blue) regions; dashed lines indicate T–H boundaries. e, Quantification of hepatocyte PHGDH staining intensity in d. The proportion of samples with strong hepatocyte PHGDH staining (strong vs negative/weak) was compared between normal livers and livers bearing sgPhgdh metastases (two-sided Fisher’s exact test). Numbers in bars indicate the number of mice. f, Representative immunofluorescence of albumin (green) and PHGDH (red) in mouse PDAC liver metastases from d. Dashed lines indicate T–H boundaries. Scale bar, 20 μm. g, Phgdh mRNA expression in hepatocytes isolated from normal livers or livers bearing sgTOM or sgPhgdh HY19636 metastases (n = 2 mice per group), quantified by qPCR.
Source data
Extended Data Fig. 2 Hepatocyte-specific knockout of Phgdh suppresses the metastatic growth of exSer-dependent PDAC cells.
a, Schematic of hepatocyte-specific in vivo CRISPR editing in Rosa-LSL-Cas9-EGFP mice using AAV2/8-U6-sgRNA-TBG-Cre (AAV-sgRNA). The TBG promoter drives hepatocyte-specific Cre expression, enabling Cas9 activation and EGFP labelling in hepatocytes. EGFP+ hepatocytes were isolated for validation of target-gene disruption. b, Immunoblot showing efficient PHGDH ablation in EGFP+ hepatocytes isolated from mice injected with AAV-sgScr or AAV-sgPhgdh. β-actin was used as a loading control. c, Plasma biochemistry in Cas9 mice maintained on a −SG diet after injection with AAV-sgScr or AAV-sgPhgdh (n = 4 mice per group). Each symbol represents an individual mouse; bars indicate mean ± s.e.m. P values were determined by two-sided unpaired Welch’s t-test. d–f, Whole liver weight and representative gross liver images from Cas9 mice injected with AAV-sgScr or AAV-sgPhgdh and subsequently inoculated intrasplenically with sgTOM or sgPhgdh HY19636 cells. d, Female mice on −SG diet (n = 6, 7, 7 and 8 mice, respectively). e, Male mice on −SG diet (n = 6, 5, 7 and 7 mice, respectively). f, Female mice on control AA diet (n = 3, 3, 3 and 5 mice, respectively). Each symbol represents an individual mouse; bars indicate mean ± s.e.m. Statistical significance was assessed using one-way ANOVA with Holm–Šidák’s multiple-comparisons test. g, Representative PHGDH IHC images from livers in Fig. 1e. h, Hepatocyte PHGDH staining intensity in g, scored as negative/weak versus strong (cut-offs as indicated). Numbers in bars indicate the number of mice. Exact P values are shown in the relevant panels.
Source data
Extended Data Fig. 3 Hepatocyte-derived serine supports exSer-dependent PDAC cells under serine and glycine deprivation.
a, Immunoblot validation of PHGDH expression in HY19636 PDAC cells after CRISPR-mediated Phgdh knockout (sgPhgdh) and reconstitution with human PHGDH(WT) or catalytically inactive PHGDH (PHGDH(CD)). β-actin was used as a loading control. b, Crystal violet–based viability assay (absorbance at 590 nm) of sgPhgdh HY19636 cells reconstituted with PHGDH(WT) (left) or PHGDH(CD) (right), cultured under +SG or −SG conditions, with or without transwell co-culture with primary mouse hepatocytes. n = 3 biological replicates. Bars indicate mean ± s.d. Statistical significance was assessed using one-way ANOVA with Tukey’s multiple-comparisons test. c, Representative immunofluorescence images of liver metastases formed by sgPhgdh PDAC cells in mice that were fed a −SG diet stained for CK19 (green), Ki67 (red) and DAPI (blue), showing a field near the tumour–liver interface (periphery) and an adjacent field toward the tumour core. Asterisks mark the same landmark in adjacent fields; dashed lines indicate tumour (T)–hepatocyte (H) boundaries. d, Quantification of proliferating tumour cells in c, shown as the percentage of Ki67+ cells among CK19+ tumour cells across concentric regions from the tumour periphery toward the core (regions 1–3). Bars indicate mean ± s.e.m. Statistical significance was assessed using one-way ANOVA with Tukey’s multiple-comparisons test. Each symbol represents one microscopic field; two fields were analysed per mouse from four mice (n = 8 fields per region). e, 13C6-glucose tracing showing fractional isotopologue distribution of intracellular serine in sgPhgdh HY19636 cells reconstituted with PHGDH(WT) or PHGDH(CD), cultured with or without hepatocytes under −SG conditions, measured by GC–MS. Stacked bars indicate mean ± s.d. from n = 3 biological replicates. f, Time course of 13C6-glucose incorporation into intracellular serine (M + 0 to M + 3) in sgTOM or sgPhgdh HY19636 PDAC cells and AML12 hepatocytes, measured by GC–MS. Stacked bars indicate mean ± s.d. from n = 3 biological replicates. g, GSEA of RNA-seq data from hepatocytes cultured alone or co-cultured with sgPhgdh HY19636 cells under −SG conditions, showing a trend toward enrichment of glycine–serine–threonine metabolism. NES, normalized enrichment score. h, Crystal violet–based viability assay testing the serine concentration required to support growth of sgPhgdh HY19636 PDAC cells. sgTOM and sgPhgdh HY19636 cells were cultured in serine- and glycine-free medium supplemented with the indicated concentrations of serine, and viability was quantified by absorbance at 590 nm after crystal violet staining. n = 6 biological replicates. Bars indicate mean ± s.d.; statistical significance was assessed using one-way ANOVA with Tukey’s multiple-comparisons test. Exact P values are shown in the relevant panels.
Source data
Extended Data Fig. 4 exSer-dependent PDAC upregulates CXCL5 in the presence of hepatocytes, which is required to induce PHGDH expression in hepatocytes.
a, CXCL5 concentrations in culture media from sgTOM or sgPhgdh HY19636 cells co-cultured with hepatocytes under the indicated serine and glycine conditions, quantified by ELISA. n = 2 biological replicates. b, Cxcl5 mRNA expression in PDAC cells isolated from metastasis-bearing mouse livers and cultured in vitro for 2 weeks in +SG or −SG medium, quantified by qPCR and expressed relative to sgTOM cells cultured in −SG medium. n = 2 biological replicates. c, CXCL5 concentrations in culture media from PDAC cells in b, quantified by ELISA. n = 2 biological replicates. d, NGF concentrations in culture media from sgTOM or sgPhgdh HY19636 cells cultured alone or co-cultured with hepatocytes under the indicated serine and glycine conditions, quantified by ELISA. n = 6 biological replicates. e, NGF concentrations in TIFs from mouse livers bearing sgTOM or sgPhgdh HY19636 liver metastases under +SG or −SG diet conditions (related to Fig. 3b), quantified by ELISA. n = 3 mice per group. f, PGP9.5 immunohistochemistry (IHC) as a readout of tumour innervation. Left, positive-control staining in brain, skin and periportal liver. Middle, representative PGP9.5 IHC in sgTOM or sgPhgdh HY19636 liver metastases under +SG or −SG diet conditions. Right, quantification of PGP9.5+ area. Each symbol represents the mean of analysed fields from one mouse. n = 3 mice per group. g, Cxcl5 mRNA expression in orthotopic pancreatic tumours formed by sgPhgdh HY19636 cells in mice fed a control or −SG diet, quantified by qPCR and expressed relative to tumours from control diet–fed mice (set to 1). n = 5 mice per group. h, Immunoblot validation of CXCL5 knockout in sgTOM and sgPhgdh HY19636 cells under reducing and non-reducing conditions. i, Immunoblot of PHGDH in primary mouse hepatocytes cultured alone or co-cultured with HY19636 cells, with or without rCXCL5. j, Phgdh mRNA expression in primary mouse hepatocytes cultured under −SG conditions and treated as indicated, quantified by qPCR and normalized to monocultured hepatocytes (dashed line). n = 2 biological replicates. k, Heat map of differentially expressed genes in hepatocytes co-cultured with sgPhgdh HY19636 cells and treated with DMSO or a CXCR2 inhibitor. The inset highlights SSP genes (Phgdh, Psat1 and Psph). l, Representative immunofluorescence images of CK19 (green), Ki67 (red) and DAPI (blue) in mouse liver metastases formed by Phgdh-KO single-knockout (SKO) or Phgdh/Cxcl5 double-knockout (DKO) HY19636 cells. m, Quantification of proliferating tumour cells in l, shown as the percentage of Ki67+ cells among CK19+ tumour cells (n = 6 mice per group). n, (Related to Fig. 3j,k) Immunoblot showing efficient CXCR2 ablation in EGFP+ hepatocytes isolated from mice injected with AAV-sgScr or AAV-sgCxcr2 and subsequently subjected to intrasplenic injection of sgTOM or sgPhgdh HY19636 cells. β-actin was used as a loading control. For d–g,m, bars indicate mean ± s.e.m. Statistical significance was assessed by one-way ANOVA with Tukey’s multiple-comparisons test for d,e, one-way ANOVA with Holm–Šidák’s multiple-comparisons test for f and two-sided unpaired Welch’s t-test for g,m. Exact P values are shown in the relevant panels.
Source data
Extended Data Fig. 5 CXCR2–PI3K–AKT signalling suppresses FOXO3 activity to induce hepatocyte PHGDH and serine production.
a, Immunoblot analysis of primary mouse hepatocytes cultured alone or co-cultured with sgTOM or sgPhgdh HY19636 cells under −SG conditions and treated with CXCR2i (SB225002, 100 nM) or PI3Ki (BKM120, 50 nM) for 24 h. b, Peak areas of intracellular serine isotopologues corresponding to Fig. 4d. Stacked bars indicate mean ± s.d. from n = 3 biological replicates. c, Crystal violet–based viability assay of sgPhgdh HY19636 cells cultured alone or in transwell co-culture with primary mouse hepatocytes, with CXCL5-neutralizing antibody, CXCR2 inhibitor (SB225002, 100 nM) or PI3K inhibitor (BKM120, 50 nM) as indicated. High SG, 400 µM; low SG, 150 µM. n = 3 biological replicates. d, Crystal violet–based viability assay of sgPhgdh HY19636 cells reconstituted with PHGDH(WT) or catalytically inactive PHGDH (PHGDH(CD)) under the conditions in c. n = 3 biological replicates. e, Schematic of RNA-seq comparisons used for f,g: primary mouse hepatocytes cultured alone versus co-cultured with sgPhgdh HY19636 cells (top), and hepatocytes co-cultured with sgPhgdh HY19636 cells treated with CXCR2 inhibitor (SB225002, 100 nM) versus vehicle (bottom), under −SG conditions. f, GSEA showing negative enrichment of the FoxO signalling pathway in hepatocytes co-cultured with sgPhgdh HY19636 cells compared with monocultured hepatocytes (top), and positive enrichment of the FoxO signalling pathway in hepatocytes co-cultured with sgPhgdh HY19636 cells treated with CXCR2 inhibitor compared with untreated controls (bottom). g, GSEA showing negative enrichment of the FOXO3 target-gene set in hepatocytes co-cultured with sgPhgdh HY19636 cells compared with monocultured hepatocytes (top), and positive enrichment of the FOXO3 target-gene set in hepatocytes co-cultured with sgPhgdh HY19636 cells treated with CXCR2 inhibitor compared with untreated controls (bottom). NES, normalized enrichment score. h, Immunofluorescence analysis of human PDAC liver metastases from Fig. 1a. Representative images are shown for regions with strong or weak hepatocyte PHGDH staining. In regions with strong hepatocyte PHGDH staining, hepatocyte FOXO3A was predominantly non-nuclear (Pearson’s R = 0.18; Spearman’s ρ = 0.19). By contrast, in regions with weak hepatocyte PHGDH staining, FOXO3A showed predominantly nuclear localization (Pearson’s R = 0.52; Spearman’s ρ = 0.54). i, Representative immunofluorescence images of mouse liver metastases formed by sgPhgdh HY19636 cells, with or without CXCR2 inhibitor treatment. Pearson’s correlation coefficients between DAPI and FOXO3A signals are indicated. Similar FOXO3A localization patterns were observed across the indicated biologically independent samples (h, eight and seven patient specimens, respectively; i, six mice per group). For h,i, the reported correlation coefficients were calculated from the representative images shown. For c,d, bars indicate mean ± s.d.; statistical significance was assessed using one-way ANOVA with Tukey’s multiple-comparisons test. Exact P values are shown in the relevant panels.
Source data
Extended Data Fig. 6 Phospho-mutant FOXO3A phenocopies the effects of abrogating CXCR2 and PI3K axes on hepatocyte serine biosynthesis.
a, Immunoblot showing FOXO3A expression in AML12 mouse hepatocytes, including parental cells, Foxo3a-knockout cells (Foxo3a-KO), and Foxo3a-KO cells reconstituted with FOXO3A(WT) or non-phosphorylatable FOXO3A(AAA) (T32A/S253A/S315A). β-actin was used as a loading control. b, Representative immunofluorescence images of FOXO3A (green) and DAPI (blue) in AML12 (Foxo3a-KO) cells reconstituted with FOXO3A(WT) or FOXO3A(AAA) and treated with DMSO or a PI3K inhibitor (PI3Ki). c, Quantification of FOXO3A nuclear localization in b, shown as the correlation coefficient between DAPI and FOXO3A signals. n = 3 biological replicates. d, Dual-reporter assay of PHGDH promoter activity in AML12 (Foxo3a-KO) cells reconstituted with FOXO3A(WT) or FOXO3A(AAA), expressed as the ratio of secreted Gaussia luciferase (GLuc; driven by the PHGDH promoter) to secreted alkaline phosphatase (SEAP) internal control. n = 4 biological replicates. e, Schematic model illustrating how non-phosphorylatable FOXO3A(AAA) (T32A/S253A/S315A) constitutively localizes to the nucleus and represses Phgdh transcription in hepatocytes, compared with FOXO3A(WT). f, FOXO3A ChIP–qPCR at the Phgdh promoter (P2; Fig. 4h) in AML12 (Foxo3a-KO) cells reconstituted with FOXO3A(WT) or FOXO3A(AAA) after co-culture with sgPhgdh HY19636 cells under −SG conditions. Data are shown as fold change (FOXO3A/IgG) normalized to 5% input. n = 3 biological replicates. g, Schematic of the 13C6-glucose tracing experiment in sgPhgdh HY19636 cells cultured alone or in transwell co-culture with AML12 cells reconstituted with FOXO3A(WT) or FOXO3A(AAA), under −SG conditions, followed by GC–MS analysis. h, Fractional isotopologue distribution (M+0 to M+3) of intracellular serine in sgPhgdh HY19636 cells from g, shown as percentage serine labelling from 13C6-glucose. Stacked bars indicate mean ± s.d. from n = 3 biological replicates. i, Immunoblot of CXCR2 in AML12 cells (wild type; WT or Cxcr2-KO) used for co-culture with sgTOM or sgPhgdh HY19636 cells. β-actin was used as a loading control. j, Fractional isotopologue distribution (M + 0 to M + 3) of intracellular serine in sgTOM or sgPhgdh HY19636 cells cultured alone or in transwell co-culture with AML12 cells (WT, Cxcr2-KO, or FOXO3A(AAA)) in −SG medium, shown as percentage serine labelling from 13C6-glucose. Stacked bars indicate mean ± s.d. from n = 3 biological replicates. k, Crystal-violet-based proliferation assay (absorbance at 590 nm) of HY19636 sgTOM or sgPhgdh PDAC cells cultured alone or in transwell co-culture with AML12 cells (WT, Cxcr2-KO, or FOXO3A(AAA)) under low SG conditions for 5 days. n = 3 biological replicates. Bars indicate mean ± s.d. for c,k and mean ± s.e.m. for d,f. Each symbol in graphs c,d,f,k represents an independent sample or biological replicate, as indicated. P values were calculated using one-way ANOVA followed by Holm–Šidák’s multiple-comparisons test for k, two-way ANOVA with Tukey’s multiple-comparisons test for c,h,j, and two-sided unpaired Student’s t-test for d,f. Exact P values are shown in the relevant panels.
Source data
Extended Data Fig. 7 CXCR2 inhibition reduces hepatocyte PHGDH induction and does not have any overt effects on liver weight.
a, Experimental scheme for b. b, Whole liver weights of sham-operated mice (n = 5, 4, 4 and 4 mice, respectively). The mean liver weight across all sham-operated mice was 958.2 mg and is shown as a blue dotted line in Fig. 5b. Each symbol represents an individual mouse. Bars indicate mean ± s.e.m. Statistical significance was assessed using one-way ANOVA followed by Holm–Šidák’s multiple-comparisons test. Exact P values are shown in the panel. c,d, Representative images (c) and quantification (d) of PHGDH staining in mouse livers bearing sgTOM or sgPhgdh HY19636 metastases under a −SG diet and treated with vehicle or a CXCR2 inhibitor. Numbers in bars indicate the number of mice.
Source data
Supplementary information
Source data
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Reprints and permissions
About this article
Cite this article
Ganguly, K., Yamamoto, K., Rodencal, J. et al. Hepatocytes promote liver metastasis of pancreatic cancer by providing serine. Nature (2026). https://doi.org/10.1038/s41586-026-11051-z
Download citation
Received:
Accepted:
Published:
Version of record:
DOI: https://doi.org/10.1038/s41586-026-11051-z