Isotope tracing reveals distinct substrate preference in murine melanoma subtypes with differing anti-tumor immunity
Cancer & Metabolism volume 10, Article number: 21 (2022)
Research about tumor “metabolic flexibility”—the ability of cells to toggle between preferred nutrients depending on the metabolic context—has largely focused on obesity-associated cancers. However, increasing evidence for a key role for nutrient competition in the tumor microenvironment, as well as for substrate regulation of immune function, suggests that substrate metabolism deserves reconsideration in immunogenic tumors that are not strongly associated with obesity.
We compare two murine models: immunologically cold YUMM1.7 and immunologically-hot YUMMER1.7. We utilize stable isotope and radioisotope tracer-based metabolic flux studies as well as gas and liquid chromatography-based metabolomics analyses to comprehensively probe substrate preference in YUMM1.7 and YUMMER1.7 cells, with a subset of studies on the impact of available metabolites across a panel of five additional melanoma cell lines. We analyze bulk RNA-seq data and identify increased expression of amino acid and glucose metabolism genes in YUMMER1.7. Finally, we analyze melanoma patient RNA-seq data to identify potential prognostic predictors rooted in metabolism.
We demonstrate using stable isotope tracer-based metabolic flux studies as well as gas and liquid chromatography-based metabolomics that immunologically-hot melanoma utilizes more glutamine than immunologically-cold melanoma in vivo and in vitro. Analyses of human melanoma RNA-seq data demonstrate that glutamine transporter and other anaplerotic gene expression positively correlates with lymphocyte infiltration and function.
Here, we highlight the importance of understanding metabolism in non-obesity-associated cancers, such as melanoma. This work advances the understanding of the correlation between metabolism and immunogenicity in the tumor microenvironment and provides evidence supporting metabolic gene expression as potential prognostic factors of melanoma progression and may inform investigations of adjunctive metabolic therapy in melanoma.
Deidentified data from The Cancer Genome Atlas were analyzed.
Melanoma is the most lethal form of skin cancer . Checkpoint-blockade immunotherapy has significantly decreased mortality from advanced-stage melanoma. Long-term outcomes in melanoma patients treated with checkpoint inhibitors are limited by the development of immunotherapy resistance, especially in late-stage melanoma patients with metastasis. This limitation means that there is an urgent need to develop better (or adjuvant) therapies for melanoma treatment. Metabolic therapies, which may decrease tumor cell fitness relative to healthy cells, could hold therapeutic promise in this area. Our study begins to investigate this possibility by seeking to identify which tumor type—those that are responsive to immunotherapy (“immune hot”) and those that are not (“immune cold”)—might respond best to metabolic therapies.
For more than a century, researchers have recognized metabolism as a mediator of tumor cell fitness. The Warburg effect describes tumor cells’ pathognomonic reliance on glycolysis to produce ATP and building blocks for rapid cell division [2,3,4]. Investigators have taken multiple approaches to exploit the Warburg effect. Some of these approaches have been effective preclinically. For example, dichloroacetate (DCA) inhibits pyruvate dehydrogenase kinase, activating pyruvate dehydrogenase and driving pyruvate into the tricarboxylic acid (TCA) cycle . In vitro DCA treatment increases apoptosis, slows tumor cell division, and increases chemotherapy efficacy [6,7,8,9], likely by diverting substrates away from biosynthetic pathways required for cell division.
Recent work has also highlighted the importance of metabolism in immune cell function. During the anti-tumor immune response, tumor-infiltrating lymphocytes (TILs) undergo a metabolic switch. Following activation, TILs change their primary energy source from oxidative phosphorylation to glycolysis to meet increased energy requirements and generate necessary material for cell division and differentiation [10, 11]. In short, high-functioning TILs rely on glycolysis. Therefore, researchers have sought to improve responsiveness to immunotherapy by modulating glycolysis using small molecules [12,13,14,15,16]. Given their expected impact to disturb normal metabolic processes in vivo, however, the translational potential of such systemic metabolic strategies may be limited.
Thus, researchers have turned to pathways that may enhance T cell activation and/or impede exhaustion while exerting a less profound effect on metabolic homeostasis. In examinations of other immune cell types, researchers have described glutamine as the “fuel of the immune system” . Lymphocytes increase the uptake of glutamine during the anti-tumor response. Researchers have also suggested that glutamine may be essential for lymphocyte proliferation and function . Supplementing lymphocytes and macrophages with glutamine, for instance, enhances function [19, 20]. Recent studies, however, indicate that the direction of this relationship is unclear. For example, Nabe et al. recently showed that glutamine deprivation during in vitro activation improves T cell tumor infiltration and responsiveness to anti-PD1 treatment in tumor-bearing mice . Most investigations into the role of glutamine have focused on the impact of glutamine metabolism in TILs in vitro, while largely ignoring the impact of glutamine on tumor cell division and function in the in vivo tumor microenvironment (TME). However, pharmacologic interventions to manipulate metabolism in vivo are difficult to administer in a cell-specific manner; most small molecules administered systemically will affect immune cell metabolism and, simultaneously, metabolism in tumor cells and other tissues. Our study seeks to combine in vivo and in vitro studies implanting [13C] tracer methodology combined with gas and liquid chromatography-based analysis to address these gaps.
In this report, we compare two common murine melanoma models: YUMM1.7 and YUMMER1.7 , adding five additional cell lines with a variety of immunogenicity for confirmatory studies. YUMM1.7 is “immunologically cold,” meaning that lymphocytes minimally infiltrate YUMM1.7 tumors and that YUMM1.7 does not respond to immune checkpoint treatments. YUMMER1.7 is derived from YUMM1.7, but in contrast to YUMM1.7, YUMMER1.7 is “immunologically hot.” YUMMER1.7 therefore has high T cell infiltration and responds to immune checkpoint treatments . However, the interplay between metabolic changes and immunogenicity is not fully understood.
In this study, we demonstrate an association between immune “hot” or “cold” status and substrate utilization in the TME by comparing gene expression, metabolic flux, and cell survival under nutrient-depleted conditions in two murine melanoma models, with confirmatory studies in five more. We demonstrate that glutamine transporter expression is positively correlated with lymphocyte infiltration and anti-tumor effector function. This work advances the understanding of the relationship between tumor substrate preference/utilization and tumor cell division in the TME. This work may also influence the design of clinical trials investigating metabolic adjuvants to immunotherapy in melanoma, with a particular emphasis on targeting tumor metabolism to, in turn, optimize the ability of the immune system to fight melanoma.
We purchased YUMM1.7, B2905, and B16-F10 cells from ATCC (CRL-3362, Cellosaurus Accession Number CVCL_JK16; and CRL-3476, Cellosaurus Accession Number CVCL_B0CG; and CRL-6475, Cellosaurus Accession Number CVCL_0159, respectively) and D4M-3A cells from Sigma (SCC428, Cellosaurus Accession Number CVCL_0P27) and obtained YUMMER1.7 cells from the Bosenberg laboratory, where they are regularly authenticated. We did not authenticate the cells purchased commercially, because this is routinely done by the suppliers. We maintained cells from both lines in the media shown in Table 1.
In all cases, the FBS was not dialyzed, and 1% penicillin/streptomycin was also included. For the irradiation studies, B2905 and D4M-3A were thrice UVB irradiated (600 J/m2)  within 24 h prior to injection. In all cases, we injected cells into mice after passage <10. In in vitro flux studies, we serum-starved cells for 14 h prior to incubation in a tracer. We supplied tracer media containing [13C6] glucose (5 mM), [13C] palmitate/oleate (2.5 mM of each) or [13C5] glutamine (5 mM). We also provided unlabeled substrates at these concentrations to ensure that substrates were not limiting, regardless of 13C label (in other words, in cell survival experiments and all other studies, “full media” refers to 5 mM glucose/5 mM fatty acid (1:1 palmitate to oleate)/5 mM glutamine). We collected cells in 50% ethanol after the incubation times indicated in the figures. For cell survival experiments, we plated cells in recommended media in 24-well plates, with three wells per condition (i.e., media containing glucose and glutamine, but no palmitate; media containing glucose and palmitate, but no glutamine, etc.). At 24 h, 48 h, and 72 h after plating, cells are stained with Trypan blue and counted (3 measurements per condition) by a LUNA-II automated cell counter.
Lactate production was measured in cells incubated in lactate-free but otherwise complete media for 6 h. The rate of lactate production was assumed to be linear over the 6-h period and was normalized to total protein content measured using the BCA protein assay.
Mitochondrial stress test
To examine mitochondrial and non-mitochondrial respiration, we performed a mitochondrial stress test using the Seahorse XF Flux Analyzer. Studies were performed with cells at 70–80% confluency, and data normalized to total protein. The protocol used for the mitochondrial stress test has been published previously .
The Yale Institutional Animal Care and Use Committee approved all animal studies. We purchased wild-type male C57bl/6J mice from Jackson Laboratories (Bar Harbor, ME) at 8 weeks of age. After 1 week of acclimation, we injected mice subcutaneously in the right chest with 5x105 YUMM1.7 or YUMMER1.7 tumor cells after confirmation that they were pathogen-free (including mycoplasma) by the Yale Comparative Pathology Research Core. Three weeks after tumor cell injection, we performed surgery under isoflurane anesthesia to implant catheters in the jugular vein for tracer infusion (described in the “Flux studies” section). Tumor interstitial fluid was isolated from a separate set of mice not infused with tracer. Tumors were minced in PBS (2:1 v/w PBS to tumor), transferred to a centrifuge tube, and 8:1 PBS to tumor weight added. Samples were then incubated at 37°C for 60 min, centrifuged at 2000 rpm for 10 min, and the supernatant (tumor interstitial fluid) was collected.
We measured glucose and lactate concentrations enzymatically using the YSI 2500 Glucose/Lactate Analyzer (Yellow Springs, OH, USA). We measured glutamine concentration and enrichment by LC-MS/MS (AbSCIEX 6500 QTRAP with a Shimadzu ultrafast liquid chromatography system in negative ion mode) as we have previously reported  after spiking with [2H4] taurine as an internal standard. We measured citrate, succinate, and malate concentration and enrichment using the same LC-MS/MS method. Mass/charge ratios for the unlabeled metabolites are shown in Table 2. Due to space constraints, the labeled intermediate mass/charge ratios are not shown.
We measured amino acid concentrations with gas chromatography/mass spectrometry after derivatization as described previously . We measured non-esterified fatty acid (NEFA) concentrations using the Wako NEFA-HR assay (Mountain View, CA, USA).
Mice underwent a primed (3X)-continuous infusion of [U-13C5] glutamine, [U-13C6] glucose, or [U-13C16] palmitate (continuous infusion rate for each tracer, all 1 mg/kg body weight/min), or [14C] glutamine, [3H] glucose, or [14C] palmitate (continuous infusion rate 1 μCi/kg/min) for 60 min, unless otherwise stated. All stable isotopes were obtained from Cambridge Isotopes (Tewksbury, MA, USA), and all radioisotopes from PerkinElmer (Waltham, MA, USA). Within 30 s of euthanizing animals with IV pentobarbital (Covetrus; Portland, ME, USA), we isolated and freeze-clamped tumors in liquid nitrogen.
In method-developing trials, we observed that regardless of tracer identity ([13C] glutamine, glucose, or palmitate), accumulation of [13C] label was linear between 0 and 30 min and by comparing citrate enrichment in tumors of mice infused for 60 and 150 min, ascertained that 60 min of infusion was sufficient to reach steady-state under these conditions. Therefore, we infused mice with [13C] isotope for 60 min in the steady-state experiments. We calculated relative VCS as the slope of [13C] citrate enrichment measured at 60 min vs. time and normalized to YUMM1.7 tumors infused with [13C6] glucose. In mice infused with [U-13C6] glucose, the ratio of [13C3] pyruvate/[13C6] glucose was used as a readout of glycolysis, although we acknowledge that this ratio is not a direct readout of glycolysis. In in vivo radioisotope infusion studies, we assumed that [14C] or [3H] label would reach a steady state at a similar timepoint as [13C] label and thus terminated radioisotope infusion studies at 60 min.
We measured in vitro cell viability and PD-L1 expression using the BD LSR II flow cytometer. Cells were digested into single cells with trypsin EDTA (0.25%), followed by cell viability staining (ThermoFisher Scientific, #L34964) for 20 min at room temperature. After washing with PBS two times, cells were either ready for viability analysis or stained with PD-L1 antibody (ThermoFisher, #12-5982-82) for 20 min under room temperature and washed with PBS three times. In both cases, cells were fixed by 4% polymeric formaldehyde (PFA) after staining and analyzed within 24 h.
All RNA-seq analysis follows standard DESeq workflow. YUMM1.7 and YUMMER1.7 bulk RNA-Seq data are available upon request from author Marcus W. Bosenberg and the remainder of the data from corresponding author Rachel J. Perry. The human PRECOG database, a meta-analysis of gene expression, and overall survival data are available at https://precog.stanford.edu/, and the TCGA human tumor gene expression data (TCGA Melanoma (SKCM)) can be downloaded pre-filtered and normalized at https://xenabrowser.net/datapages/.
The principal component analysis, volcano plot, and correlation maps were generated in R. Statistical significance for genes revealed by differential expression analysis was defined as a Log2(Fold Change) >2 or ≤2 and a Benjamini-Hochberg corrected p value of <0.01. PRECOG datasets were downloaded, and the melanoma and metastatic melanoma meta-analyses were filtered to metabolic gene sets downloaded from the Kyoto Encyclopedia of Genes and Genomes (KEGG). Z scores in the PRECOG dataset reflect the effect of high gene expression on overall survival in publicly available tumor gene expression datasets, where a Z score > 3.09, which correlates to a p value 0.001, corresponds to a detrimental effect on survival, while a Z score < −3.09 corresponds to a beneficial effect on survival with a p value of 0.001.
The heatmap of gene expression was generated in Python, with an automated clustering algorithm was employed using Euclidean distances between genes, and with a Z score applied to each row, using Seaborn’s Clustermap function.
All code used to generate figures is available online (https://github.com/xz710/YUMM-YUMMER-comparison.git).
We used the 2-tailed unpaired Student’s t test to compare 2 groups, and ANOVA with Tukey’s multiple comparisons test to compare 3 or more groups, in GraphPad Prism Version 9. In all cases, we verified that the data met the assumptions of the statistical test employed.
YUMMER1.7 cells utilize more glucose and glutamine than YUMM1.7 cells in vitro
We first sought to identify differences between the RNA expression of metabolic genes in YUMM1.7 and YUMMER1.7 murine melanoma models. We hoped such differences might indicate a mechanism contributing to the immunogenic capacity of the tumor. To avoid confounding by the presence or absence of TILs, we compare bulk RNA sequencing data from YUMM1.7 and YUMMER1.7 cell lines. Principal component analysis (PCA) revealed distinct separation between the two cell lines, with 36.4% of the variance in the datasets explained by the first principal component, along which the samples appeared to separate (Fig. 1A). We identified 862 genes expressed differently between these models (Fig. 1B). Some of these genes were key metabolic regulators. For instance, Slc43a1 and Slc43a3 encode neutral amino acid transporters, Cd36 encodes a glycoprotein functioning as a fatty acid transporter in many tissues, and Gfpt2 encodes glutamine-fructose-6-phosphate transaminase 2, which controls glucose flux into hexosamine pathway (Supplementary Table S1). We then tested whether the observed differences in gene expression correlated to substrate flux through these enzymes.
To do so, we used stable isotope tracers to compare oxidation of glucose, glutamine, and fatty acids between YUMM1.7 and YUMMER1.7 cells (Fig. 1C–F). We chose these metabolites as they are the three most important energy resources in most cells. We calculated the TCA cycle flux driven by each nutrient by measuring the enrichment of 13C-labeled metabolites using gas and liquid chromatography-mass spectrometry (GC-MS and LC-MS/MS, respectively). We first confirmed that cells reached an isotopic steady state in all isotopologues 13C citrate after incubation in media containing either [U-13C6] glucose, [U-13C5] glutamine, or [U-13C16] palmitic acid within 30 min (Supplementary Figure S1A-S1R).
We measured [13C2] citrate derived from [13C6] glucose as a readout of the contribution of glucose to the TCA cycle. Relative citrate synthase flux (Vcs), defined as the slope of [13C] citrate enrichment measured at 30 min, represents the metabolic flux entering the TCA cycle from glycolysis. We observed a higher contribution of glucose to citrate in YUMMER1.7 as compared to YUMM1.7 cells, indicating higher glucose oxidation in immune “hot” YUMMER1.7 cells (Fig. 2A). Moreover, we observed a higher dilution of [13C] glucose label throughout the TCA cycle in YUMMER1.7, indicating higher anaplerotic influx into the TCA cycle from other unlabeled substrates (Fig. 2B, Supplementary Figure S2A-E). Glutamine contribution to VCS was also higher in YUMMER1.7 cells, as was the dilution of the label through the TCA cycle, reflecting anaplerosis (Fig. 2C, D and Supplementary Figure S2F-J). The contribution of fatty acids to the TCA cycle was similar between the two models, but anaplerosis, as expected, remained greater in YUMMER1.7 cells (Fig. 2E, F and Supplementary Figure S2F). These data were expected because rates of anaplerosis (i.e., unlabeled substrate entry into the TCA cycle) are unaffected by which tracer is used.
We also measured tumor cell proliferation while depriving cells of key nutrients in vitro. Our data demonstrate that glucose deprivation had no impact on either YUMM1.7 or YUMMER1.7 cell proliferation in vitro, whereas glutamine deprivation inhibited cell proliferation in both models (Fig. 2G, I). Moreover, incubation in free fatty acid had a deleterious effect in vitro, which further suggests that fatty acid is unlikely to promote melanoma anti-tumor immune cell function in the melanoma TME. This effect may be attributable in part to palmitate toxicity, which has been observed in multiple cell types including adipocytes , hepatocytes , and cardiomyocytes , among others; however, we anticipate that the inclusion of the unsaturated fatty acid oleate may mitigate this toxicity, as has been demonstrated recently in additional cell types [30, 31]. The fact that the tumor cell number curves are not exponential growth suggests that cells may be approaching a horizontal asymptote in the division at the 72-h time point. This effect is unlikely to result from substrate limitation considering the abundant, substrate-replete, media provided, but is more likely a result of approaching confluency and resulting physical changes. It is unlikely that these effects were the result of increased cell death with substrate deprivation, while glucose deprivation had a small (2.5%) effect to reduce the fraction of live YUMM1.7 cells, glutamine deprivation had no effect, and neither substrate affected the fraction of live YUMMER1.7 cells (Fig. 2J, K).
We then expanded the scope of our study of the impact of available substrates on melanoma cells by examining the impact of substrate depletion on five additional melanoma cell lines, with a range of immunogenicity: UVB-irradiated B2905 , D4M-3A [32, 33], B16-F10 , not irradiated B2905 , and irradiated D4M-3A. Our data reveal an increasing dependence upon glutamine in the more immunogenic cell lines: incubation in media lacking glutamine had a greater effect to suppress cell division in B2905 and irradiated D4M-3A cells, both of which are considered more immunogenic melanoma cell lines, as compared to the others (Fig. 2L).
Melanoma cells shift their metabolism, to varying effects, in response to substrate deprivation. YUMMER1.7 cells showed a greater reliance upon glucose as reflected by an increased conversion of 13C glucose to 13C pyruvate and increased production of lactate in vitro, both resulting from glycolysis (Fig. 3A, B) in addition to basal respiration, ATP production, and maximal respiration during a mitochondrial stress test (Fig. 3C). To directly understand how substrate deprivation affects tumor cell metabolism, we performed a series of in vitro studies incubating YUMM1.7 and YUMMER1.7 cells in glucose-, fatty acid-, and glutamine-deprived media. We observed that upon incubation in glucose- or glutamine-free media, cells increased their oxidation of the other substrate, whereas incubation in fatty acid-free media did not alter the oxidation of glucose or glutamine (Fig. 3D–F).
Glucose and glutamine deprivation increases tumor PD-L1 expression. Next, to generate a direct readout of the impact of tumor cell substrate metabolism on tumor immunogenicity, we incubated YUMM1.7 and YUMMER1.7 cells in media lacking glucose or glutamine for 48 h and used flow cytometry to demonstrate how these perturbations affected cancer cell PD-L1 expression. As expected, PD-L1 expression was lower in “immune cold” YUMM1.7 than in “immune hot” YUMMER1.7 (Fig. 3G, H). Glucose deprivation, and to a greater extent glutamine deprivation, increased PD-L1 expression in both cell lines. Checkpoint inhibitors may have greater potential efficacy in tumors with higher PD-L1 expression, so while these data do not prove the impact of substrate metabolism on immunogenicity in melanoma cells, they do directly link tumor cell substrate metabolism to immunogenicity status.
Next, we employed LC-MS/MS to measure metabolite concentrations in tumor interstitial fluid. We observed a decrease in concentrations of glucose, glutamine, leucine, phenylalanine, and proline in the YUMMER1.7 TME (Fig. 4A, B). In contrast, the concentrations of lactate and lysine were significantly higher in the YUMMER1.7 TME.
YUMMER1.7 tumors utilize more glucose and glutamine than YUMM1.7 tumors in vivo
We performed steady-state (Supplementary Figure S3A-C) stable and radioisotope infusions to study substrate uptake in melanoma tumors in awake mice in vivo. Our data demonstrate that despite showing no differences in size, YUMMER1.7 tumors had higher uptake of glucose (Fig. 5A) and glutamine (Fig. 5C), while palmitate uptake remained the same between the two melanoma models (Fig. 5E). We used stable isotope infusions to examine the relative contributions of glucose, glutamine, and palmitate to VCS in tumor tissues and observed a similar result to the in vitro study. Glucose (Fig. 5B) and glutamine (Fig. 5D) had higher oxidation rates in YUMMER1.7 relative to YUMM1.7 tumors, whereas palmitic acid oxidation did not significantly differ between the two models (Fig. 5F).
Glutamine transporter expression is a prognostic factor in melanoma patients
Finally, we examined the prognostic impact of glutamine metabolism in melanoma patients, utilizing melanoma tumor gene expression from the PRECOG database . A positive Z score indicates a detrimental effect on overall survival, and a negative Z score indicates improved survival, with (−)3.09 correlating to a p value of 0.001, which we deemed as significant as described previously . Glucose and fatty acid metabolic genes mostly had a detrimental effect in both early- and late-stage melanoma patients (Supplementary Figure S4A and S4B). By contrast, higher expression of metabolic genes, as well as GLUL, correlated with an improved effect on prognosis, especially for late-stage melanoma patients with metastasis (Fig. 6A). Furthermore, although high expression of most amino acid transporters was correlated with increased mortality, elevated expression of SLC38A1 (glutamine transporter) and SLC7A7 (lysine transporter) had strong benefits for patients with metastatic melanoma (Fig. 6B) . Though there was no obvious subsetting of melanoma patients by automated clustering, there were distinct amino acid gene clusters that varied in high vs. low expression across all melanoma patients (Fig. 6C).
Although our study is primarily focused on the impact of tumor metabolism on tumor-immune cell metabolic crosstalk and outcomes, we also examined the connection between tumor metabolism and anti-tumor lymphocyte infiltration and function. We calculated the correlation coefficient of leukocyte fraction (identified from pathology images using deep learning ) to amino acid metabolism-related gene expression in 369 melanoma patients from the TCGA database (Supplementary Figure S5B). We identified 14 genes (SLC7A7, IL41, SLC43A2, SLC38A5, ACY3, SLC6A12, SLC38A1, DDO, GLUL, ASS1, ABAT, GFPT2, SLC6A6, SLC36A1, SLC38A4) that correlate positively and 13 genes (SLC7A5, SLC16A10, SLC3A1, GLUD2, GAD2, SLC3A2, SLC43A1, PPAT, ADSL, CAD, SLC6A15, SLC7A1) that correlate negatively with leukocyte fraction in tumor tissue (Fig. 6D), with similar patterns observed with gene expression correlations with cytolytic activity (Fig. 6E). In contrast, we found that patients’ overall survival time did not correlate with leukocyte fraction within tumor tissue (Supplementary Figure S5A), implying that other factors, including anti-tumor function of lymphocytes, are important than strict lymphocyte fraction for clinical outcomes.
We approximated this lymphocyte effector function by estimating immune cytolytic activity (CYT). CYT is the geometric mean of the expression level of GZMA—which encodes serine protease—and PRF1—which encodes perforin. GZMA and PRF1 are the key enzymes involved in cytotoxic T lymphocyte effector function . Higher CYT indicates greater cytotoxic T cell activity. By studying 159 patients from the TCGA database (Supplementary Figure S3C), we identified 10 genes (SLC7A7, IL41, SLC43A2, GLUL, SLC6A12, SLC38A5, SLC38A1, DDO, SLC38A4) that correlate positively and 11 genes (ADSL, SLC7A5, NIT2, ASL, SLC38A2, GOT2, CAD, SLC16A10, SLC6A15, SLC7A1) that correlate negatively with CYT (Fig. 5E).
Substrate preference is a teleological concept based on correlations between nutrient utilization and outcomes, rather than a single measurable variable. For this reason, we take a “multi-omics” approach to define how differences in metabolism may correlate with tumor immunogenicity in melanoma cells. We combine several analyses, including gene expression, metabolic flux, and survival in environments with restricted substrate availability. We use these methods to examine how the immunologically “hot” or “cold” identity of melanoma tumors reciprocally affects tumor cell substrate preference within the TME.
Researchers have long considered glutamine an essential nutrient for lymphocyte proliferation and effector function [38, 39], although a recent study on melanoma-bearing mice paradoxically suggested that glutamine deprivation enhances anti-tumor function and responsiveness to anti-PD1 treatment . Most of these studies, however, were conducted in an in vitro T cell culture system and lack in vivo evidence to confirm the physiologic impact of T cell metabolism on tumor cells in the TME. Our data demonstrate that glutamine is a key substrate for tumor cells themselves in immunologically “hot” melanoma, both in vivo and in vitro.
Our study reveals divergent roles for glutamine in the TME between two genetically similar, immunologically “cold” (YUMM1.7) and “hot” (YUMMER1.7) murine tumor models. We compared RNA expression between YUMM1.7 and YUMMER1.7 to avoid confounding by the inclusion of TIL gene expression. We found that genes belonging to glucose and amino acid metabolism pathways are enriched in YUMMER1.7 compared to YUMM1.7. We then used isotope tracing methods to compare metabolic fluxes in both in vitro cell culture and in vivo tumor tissue. We discover that the uptake of glutamine and glucose, as well as the relative oxidation of both glucose and glutamine, is higher in YUMMER1.7 than in YUMM1.7. These data are consistent with our RNA-seq analysis, which indicated that glucose and glutamine metabolism pathways are enriched in YUMMER1.7. The uptake and oxidation rate of palmitic acid, by contrast, was similar between the two models.
We recognize that there are potential limitations to any tracer method. These limitations include (1) the invalidity of data when intracellular tracer enrichment and concentration are not at steady-state, (2) the possibility of entry of unaccounted anaplerotic substrates or metabolic pathways (for example, dilution of [13C] acetyl CoA by unlabeled citrate through ATP citrate lyase), and (3) incomplete equilibration between cytosolic and mitochondrial substrate pools. Additionally, we cannot rule out a contribution of altered metabolism in other cell populations within the tumor, such as immune and stromal cells, particularly if the fraction of these cell types differs between YUMM1.7 and YUMMER1.7 tumors. However, the results of our in vitro studies demonstrating differing effects of glucose and glutamine in melanoma cell lines with varying immunogenicity indicate that there are cell-specific effects of metabolic changes in melanoma. Altogether, we suggest that studies should use isotope tracer methods in combination with independent functional readouts, such as the in vitro substrate incubation studies included herein.
We also measured the concentrations of metabolites in tumor interstitial fluid. We identified a decrease in key nutrients, including glucose, alanine, proline, phenylalanine, and glutamine, in the YUMMER1.7 TME. These data indicate that a smaller quantity of key nutrients are available to CD8+ T cells in the immunotherapy-responding melanoma model (that is, in YUMMER1.7, the “hot” tumor model). Considering the increased anaplerotic flux in YUMMER1.7 cells, the lower substrate concentrations in the YUMMER1.7 TME likely represent increased substrate uptake. This highlights the importance of using tracers to study metabolism, rather than simply relying on nutrient concentrations.
In contrast, the concentration of lysine was higher in the YUMMER1.7 TME. Although lysine is an essential amino acid that enters the TCA cycle , the relationship between lysine availability and immune cell function is not well understood. Further work is necessary to explore whether therapies seeking to improve melanoma immunotherapy efficacy should target lysine metabolism.
We also found an increase in the lactate concentration in the YUMM1.7 TME. This might be due to increased rates of glycolysis in YUMM1.7. Glycolysis is elevated in YUMM1.7 to compensate for the relative decrease in anaplerotic substrate supply and, consequently, the expected decrease in oxidative phosphorylation.
Finally, and most importantly with regard to clinical implications, we studied the impact of amino acid uptake and amino acid metabolism on human melanoma patients. Unlike elevated glucose and fatty acid metabolic gene expression, which is generally deleterious for patients, heightened anaplerotic gene expression (e.g., D-aspartate oxidase (DDO) or glutaminase 2 (GLS2)) is a prognostic indicator of a beneficial effect. This effect is especially apparent in late-stage melanoma patients. We also used the TCGA database to study the impact of amino acid metabolism gene expression on TIL infiltration and cytolytic activity (CYT) in melanoma patients. We found that glutamine transporter expression in tumor tissue is positively correlated with both leukocyte infiltration and CYT. Based on our and others’ data demonstrating that tumor cells actively take up and oxidize glutamine, tumor tissue with higher glutamine transporter expression have a higher capacity to “soak up” glutamine, leading to glutamine deprivation in the TME. This again suggests that lower glutamine concentration in the TME is beneficial for CD8+ T cell anti-tumor function.
This study has the inevitable caveat that human RNA-seq data is obtained from tumor samples that may or may not contain TILs. Thus, one can argue the positive association between anaplerotic gene expression and lymphocyte activation is due to increased activated lymphocytes in sequencing tissue rather than due to increased anaplerotic gene expression in tumor tissue itself. Regardless, our results validate anaplerotic gene expression as a prognostic factor for lymphocyte effector function and potentially predict responsiveness to immune checkpoint inhibitors.
While our study lays a foundation for further studies of competition for substrates between tumor and immune cells in the TME in “immunologically hot” melanoma, we do not test competition with immune cells. Further experiments will be needed to directly test the “substrate competition” hypothesis. These studies must also understand that the idea of substrate competition by definition is teleological and will be almost impossible to conclusively prove experimentally. In other words, we can measure substrate uptake and use, but without the ability to interview tumor cells and immune cells, it is difficult to prove which substrate the cells “seek.” Even so, the use of isotope tracers enables a unique window into substrate utilization in tumor cells in vivo and in vitro.
In sum, our study provides insight into how different tumor subtypes vary in metabolic flexibility and substrate preference, suggesting that human tumor subtypes, too, may vary along such dimensions. Unlike many prior studies, the focus of the current work is on how tumor metabolism predicts immunogenicity, as opposed to primarily centering immunometabolism. These differences in metabolism between immunologically hot and cold melanoma cell lines may be exploited to design metabolic therapies to enhance immunotherapy efficacy and overcome immunotherapy resistance in patients with melanoma and potentially other cancers.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Non-esterified fatty acid
Principal Component Analysis
TriCarboxylic acid cycle
Bolick NL, Geller AC. Epidemiology of melanoma. Hematol Oncol Clin North Am. 2021;35(1):57–72.
Vander Heiden MG, Cantley LC, Thompson CB. Understanding the Warburg effect: the metabolic requirements of cell proliferation. Science. 2009;324(5930):1029–33.
Warburg O. On the origin of cancer cells. Science. 1956;123(3191):309–14.
Warburg O, Wind F, Negelein E. The metabolism of tumors in the body. J Gen Physiol. 1927;8(6):519–30.
do Nascimento RS, et al. Sodium dichloroacetate attenuates the growth of B16-F10 melanoma in vitro and in vivo: an opportunity for drug repurposing. Anticancer Drugs. 2021;32(2):111–6.
Tataranni T, et al. Dichloroacetate affects mitochondrial function and stemness-associated properties in pancreatic cancer cell lines. Cells. 2019;8(5).
Kim TS, et al. Metformin and dichloroacetate suppress proliferation of liver cancer cells by inhibiting mTOR complex 1. Int J Mol Sci. 2021;22(18).
Woo SH, et al. Dichloroacetate potentiates tamoxifen-induced cell death in breast cancer cells via downregulation of the epidermal growth factor receptor. Oncotarget. 2016;7(37):59809–19.
Kinnaird A, et al. Metabolic modulation of clear-cell renal cell carcinoma with dichloroacetate, an inhibitor of pyruvate dehydrogenase kinase. Eur Urol. 2016;69(4):734–44.
van der Windt GJ, Pearce EL. Metabolic switching and fuel choice during T-cell differentiation and memory development. Immunol Rev. 2012;249(1):27–42.
Chang CH, et al. Posttranscriptional control of T cell effector function by aerobic glycolysis. Cell. 2013;153(6):1239–51.
Meng G, et al. Targeting aerobic glycolysis by dichloroacetate improves Newcastle disease virus-mediated viro-immunotherapy in hepatocellular carcinoma. Br J Cancer. 2020;122(1):111–20.
Ohashi T, et al. Dichloroacetate improves immune dysfunction caused by tumor-secreted lactic acid and increases antitumor immunoreactivity. Int J Cancer. 2013;133(5):1107–18.
Liu W, et al. Pleiotropic effects of metformin on the antitumor efficiency of immune checkpoint inhibitors. Front Immunol. 2020;11:586760.
Afzal MZ, Mercado RR, Shirai K. Efficacy of metformin in combination with immune checkpoint inhibitors (anti-PD-1/anti-CTLA-4) in metastatic malignant melanoma. J Immunother Cancer. 2018;6(1):64.
Ma G, et al. Targeted glucose or glutamine metabolic therapy combined with PD-1/PD-L1 checkpoint blockade immunotherapy for the treatment of tumors - mechanisms and strategies. Front Oncol. 2021;11:697894.
Rodas PC, et al. Glutamine and glutathione at ICU admission in relation to outcome. Clin Sci (Lond). 2012;122(12):591–7.
Carr EL, et al. Glutamine uptake and metabolism are coordinately regulated by ERK/MAPK during T lymphocyte activation. J Immunol. 2010;185(2):1037–44.
Oh MH, et al. Targeting glutamine metabolism enhances tumor-specific immunity by modulating suppressive myeloid cells. J Clin Invest. 2020;130(7):3865–84.
Wei J, et al. Nutrient and metabolic sensing in T cell responses. Front Immunol. 2017;8:247.
Nabe S, et al. Reinforce the antitumor activity of CD8(+) T cells via glutamine restriction. Cancer Sci. 2018;109(12):3737–50.
Wang J, et al. UV-induced somatic mutations elicit a functional T cell response in the YUMMER1.7 mouse melanoma model. Pigment Cell & Melanoma Research. 2017;30(4):428–35.
Wolf Y, et al. UVB-induced tumor heterogeneity diminishes immune response in melanoma. Cell. 2019;179(1):219–235 e21.
Camporez JP, et al. Cellular mechanisms by which FGF21 improves insulin sensitivity in male mice. Endocrinology. 2013;154(9):3099–109.
Perry RJ, et al. Leptin mediates a glucose-fatty acid cycle to maintain glucose homeostasis in starvation. Cell. 2018;172(1-2):234–248 e17.
Leimer KR, Rice RH, Gehrke CW. Complete mass spectra of N-trifluoroacetyl-n-butyl esters of amino acids. J Chromatogr. 1977;141(2):121–44.
Davis JE, et al. The c-Jun N-terminal kinase mediates the induction of oxidative stress and insulin resistance by palmitate and toll-like receptor 2 and 4 ligands in 3T3-L1 adipocytes. Horm Metab Res. 2009;41(7):523–30.
Karim N, et al. Pelargonidin-3-O-glucoside encapsulated pectin-chitosan-nanoliposomes recovers palmitic acid-induced hepatocytes injury. Antioxidants (Basel). 2022;11(4).
Joseph LC, et al. Inhibition of NAPDH Oxidase 2 (NOX2) Prevents oxidative stress and mitochondrial abnormalities caused by saturated fat in cardiomyocytes. PLoS One. 2016;11(1):e0145750.
Urso CJ, Zhou H. Palmitic acid lipotoxicity in microglia cells is ameliorated by unsaturated fatty acids. Int J Mol Sci. 2021;22(16).
Tumova J, et al. Protective effect of unsaturated fatty acids on palmitic acid-induced toxicity in skeletal muscle cells is not mediated by PPARdelta activation. Lipids. 2015;50(10):955–64.
Jenkins MH, et al. Multiple murine BRaf(V600E) melanoma cell lines with sensitivity to PLX4032. Pigment Cell Melanoma Res. 2014;27(3):495–501.
Lo JA, et al. Epitope spreading toward wild-type melanocyte-lineage antigens rescues suboptimal immune checkpoint blockade responses. Sci Transl Med. 2021;13(581).
Zoller M. IFN-treatment of B16-F1 versus B16-F10: relative impact on non-adaptive and T-cell-mediated immune defense in metastatic spread. Clin Exp Metastasis. 1988;6(5):411–29.
Gentles AJ, et al. The prognostic landscape of genes and infiltrating immune cells across human cancers. Nat Med. 2015;21(8):938–45.
Saltz J, et al. Spatial organization and molecular correlation of tumor-infiltrating lymphocytes using deep learning on pathology images. Cell Rep. 2018;23(1):181–193 e7.
Rooney MS, et al. Molecular and genetic properties of tumors associated with local immune cytolytic activity. Cell. 2015;160(1-2):48–61.
Johnson MO, et al. Nutrients and the microenvironment to feed a T cell army. Semin Immunol. 2016;28(5):505–13.
O'Neill LA, Kishton RJ, Rathmell J. A guide to immunometabolism for immunologists. Nat Rev Immunol. 2016;16(9):553–65.
Reimer LC, Will SE, Schomburg D. The fate of lysine: non-targeted stable isotope analysis reveals parallel ways for lysine catabolization in Phaeobacter inhibens. PLoS One. 2017;12(10):e0186395.
The authors thank members of the Perry and Bosenberg labs for helpful discussions. The results shown here are in part based upon data generated by the TCGA Research Network: https://www.cancer.gov/tcga. The authors thank Yale Flow Cytometry for their assistance with flow cytometry. The core is supported in part by an NCI Cancer Center Support Grant # NIH P30 CA016359. The BD symphony was funded by shared instrument grant # NIH S10 OD026996.
This study was funded by a Young Investigator Award from the Melanoma Research Alliance (to R.J.P.), the National Institutes of Health Medical Scientist Training Program Training Grant T32GM007205 (supporting B.P.L), a Timothy Dwight Richter Fellowship (to A.A.H.), and an award from the Yale Westchester Alumni Association (to A.A.H.).
Ethics approval and consent to participate
Animal studies were approved by the Yale University Institutional Animal Care and Use Committee (protocol 2019-20290). All samples in TCGA have been collected and utilized following strict human subject protection guidelines, informed consent, and IRB review of protocols.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Table S1. Differentially expressed genes comparing YUMMER1.7 to YUMM1.7. This 650 page table is uploaded as supplementary material exclusively, due to its length.
Mass isotopomer distribution in YUMM1.7 and YUMMER1.7 cells incubated in [U-13C6] glucose, [U-13C5] glutamine, or [U-13C16] palmitate. *P<0.05, **P<0.01, ***P<0.001, ****P<0.0001 by the 2-tailed unpaired Student’s t-test.
YUMMER1.7 cells exhibit more anaplerosis in vitro than YUMM1.7 cells. (A)-(E) 13C enrichment of TCA cycle intermediates in cells incubated in [U-13C6] glucose. (F)-(J) 13C enrichment of TCA cycle intermediates in cells incubated in [U-13C5] glutamine. VCS from [U-13C5] glutamine. (F) VCS from [U-13C16] palmitate. In panels (E) and (F), data were normalized to VCS data from YUMM1.7 cells incubated in the same tracer. (F) VCS from [U-13C16] palmitate. In all panels, data were compared by the 2-tailed unpaired Student’s t-test.
Supplementary Figure S3. Mice reach steady-state in the contribution of precursors to citrate within 60 minutes. (A) [U-13C6] glucose, (B) [U-13C5] glutamine, and (C) [U-13C16] palmitate tracers were used. The solid lines denote the mean of n=5 per cell line. No significant differences between 60 and 150 min were observed using the 2-tailed unpaired Student’s t-test.
Supplementary Figure S4. PRECOG data of glucose and fatty acid metabolic enzymes in melanoma patients. (A) Z scores of survival effect due to mRNA expression of glucose metabolic enzymes in melanoma patients from PRECOG database. (B) Z scores of survival effect due to mRNA expression of fatty acid metabolic enzymes in melanoma patients from PRECOG database.
Supplementary Figure S5. Analysis of RNA seq data in melanoma patients from TCGA database. (A) Dot plot of overall survival time and leukocyte infiltration fraction in 369 melanoma patients. Adjusted R2=-0.001507, P value=0.5046. (B) Correlation matrix of leukocyte fraction to amino acid metabolic enzyme and transporter gene expression in 369 melanoma patients from TCGA database. (C) Correlation matrix of cytolytic activity to amino acid metabolic enzyme and transporter gene expression in 159 melanoma patients from TCGA database.
About this article
Cite this article
Zhang, X., Halberstam, A.A., Zhu, W. et al. Isotope tracing reveals distinct substrate preference in murine melanoma subtypes with differing anti-tumor immunity. Cancer Metab 10, 21 (2022). https://doi.org/10.1186/s40170-022-00296-7