Article Text
PDF
Abstract
Background Pancreatic ductal adenocarcinoma (PDAC) is one of the deadliest cancers. Despite the successful application of immune checkpoint blockade in a range of human cancers, immunotherapy in PDAC remains unsuccessful. PDAC is characterized by a desmoplastic, hypoxic and highly immunosuppressive tumor microenvironment (TME), where T-cell infiltration is often lacking (immune desert), or where T cells are located distant from the tumor islands (immune excluded). Converting the TME to an immune-inflamed state, allowing T-cell infiltration, could increase the success of immunotherapy in PDAC.
Method In this study, we use the KPC3 subcutaneous PDAC mouse model to investigate the role of tumor-derived sialic acids in shaping the tumor immune landscape. A sialic acid deficient KPC3 line was generated by genetic knock-out of the CMAS (cytidine monophosphate N-acetylneuraminic acid synthetase) enzyme, a critical enzyme in the synthesis of sialic acid-containing glycans. The effect of sialic acid-deficiency on immunotherapy efficacy was assessed by treatment with anti-programmed cell death protein 1 (PD-1) and agonistic CD40.
Result The absence of sialic acids in KPC3 tumors resulted in increased numbers of CD4+ and CD8+ T cells in the TME, and reduced frequencies of CD4+ regulatory T cells (Tregs) within the T-cell population. Importantly, CD8+ T cells were able to infiltrate the tumor islands in sialic acid-deficient tumors. These favorable alterations in the immune landscape sensitized sialic acid-deficient tumors to immunotherapy, which was ineffective in sialic acid-expressing KPC3 tumors. In addition, high expression of sialylation-related genes in human pancreatic cancer correlated with decreased CD8+ T-cell infiltration, increased presence of Tregs, and poorer survival probability.
Conclusion Our results demonstrate that tumor-derived sialic acids mediate T-cell exclusion within the PDAC TME, thereby impairing immunotherapy efficacy. Targeting sialic acids represents a potential strategy to enhance T-cell infiltration and improve immunotherapy outcomes in PDAC.
- Tumor Microenvironment
- Immunotherapy
- Lymphocytes, Tumor-Infiltrating
Data availability statement
Data are available upon reasonable request.
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See http://creativecommons.org/licenses/by-nc/4.0/.
Statistics from Altmetric.com
WHAT IS ALREADY KNOWN ON THIS TOPIC
Aberrant glycosylation, particularly overexpression of sialic acids (a process called sialylation), is a common feature in various cancer types contributing to immune evasion and tumor progression. Preclinical studies in other cancer types have shown that targeting sialic acids can improve survival and synergize with immunotherapy. However, the impact of tumor sialylation in pancreatic ductal adenocarcinoma (PDAC) progression, immune evasion and immunotherapy efficacy in vivo has remained unexplored until now.
WHAT THIS STUDY ADDS
Our study provides valuable insights into the influence of tumor sialylation on PDAC progression, immune evasion, and immunotherapy efficacy in an in vivo setting. We show that sialic acids in PDAC promote T-cell exclusion and hamper immunotherapy efficacy in vivo. Furthermore, we show that high expression sialylation-related genes in human PDAC is associated with poorer survival probability, decreased CD8+T cell numbers, and increased regulatory T cell frequencies. These results elucidate a novel mechanism underlying T-cell exclusion in the PDAC tumor microenvironment (TME).
HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY
By understanding the role of sialic acids in shaping the immune microenvironment and promoting T-cell exclusion, this research opens up new avenues for the development of targeted immunotherapeutic strategies in PDAC. Targeting tumor sialylation presents a promising approach to enhance the efficacy of immunotherapy in PDAC. Additionally, the results suggest that expression levels of sialylation-related genes could serve as potential biomarkers for predicting immune responses and treatment outcomes in PDAC.
Background
Pancreatic ductal adenocarcinoma (PDAC) is one of the most aggressive cancer types, often diagnosed at an advanced stage with limited treatment options, resulting in a 5-year overall survival of only 11%.1 While immune checkpoint blockade (ICB) has revolutionized cancer therapy in various cancer types, its effectiveness in PDAC remains limited.2 The limited effectiveness of ICB in PDAC is attributed to the relatively low mutational burden and the complex immunosuppressive tumor microenvironment (TME).3 Understanding the underlying mechanisms responsible for the immunosuppressive TME is crucial to overcome the immune inhibition and to develop novel therapeutic strategies that improve the survival of patients with PDAC.
The TME of PDAC is predominantly composed of desmoplastic stroma, which harbors suppressive immune cells like tumor-associated macrophages (TAMs) and regulatory T cells (Tregs), that hinder effective antitumor immune responses.4 Both TAMs and Tregs dampen antitumor immunity through the production of tolerogenic cytokines and by inhibiting the activation and function of CD8+ T cells, natural killer (NK) cells and dendritic cells (DCs).5 6 Consequently, the abundance of Tregs and TAMs is associated with bad prognosis in patients with PDAC.7–10 The immunosuppressive nature of the TME is further shaped by cancer-associated fibroblasts (CAFs), which generate dense fibrotic stroma acting as a physical barrier that restricts immune cell infiltration.4 11 Cytotoxic T cells, the crucial players in antitumor immunity and the effectiveness of ICB, are scarce within the PDAC TME.12 Even in patients with relatively high T cell numbers, T cells are mainly located within the stroma, distant from the tumor islands.3 13 14 As a result, many patients exhibit an immune desert phenotype characterized by limited T cells or an immune-excluded phenotype where T cells are present but excluded from the tumor beds.3 Importantly, studies have shown that the amount of CD8+ T cells and the proximity of T cells to PDAC tumor cells correlates with increased patient survival.9 14–16 Therefore, reversing the immunological “cold” PDAC microenvironment to an “inflamed” state, allowing CD8+ T cells to infiltrate, holds potential to improve the effectiveness of ICB in PDAC.
Aberrant glycosylation is a hallmark of cancer.17–19 Specific alterations in glycan structures contribute to immune evasion.19 20 Sialic acids, located at the terminal side of glycans, are commonly overexpressed on various glycoproteins, including mucins like MUC1, in different cancer types.19 21 Sialic acid on these glycan carriers can interact with sialic acid-binding immunoglobulin-like lectin (Siglecs) receptors expressed on innate and adaptive immune cells, modulating immune responses.22 23 Numerous studies have demonstrated that tumor sialylation inhibits CD8+ T and NK cell cytotoxicity.24–30 Furthermore, activation of the sialic acid-Siglec axis on DCs hampers DC maturation and subsequent priming of T cells and promotes Treg induction.31–33 Additionally, tumor sialylation polarizes TAMs to a more M2-like phenotype with tumor promoting properties.34–38 Notably, in mouse models of melanoma, colorectal, and breast cancer, removal or reduction of tumor sialylation has been shown to result in decreased tumor growth and improved survival,29 34 39 40 which in melanoma is associated with a favorable shift in the CD8+/Treg balance.29 39 Given their immunosuppressive properties, sialic acids are now recognized as a potential target for immunotherapy.20
Preclinical studies targeting tumor sialylation through sialic acid mimetics, sialidase-coupled antibodies, or interference with Siglecs have demonstrated improved survival in mouse models and synergistic effects with ICB in melanoma, colorectal, and breast cancer.34 39–41 While previous research has highlighted the impact of sialylation on tumor progression and immune modulation in other cancer types, its specific implications in PDAC remain largely unexplored. We have previously shown that sialic acids are overexpressed in patients with PDAC and that high expression of sialylation-related genes is associated with shorter survival probability.42 Additionally, tumor sialylation promoted the differentiation of monocytes to M2-like TAMs in vitro.42 Importantly, the significance of targeting sialic acids in cancer therapy is gaining recognition, as evidenced by the ongoing phase 1 trial (NCT05259696) evaluating the first agent specifically designed to target sialic acids, which has remarkably enrolled several patients with PDAC (website: https://palleonpharma.com/education-hub). However, despite these advancements, there is a notable gap in our understanding of the in vivo effects of sialic acids on tumor growth, the immune landscape, and immunotherapy response in PDAC.
The aim of this study was to investigate the influence of tumor sialylation on the in vivo dynamics of PDAC, including tumor growth, the immune cell landscape of the TME and draining lymph nodes, as well as immunotherapy response. By employing a sialic acid-deficient model, our objective was to unravel the interplay between sialylation and the immune landscape of PDAC, ultimately providing valuable knowledge for the application of sialic acid targeted therapeutic strategies in this challenging disease.
Method
Cell lines and cell culture
The murine pancreatic cancer cell line KPC3 was obtained from a primary KPC tumor with mutant p53 and K-ras43 from a female C57BL/6 mouse (RRID:CVCL_A9ZK).44 To generate a sialic acid knock-out (KO) line, KPC3 cells were transfected with a Mock CRISPR/Cas9 construct or a CRISPR/Cas9 plasmid including a gRNA for the Cmas gene. Detailed procedure can be found in the online supplemental methods.
Supplemental material
Glycan profiling
Measurement of glycosphingolipids, N-glycans and O-glycans was done on 10 million cells per KPC3 cell line by PGC nano-LC-ESI-MS/MS. Additional details are included in the online supplemental methods.
Cell titer blue viability assay
The viability of the KPC3 cells was analyzed with CellTiter-Blue Cell Viability assay (Promega) as an indicator of proliferation after 48 hours of culture, and measured on the FLUOstar Galaxy (MTX Lab systems, excitation 560 nm, emission 590 nm).
Immunohistochemistry and immunofluorescent stainings
Fresh frozen tumor tissues were cryosectioned at 10 µm thickness and stored at −20 degrees. Tissue slides were rehydrated and blocked with Carbo-free blocking solution (Vector Labs) before incubation with 2 µg/mL pan-Lectenz (Lectenz Bio) pre-complexed with streptavidin-horseradish peroxidase (HRP). After incubation with pan-Lectenz for 2 hours at room temperature, slides were washed and developed with DAB (DAKO). Slides were dehydrated and mounted with Entallan before brightfield imaging on Vectra Polaris (PerkinElmer).
For immunofluorescent staining, slides were fixed with 1% paraformaldehyde after staining with pan-Lectenz pre-complexed with streptavidin AF488. Slides were incubated overnight at 4°C with anti-CD8 AF555 and anti-EPCAM AF647 in the presence of anti-CD16/CD32 for Fc blocking, washed and mounted with Fluoramount-G (Thermo Fisher Scientific) containing DAPI (1:2000; Invitrogen). All slides were imaged by Vectra Polaris (PerkinElmer).
Flow cytometry
Staining with 2 µg/mL pan-Lectenz (Lectenz Bio) was performed in Hanks' Balanced Salt Solution (HBSS) containing magnesium and calcium (Gibco) supplemented with 0.5% fatty-acid free bovine serum albumin (BSA) (Sigma). The pan-Lectenz reagent was pre-incubated with 1 µg/mL streptavidin-APC, and incubated with the cells after life dead staining with Fixable Viability Dye eFluor 450 (eBioscience). Other antibody stainings for flow cytometry were done in phosphate-buffered saline (PBS)/1% BSA, after life dead staining with Fixable Viability Dye eFluor 450. Fluorescence intensities were measured using the Fortessa X‐20 and analyzed with FlowJo V.10.
In vivo tumor experiments
For tumor experiments, 1×105 KPC3 tumor cells were inoculated subcutaneously (s.c.) in the flank of mice under 2–3% isoflurane. Tumor size was measured three times a week and mice were terminated at day 21 for KPC3 wild-type (WT) tumors, and terminated at day 27 for KPC3 Mock or CMAS KO tumors.
In indicated experiments, mice were treated by intraperitoneal (i.p.) injection of 250 µg antagonistic programmed cell death protein 1 (PD-1) antibody (clone RMP1–14, in house, endotoxin-free) two times a week and with 40 µg agonistic CD40 antibody (clone 1C10) one time per week, starting from day 16 when tumors were at an average size of approximately 50–100 mm3. The mice were sacrificed individually when tumors reached >800 mm3, or when tumors showed beginning ulcer formation to prevent severe discomfort of mice. Untreated mice were i.p. injected with PBS. Additional details have been included in the online supplemental methods.
Tissue processing and immune cell profiling
Tissues were digested and stained with multicolor antibody panels (online supplemental table 1). Samples were acquired on Cytek Aurora and data was unmixed using the SpectroFlo software. The data was analyzed using FlowJo V.10 and Omiq. Detailed description of tissue digestion, antibody staining and immune cell profiling can be found in the online supplemental methods.
Supplemental material
Transcriptomic analysis
The single-cell RNA sequencing data of patients with PDAC published by Peng et al was downloaded from the Genome Sequence Archive project PRJCA001063 as preprocessed row data and imported into the package Seurat (v3) for downstream analysis as described previously.42 45 The function AddModuleScore was used to generate a gene score for the donor synthesis gene set (GNE, NANS, NANP, CMAS, SLC35A1). The data analysis on bulk-sequencing data was performed on 150 pancreatic cancer samples from the PanCanAtlas Publications (https://gdc.cancer.gov/about-data/publications/pancanatlas) using R V.4.2.2. Samples classified as normal, pseudo-normal, or non-PDAC were excluded from the analysis.46 Cell types were estimated by bulk RNA-sequencing deconvolution (Cibersort) in package IOBR V.0.99.9. Enrichment score for the donor synthesis gene set per individual sample was calculated using the gsva method of the GSVA package. Patients were classified between high and low expression based on the median, and survival analysis was performed using the package survival. This analysis was validated in a second data set from Bailey et al.47
Statistical analysis
Statistical analysis was performed in Graphpad Prism V.9.3.1. Comparison of KPC3 Mock tumors with KPC3 CMAS KO tumors was done using two-tailed unpaired Student t-test unless stated otherwise in the legend. The bars in graphs represent the SD and statistical significance was determined by a value of<0.05 (*p≤0.05, **p≤0.01, ***p≤0.001, ****p≤0.0001; ns, not significant).
Results
Establishing the mouse model to study tumor sialylation in vivo
To evaluate the impact of tumor sialylation on various dynamics of PDAC in vivo, we employed murine KPC3 cells, which is a PDAC model induced by KRAS and P53 mutations.43 44 The KPC3 model recapitulates human PDAC features as the tumors contain an abundance of stroma, and tumor islands expressed sialic acids as evaluated by immunohistochemistry using pan-Lectenz staining (figure 1A).To specifically investigate the influence of tumor sialylation, we generated a sialic acid-deficient variant of the KPC3 cells using CRISPR/Cas9 gene editing. This was achieved by knocking-out cytidine monophosphate N-acetylneuraminic acid synthetase (CMAS), an essential enzyme in the sialic acid metabolism pathway (figure 1B). Consequently, the CMAS KO resulted in removal of sialic acids on the membrane of the KPC3 cells both in vitro and in vivo (figure 1C,D). The absence of sialic acids on the cell surface did not impact proliferation of the KPC3 cells, nor did it affect major histocompatibility complex (MHC) class I expression (figure 1E, online supplemental figure S1A).
Supplemental material
CMAS KO in the KPC3 cells eliminates sialic acids from glycoproteins. (A) Sialic acid staining using pan-Lectenz in s.c. KPC3 WT tumor. Scale bar represents 100 um, representative image of n=5. (B) Schematic representation on the role of the CMAS enzyme in the sialylation pathway. (C) Sialic acid expression as evaluated by pan-Lectenz in KPC3 cell lines transfected with Mock construct (KPC3 Mock), or CRISPR/Cas9 construct for CMAS KO (KPC3 CMAS KO). (D) Sialic acid levels on EPCAM+ tumor cells from s.c. KPC3 Mock and KPC3 CMAS KO tumors. Statistical analysis with Mann-Whitney test, n=9 for KPC3 Mock and n=10 for KPC3 CMAS KO. (E) Viability of KPC3 variants after 48 hours of culture as an indicator of proliferation of these cell lines. Data combined from three independent experiments. (F) N-glycan profiles of KPC3 variants determined by mass spectrometry. (G) Glycosphingolipid glycan profiles of KPC3 variants determined by mass spectrometry. (H) O-glycan profiles of KPC3 variants determined by mass spectrometry. Data presented as mean±SD P value represented as ****p≤0.0001. CMAS, cytidine monophosphate N-acetylneuraminic acid synthetase; KO, knock-out; s.c., subcutaneously; WT, wild-type.
In order to validate the KO of sialic acid and to gain insight into the specific glycan structures expressed by the KPC3 cells, a comprehensive glycan profiling analysis was conducted, including the O-glycosylation, N-glycosylation and glycosphingolipid glycosylation profiles. The KPC3 Mock transfected cells exhibited similar glycan profiles to the KPC3 WT cells, indicating that the transfection process did not alter the overall glycan profile of the cell lines (figure 1F–H, online supplemental figure S1B, online supplemental table 2). The analysis revealed that nearly all O-glycans and glycosphingolipid glycans were sialylated in the KPC3 WT and KPC3 Mock cells. The N-glycan profile showed that 20–30% of N-glycans were sialylated in KPC3 WT and Mock cells. Regarding the overall glycan profile of KPC3 WT and Mock cells, we observed that the cells expressed relatively high levels of sialyl-T (ST) and disialyl-T (di-ST) antigens (figure 1H), which are specific glycan structures often found on mucins. The presence of these glycans is associated with tumor progression and poor prognosis in several cancer types.48–51
Supplemental material
In the context of the KPC3 CMAS KO cells, sialylation was significantly reduced or even absent, particularly on N-glycans and O-glycans (figure 1F–H, online supplemental figure S1B). Surprisingly, the glycosphingolipid profile showed distinct features compared with the KPC3 WT and Mock, and although sialylation was reduced on glycosphingolipid glycans, the majority of glycosphingolipids remained sialylated in the KPC3 CMAS KO cells (figure 1G, online supplemental figure S1B). Additionally, the absence of sialic acid led to elongation of O-glycans, a pattern not observed for N-glycans (figure 1F,H).
Together, these data show that KO of CMAS in KPC3 cells results in a sustained lack of sialic acids on tumor cells in vitro and in vivo.
Sialic acid-deficiency in KPC3 leads to a favorable T-cell landscape in the TME
To evaluate the effect of tumor sialylation on KPC3 tumor growth, mice were injected s.c. with the KPC3 cell lines. Loss of tumor sialylation was verified by immunohistochemistry (online supplemental figure S2A). Given the previously reported influence of sialic acid loss on tumor growth in other cancer types,29 34 39 40 we hypothesized that sialic acid-deficiency would result in reduced KPC3 tumor growth as compared with sialic acid-expressing KPC3 tumors. Unexpectedly, the KPC3 CMAS KO cell line, that lack sialylated N-glycans and O-glycans, showed comparable tumor growth to the KPC3 Mock cells in vivo, indicating that tumor sialylation does not influence tumor growth (figure 2A,B, online supplemental figure S2B). Subsequently, we analyzed the cellular composition of the TME, including immune cells, CAFs, and endothelial cells. Interestingly, the KPC3 CMAS KO tumors demonstrated an increased proportion of immune cells, based on CD45 expression within in the TME (figure 2C). Notably, the sialylation KO had no effect on the relative abundance of endothelial cells (CD31+) and CAFs (PDPN+) within the tumor (figure 2C).
Supplemental material
Ablation of tumor sialylation alters the immune composition of the TME in murine PDAC. (A) Tumor size over time of s.c. injected KPC3 Mock and CMAS KO tumors. (B) KPC3 tumor weight at day 27 after injection. (C) Percentage of CD45+ immune cells, CD31+ endothelial cells and PDPN+ CAFs in KPC3 Mock and CMAS KO tumors of A–B analyzed by spectral flow cytometry. (D) tSNE of unsupervised clustering analysis, showing the dimensional reduction and clustering of spectral flow cytometry data on alive cells in KPC3 tumors. (E) Heatmap of cell marker expression across cell clusters shown in D. (F) Percentage of different immune cell types in the different KPC3 tumors gated over alive cells. (G) Relative cell frequencies within the total CD3+ T cells. Statistical analysis performed using two-way ANOVA with Bonferroni multiple comparison test comparing Mock and CMAS within the cell subsets. (H) Percentage of Treg cells over total CD3+ T cells. (I) Representative image of immunofluorescence staining of EPCAM (red) and CD8 (green) in KPC3 tumors. Cell nuclei are stained with DAPI (blue). Data presented as mean±SD. For KPC3 Mock n=9 and for KPC3 CMAS KO n=10. P value represented as *p≤0.05, **p≤0.01, ***p≤0.001, ****p≤0.0001. ANOVA, analysis of variance; CAFs, cancer-associated fibroblasts; cDCs, classical DCs; CMAS, cytidine monophosphate N-acetylneuraminic acid synthetase; DCs, dendritic cells; KO, knock-out; NK, natural killer; PDAC, pancreatic ductal adenocarcinoma; pDCs, plasmacytoid DCs; s.c., subcutaneously; TAMs, tumor-associated macrophages; TME, tumor microenvironment; Tregs, regulatory T cells; t-SNE, t-Distributed Stochastic Neighbor Embedding.
To gain further insight into the distinct immune cell subsets within these tumors, we employed a 24-color panel for spectral flow cytometry, enabling the identification of 19 distinct cell populations (figure 2D,E, online supplemental figure S2C). The majority of immune cells in KPC3 tumors belonged to the myeloid lineage, with the predominant cell type being neutrophils (figure 2D, online supplemental figure S2C,D). In contrast, lymphocytes were present in limited numbers, accounting for less than 10% of all immune cells (figure 2D, online supplemental figure S2D). Remarkably, the CD8+ T-cell population was clustered in two subsets, distinguished by the expression of Ly6C (figure 2D,E, online supplemental figure S2C). While the precise functional distinctions conferred by Ly6C expression on CD8+ T cells remain unclear, Ly6C-positive cells are associated with T-cell homing and have been linked to central memory and effector T-cell subsets.52–54 To determine which cell clusters exhibited significant differences between KPC3 Mock and CMAS KO, we manually gated the identified populations (online supplemental figure S3). Interestingly, the 3 T-cell subsets, CD8+ T, CD8+ Ly6C+ T and CD4+ T cells, showed a significant increase in KPC3 CMAS KO tumors compared with KPC3 Mock tumors (figure 2F). Next, we examined the immunosuppressive cell populations, specifically TAMs and Tregs. Surprisingly, we did not observe significant differences in TAM frequencies or MHC class II expression on TAMs (online supplemental figure S2E). However, a notable decrease in Tregs (CD4+ Foxp3+ T cells) was observed in the sialic acid-deficient tumors when analyzed within the T-cell population, indicating a shift in the CD8+T cell/Treg balance in the KPC3 CMAS KO tumors (figure 2G,H).
Supplemental material
Collectively, these data demonstrate that tumor sialylation does not impact tumor growth, but exerts profound effects on the lymphoid immune landscape within the TME leading to a favorable T-cell profile.
Given the prevalence of T-cell exclusion in PDAC, we aimed to investigate the spatial distribution of the increased number of CD8+ T cells within the tumor. To accomplish this, we performed immunofluorescent staining on KPC3 Mock and KPC3 CMAS KO tumors. Interestingly, CD8+ T cells in the KPC3 CMAS KO tumors were localized in close proximity to the tumor islands in comparison to KPC3 Mock tumors (figure 2I). These findings provide evidence that absence of tumor sialylation can have the potential to reverse the immune-excluded phenotype found in PDAC.
Tumor sialylation modulates the draining lymph node
Tumor cells can secrete sialylated carbohydrate antigens, such as the CA19-9 antigen (also known as sialyl Lewis A), which serves as a well-known biomarker for pancreatic cancer.55 These sialylated antigens have the potential to drain to the lymph node and modulate immune function. Previous research has demonstrated that a reduction in tumor sialylation in melanoma can alter immune cell frequencies in the draining lymph node.29 Therefore, we investigated the possible contribution of KPC3-derived sialic acids to immune suppression beyond the local tumor site, by examining immune subsets and co-stimulatory markers in both the tumor draining and non-draining lymph nodes (figure 3A,B, online supplemental figure S4A,B). The non-draining lymph nodes did not exhibit any differences in cell populations or activation states between KPC3 Mock and CMAS KO bearing mice (figure 3A–E, online supplemental figure S4C). However, within the tumor-draining lymph node, we observed an increase in B cells in the KPC3 CMAS KO tumor-bearing mice compared with KPC3 Mock tumor-bearing mice, implying an enhanced B-cell response (figure 3C). Furthermore, while the decrease in T-cell percentages did not reach statistical significance, it suggests a possible migration of T cells from the lymph node to the tumor site (figure 3D). Moreover, a reduction in CD4+ CD25+ T cells, potentially corresponding to Tregs, was observed, indicating a diminished suppressive immune response within the draining lymph node of KPC3 CMAS KO bearing mice (figure 3D). Analysis of the antigen-presenting cells revealed higher expression levels of the co-stimulatory molecule CD80 and the inhibitory molecule programmed death-ligand 1 (PD-L1) in KPC3 CMAS KO bearing mice, specifically in the CD11bmidSirpa+ subset (figure 3E). Percentages of other antigen presenting cell subsets were not affected by the CMAS KO (online supplemental figure S4C). Taken together, these results indicate that the absence of tumor-derived sialic acids selectively affects immune responses within the draining lymph node, suggesting a localized rather than systemic impact on immune regulation by tumor-derived sialic acids.
Supplemental material
Tumor sialylation alters immune cell populations in tumor-draining lymph nodes. (A) Percentage of CD45+ immune cells in draining and non-draining lymph nodes of KPC3 tumor-bearing mice. (B) Dimensional reduction of CD45+ immune cells in the lymph node based on spectral flow cytometry data. (C) Percentage B cells in the lymph nodes of KPC3 tumor-bearing mice. (D) Percentage of lymphoid subsets cells in the lymph nodes of KPC3 tumor-bearing mice. (E) Expression of CD80 and PD-L1 on antigen-presenting cell subsets in lymph nodes of KPC3 tumor-bearing mice. Data presented as mean±SD. Statistical analysis in C–E was done using two-way ANOVA with Bonferroni multiple comparison test comparing Mock and CMAS within the draining and non-draining lymph nodes. P value represented as *p≤0.05, **p≤0.01, ***p≤0.001, ****p≤0.0001. Total n=10 per group, lymph nodes of two mice were pooled after harvesting. ANOVA, analysis of variance; APC, antigen-presenting cell, cDCs, classical dendritic cells; CMAS, cytidine monophosphate N-acetylneuraminic acid synthetase; KO, knock-out; MFI, mean fluorescent intensity; MHC, major histocompatibility complex; PD-L1, programmed death-ligand 1.
Sialic acid-deficiency synergizes with immunotherapy
Sialic acid deficient tumors showed a favorable T-cell landscape characterized by enhanced CD8+ T cells and decreased Tregs. Despite this, their tumor growth remained comparable to sialic acid-expressing tumors, suggesting the presence of additional local immune cell constraints within the TME. Consequently, we hypothesized that immunotherapy could potentially augment T-cell activity in the KPC3 CMAS KO tumor, resulting in reduced tumor growth. Supporting this hypothesis, we observed the expression of PD-1 on both CD8+ and CD4+ T cells, with elevated expression on CD4+ T cells in the KPC3 CMAS KO tumors (figure 4A). Given the ineffectiveness of ICB alone in PDAC, and the key role of myeloid cells in immune suppression, we aimed to target both T cells and myeloid cells using combination treatment of anti-PD-1 and agonistic CD40 antibodies.56 57 Remarkably, the application of immunotherapy effectively reduced tumor growth and significantly improved the survival in the KPC3 CMAS KO bearing mice, while no discernible effect was observed in KPC3 Mock tumors (figure 4B,C). These findings illustrate that the lack of sialic acids synergizes with immunotherapy in murine PDAC, overcoming the intrinsic resistance of PDAC to ICB.
Sialic acid-deficiency synergizes with immunotherapy. (A) Percentage of PD-1 expressing cells within lymphoid populations in KPC3 tumors. Data presented as mean±SD. For KPC3 Mock n=9 and for KPC3 CMAS KO n=10. P value represented as ***p≤0.001. (B) Individual tumor growth curves and tumor growth factor per mouse, treated or not with anti-PD-1 and agonistic CD40 antibodies (aPD-1/aCD40). KPC3 Mock or KPC3 CMAS KO tumors were established and treated at day 16. Tumor growth factor was calculated by fitting a linear curve (γ=αx+β) to the tumor growth measurements from day 16 until the day of sacrifice, and plotting the α (ie, average growth per day). Data presented as mean±SD, with n=16–22 per group (Mock n=17, Mock – aPD-1/aCD40 n=16, CMAS KO n=22, CMAS KO – aPD-1/aCD40 n=19), and statistical analysis was performed with Mann-Whitney test. P value represented as **p≤0.01. (C) Survival curves of the mice from B. Statistical analysis was performed with log-rank (Mantel-Cox) test. Statistical significance was reached when p≤0.0167 (Bonferroni corrected) as marked by *. CMAS, cytidine monophosphate N-acetylneuraminic acid synthetase; KO, knock-out; PBS, phosphate-buffered saline; PD-1, programmed cell death protein 1.
High expression of sialylation genes correlates with decreased CD8+ T cells and increased Treg cell numbers in human PDAC
In order to extend the relevance of our findings from murine PDAC to human PDAC, we analyzed the expression of sialylation-related genes and their correlation with immune cell abundance and survival probability. Specifically we examined the genes involved in the donor synthesis pathway of sialylation. This pathway is critical for the biosynthesis of sialic acid precursors which can be attached to glycoproteins and glycolipids, and showed to be highly expressed in tumor cells (figure 5A), thus serving as a potential surrogate marker for tumor sialylation. Interestingly, the expression of donor synthesis genes showed the most significant correlation with CD8+ T-cell abundance, with high expression of donor genes associated with low amounts of CD8+ T cells (figure 5B, online supplemental figure S5A). In addition, patients with enriched expression of donor synthesis genes showed worse overall survival probability and disease-free survival probability (figure 5C, online supplemental figure S5B). Among the analyzed genes, CMAS exhibited the most pronounced expression in tumor cells compared with other cell types (figure 5D). Importantly, CMAS expression showed a significant positive correlation with the abundance of Tregs (figure 5E, online supplemental figure S5C). This observation suggests that elevated tumor sialylation in human PDAC may be associated not only with reduced CD8+ T-cell infiltration but also with an increased number of Tregs. Collectively, these findings show the potential impact of tumor sialylation on survival probability, lymphocyte abundance and immune regulation in the context of human PDAC.
Supplemental material
Sialylation gene expression correlates with T-cell abundance in human PDAC. (A) Expression of donor synthesis genes, based on a gene score, in different cell subsets in PDAC, analyzed using the single-cell RNA sequencing data set of PDAC patients from Peng et al.45 (B) Spearman correlation between donor synthesis gene score and CD8+ T-cell abundance as identified by tumor microenvironment deconvolution (Cibersort) in human PDAC tumors, analyzed in bulk sequencing data (TCGA). (C) Kaplan-Meier curve showing the survival probability of patients with PDAC from TCGA data set, stratified by high versus low expression of donor synthesis genes. Statistical analysis performed using Cox Regression. (D) Expression of CMAS in different cell subsets in PDAC, analyzed using the single-cell RNA sequencing data set of patients with PDAC from Peng et al.45 (E) Spearman correlation between the expression of CMAS enzyme and the estimated abundance of Treg cells as identified by tumor microenvironment deconvolution (Cibersort) in human PDAC tumors, analyzed in bulk sequencing data (TCGA). CMAS, cytidine monophosphate N-acetylneuraminic acid synthetase; PDAC, pancreatic ductal adenocarcinoma; TCGA, The Cancer Genome Atlas; Tregs, regulatory T cells.
Discussion
The application of immunotherapy in PDAC has encountered substantial challenges due to the immunosuppressive TME, limited tumor immunogenicity, and the development of resistance mechanisms that hinder the efficacy of immune checkpoint inhibitors. CD8+ T-cell infiltration, which is crucial for the effectiveness of ICB, is typically limited in PDAC tumors. In this study, we aimed to investigate the role of tumor sialylation, a novel potential target for immunotherapy, in tumor growth and modulation of the immune TME. Our data demonstrates that sialic acid deficiency can convert the TME to a state that allows for enhanced numbers of CD8+ T cells infiltrating the tumor. Consequently, by reversing the immune excluded phenotype, we observed a significant improvement in the survival of tumor-bearing mice on treatment with immunotherapy.
Previous preclinical studies in other cancer types consistently demonstrated that reducing sialic acid expression on tumor cells decreased tumor growth and improved survival.29 34 39 40 In contrast to these reports, our findings revealed that loss of sialic acids by itself did not affect tumor growth in murine PDAC (figure 2A,B). This discrepancy can be attributed to several factors specific to PDAC. First, PDAC exhibits low levels of tumor antigens as compared with other tumor types, limiting the recognition of tumor cells by the immune system.3 Second, the immunosuppressive environment creates formidable barriers to reduce tumor growth through immune modulation.4 Moreover, a notable distinction between PDAC and other tumor types is the presence of dense fibrotic stroma, which profoundly shapes the TME and contributes to disease aggressiveness.4 Targeting both the tumor and stroma may be necessary for therapeutic success in PDAC, and this concept can extend to the manipulation of sialic acids. Notably, stromal cells themselves have been reported to express sialic acids,58 indicating that solely targeting tumor cell sialylation may not be sufficient to inhibit tumor growth in PDAC, necessitating the consideration of stromal sialylation as well.
Although the lack of tumor sialylation did not directly impact tumor growth, our study revealed significant alterations in local immune compositions within the TME and tumor-draining lymph nodes. In the TME, we observed an augmentation of tumor-infiltrating lymphocytes and a decrease in Tregs, indicating a potential shift towards an immune-inflamed state. Notably, classical DCs (cDC), which are essential for an effective antitumor T-cell response and immunotherapy response,59 60 did not exhibit significant alterations in their subsets or co-stimulatory molecule expression in the draining lymph node. This suggests that the increased lymphocyte presence in the KPC3 CMAS KO tumors may involve mechanisms independent of priming by cDCs in the draining lymph node.
Sialic acids have been widely investigated for their involvement in driving macrophage polarization towards a tumor-promoting phenotype.34–38 41 42 61 However, in contrast to these findings, our study did not observe significant alterations in frequencies or MHC class II expression of TAMs in sialic acid-deficient tumors. This finding was unexpected, considering that KPC3 cells express sialylated glycans such as the ST-antigen (figure 1H), which has been reported to influence myeloid differentiation to a TAM phenotype in the context of breast cancer.36 Furthermore, the sialylated glycans on KPC3 resemble the glycan profile of human PDAC cell lines which also express ST and di-ST,17 and we have previously demonstrated that PDAC cell lines modulate TAM differentiation through sialic acids in vitro.42 Besides, the effects of sialic acids on tumor growth and immune modulation have been shown to be dependent on mouse Siglec-E, primarily expressed on myeloid cells.34 Although T cells can express Siglec-E in certain settings,27 28 we did not observe its expression on T cells in our study (data not shown). Thus, while the frequencies of myeloid populations and MHC class II expression on TAMs remained unchanged, it is plausible that functional alterations in myeloid cells, potentially mediated through sialic acid-Siglec interactions, contributed to the establishment of a more inflamed TME that influenced T-cell infiltration.
Interestingly, sialic acid-deficiency made our KPC3 tumor responsive to immunotherapy (figure 4B). Surprisingly, the treatment with anti-PD-1 and agonistic CD40 antibodies was shown to be ineffective in the KPC3 Mock tumors, despite its reported therapeutic efficacy in other PDAC studies.57 62 63 The lack of response to immunotherapy is often associated with the limited neo-antigen expression in PDAC or inadequate immune cell infiltration.3 The inadequate immune cell infiltration, particularly of CD8+ T cells, was reversed by the absence of tumor sialylation in the KPC3 CMAS KO, which could underlie the responsiveness to immunotherapy in these tumors. Nevertheless, although immunotherapy significantly improved survival and delayed tumor growth in the KPC3 CMAS KO tumor-bearing mice, the magnitude of this effect remained small and no mice were cured. Combining immunotherapy with sialic acid removal in a setting with sufficient tumor-antigen specific T cells, for example, by vaccination strategies, may further enhance the therapeutic efficacy.
The amount of T cells and their proximity to tumor cells is known to be associated with better survival outcomes in patients with PDAC.9 14–16 Consequently, targeting sialic acids in patients with PDAC offers the potential to augment T cell numbers and overcome the immune excluded phenotype of the TME. In recent years, several strategies have emerged for targeting sialic acids, including tumor-targeted degradation of sialic acids and sialic acid mimetics that inhibit sialyltransferase activity.39 40 64 A phase 1 clinical trial, which includes patients with PDAC, is currently investigating the safety and potential of sialidase treatment (NCT05259696). Initial results of this trial with systemic administration of the bi-sialidases show that the compound was well tolerated with no dose limiting toxicities. These therapeutic approaches represent promising avenues for combination therapy in PDAC, particularly when combined with agents that enhance T-cell function.
In this study we use a complete and durable KO of sialic acids to study the role of sialic acids in PDAC. However, in patients, we observed heterogeneity in donor synthesis gene expression including CMAS, and interestingly, we found a correlation between these gene levels and the abundance of CD8+ T cells and Tregs. This suggests that CMAS expression levels could potentially serve as a biomarker for a favorable T-cell landscape. Taking into account the absolute deficiency of sialic acids compared with the varying levels in humans, it is crucial to consider the consequences of sialylation loss on underlying glycan structures. This loss exposes terminal galactose residues (figure 1H), which can serve as binding sites for galectins, introducing an additional layer of complexity. For instance, it has been described that silencing of ST6GALNAC2 decreases di-ST expression and creates binding sides for galectin-3, which in turn leads to increased metastatic burden in breast cancer.65 Modulating sialylation levels in patients to decrease their sialic acid expression may offer a more favorable approach compared with the complete removal of sialic acids. Next to potential binding of galectins, loss of tumor sialylation may also affect tumor cell survival, invasion or metastasis.66
In conclusion, our study demonstrates that the absence of sialic acids, achieved through genetic depletion of CMAS, leads to favorable alterations in the immune landscape of pancreatic tumors and the draining lymph node. Specifically, we observed increased frequencies of CD4+ helper and CD8+ cytotoxic T cells, accompanied by a reduction in Tregs within the T-cell population in the TME. Notably, mice with sialic acid-deficient tumors showed improved survival in response to immunotherapy, while those with sialylated tumors did not exhibit survival benefits. Importantly, our findings extend to the clinical context, as higher expression of sialylation genes in human pancreatic cancer is associated with decreased CD8+ T-cell infiltration and increased Treg presence, and worse patient outcome. These results highlight the critical role of tumor sialylation in shaping the tumor immune microenvironment and provide insights for the development of targeted immunotherapeutic strategies in pancreatic cancer.
Supplemental material
Data availability statement
Data are available upon reasonable request.
Ethics statements
Patient consent for publication
Ethics approval
All animal experiments were approved by the Animal Experiments Committee of the VU University (DEC201512, protocols 1024-MCB20-45; 1024-MCB20-58; 1024-MCB20-59; 1024-MCB21-66; 1024-MCB21-67) and performed in accordance with national and international guidelines and regulations
Acknowledgments
We acknowledge the Microscopy and Cytometry Core Facility at the Amsterdam UMC – Location VUmc for providing assistance in our microscopy and cytometry work. In particular, we would like to acknowledge Liesbeth Paul for the design and advice on Spectral flow cytometry panels.
References
Supplementary materials
Supplementary Data
This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.
Footnotes
Contributors KB and YVK conceived the study. KB, LK, DW, CMdW, IvdHA, AdH, LB, REM, NvM, MW, SVV, YVK designed experiments or provided essential technical support. KB, LK, BS, LW acquired experimental data. KB, CMdW, BS analyzed data. DW, MW performed glycan profiling. KB, DL, ER performed transcriptomic analysis. KB and YVK drafted the manuscript. KB, LK, DW, CMdW, DL, ER, BS, IvdHA, AdH, LB, REM, NvM, MW, JMMDH, SVV, YVK provided critical intellectual content. YVK supervised the study and is responsible for the overall content as guarantor.
Funding This work is financially supported by KWF VU2014-7200 to KB; SPINOZANWO SPI-93-538 to KB, LK, LW and YVK; LSH-TKI project DC4Balance LSM1806-SGFLSH-TKI to AdH; CMdW is supported by Cancer Center Amsterdam (grant no. CCA2019-9-57 and CCA2020-9-73).
Competing interests None declared.
Provenance and peer review Not commissioned; externally peer reviewed.
Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.