- 1Senior Dept of Oncology, The Fifth Medical Center of People's Liberation Army (PLA) General Hospital, Beijing, China
- 2Dept of Oncology, The Forth Medical Center of People's Liberation Army (PLA) General Hospital, Beijing, China
- 3Department of Research and Development (R&D), Hangzhou Repugene Technology Co., Ltd., Hangzhou, China
Introduction: Tumor-infiltrating myeloid cells (TIMs) are key regulators in tumor progression, but the similarity and distinction of their fundamental properties in pancreatic ductal adenocarcinoma (PDAC) remain elusive.
Method: In this study, we conducted scRNA-seq data analysis of cells from 12 primary tumor (PT) tissues, 4 metastatic (Met) tumor tissues, 3 adjacent normal pancreas tissues (Para), and PBMC samples across 16 PDAC patients, and revealed a heterogeneous TIMs environment in PDAC.
Result: Systematic comparisons between tumor and non-tumor samples of myeloid lineages identified 10 necroptosis-associated genes upregulated in PDAC tumors compared to 5 upregulated in paratumor or healthy peripheral blood. A novel RTM (resident tissue macrophages), GLUL-SQSTM1- RTM, was found to act as a positive regulator of immunity. Additionally, HSP90AA1+HSP90AB1+ mast cells exhibited pro-immune characteristics, and JAK3+TLR4+ CD16 monocytes were found to be anti-immune. The findings were validated through clinical outcomes and cytokines analyses. Lastly, intercellular network reconstruction supported the associations between the identified novel clusters, cancer cells, and immune cell populations.
Conclusion: Our analysis comprehensively characterized major myeloid cell lineages and identified three subsets of myeloid-derived cells associated with necroptosis. These findings not only provide a valuable resource for understanding the multi-dimensional characterization of the tumor microenvironment in PDAC but also offer valuable mechanistic insights that can guide the design of effective immuno-oncology treatment strategies.
Highlights
• ScRNA-seq revealed a heterogeneous tumor-infiltrating myeloid (TIMs) environment in PDAC.
• 10 necroptosis-associated genes were found to be upregulated in PDAC tumors, while 5 genes were upregulated in paratumor or healthy peripheral blood.
• A novel RTM subset, GLUL-SQSTM1- RTM, was found to act as a positive regulator of immunity.
• HSP90AA1+HSP90AB1+ mast cells were identified to be pro-immune, while HSP90AA1-HSP90AB- mast cells were anti-immune.
• JAK3+TLR4+ CD16 monocytes were found to have anti-immune characteristics, while JAK3-TLR4- CD16 monocytes displayed pro-immune properties.
Introduction
Pancreatic cancer is a highly malignant tumor of the digestive system, with the most common subtype being pancreatic duct adenocarcinoma (PDAC), and its morbidity and mortality are increasing year by year worldwide (1). The insidious primary lesion determines that more than 80% of cases cannot undergo surgical resection due to regional or distant metastasis, and the postoperative recurrence rate is as high as 85% for resectable patients (2). For radiotherapy and chemotherapy, the mainstay therapeutic strategy, resistances are usually inevitable, leading to limited clinical benefits, especially for advanced patients. Recently, researchers reported the therapeutic potential of immunotherapy by recruiting and activating the host’s T cells to recognize tumor-specific antigens, however, cancer cells developed mechanisms to escape the cytotoxicity effect of T cells. Moreover, the unique immunosuppressive microenvironment of pancreatic cancer hinders the promotion of anti-tumor immune responses through immune checkpoint manipulation (3–5). Therefore, it is urgent to explore novel therapeutic strategies that can significantly improve patient survival and prognosis (6).
Investigations of converting immune “cold” tumors into “hot” tumors are ongoing in immunotherapy. Regulated cell death (RCD), one of the hallmarks of cancer, has been identified as a potential therapeutic target due to its association with anti-tumor immunity. Necroptosis, a recently proposed form of cell death proposed (7), can be inhibited by Necrostain-1 and mediated by receptor-interacting serine/threonine protein kinase 1 (RIPK1) (8). Necroptosis often triggers a robust inflammatory response by releasing cellular contents into the extracellular environment, and this process has been implicated in the pathogenesis and progression of various diseases (9). It is reported that the impact of necroptosis on cancer development, whether inhibitory or promotive, often varies depending on the specific tumor type and stage (10).
The involvement of necroptosis in dysregulated tumor immune microenvironment (TIME) has been demonstrated, especially for myeloid lineages. For instance, damage-associated molecular patterns (DAMPs) were released by tumor cells through necroptosis to stimulate the antigen presentation by dendritic cells (DCs), further enhancing the cytotoxicity of CD8+ T cells (11, 12). RIPK3, the effector of necroptosis, contributes to NF-κB activation, tissue repair of DC cells, and infiltration of CD8+ T cell (13, 14). Previous studies have reported that the anti-tumor immune response can be activated by NF-κB signaling via necroptosis of fibroblasts (15, 16). However, the regulatory mechanisms of necroptosis in tumor progression in other studies seem to differ from the aforementioned processes. RIPK1 was found to be upregulated in tumor-associated macrophages (TAMs) during M2 Macrophages polarization in a PDAC mice model (17). Necroptosis mediated by RIPK3 promoted the accumulation of immunosuppressive myeloid-derived suppressor cells (MDSCs) in tumor microenvironment (TME) of pancreatic cancer through producing C-x-c motif chemokine ligand 1 (CXCL1) and CXCL5 (18, 19). In an intestinal tumor model, RIPK3 in intermediate MDSC subpopulation was found to increase tumor size (20). Taken together, these findings suggest that myeloid clusters and associated necroptosis may play critical roles in tumor progression and immune evasion.
The exploration of necroptosis-associated myeloid subpopulation can provide a better understanding of the mechanisms underlying immune evasion and therapy resistance in PDAC. Recently, single-cell transcriptomic has made remarkable breakthroughs in deciphering the heterogeneity at the individual cell level. Accumulating evidence has demonstrated the abundance of myeloid cells in tumor immune microenvironment (TIME) of PDAC, serving as key regulators in immune response and treatment resistance (21). By refining the clustering of tumor-associated macrophages (TAMs) in human and mouse samples, researchers have identified significant upregulation of proliferating tissue-resident macrophages and inflammatory macrophages in PDAC TIME received chemotherapy. Conversely, monocyte-derived antigen-presenting cells (APCs) and Marco+ macrophages highly expressed the scavenger receptor MARCO, showed decreased expression. Results from multiplex immunohistochemistry (mIHC) further supported the chemotherapy resistance of proliferating tissue-resident macrophages (22). The deficiency of DCs has been linked to dysfunctional T cell-mediated immunity in early-stage PDAC, indicating their vital role in immune escape and tumor progression (23). Although some progress have been made, further detailed characterization of myeloid cell lineage is needed, and the therapeutic application of myeloid cells in pancreatic cancer remains limited.
This study aims to further elucidate the unique microenvironment of PDAC, explore its intrinsic mechanisms in the tumor occurrence and progression, and provide a potential novel approach for the treatment of PDAC patients. Leveraging a publicly available scRNA-seq resource (24), we revealed a tumor-associated myeloid environment in PDAC. Specifically, we identified upregulated necroptosis genes and immune-related novel clusters in PDAC. Furthermore, we discovered cell-specific signaling pathways and receptor-ligand pairs within these new clusters, which have the potentially to either promote or suppress tumor development. In general, utilizing this unique resource, we analyzed myeloid cell lineages, necroptosis-associated networks, and cell-cell crosstalk in PDAC. This sheds light on the myeloid ecosystems underlying PDAC initiation and progression, and may provide a myeloid-modulating therapeutic strategies from pre-clinical models to pancreatic cancer treatment.
Methods
Data source and preprocessing
The PDAC dataset GSE155698 (24) was downloaded from the GEO database, including 12 primary tumor (PT) patients, 4 metastatic (Met) patients, 3 adjacent normal pancreas tissue (Para) patients, and all samples were coupled with peripheral blood (Figure 1A). The original dataset contained a total of 25,236 genes and 142,353 cells. The raw UMI count matrices were processed using the R package Seurat (version 3.2.3) (25). The data underwent several filtering steps: 1) cells with a low number of unique detected genes (< 200) and a high number of 5000 were removed; 2) cells with more than 30,000 UMIs were discarded; 3) cells with mitochondrial content higher than 30% were removed; 4) cell cycle genes were regressed out. After excluding low-quality cells, 124,575 single cells remained for downstream analysis. Additionally, another publicly available scRNA-seq data from CRA001160 (26), including a total of 57,539 cells from 24 primary PDAC tumors and 11 control pancreases, was utilized to validate the findings.
Figure 1 Dissection of the tumor microenvironment in PDAC. (A) Samples were collected from GSE155698. (B) Visualization of single-cell RNA-seq data of 124,575 cells by t-SNE. (C) Single-cell resolution heatmap of top expressed genes for each cell type. (D) Proportions of 11 cell types among 39 samples including 16 primary or metastatic tumor tissues coupled with 16 PDAC PBMC samples, 3 paratumor tissues and 4 healthy PBMC samples.
Public clinical data and gene expression information were retrieved from the TCGA database (https://www.cbioportal.org/). A total of 178 samples from the TCGA-PAAD (Pancreatic adenocarcinoma) cohorts were included for further analysis.
Sing-cell RNA-seq data clustering and dimensional reduction
First, we performed data normalization on the merged data using the NormalizeData function and identified the first 2000 highly variable genes through the FindVariableFeatures function, which is based on the variance stabilization transformation (“vst”). Simultaneously, all genes were scaled using the ScaleData function, and the RunPCA function was applied to reduce the dimensionality of the data using PCA for previously identified highly variable genes. We selected a dimensionality reduction of 30 (dim = 30) and clustered the cells using the FindNeighbors and FindClusters functions with a resolution of 1.2, enabling the identification of distinct cell clusters. To further reduce dimensionality and visualize the data, we employed the UMAP and tSNE methods using the top 30 principal components. Specifically focusing on myeloid cell types, we repeated the clustering protocol to identify clusters within the aforementioned myeloid lineages. To address batch effects, we utilized the runHarmony function from the Harmony package (version 0.1.0) (27). Finally, we employed the FindAllMarkers function to screen the marker genes of 40 subgroups, considering a log-fold change (logfc) threshold of 0.25 for differential expression and a minimum percentage (min.pct) of 0.25 for the expression ratio of the least differentially expressed genes. We applied a corrected p threshold of less than 0.05 to screen the significant marker genes.
Cell type identification
The annotation of each cell cluster was confirmed by the expression of canonical marker genes. Epithelial cells were identified using the higher expression of EPCAM, ACTA2, KRT7, KRT8, KRT18, KRT19, CDH1, PRSS1, CTRB2, REG1A, CLU, MKI67, SPINK1, TFF1, and MUC1, and other cell types were annotated using: T cells (CD3D, CD3E, CD3G, CD4, CD8A, IL7R, and LEF1), B cells (MS4A1, CD79A, CD79B, CD52, CD19, SDC1, IGJ, IGLL5, CXCR4, KIT, CD27, and HLA-DRA), NK (natural killer) cells (NCR3, FCGR3A, NCAM1, KLRF1, KLRC1, and CD38), acinar cells (PRSS1, CTRB1, CTRB2, REG1B, SPINK1, and AMY2A), mast cells (TPSAB1, and CPA3), fibroblast (LUM, DCN, COL1A1, ACTA2, SPARC, CDH11, PDGFRA, PDGFRB, COL3A1, RGS5, IGFBP7, PDPN, MCAM, IL6, APOE, GLI1, GLI2, GLI3, and PDGFA), myeloid cells (CD14, ITGAM, MNDA, MPEG1, ITGAX, FCGR3A, FCGR3B, APOE, C1QA, MARCO, LYZ, and HLA-DRA), stellate cells (RGS5, ACTA2, PDGFRB, and ADIRF), and endothelial cells (CDH5, PLVAP, VMF, VLDN5, KDR, and PECAM1).
Among myeloid cells, cell clusters were identified using genes previously reported. Mast cells were identified by the high expression of KIT, CPA3, and TPSAB1, and other myeloid lineages were annotated using: granulocytes (CXCR2, FCGR3B, IFTIM2, SLC25A37, IL1R2, CXCR1, SIRPA, and S100A8), macrophages (C1QC, C1QA, APOE, CCL4, PLTP, and IL1B), CD14 monocyte (CD14, FCN1, S100A8, and S100A9), CD16 monocyte (FCGR3A, LST1, and LILRB2), and DCs (dendritic cells) (IL7R, CCR7, GZMB, LYZ, IL3RA, and IL32).
Tissue distribution of clusters
We quantified the tissue preference of each cluster by calculating the ratio of observed to expected cell numbers (Ro/e) in different tissue (28, 29). The expected cell numbers for each combination of cell clusters and tissues were determined using the chi-square test. A cluster was considered enriched in a specific tissue if Ro/e value was greater than 1.
Differential expression and pathway analysis
To identify differentially expressed genes between two groups of clusters, we used the Wilcox method in the FindMarkers function in Seurat to evaluate the significance of each gene, with multiple hypothesis correction using the Benjamini-Hochberg procedure. Genes with adjusted P-values less than 0.05 were considered as differentially expressed genes. In addition, the log2 fold change (log2FC) for each gene was calculated by subtracting the log2 transformed mean count in each group.
KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment was performed using the clusterProfiler package (version 3.14.3) (30), with a Benjamini-Hochberg multiple testing adjustment. Gene sets with FDR-corrected P < 0.01 were considered to be significantly enriched.
TCGA data analysis
Two endpoints (overall survival (OS) and disease-free interval (DFI) from the TCGA-PAAD) were used to analyze patients’ clinical outcomes. We employed the Cox proportional hazards model implemented in the survival package to assess the correlation between selected genes and patients’ survival. Kaplan-Meier survival curves were generated using the R function ggsurvplot from the survminer package.
Specifically, to examine the relationship between clusters and patients’ survival, we utilized their signature genes. The signature genes including the following categories: macrophage (C1QC, C1QA, APOE, MACRO, INHBA, IL1RN, CCL4, NLRP3, EREG, IL1B, LYVE1, PLTP, SEPP1), granulocyte (FCGR3B, IFITM2, CXCR2, S100A8, SLC25A37, CXCR1, IL1R2), CD14 monocyte (CD14, FCN1, S100A8, S100A9, S100A12, VCAN, CD36), CD16 monocyte (FCGR3A, LST1, LILRB2, IFITM2, SIGLEC10, CX3CR1, LILRB1, LIBRA1, TCF7L2, MTSS1, RHOC), DCs (GZMB, JCHAIN, MZB1, CLIC3, CXCL8, IL7R, CCR7, MMP7, and IL32), and mast cells (KIT, CPA3, TPSAB1, HDC, GATA2, HPGDS, TPSD1, SLC18A2, MS4A2, IL1RL1, and VWA5A). The mean expression of the signature genes was used to classify samples into a high and low groups based on risk score (high: risk score > 0; low: risk score ≤ 0). The Cox model was employed to adjust for OS and status in the survival analysis.
Cell-cell interaction analysis
To investigate the potential interactions between different cell types in the TME of PDAC, we conducted cell-cell interaction analysis using CellChat (v1.1.3), which integrates a curated repository of ligand-receptor (L-R) pairs and employs a statistical framework (31). We combined CD4+T cells, CD8+T cells, B cells, NK cells, and epithelial cells with the newly identified clusters, including GLUL-SQSTM1- RTM, GLUL+SQSTM1+ Macro., HSP90AA1+HSP90AB1+ Mast and JAK3+TLR4+ Mono. Interactions networks between cell clusters were investigated.
RNA fluorescence in situ hybridization
The samples used for RNA FISH were obtained from tumor paraffin sections of patients diagnosed with PDAC. Isolated cancer-associated cells were adhered onto laminin coated #1 coverslips (ThermoScientific) were fixed for 10 min at room temperature with Fixation Buffer (3.7% formaldehyde in PBS), washed twice in 1x PBS and permeabilized with 70% EtOH at 4°C for at least an hour. RNA FISH was performed using 20-mer Stellaris Biosearch Probes for LINCMs and core gene conjugated to Quasar 670 or CAL Fluor Red 610. Briefly, cells were washed with Wash Buffer (10% formamide in 2x SSC) prior to overnight 37°C hybridization with target probes (125 nM) in Hybridization buffer (100 mg/ml Dextran Sulfate, 10% Formamide in 2x SSC). After hybridization, cells were washed in Wash Buffer for 30 min at 37°C, counterstained with DAPI (5 ng/ml in Wash Buffer) for 30 min at 37°C, and washed in 2x SSC at room temperature. Coverslips were transferred onto glass slides with mounting medium (Vectashield) and imaging was performed immediately on upright microscope (Nikon, Ni-E) with 100x Objective (Nikon) on a cooled CCD/CMOS camera (Qi-1, Qi-2, Nikon).
For the notable exception of S100A4, SQSTM1 and GLUL RNA FISH co-staining, RNA FISH was performed using 50-mer ZZ ACD RNAScope probes due to the short unique sequence of the antibodies available for probe design. Cells were fixed and permeabilized as described above in 70% EtOH, washed in 1x PBS and 1x Hybwash buffer for 10 and 30 min, respectively. They were then incubated with 1x Target Probe Mix at 40°C for 3 hours. Cells were washed thrice in 1x Hybwash at room temperature, incubated in 1x Pre Amp Mix for 40 min at 40°C, washed thrice in 1x Hybwash at r.t.p, incubated in 1x Amp Mix for 30 min at 40°C, washed twice in 1x Hybwash before incubation in 1x Label Probe Mix (Alexa Fluo 488, ATTO0550) at 40°C for 25 min. Cells were washed thrice in 1x Hybwash in dark at r.t.p, counterstained with DAPI (5ng/ml) prior to mount and imaging.
Statistical analysis
All statistical analyses were conducted using R software. Comparisons between two groups of samples were evaluated using Wilcoxon rank-sum test (Mann-Whitney U-test) for statistical analysis. Statistical significance was denoted as *P < 0.05, **P < 0.01, ***P < 0.001.
Results
Overall characteristics of the cell cluster composition in pancreatic cancer
To gain a comprehensive understanding of TME, and explore its heterogeneity between PDAC and normal tissues, we investigated 39 PDAC samples consisting of primary tumors, metastatic tumors, adjacent normal tissues and paired peripheral blood (Figure 1A) from GSE155698 (24). Based on canonical cell markers mentioned in Methods (Table S1), a total of 124,575 cells were classified into distinct cell populations, including T cells (39,372 cells, 31.61%), myeloid cells (48,054, 38.57%), epithelial cells (14,998, 12.04%), NK cells (7,758, 6.23%), fibroblast (2,617, 2.10%), B cells (4,548, 3.65%), acinar cells (2,488, 2.00%), mast cells (2,086, 1.67%), stellate cells (1,324, 1.06%), endothelial cells (1,238, 0.99%) and minor unknown cells (92, 0.74%) (Figures 1B, C, S1A).
Compared to primary tumors and paracancerous samples, metastatic tumors exhibited higher composition of T cells (34.52%, 27.55%, 23.59%, respectively) and epithelial cells (15.17%, 12.54%, 9.47%) (Figure S1B, Table S2). Conversely, compared to primary and metastatic tumor samples, adjacent normal tissues had higher proportions of stromal cells, including acinar cells (14.27%, 1.11%, 1.21, respectively), stellate cells (6.26%, 0.95%, 0.57%) and endothelial cells (3.90%, 0.92%, 0.89%) (Figure S1B, Table S2). As for PDAC samples across clinical stage I to IV, the proportions of epithelial cells continued to rise, from 5.58% to 15.17%, while myeloid cell compositions decreased from 71.63% to 35.13% (Figure S1C, Table S2). Among the PDAC samples, T cells (28.68%), myeloid cells (40.25%), and epithelial cells (12.86%) were the most abundant populations. In normal samples, T cells (51.88%), myeloid cells (27.00%), and NK cells (9.26%) accounted for 88.14% of the cell population (Table S2, Figure S1D). The peripheral blood from PDAC patients had a higher proportion of T cells (36.87% vs. 22.41%), myeloid cells (41.46% vs. 33.53%), and NK cells (7.53% vs. 3.95%) compared to solid tissue, whereas the proportion of epithelial cell was higher in solid tissue (20.19%) than peripheral blood (7.37%) (Figure S1E, Table S2).
Similar to previous studies (26, 32), there was significant variation in the portions of epithelial, stromal, and immune cells among the samples, which could be attributed to intrinsic differences in tumor stages or specific locations within tumor where biopsies were taken (Figure 1D). For example, PDAC patients 15 and 16 (stage II) exhibited a highly immune-rich microenvironment, with nearly 70% T cells in peripheral blood compared to only 10% in solid tissue. Furthermore, patients 15 and 16 (stage II) had higher T cell portions (70%) than patients 2 and 3 (stage IV, approximately 20%). These findings indicated that the formation and progression of metastases in PDAC may necessitate a more immunosuppressive TME compared to primary tumors.
B cell may play a tumor-suppressive role in PDAC
Subsequently, to assess the clinical significance of these cell types in PDAC, we identified the top 20 genes that predominantly determined the identity of each cell type through ROC analysis. The correlation between the expression levels of these genes and the patient prognosis was then computed using multivariate Cox regression on TCGA-PAAD data (Table S3). Our analysis revealed that genes exclusively expressed in C3 (epithelial cell) (ave.cox = 0.182), C13 (epithelial cell) (ave.cox = 0.182), C21 (epithelial cell) (ave.cor = 0.193) were associated with poor prognosis (Figure S2A), where ave.cox represents the average Pearson correlation coefficient. On the other hand, genes expressed in C0 (T cells) (ave.cox = -0.076) and C14 (B cells) (ave.cox = -0.062) were correlated with a favorable prognosis in PDAC, suggesting potential tumor-suppressive functions of these cells. It is worth noting that B cells are prominent features of PDAC tumors, although their roles in this disease remain controversial (33). Notably, higher expression levels of genes exclusively expressed in the C14 and C38 B cell types (such as BCL11A and DNASE1L3) were positively associated with favorable prognoses(Figure S2B, Table S3), indicating the tumor-suppressive functions of C14 and C38 cells in the PDAC microenvironment. Correlation analysis revealed that BCL11A (R = 0.419, p = 3.52e-09) and DNASE1L3 (R = 0.689, p < 2.2e-16) were positively correlated with CD8A (Figure S2C). Clinical outcomes demonstrated that higher expressions of BCLAA1 and DNASE1L3 were significantly associated with improved survival (Figures S2E, F). Additionally, these genes exhibited higher expression values in PDAC tumors compared to normal samples (Figure S2D), suggesting B cells may exert tumor-suppressive roles as tumor-infiltrating B cells. Consistent with our findings, previous studies have demonstrated that tumor-infiltrating B cells are a positive prognosis factor, both in PDAC and other cancers (34, 35).
Myeloid cells exert immune-suppressive potential
The presence and functional activities of myeloid cells in tumors have garnered increasing interest due to their relevance as modulators of anticancer therapies and potential targets for specific treatment. In this study, we focused on unraveling the potential roles of myeloid cells in PDAC (Figure 2). Correlation analysis showed that the genes exclusively expressed in C1 (ave.cox = 0.023), C4 (ave.cox = 0.019), C7 (ave.cox = 0.026), C10 (ave.cox = 0,023), C17 (ave.cox = 0.034), C22 (ave.cox = 0.042) (myeloid cell) were associated with poor prognosis (Figure 2A), suggesting the immune-suppressive functions of myeloid cells in the microenvironment of PDAC. Besides, Immune-suppressive markers as previously reported (37), SPP1, MACRO, APOE, CD68, and SIPRA, were exclusively expressed in myeloid cells (Figure 2B). Additionally, myeloid cells had a relatively higher stemness score compared to other stromal cells (Figure 2C).Previous studies have demonstrated that cancer progression involves a gradual loss of differentiated phenotype and the acquisition of progenitor-like, stem cell-like features (38). Furthermore, myeloid cells exhibited heterogeneous expression of immune checkpoint receptors (CD86, HAVCR2, CD48, and VSIR) (Figure 2D). Collectively, these findings suggest that myeloid cells may play an immune-suppressive role in the PDAC tumor environment, consistent with previous observations and supporting the notion that myeloid cells are a key immunosuppressive component in TME (39).
Figure 2 Myeloid cells exert immune-suppressive potentials. (A) Bar Graphs illustrating the average coefficients (ave.cor) for the prognostic effect of genes exclusively expressed in each clusters. Positive values of ave.cor indicate the associations with anti-tumor immunity and good clinical outcomes, opposite from negative ave.cor values. (B) Expression landscapes of immunosuppressive markers of SPP1, MARCO, APOE, CD58, and SIRPA. (C) The high stemness score profile of myeloid cells calculated by scCancer (36). (D) Heatmap of immune checkpoint inhibitors in 11 cell clusters.
Characterization of major myeloid cell lineages
To investigate tumor-infiltrating myeloid cells (TIMs) in PDAC, we firstly excluded all cells from healthy samples, resulting in 45,859 myeloid cells for further analysis. Subsequently, we performed unsupervised clustering and cell annotation of myeloid cells using canonical markers (Figures 3A, C), as described in the Methods section. It revealed the presence of 6 distinct subclusters within the myeloid lineage, including granulocytes (24,786, 54.0%), CD14 monocytes (8,713, 19.0%), macrophages (7,767, 16.9%), mast cells (1,999, 4.4%), CD16 monocytes (1,270, 2.8%), and dendritic cells (DCs) (1,324, 2.9%) (Figure 3B).
Figure 3 Myeloid-derived cell components in PDAC. (A) t-SNE plot showing 6 myeloid clusters of samples from PDAC patients. (B) Proportion of each myeloid cell lineage from the primary tumor, metastatic tumor, paracancerous tissue, and peripheral blood. (C) Bubble plot showing selected cell type-specific markers across all clusters. The size of dots represents the fraction of cells expressing a particular marker, and the intensity of the color indicates the levels of average mean expression. (D) t-SNE plots showing the expression of specific immunosuppressive markers, which were denoted in Figure B, in the myeloid subclusters. (E–I). Kaplan-Meier survival analyses of some myeloid subclusters markers, including macrophage, mast cells, CD14 monocyte, CD16 monocyte, and DCs.
To visualize the distribution of cell populations of myeloid subclusters across different histologic types, we utilized unbiased hierarchical clustering algorithms and supervised annotation on PBMC samples(Figures S3A, B). Compared to adjacent normal tissues, tumor tissues exhibited overall increases in DCs, CD14 monocytes, and CD16 monocytes, indicating a redirected immune response (Figure S3C). Besides, PBMC samples predominantly consisted of granulocytes, while tumor and paratumor tissues exhibited abundant macrophages (Figure S3C). The distict cellular compositions suggested a heterogeneous TIMs environment in tumor.
We subsequently investigated the expression patterns of immune-suppressive markers mentioned above. The results revealed elevated expression of these markers in macrophages (Figure 3D), suggesting that macrophages may contribute to the immunosuppressive effects of TIMs in PDAC. To further assess the clinical impact of the signature (Figure 3C) for myeloid subclusters, including macrophages, mast cells, CD14 monocytes, CD16 monocytes and DCs, we utilized an independent PAAD cohort from TCGA. Patients with high expression of signature genes exhibited worse OS compared to those with low expression(two-sided log-rank test p < 0.01) (Figures 3E–I). These findings further underscored the immunosuppressive effects of myeloid lineages at the bulk level.
PDAC tumor progression is associated with necroptosis
In order to comprehensively investigate the diverse range of myeloid cell populations in PDAC, we first performed tissue prevalence analysis. Compared to paratumor samples, macrophages and mast cells were highly enriched in tumor tissues, suggesting the coexistence of host immune response and tumor escape in the PDAC milieu. Moreover, macrophages exhibited a higher enrichment in paratumor samples compared to primary and metastatic tumor samples, in contrast to CD14 monocytes and granulocytes (Figure S3C, left). Moreover, macrophages and mast cells exhibited a preferential enrichment in PDAC tissues rather than peripheral blood samples (Figure S3C, middle). Subsequently, we performed differential gene expression (DGE) analysis between tumor and paratumor samples (Figure 4A) and gene set enrichment analysis of each cluster’s upregulated genes (Figures 4B-F). Strikingly, gene ontology (GO) characteristics related to necroptosis were detected across all myeloid lineages in tumor samples. Necroptosis can either elicit robust adaptive immune responses that may impede tumor progression, or it can recruit inflammatory responses that may potentially facilitate tumorigenesis, cancer metastasis and the generation of an immunosuppressive tumor microenvironment (40).
Figure 4 Differential gene analyses of PDAC tumors and paracancerous tissues. (A) Differential gene expression analysis showing up- and down-regulated genes across all cell types between tumor and paratumor samples from PDAC patients. The top 10 DE genes were shown, and the points dotted in red indicate significant genes. An adjusted p-value < 0.01 is indicated in red, while an adjusted p-value ≥ 0.01 is indicated in black. (B–F). Differential pathways enriched in tumor and paratumor for each cell type.
To investigate the potential enrichment of necroptosis in peripheral blood, we only kept myeloid cells from peripheral blood samples (Figure S3A). Compared with normal blood samples, blood samples from PDAC patients showed higher proportions of CD14 monocytes and lower proportions of DCs, and CD16 monocytes (Figure S3B), indicating a heterogeneous myeloid environment in peripheral blood. However, myeloid lineages showed comparable enrichment in peripheral blood except for DCs, CD16 monocytes, and macrophages (Figure S3C, right). We performed DGE and enrichment analyses of PBMC samples from PDAC or healthy individuals, and results demonstrated the presence of necroptosis in PBMC samples from PDAC (Figures S3D-I), indicating that the necroptosis event was not tissue-specific. To investigate the differences between tumor and paratumor tissues, we excluded myeloid cells from peripheral blood samples. DGE and enrichment analyses showed that necroptosis remained specifically enriched in all myeloid lineages within tumor tissues (Figures S4A-F). Moreover, we removed myeloid cells from adjacent normal samples to compare the differences between tumor tissue and peripheral blood from PDAC. The results showed that necroptosis was specifically enriched in all myeloid lineages from tumor tissues, rather than peripheral blood (Figures S5A-F), suggesting a propensity for necroptosis events to occur in solid tumor tissues. Detailed information on all DEGs and necroptosis-associated DEGs could be found in Tables S4,S5, respectively.
We next performed overlapping analyses to find key necroptosis-pathway-associated (NPA) DEGs that exhibited significant up- and down-regulation within each myeloid lineage (Figure 5; Figure S6). The results revealed specific patterns in the expression of NPA genes in different myeloid cell types within tumor tissues. In macrophages from tumor samples, two NPA genes, GLUL and SQSTM1, were found to be up-regulated (Figure 5A), while SLC25A6 exhibited down-regulation (Figure 5B). Besides, in mast cells, two NPA genes HSP90AA1 and HSP90AB1 were up-regulated in tumor tissues (Figure 5A), whereas BIRC3 was down-regulated (Figure 5B). Interestingly, BIRC3 displayed an opposite expression pattern, being up-regulated in CD14 monocytes and CD16 monocytes but down-regulated in mast cells (Figure 5B), indicating potential distinct roles of this NPA gene in different cell types. Additionally, in CD16 monocytes, the remaining NPA genes, JAK3, PPIA, and TLR4 were up-regulated in tumor tissues (Figure 5A), while IFNGR1 was down-regulated in tumor samples (Figure 5B). In the case of DCs, the NPA gene, CHMP1B was up-regulated in tumor tissues (Figure 5A), while PARP1 was down-regulated in tumor samples (Figure 5B).
Figure 5 Overlapped differentially expressed genes were associated with necroptosis. Analysis of overlapped necroptosis-associated genes up-regulated (A) or down-regulated (B) in tumor for each cell type. (C) TCGA survival analysis of SLC25A6. (D) The predicted regulatory pathways of overlapped necroptosis-associated genes we identified based on public Necroptosis pathway network. Genes in dashed-circles were not identified in overlapped DEGs. ‘pdac.tumor.tissues&pbmc vs pdac.paratumor’, primary tumor tissues, metastatic tumor tissues and PBMCs samples from PDAC patients vs paratumor samples from PDAC patients.; ‘pdac.tumor.tissues vs pdac.paratumor’, primary tumors and metastatic tumors from PDAC patient vs paracarcinoma tissues from PDAC patient; ‘pdac.pbmc vs healthy.pbmc’, PBMCs of PDAC Patients vs PBMCs from healthy controls; ‘pdac.tumor.tissues vs pdac.pbmc’, primary tumor tissues and metastatic tumor tissues of PDAC patients vs PBMCs of PDAC patients; LMP, lysosome membrane permeabilization; ΔΨm, mitochondrial membrane potential.
Based on the Necroptosis pathway network (https://www.kegg.jp/pathway/map04217), several NPA DEGs that identified within PDAC myeloid cells were involved in this network, thereby the interactive relationship was mapped. This approach allowed us to gain insights into the underlying contribution of necroptosis in the context of myeloid cell-mediated immune responses within TME of PDAC. Of note, we designated this mapping as the representation of the necroptosis pathway associated with PDAC myeloid cells, while experimental studies are necessary to confirm their precise roles and interactions(Figure 5D). In this predicted simplified model, TLR4 acts as an upstream regulator that promotes the phosphorylation of RIPK3, a key regulator of necroptosis (41). This phosphorylation event leads to the subsequent phosphorylation of GLUL, which contributes to increased lysosome membrane permeabilization (LMP), a common phenomenon in cancer cells (42). Besides, JAK3, a downstream signaling molecule of IFN, BIRC3, a downstream molecule of TNF, and SQSTM1, a dissociated molecule, work together to promote the phosphorylation of RIPK1, another core regulator of necroptosis (41). The phosphorylated RIPK1, in turn, represses the expression of SLC25A6, result in a transient increases in mitochondrial transmembrane potential (ΔΨm), which is highly related to cancer malignancy (43). Moreover, HSP90AA1 and HSP90AB1 can simultaneously promote the phosphorylation of necrosome, including RIPK1, RIPK3, and MLKL. This leads to the activation of various necroptosis pathways, such as MLP, ΔΨm and mitochondrial fission (44). Mitochondrial fission facilitates the proliferation, metastasis, and drug resistance of cancer cells (45). Despite their necroptosis-promoting function, HSP90AA1 and HSP90AB1 can also act as upstream regulator of ESCRT-III, which helps maintain membrane integrity during the initiation of necroptosis, thereby promoting cell survival (46).
To explore the clinical relevance of necroptosis-associated genes that were down-regulated in tumor (Figure 5B), we conducted survival analysis. The results revealed that higher expression of SLC25A6 was correlated with improved survival outcomes (Figure 5C). Interestingly, we observed a significantly higher expression of SLC25A6 in HSP90AA1+HSP90AB1+ mast cells compared to HSP90AA1-HSP90AB1- mast cells (Figure 6E). These findings indicated that HSP90AA1 and HSP90AB1 might have an unknown mechanism of targeting SLC25A6 in necroptosis pathway (Figure 5D).
Figure 6 Schema to identify HSP90AA1+HSP90AB1+ mast cells. (A) Schema showing the procedures to distinguish mast cells based on HSP90AA1 and HSP90AB1. (B) Violin plots displaying the expression of HSP90AA1+HSP90AB1+ mast cluster-specific genes. (C) Scatterplot showing the Spearman correlation between the HSP90AA1, HSP90AB1, and CD8A in this study. (D) Scatterplot illustrating the correlations between the top 6 HSP90AA1+HSP90AB1+ mast cluster-specific genes and CD8A in this study. (E) Expression patterns of cytokines (AVP, CTSG, NAMPT, SLC25A6) in mast cells.
A novel immunological RTM population is specific to paratumor tissue
Based on the identified NPA genes mentioned above, we proceeded to investigate the relations within myeloid lineages. Firstly, myeloid cells from tumor and paratumor samples were selected for subsequent analysis (Figure 7A, left). Subsequently, an unsupervised clustering analysis was performed on macrophage subsets (Figure 7A, middle). Cluster 3, characterized by high expression of ITGAX, CD86, HLA-DRA, and HLA-DRB1, was identified as M1 macrophages (Figure 7B). Clusters 0, 2, 4, and 6 were classified as M2 macrophages as they highly expressed SPP1, MACRO, APOE, FABP5, and LAMP1. Clusters 1 and 5, which displayed elevated expression profiles of S100A4, RGS1, CD74, and CSF1R,were designated as RTM (resident tissue macrophage) subset.
Figure 7 Pipeline to identify novel immunological RTM (resident tissue macrophage) population. (A) Schema illustrating the procedures for identifying the subclusters of macrophages. (B) Bubble plot displaying selected cell type-specific markers across all clusters. The size of dots represents the fraction of cells expressing a particular marker, and the color intensity indicates the level of average mean expression. (C) t-SNE plots showing the expression of GLUL and SQSTM1. (D) t-SNE plots showing the subclusters of macrophages between tumor and paratumor samples. (E) Violin plots presenting RTM (GLUL-SQSTM-) cluster-specific markers. (F) Correlations between the GLUL-SQSTM1- RTM cluster-specific genes and CD8A. (G) TCGA validation of the genes illustrated in (F). (H) Expression and distribution patterns of cytokines in macrophages. The left figure represents macrophage subsets, while the right figures are expression patterns of some cytokines.
Focusing on the expression profiles of NPA genes, including GLUL and SQSTM1 in macrophages, strikingly, we found that GLUL and SQSTM1 were elevated in all macrophage subsets except for a minor RTM cell population (Figure 7C). Therefore, we named this RMT subset as GLUL-SQSTM1- RTM. Unexpectedly, this subset was deficient in tumor samples compared to paratumor samples (Figure 7D), indicating that GLUL-SQSTM1- RTM may act as a positive regulator of immunity. To validate this hypothesis, we performed a correlation analysis of GLUL-SQSTM1- RTM and CD8+T cells using GLUL-SQSTM1- RTM markers, including PRSS1, CTRB1, CLPS, PLA2G1B, PNLIP and CPA1 (Figure 7E). The results showed positive correlations between these markers with CD8+T cells both in this cohort (Figure 7F) and TCGA-PAAD cohort (Figure 7G), which validated that GLUL-SQSTM1- RTM represents an immunological cell population. Furthermore, to investigate whether this subpopulation directly interacts with epithelial cells, we performed correlation analyses. The results indicated that this subpopulation had no direct correlations with epithelial cells, neither in this cohort (Figure S7A) nor in the TCGA-PAAD cohort (Figure S7B). This suggests that this subpopulation may not directly exert immunological functions on epithelial cells.
Cytokines, which are small proteins crucial in controlling the growth and activity of the immune system, play a significant role in the immune and inflammatory responses of all cells in the body (47). Alternatively, cancers can respond to host-derived cytokines that promote growth, inhibit apoptosis and facilitate invasion and metastasis (48). In this study, we investigated the expression profiles of several cytokines in macrophages. The results showed that ADM (49), CCL2 (50), CCL4 (51), CXCL3 (52), MIF (53), SPP1 (54), VEGFA (55) and VEGFB (56), which have previously been reported to promote tumor progression and metastasis, were specifically deficient in GLUL-SQSTM1- RTM (Figure 7H). This further emphasizes the immunological role of GLUL-SQSTM1- RTM in TME.
HSP90AA1+HSP90AB1+ mast cells are pro-immune
After investigating macrophages, our focus shifted to mast cells (Figure 6). We extracted mast cells from the myeloid cell population to perform unsupervised clustering. A total of 6 clusters were identified, and cluster 3 lacked HSP90AA1 and HSP90AB1 expression. Consequently, we termed cluster 3 as HSP90AA1-HSP90AB1- mast cells, while the remaining mast cells were classified as HSP90AA1+HSP90AB1+ mast cells (Figure 6A). To determine whether HSP90AA1+HSP90AB1+ mast cells exhibited pro-immune or anti-immune characteristics, we evaluated the correlations between HSP90AA1+HSP90AB1+ mast cell markers (Figure 6B) and CD8+T cells. The results showed that all of the HSP90AA1+HSP90AB1+ mast cell markers, including HSP90AA1, HSP90AB1, TPSAB1, AREG, CPA3, JUN, LTC4S, CLU, KIT, FAU, etc., exhibited positive correlations with CD8A (Figures 6C, D, S8A), demonstrating that HSP90AA1+HSP90AB1+ mast cell may act as positive regulator of immunity. In contrast, HSP90AA1-HSP90AB1- mast cells were anti-immune, and their markers, including CLC, RUNX1, FAM101B, SORL1, PIM1, CSF3R, ATP100, MAF, MYO1F, etc. (Figure S8B), showed negative correlations with CD8+ T cells (Figure S8C). To further investigate the novel mast clusters, we evaluated the expression patterns of certain cytokines (Figure 6E), including AVP, CTSG, and NAMPT. AVP (57) and CTSG (58), known to play important roles in inflammation and immune responses, were exclusively sufficient in HSP90AA1+HSP90AB1+ mast cells, while NAMPT (59), previously reported to be associated with maintaining cancer stemness, was highly expressed in HSP90AA1-HSP90AB1- mast cells. Additionally, the expression of CD8A was significantly higher in HSP90AA1+HSP90AB1+ mast cells (Figure S8D), providing further supporting for the notion that HSP90AA1+HSP90AB1+ mast cells were immune-promoting.
JAK3+TLR4+ CD16 monocytes are anti-immune
Furthermore, we investigated CD16 monocytes (Figure 8). Firstly, we extracted CD16 monocytes from myeloid cells and performed unsupervised clustering. As a result, a total of 7 clusters were identified, in which clusters 1, 4, and 5 were all JAK3-TLR4- (Figure 8A). Therefore, we classified CD16 monocytes into JAK3+TLR4+ CD16 monocytes and JAK3-TLR4- CD16 monocytes (Figure 8A). We then proceeded to examine the distinctive features of these two clusters. In contrast to HSP90AA1+HSP90AB1+ mast cells, JAK3+TLR4+ CD16 monocytes were immunosuppressive, as indicated by the negative correlation between their markers (JAK3, TLR4, CRIP1, IFI6, ZBTB7A, ZYX, etc.) (Figure 8B) and CD8+T cells at a significant levels (Figures 8C, D). In contrast, JAK3-TLR4- CD16 monocytes displayed upregulation of EEF1D, MS4A4A, TMEM66, and TNF (Figure S9A), all of which were positively correlated with CD8+ T cells at significant levels (Figure S9B), suggesting a pro-immune role of JAK3-TLR4- CD16 monocytes. Cytokine analysis further confirmed these results. CAT (60), CECR1 (61), GPI (62), HDGF (63), and MIF (53), previously reported to promote tumor development and progression, were specifically abundant in JAK3+TLR4+ CD16 monocytes, rather than JAK3-TLR4- CD16 monocytes (Figure 8E).
Figure 8 Pipeline to identify JAK3+TLR4+ CD16 monocytes. (A) Workflow showing the procedures to identify JAK3+TLR4+ mast cells. (B) Violin plots displaying the expression of JAK3+TLR4+ mast cluster-specific genes. (C) Scatterplot showing the correlations between JAK3, TLR4, and CD8A in this study. (D) Scatterplot illustrating the correlation between JAK3+TLR4+ mast cluster-specific genes and CD8A in this study. (E) Violin plot showing the expression profiles of some cytokines in CD16 monocytes. **** Represents a statistical significance level of P<0.0001.
Furthermore, we explored the NPA genes in other myeloid cells, however, the correlations between these genes and CD8A were insignificant. For example, in BIRC3+ CD14 monocytes, the marker BIRC3, did not show a significantly correlation with CD8+ T cells (Figure S10A). Moreover, in BIRC3+ CD16 monocytes, although BIRC3 exhibited a significant positive correlation with CD8A, other markers such as APRT, C1QB, GABARAP, and IFITM1 were not correlated with CD8+ T cells (Figure S10B). Moreover, in CHMP1B+ DCs, the marker CHMP1B was not correlated with CD8+ T cells (Figure S10C). These results indicated that these NPA genes may function as immune mediators in an unknown manner.
Cluster-specific cellular interaction networks that mediate immunity
To gain a comprehensive understanding of the interactions among the novel clusters and the T/B/NK/epithelial cell populations, as well as their collective contribution to the PDAC tumor microenvironment, we inferred a putative cellular interaction network based on the receptor-ligand database (31). Our findings unveiled specific interactions within various signaling pathways. Specifically, we observed that EPO/EPOR interaction of EPO signaling pathway was unique to the GLUL-SQSTM1- RTM cluster. In this cluster, the ligand EPO was predominantly expressed in GLUL-SQSTM1- RTM, while EPOR receptor was present in CD4+ T cells, CD8+ T cells, B cells, NK cells and epithelial cells (Figure 9A). The EPO/EPOR interaction has been reported to initiate a signaling cascade that activated and recruited a variety of Src homology-2 (SH2) domain-containing proteins, subsequently triggering downstream signaling pathways such as ERK-1/2 and JAK-2 (64). Interestingly, a recent study has demonstrated that EPO/EPOR could reduce the variability of myeloma cell lines and malignant primary plasma cells (65). And ANXA1/FPR1 interaction within the ANNEXIN signaling pathway was specific to GLUL+SQSTM1+ macrophage cells, and ligand ANXA1 was specific to GLUL+SQSTM1+ macrophage, while receptor FPR1 was in CD4+ T, CD8+ T, B, NK and epithelial cells (Figure 9B). FPR1, previously reported to promote chemotherapy-induced antitumor immune response (66), was demonstrated to act as a receptor for ANXA1, promoting cell death through the necroptosis pathway (67). Furthermore, we identified the IL10/IL10RA interaction within the IL10 signaling pathway, which was specific to HSP90AA1+HSP90AB1+ mast cell, in which ligand IL10 was expressed in CD4+ T, CD8+ T, B, NK, and epithelial cells, while the receptor IL10RA was specific to HSP90AA1+HSP90AB1+ mast cell (Figure 9C). IL10- and IL10R-dependent signaling have been reported to play critical roles in controlling immune responses in both innate and adaptive immune systems (68). Finally, TNF-α (TNF)/TNFR1 (TNFRSF1A) interaction was found to be specific to JAK3+TLR4+ CD16 monocytes. In this interaction, the ligand TNF was specific to JAK3+TLR4+ CD16 monocytes, while receptor TNFR1 was expressed in CD4+ T, CD8+ T, B, NK and epithelial cells (Figure 9D). TNF-α signaling meditated by TNFR1 in the TME has been reported to promote gastric tumor development and maintain tumor cells in an undifferentiated state (69). Taken together, the intercellular interactions revealed a close relationship between immune cell and cancer cell dynamics, as well as the molecular features of novel clusters (Figure 10). These interactions may play a crucial role in determining the prognostic and therapeutic response in PDAC.
Figure 9 Cluster-specific interaction network. Signaling pathway network of ligand-receptor pairs that were specific in (A) GLUL-SQSTM1- RTM, (B) GLUL+SQSTM1+ macro, (C) HSP90AA1+HSP90AB1+ mast, (D) JAK3+TLR4+ CD16 Mono.
Figure 10 The crosstalk between novel myeloid-derived novel clusters, CD8+T cells Macrophage (A), GLUL-SQSTM1-Macrophage (B), HSP90AA1+ HSP90AB1+ Mast (C) and JAK3+ TLR4+ CD16+ Monocyte (D), and CD8+T cells that mediate immunity. Schematic for cellular crosstalk and corresponding signaling pathways in PDAC TIMs that contribute to immunity or immune suppression. The novel clusters secrete various cytokines and ligands that signal to their respective receptors, thus activating the corresponding signaling, respectively. Red arrows indicate up-regulated cytokines while green arrow indicates down-regulated cytokines.
Validation of the existence of the novel clusters
We proceeded to validate the presence of GLUL-SQSTM1- RTM, HSP90AA1+HSP90AB1+ mast cell, and JAK3+TLR4+ CD16 monocyte cellular clusters in other PDAC cohorts. To accomplish this, we analyzed publicly available scRNA-seq data from the CRA001160 dataset (26). In this dataset, 1,047, 3,098, and 1,464 cells were annotated as RTM, mast cells, and CD16 monocytes, respectively (Figure S1). To focus specifically on RTM, mast cells, and CD16 monocytes, we distinguished these cell types based on the expression of GLUL/SQSTM1 (RTM), HSP90AA1/HSP90AB1 (mast cell), and JAK3/TLR4 (CD16 monocytes), respectively (Figure 11A). Consequently, we obtained 49 GLUL-SQSMT1- RTM, 2,151 HSP90AA1+HSP90AB1+ mast cells, and 1,234 JAK3+TLR4+ CD16 monocytes (Figure 11A).
Figure 11 Detection and validation of the cellular clusters. (A) t-SNE plot combing GLUL+SQSTM1+ RTM (n = 998, in purple), HSP90AA1+HSP90AB1+ mast cells (n = 2,151, in blue), and JAK3+TLR4+ CD16 monocyte (n = 1,234, in orange) from PDAC (Peng et al., 2019). (B) RNA FISH staining in the PDAC tissues. One representative image for each gene is shown. Spectrum orange dots indicate GLUL, spectrum red dots indicate SQSTM1, spectrum green dots indicate S100A4, spectrum gold indicate CD8A. Scale bars, 90 μm.
Although GLUL-SQSMT1- RTM was specifically present in PDAC tumor tissues, it was noteworthy that the majority of RTM in the TME are characterized by the expression of GLUL and SQSTM1. Given the prevalence of GLUL+SQSTM1+RTM cells in PDAC tumor tissues, and their interaction with CD8+T cells revealed by the cell chat analysis, it becomes apparent that a deeper exploration is warranted. To validate the expression of these genes at single-cell level, visualize their spatial distribution within complex tissue structures and validate co-expression patterns, functional gene interactions, we performed RNA fluorescence in situ hybridization (FISH) on paraffin sections of PDAC tumor tissue (Figure 11B).Our analysis revealed frequent overlap of the GLUL, SQSTM1 and S100A4, providing evidence for the existence of the GLUL+SQSTM1+ RTM. Additionally, we conducted further investigations and confirmed the adjacent spatial relationship between the GLUL+SQSTM1+ RTM and CD8+ T cells in PDAC solid tumors. This observation suggests potential functional interactions between these cell populations within the tumor microenvironment.
Discussion
It is well-established that necroptosis, a programmed form of necrosis or inflammatory cell death (70), has gained significant attention in cancer research due to its implications in pathogenesis and therapy (40, 71, 72). The involvement of necroptosis in recruiting immune cells, regulating pro- or anti-tumor components in TME, and modulating immune responses through the release of DAMPs, chemokines and other cytokines. However, the specific roles of necroptosis in PDAC have not been systematically described and remain to be deciphered. In this study, public scRNA-seq data from Gene Expression Omnibus (GEO) database was downloaded, which covered more than 124,000 cells of 16 PDAC patients across multiple immune-relevant tissue sites (24). Our analysis focused on charactering major myeloid cell lineages, and identifying three necroptosis-associated subsets of myeloid-derived cells. These findings provide a valuable resource for comprehensively understanding multi-dimensional characterization of the tumor microenvironment in PDAC.
In this study, we identified HSP90AA1+HSP90AB1+ mast cells that exerting anti-tumorigenic effects in PDAC, contributing to better clinical outcomes. the Heat shock protein 90 (HSP90) protein family, including HSP90AA1 and HSP90AB1, plays prominent roles in various biological processes such as protein folding (73), apoptosis (74), cell-cycle regulation (75), as well as signal transduction (76). Of note, previous studies have suggested that HSP90 can influence the activation and stability of crucial regulators involved in the necroptosis process, such asRIPK1, RIPK3 and MLKL, thereby contributing to immune cell recruitment and immunogenic cell death of tumor cells (77–79). In addition, necroptosis is deemed to trigger an adaptive immune response by releasing cytokines (80). The novel mast cells were predicted to interact with CD8+ T cells via IL10 signaling pathway, and highly expressed AVP and CTSG. Consistently, AVP (57) and CTSG (58) have been reported to play important roles in inflammation and immune response, and IL10 has been shown to potentiate IFN-γ and induct the cytotoxicity of CD8+ T cells (68, 81–83), thereby triggering anti-tumor immune responses. These findings collectively suggested that HSP90AA1+HSP90AB1+ mast cells are functionally important in necroptosis process and are involved in immune cell recruitment through the IL10 signaling pathway. Future research should focus on designing effective drugs that modulate HSP90 activity and developing diagnostic tools for accurate patient stratification for therapy with HSP90 agonists or HSP90 antagonist.
Another cluster of special interest is the GLUL-SQSTM1- RTM subpopulation, which is enriched in PDAC tumor tissues and acts as a positive regulator of immunity. Glutamine synthetase (GLUL) has been associated with RIP3-mediated necroptosis (41), and the p62/SQSTM1 complex binding to necroptosis-related proteins RIP1 and RIP3 facilitates the transition from autophagy to necroptosis (84). Moreover, it has been reported that both GLUL and p62/SQSTM1 may influence the recruitment, activation and polarization of macrophage. GLUL is known to be associated with the differentiation and function of macrophage, more specifically, enhancing M2- polarization (85, 86). Exogenous p62/SQSTM1 has been shown to induce M1 polarization of macrophage through activation of the NF-κB pathway (87). Given the potential roles of GLUL and p62/SQSTM1 in macrophage polarization, we evaluating the expression levels of molecular markers in macrophages. For example, CCL2/CCR2 axis is a major player in macrophage polarizing towards to M2 phenotype (88, 89). Stimulation of CXCL13 may activate Akt pathway, suggesting an increase in M2 macrophage in renal cell carcinoma (90). The vascular endothelial growth factor (VEGF) family, including VEGFA and VEGFB, can contribute to M2 polarization in the decidua (91). MIF (92), SPP1 (93) and ADM (94) have all been previously reported to be associated with a dominant M2 polarization and a loss of M1 function. In our study, the deficiency of these chemokines in GLUL-SQSMT1- RTM cells may suggest the pro-inflammatory M1-polarized phenotype, thereby modulating the antitumor response. These results indicated that the GLUL-SQSMT1- RTM cells with a pro-immune profile may evade programmed necroptosis and abundantly infiltrate in the PDAC TME of patients with superior efficacy. The combination of GLUL and SQSTM1 inhibitors in precisely characterized patients may have superior effects against cancer compared to immunotherapy alone.
The novel JAK3+TLR4+ CD16 monocyte subset exhibits anti-immune properties and is associated with unfavorable clinical outcomes. Janus Kinase 3 (JAK3) is a tyrosine kinase that belongs to the Janus family of kinases. Hyper-activation of the JAK3-STAT signaling pathway has been linked to tumor development and progression by inducing factors associated with suppressive immune cell recruitment, angiogenesis and neo-vascularization (95–99). TLR4, a member of the toll-like receptor (TLR) family, can lead to the activation of NF-κB pathway, which is essential for necroptosis signaling, as well as the production of pro-inflammatory cytokines and angiogenetic factors (100, 101). In addition, we observed that this novel monocyte subset was predicted to interact with CD8+ T cells via TNF-TNFR1 signaling pathway and expressed high levels of pro-tumor cytokines CAT, HDGF, CECER1, GPI and MIF. NFR1-dependent TNF signaling has been reported to impair the accumulation of tumor-infiltrating lymphocyte (TILs) and induce significant death of activated CD8+ T cells (102, 103). HDGF is considered as an angiogenic and anti-apoptotic factor, contributing to tumorigenesis in several malignant diseases (104–106). MIF plays an essential role in inhibiting cytotoxic T lymphocytes (CTLs) and regulating lymphocyte transmigration (107, 108). CAT, CECR1 and GPI have all been previously associated with promoting tumor progression (60–62). Altogether, these results align with our observations. We hypothesized that this monocyte subset (JAK3+TLR4+ CD16 monocyte) contributes to shaping pro-tumor immunity in TME, ultimately accelerating malignant transformation and tumor progression. However, further investigation is needed to elucidate the underlying mechanisms. Compared to other JAKs, JAK3 has a more restricted expression profile, primarily confined to immune system. Therefore, selective targeting of JAK3 represents a potent immunosuppressant strategy that could minimize potential adverse effects. Inhibition of TLR4-related pathways has shown promising results in clinical trials for disease treatment with excessive immune response (109–111). Current study supports the notion that the discovery of JAK3 and TLR4 antagonists could be an ideal strategy for cancer treatment.
A major limitation of the current finding is the lack of sufficient experimental validation. For instance, the proposed NPA gene network was supposed only based on the expression profiles, and the underlying regulatory mechanism of necroptosis pathway in tumor progression remains obscure in real world. In addition, our study illustrated the indispensable roles of three novel myeloid subpopulations in tumor microenvironment and their associations with necroptosis, however, the underlying mechanisms need further investigation. And the expression profiles of GLUL and SQSTM1 of macrophages in normal samples were not explored. The implementation of advanced biological techniques and bioinformatics analysis in mammalian models of human pathological samples will be critical for gaining a better understanding of these subpopulations in the context of molecular mechanism and drug targeting.
Given the robust immunosuppressive and desmoplastic TME in PDAC, which contributes to adaptive or acquired resistance to therapy, investigating the relationship between necroptosis and tumor immunology holds promise for future treatment solutions. The identification of necroptosis-associated myeloid lineages can potentially serve as targets for therapeutic intervention, allowing for dynamically monitoring of the anti-tumor immune response and improvement of patient outcomes. Our findings provide a valuable resource for further investigation to gain deeper biological insights into the role of necroptosis in cancer. Considering the exceedingly complex and individually unique immune microenvironment of tumors, necroptosis signaling may generate a diverse array of inflammatory responses, ranging from facilitation of anti-tumor to pro-tumor signaling. The three novel necroptosis-associated myeloid subpopulations uncovered in our research may communicate with other cells to mediate ECM degradation and remodeling, signaling pathway regulation and immune cell polarization. These cells and their respective products hold potential as therapeutic targets in PDAC and other types of cancers, enabling the establishment of effective necroptosis-based cancer therapy regimens.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The studies involving humans were approved by The Forth Medical Center of PLA General Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author contributions
HXZ: Conceptualization, Methodology, Writing – original draft. WD: Conceptualization, Methodology, Writing – original draft. SX: Writing – original draft, Data curation, Software, Visualization. LZ: Data curation, Software, Visualization, Writing – original draft. TF: Data curation, Software, Visualization, Writing – original draft. LW: Writing – original draft, Investigation. HZ: Writing – review & editing. YH: Writing – review & editing. JY: Writing – review & editing. TW: Writing – review & editing, Conceptualization, Methodology, Resources, Supervision. WX: Conceptualization, Methodology, Resources, Supervision, Writing – review & editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Acknowledgments
We appreciate the assistance of the staff of the Senior Department of Oncology, The Fifth Medical Center of PLA General Hospital, who collected and providing the samples used in RNA-Fish.
Conflict of interest
Authors TW, SX, LZ, TF and LW are employed by Hangzhou Repugene Technology Co Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1263633/full#supplementary-material
Supplementary Figure 1 | Heterogeneity of TME in PDAC. (A) t-SNE plots of cells from 39 samples profiled in this study. Colored by cell types. (B) The proportion of 11 cell types among primary tumors, metastatic tumors and paratumor samples, respectively. (C) The histogram of 11 cell types among clinical-stage I-IV. (D) The percentage of cell types between PBMCs from PDAC Patients and PBMCs from healthy controls. (E) The proportion of 11 cell types compared between primary tissues and metastatic tumor tissues of PDAC patients with PBMC samples from PDAC patients.
Supplementary Figure 2 | Functional analysis of B cells in PDAC. (A) Top 100 genes were used to calculate the cox coefficient for each cluster using TCGA (PAAD, n = 183) data. (B) t-SNE plots showing the expression of BCL11A and DNASE1L3. (C) Scatterplots showing the correlation between BCL11A, DNASE1L3, and CD8A using TCGA (PAAD, n = 183) data. (D) The body maps showing the expression of BCL11A and DNASE1L3 between tumor and normal samples via GEPIA 2 (http://gepia2.cancer-pku.cn/#index). (E) Kaplan-Meier overall survival analysis of the high and low groups of BCL11A (top) and DNASE1L3 (bottom). (F) DFS (Disease-Free Survival) analysis of the high and low groups of BCL11A (top) and DNASE1L3 (bottom). The hazard ratio was calculated based on Cox PH Model, and 95% CI (Confidence Interval) was applied.
Supplementary Figure 3 | Differential gene and pathway analysis for peripheral blood. (A) t-SNE plots of myeloid lineages of PBMC from healthy individuals and PDAC patients. (B) Proportion of each myeloid cell lineage in PBMC samples from PDAC patients and healthy controls. (C) Tissue prevalence estimated by Ro/e score of primary tumor/metastatic tumor/paratumor tissues from PDAC patients (left), tissue/PBMC (middle) from PDAC patient, and PBMCs from PDAC patient/healthy controls (right). (D) Top 10 differential up- or down-regulated genes across myeloid lineages of PBMCs from PDAC patient and normal samples. Red dots indicate statistically significant genes (adjusted p-value < 0.01). (E–I) Pathways enriched by DEGs of each cell type between PBMCs from PDAC patients and PBMCs from healthy controls.
Supplementary Figure 4 | Differential genes and pathways between tumor and paratumor tissue of PDAC. (A) Differential gene expression analysis shows up- and down-regulated genes across all cell types between tumor and paratumor tissue from PDAC patients. The top 10 DE genes were shown, and the points dotted in red indicate significant genes. An adjusted p-value < 0.01 is indicated in red, while an adjusted p-value ≥ 0.01 is indicated in black. (B–F) Differential pathway enriched in tumor and paratumor from PDAC tissue for each cell type.
Supplementary Figure 5 | Differential genes and pathways between PDAC tissue and peripheral blood. (A) Differential gene expression analysis showing up- and down-regulated genes across all cell types between tumor tissue samples and PBMC from PDAC patients. The top 10 DE genes were shown, and the points dotted in red indicate significant genes. An adjusted p-value < 0.01 is indicated in red, while an adjusted p-value ≥ 0.01 is indicated in black. (B–E) Differential pathway enriched in tissue and PBMC from PDAC patients for each cell type.
Supplementary Figure 6 | Visualization of overlapped necroptosis-associated DE genes. Violin plots showing overlapped up-regulated (A) and down-regulated (B) necroptosis-associated genes (denoted in ) in tumor. *P < 0.05, **P < 0.01, ***P < 0.001.
Supplementary Figure 7 | The novel immunological RTM has no direct correlation with epithelial cells. (A) Scatterplots showing the correlation between the expression of GLUL-SQSTM1-RTM cluster-specific genes (denoted in ) and EPCAM. (B) TCGA validation of the correlations. Data are from TCGA-PAAD (n = 178).
Supplementary Figure 8 | HSP90AA1-HSP90AB1- mast cells are anti-immune. (A) Scatterplot showing the correlations between the rest of HSP90AA1+HSP90AB1+ mast cluster-specific genes (shown in ) and CD8A in this study. (B) Violin plots showing the expression of HSP90AA1-HSP90AB- mast cluster-specific genes. (C) Scatterplots showing the correlations between HSP90AA1-HSP90AB1- mast cluster-specific genes and CD8A in this study. (D) Expression profiles of CD8A in HSP90AA1+HSP90AB1+ mast cells vs. HSP90AA1-HSP90AB1- mast cells.
Supplementary Figure 9 | JAK3-TLR4- CD16 monocytes are pro-immune. (A) Violin plots showing the expression of JAK3-TLR4- mast cluster-specific high expression genes. (B) Scatterplots showing the correlations between the expression of JAK3, TLR4, and CD8A in this study.
Supplementary Figure 10 | Several myeloid subsets may act as pro/anti-immune regulators in a non-necroptosis way. Workflow showing the procedures to distinguish CD14 monocyte (A), CD16 monocyte (B), and DCs (C) by BIRC3, PPIA, and CHMP1B, respectively.
References
1. Grossberg AJ, Chu LC, Deig CR, Fishman EK, Hwang WL, Maitra A, et al. Multidisciplinary standards of care and recent progress in pancreatic ductal adenocarcinoma. CA Cancer J Clin (2020) 70(5):375–403. doi: 10.3322/caac.21626
2. Schnelldorfer T, Adams DB, Warshaw AL, Lillemoe KD, Sarr MG. Forgotten pioneers of pancreatic surgery: beyond the favorite few. Ann Surg (2008) 247(1):191–202. doi: 10.1097/SLA.0b013e3181559a97
3. Foley K, Kim V, Jaffee E, Zheng L. Current progress in immunotherapy for pancreatic cancer. Cancer Lett (2016) 381(1):244–51. doi: 10.1016/j.canlet.2015.12.020
4. Royal RE, Levy C, Turner K, Mathur A, Hughes M, Kammula US, et al. Phase 2 trial of single agent Ipilimumab (anti-CTLA-4) for locally advanced or metastatic pancreatic adenocarcinoma. J Immunother (2010) 33(8):828–33. doi: 10.1097/CJI.0b013e3181eec14c
5. Brahmer JR, Tykodi SS, Chow LQ, Hwu WJ, Topalian SL, Hwu P, et al. Safety and activity of anti-PD-L1 antibody in patients with advanced cancer. N Engl J Med (2012) 366(26):2455–65. doi: 10.1056/NEJMoa1200694
6. Gamboa AC, Lee RM, Maithel SK. The role of radiation for pancreatic adenocarcinoma. J Pancreatol (2020) 3(2):72–80. doi: 10.1097/JP9.0000000000000045
7. Degterev A, Huang Z, Boyce M, Li Y, Jagtap P, Mizushima N, et al. Chemical inhibitor of nonapoptotic cell death with therapeutic potential for ischemic brain injury. Nat Chem Biol (2005) 1(2):112–9. doi: 10.1038/nchembio711
8. Holler N, Zaru R, Micheau O, Thome M, Attinger A, Valitutti S, et al. Fas triggers an alternative, caspase-8-independent cell death pathway using the kinase RIP as effector molecule. Nat Immunol (2000) 1(6):489–95. doi: 10.1038/82732
9. Green DR. The cell's dilemma, or the story of cell death: an entertainment in three acts. FEBS J (2016) 283(14):2568–76. doi: 10.1111/febs.13658
10. Wang Z, Guo LM, Zhou HK, Qu HK, Wang SC, Liu FX, et al. Using drugs to target necroptosis: dual roles in disease therapy. Histol Histopathol (2018) 33(8):773–89.
11. Kaczmarek A, Vandenabeele P, Krysko DV. Necroptosis: the release of damage-associated molecular patterns and its physiological relevance. Immunity (2013) 38(2):209–23. doi: 10.1016/j.immuni.2013.02.003
12. Pasparakis M, Vandenabeele P. Necroptosis and its role in inflammation. Nature (2015) 517(7534):311–20. doi: 10.1038/nature14191
13. Yang H, Ma Y, Chen G, Zhou H, Yamazaki T, Klein C, et al. Contribution of RIP3 and MLKL to immunogenic cell death signaling in cancer chemotherapy. Oncoimmunology (2016) 5(6):e1149673. doi: 10.1080/2162402X.2016.1149673
14. Moriwaki K, Balaji S, McQuade T, Malhotra N, Kang J, Chan FK. The necroptosis adaptor RIPK3 promotes injury-induced cytokine expression and tissue repair. Immunity (2014) 41(4):567–78. doi: 10.1016/j.immuni.2014.09.016
15. Snyder AG, Hubbard NW, Messmer MN, Kofman SB, Hagan CE, Orozco SL, et al. Intratumoral activation of the necroptotic pathway components RIPK1 and RIPK3 potentiates antitumor immunity. Sci Immunol (2019) 4(36):eaaw2004. doi: 10.1126/sciimmunol.aaw2004
16. Yatim N, Jusforgues-Saklani H, Orozco S, Schulz O, Barreira da Silva R, Reis e Sousa C, et al. RIPK1 and NF-kappaB signaling in dying cells determines cross-priming of CD8(+) T cells. Science (2015) 350(6258):328–34. doi: 10.1126/science.aad0395
17. Wang W, Marinis JM, Beal AM, Savadkar S, Wu Y, Khan M, et al. RIP1 kinase drives macrophage-mediated adaptive immune tolerance in pancreatic cancer. Cancer Cell (2018) 34(5):757–74.e7. doi: 10.1016/j.ccell.2018.10.006
18. Seifert L, Werba G, Tiwari S, Ly NNG, Alothman S, Alqunaibit D, et al. Author Correction: The necrosome promotes pancreatic oncogenesis via CXCL1 and Mincle-induced immune suppression. Nature (2021) 591(7851):E28. doi: 10.1038/s41586-021-03322-8
19. Ando Y, Ohuchida K, Otsubo Y, Kibe S, Takesue S, Abe T, et al. Necroptosis in pancreatic cancer promotes cancer cell migration and invasion by release of CXCL5. PLoS One (2020) 15(1):e0228015. doi: 10.1371/journal.pone.0228015
20. Jayakumar A, Bothwell ALM. RIPK3-induced inflammation by I-MDSCs promotes intestinal tumors. Cancer Res (2019) 79(7):1587–99. doi: 10.1158/0008-5472.CAN-18-2153
21. Kadiyala P, Elhossiny AM, Carpenter ES. Using single cell transcriptomics to elucidate the myeloid compartment in pancreatic cancer. Front Oncol (2022) 12:881871. doi: 10.3389/fonc.2022.881871
22. Zhang J, Song J, Tang S, Zhao Y, Wang L, Luo Y, et al. Multi-omics analysis reveals the chemoresistance mechanism of proliferating tissue-resident macrophages in PDAC via metabolic adaptation. Cell Rep (2023) 42(6):112620. doi: 10.1016/j.celrep.2023.112620
23. Hegde S, Krisnawan VE, Herzog BH, Zuo C, Breden MA, Knolhoff BL, et al. Dendritic cell paucity leads to dysfunctional immune surveillance in pancreatic cancer. Cancer Cell (2020) 37(3):289–307.e9. doi: 10.1016/j.ccell.2020.02.008
24. Steele NG, Carpenter ES, Kemp SB, Sirihorachai VR, The S, Delrosario L, et al. Multimodal mapping of the tumor and peripheral blood immune landscape in human pancreatic cancer. Nat Cancer (2020) 1(11):1097–112. doi: 10.1038/s43018-020-00121-4
25. Butler A, Hoffman P, Smibert P, Papalexi E, Satija R. Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat Biotechnol (2018) 36(5):411–20. doi: 10.1038/nbt.4096
26. Peng J, Sun BF, Chen CY, Zhou JY, Chen YS, Chen H, et al. Single-cell RNA-seq highlights intra-tumoral heterogeneity and Malignant progression in pancreatic ductal adenocarcinoma. Cell Res (2019) 29(9):725–38. doi: 10.1038/s41422-019-0195-y
27. Korsunsky I, Millard N, Fan J, Slowikowski K, Zhang F, Wei K, et al. Fast, sensitive and accurate integration of single-cell data with Harmony. Nat Methods (2019) 16(12):1289–96. doi: 10.1038/s41592-019-0619-0
28. Zhang L, Yu X, Zheng L, Zhang Y, Li Y, Fang Q, et al. Lineage tracking reveals dynamic relationships of T cells in colorectal cancer. Nature (2018) 564(7735):268–72. doi: 10.1038/s41586-018-0694-x
29. Guo X, Zhang Y, Zheng L, Zheng C, Song J, Zhang Q, et al. Global characterization of T cells in non-small-cell lung cancer by single-cell sequencing. Nat Med (2018) 24(7):978–85. doi: 10.1038/s41591-018-0045-3
30. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS (2012) 16(5):284–7. doi: 10.1089/omi.2011.0118
31. Jin S, Guerrero-Juarez CF, Zhang L, Chang I, Ramos R, Kuan CH, et al. Inference and analysis of cell-cell communication using CellChat. Nat Commun (2021) 12(1):1088. doi: 10.1038/s41467-021-21246-9
32. Chen K, Wang Q, Li M, Guo H, Liu W, Wang F, et al. Single-cell RNA-seq reveals dynamic change in tumor microenvironment during pancreatic ductal adenocarcinoma Malignant progression. EBioMedicine (2021) 66:103315. doi: 10.1016/j.ebiom.2021.103315
33. Spear S, Candido JB, McDermott JR, Ghirelli C, Maniati E, Beers SA, et al. Discrepancies in the tumor microenvironment of spontaneous and orthotopic murine models of pancreatic cancer uncover a new immunostimulatory phenotype for B cells. Front Immunol (2019) 10:542. doi: 10.3389/fimmu.2019.00542
34. Brunner M, Maier K, Rümmele P, Jacobsen A, Merkel S, Benard A, et al. Upregulation of CD20 positive B-cells and B-cell aggregates in the tumor infiltration zone is associated with better survival of patients with pancreatic ductal adenocarcinoma. Int J Mol Sci (2020) 21(5):1779. doi: 10.3390/ijms21051779
35. Wouters MCA, Nelson BH. Prognostic significance of tumor-infiltrating B cells and plasma cells in human cancer. Clin Cancer Res (2018) 24(24):6125–35. doi: 10.1158/1078-0432.CCR-18-1481
36. Shan Y, Liu C, Gu J. scCancer: a package for automated processing of single-cell RNA- Q25 seq data in cancer. Brief Bioinform (2021) 22(3):bbaa127. doi: 10.1093/bib/bbaa127
37. Elyada E, Bolisetty M, Laise P, Flynn WF, Courtois ET, Burkhart RA, et al. Cross-species single-cell analysis of pancreatic ductal adenocarcinoma reveals antigen-presenting cancer-associated fibroblasts. Cancer Discovery (2019) 9(8):1102–23. doi: 10.1158/2159-8290.CD-19-0094
38. Malta TM, Sokolov A, Gentles AJ, Burzykowski T, Poisson L, Weinstein JN, et al. Machine learning identifies stemness features associated with oncogenic dedifferentiation. Cell (2018) 173(2):338–54.e15. doi: 10.1016/j.cell.2018.03.034
39. Kemp SB, Pasca di Magliano M, Crawford HC. Myeloid cell mediated immune suppression in pancreatic cancer. Cell Mol Gastroenterol Hepatol (2021) 12(5):1531–42. doi: 10.1016/j.jcmgh.2021.07.006
40. Gong Y, Fan Z, Luo G, Yang C, Huang Q, Fan K, et al. The role of necroptosis in cancer biology and therapy. Mol Cancer (2019) 18(1):100. doi: 10.1186/s12943-019-1029-8
41. Vandenabeele P, Galluzzi L, Vanden Berghe T, Kroemer G. Molecular mechanisms of necroptosis: an ordered cellular explosion. Nat Rev Mol Cell Biol (2010) 11(10):700–14. doi: 10.1038/nrm2970
42. Serrano-Puebla A, Boya P. Lysosomal membrane permeabilization as a cell death mechanism in cancer cells. Biochem Soc Trans (2018) 46(2):207–15. doi: 10.1042/BST20170130
43. Zhang BB, Wang DG, Guo FF, Xuan C. Mitochondrial membrane potential and reactive oxygen species in cancer stem cells. Fam Cancer (2015) 14(1):19–23. doi: 10.1007/s10689-014-9757-9
44. Rayamajhi M, Miao EA. The RIP1-RIP3 complex initiates mitochondrial fission to fuel NLRP3. Nat Immunol (2014) 15(12):1100–2. doi: 10.1038/ni.3030
45. Ma Y, Wang L, Jia R. The role of mitochondrial dynamics in human cancer. Am J Cancer Res (2020) 10(5):1278–93.
46. Gong YN, Guy C, Olauson H, Becker JU, Yang M, Fitzgerald P, et al. ESCRT-III acts downstream of MLKL to regulate necroptotic cell death and its consequences. Cell (2017) 169(2):286–300.e16. doi: 10.1016/j.cell.2017.03.020
47. Arai KI, Lee F, Miyajima A, Miyatake S, Arai N, Yokota T. Cytokines: coordinators of immune and inflammatory responses. Annu Rev Biochem (1990) 59:783–836. doi: 10.1146/annurev.bi.59.070190.004031
48. Dranoff G. Cytokines in cancer pathogenesis and cancer therapy. Nat Rev Cancer (2004) 4(1):11–22. doi: 10.1038/nrc1252
49. Nakamura M, Han B, Nunobiki O, Kakudo K. Adrenomedullin: a tumor progression factor via angiogenic control. Curr Cancer Drug Targets (2006) 6(7):635–43. doi: 10.2174/156800906778742442
50. Jin J, Lin J, Xu A, Lou J, Qian C, Li X, et al. CCL2: an important mediator between tumor cells and host cells in tumor microenvironment. Front Oncol (2021) 11:722916. doi: 10.3389/fonc.2021.722916
51. Mukaida N, Sasaki SI, Baba T. CCL4 signaling in the tumor microenvironment. Adv Exp Med Biol (2020) 1231:23–32. doi: 10.1007/978-3-030-36667-4_3
52. Xin H, Cao Y, Shao ML, Zhang W, Zhang CB, Wang JT, et al. Chemokine CXCL3 mediates prostate cancer cells proliferation, migration and gene expression changes in an autocrine/paracrine fashion. Int Urol Nephrol (2018) 50(5):861–8. doi: 10.1007/s11255-018-1818-9
53. Guda MR, Rashid MA, Asuthkar S, Jalasutram A, Caniglia JL, Tsung AJ, et al. Pleiotropic role of macrophage migration inhibitory factor in cancer. Am J Cancer Res (2019) 9(12):2760–73.
54. Zeng B, Zhou M, Wu H, Xiong Z, et al. SPP1 promotes ovarian cancer progression via Integrin beta1/FAK/AKT signaling pathway. Onco Targets Ther (2018) 11:1333–43. doi: 10.2147/OTT.S154215
55. Kim M, Jang K, Miller P, Picon-Ruiz M, Yeasky TM, El-Ashry D, et al. VEGFA links self-renewal and metastasis by inducing Sox2 to repress miR-452, driving Slug. Oncogene (2017) 36(36):5199–211. doi: 10.1038/onc.2017.4
56. Yang X, Zhang Y, Hosaka K, Andersson P, Wang J, Tholander F, et al. VEGF-B promotes cancer metastasis through a VEGF-A-independent mechanism and serves as a marker of poor prognosis for cancer patients. Proc Natl Acad Sci U.S.A. (2015) 112(22):E2900–9. doi: 10.1073/pnas.1503500112
57. Chikanza IC, Grossman AS. Hypothalamic-pituitary-mediated immunomodulation: arginine vasopressin is a neuroendocrine immune mediator. Br J Rheumatol (1998) 37(2):131–6. doi: 10.1093/rheumatology/37.2.131
58. Gao S, Zhu H, Zuo X, Luo H. Cathepsin G and its role in inflammation and autoimmune diseases. Arch Rheumatol (2018) 33(4):498–504. doi: 10.5606/ArchRheumatol.2018.6595
59. Lucena-Cacace A, Umeda M, Navas LE, Carnero A. NAMPT as a dedifferentiation-inducer gene: NAD(+) as core axis for glioma cancer stem-like cells maintenance. Front Oncol (2019) 9:292. doi: 10.3389/fonc.2019.00292
60. Glorieux C, Calderon PB. Catalase, a remarkable enzyme: targeting the oldest antioxidant enzyme to find a new cancer treatment approach. Biol Chem (2017) 398(10):1095–108. doi: 10.1515/hsz-2017-0131
61. Zhu C, Mustafa D, Zheng PP, van der Weiden M, Sacchetti A, Brandt M, et al. Activation of CECR1 in M2-like TAMs promotes paracrine stimulation-mediated glial tumor progression. Neuro Oncol (2017) 19(5):648–59. doi: 10.1093/neuonc/now251
62. Han J, Deng X, Sun R, Luo M, Liang M, Gu B, et al. GPI is a prognostic biomarker and correlates with immune infiltrates in lung adenocarcinoma. Front Oncol (2021) 11:752642. doi: 10.3389/fonc.2021.752642
63. Bao C, Wang J, Ma W, Wang X, Cheng Y. HDGF: a novel jack-of-all-trades in cancer. Future Oncol (2014) 10(16):2675–85. doi: 10.2217/fon.14.194
64. Chateauvieux S, Grigorakaki C, Morceau F, Dicato M, Diederich M. Erythropoietin, erythropoiesis and beyond. Biochem Pharmacol (2011) 82(10):1291–303. doi: 10.1016/j.bcp.2011.06.045
65. Vatsveen TK, Sponaas AM, Tian E, Zhang Q, Misund K, Sundan A, et al. Erythropoietin (EPO)-receptor signaling induces cell death of primary myeloma cells in vitro. J Hematol Oncol (2016) 9(1):75. doi: 10.1186/s13045-016-0306-x
66. Vacchelli E, Ma Y, Baracco EE, Sistigu A, Enot DP, Pietrocola F, et al. Chemotherapy-induced antitumor immunity requires formyl peptide receptor 1. Science (2015) 350(6263):972–8. doi: 10.1126/science.aad0779
67. Saito N, Qiao H, Yanagi T, Shinkuma S, Nishimura K, Suto A, et al. An annexin A1-FPR1 interaction contributes to necroptosis of keratinocytes in severe cutaneous adverse drug reactions. Sci Transl Med (2014) 6(245):245ra95. doi: 10.1126/scitranslmed.3008227
68. Shouval DS, Ouahed J, Biswas A, Goettel JA, Horwitz BH, Klein C, et al. Interleukin 10 receptor signaling: master regulator of intestinal mucosal homeostasis in mice and humans. Adv Immunol (2014) 122:177–210. doi: 10.1016/B978-0-12-800267-4.00005-5
69. Oshima H, Ishikawa T, Yoshida GJ, Naoi K, Maeda Y, Naka K, et al. TNF-alpha/TNFR1 signaling promotes gastric tumorigenesis through induction of Noxo1 and Gna14 in tumor cells. Oncogene (2014) 33(29):3820–9. doi: 10.1038/onc.2013.356
70. Nirmala JG, Lopus M. Cell death mechanisms in eukaryotes. Cell Biol Toxicol (2020) 36(2):145–64. doi: 10.1007/s10565-019-09496-2
71. Khoury MK, Gupta K, Franco SR, Liu B. Necroptosis in the pathophysiology of disease. Am J Pathol (2020) 190(2):272–85. doi: 10.1016/j.ajpath.2019.10.012
72. Su Z, Yang Z, Xie L, DeWitt JP, Chen Y. Cancer therapy in the necroptosis era. Cell Death Differ (2016) 23(5):748–56. doi: 10.1038/cdd.2016.8
73. Sidera K, Patsavoudi E. HSP90 inhibitors: current development and potential in cancer therapy. Recent Pat Anticancer Drug Discovery (2014) 9(1):1–20. doi: 10.2174/15748928113089990031
74. Takayama S, Reed JC, Homma S. Heat-shock proteins as regulators of apoptosis. Oncogene (2003) 22(56):9041–7. doi: 10.1038/sj.onc.1207114
75. Burrows F, Zhang H, Kamal A. Hsp90 activation and cell cycle regulation. Cell Cycle (2004) 3(12):1530–6. doi: 10.4161/cc.3.12.1277
76. Pratt WB. The role of the hsp90-based chaperone system in signal transduction by nuclear receptors and receptors signaling via MAP kinase. Annu Rev Pharmacol Toxicol (1997) 37:297–326. doi: 10.1146/annurev.pharmtox.37.1.297
77. Li D, Xu T, Cao Y, Wang H, Li L, Chen S, et al. A cytosolic heat shock protein 90 and cochaperone CDC37 complex is required for RIP3 activation during necroptosis. Proc Natl Acad Sci U S A (2015) 112(16):5017–22. doi: 10.1073/pnas.1505244112
78. Jacobsen AV, Lowes KN, Tanzer MC, Lucet IS, Hildebrand JM, Petrie EJ, et al. HSP90 activity is required for MLKL oligomerisation and membrane translocation and the induction of necroptotic cell death. Cell Death Dis (2016) 7(1):e2051. doi: 10.1038/cddis.2015.386
79. Zhao XM, Chen Z, Zhao JB, Zhang PP, Pu YF, Jiang SH, et al. Hsp90 modulates the stability of MLKL and is required for TNF-induced necroptosis. Cell Death Dis (2016) 7(2):e2089. doi: 10.1038/cddis.2015.390
80. Riegler AN, Brissac T, Gonzalez-Juarbe N, Orihuela CJ. Necroptotic cell death promotes adaptive immunity against colonizing pneumococci. Front Immunol (2019) 10:615. doi: 10.3389/fimmu.2019.00615
81. Chan IH, Wu V, Bilardello M, Mar E, Oft M, Van Vlasselaer P, et al. The potentiation of IFN-gamma and induction of cytotoxic proteins by pegylated IL-10 in human CD8 T cells. J Interferon Cytokine Res (2015) 35(12):948–55. doi: 10.1089/jir.2014.0221
82. Xi J, Xu M, Song Z, Li H, Xu S, Wang C, et al. Stimulatory role of interleukin 10 in CD8(+) T cells through STATs in gastric cancer. Tumour Biol (2017) 39(5):1010428317706209. doi: 10.1177/1010428317706209
83. Mumm JB, Emmerich J, Zhang X, Chan I, Wu L, Mauze S, et al. IL-10 elicits IFNgamma-dependent tumor immune surveillance. Cancer Cell (2011) 20(6):781–96. doi: 10.1016/j.ccr.2011.11.003
84. Liu X, Liu L, Wang X, Jin Y, Wang S, Xie Q, et al. Necroptosis inhibits autophagy by regulating the formation of RIP3/p62/Keap1 complex in shikonin-induced ROS dependent cell death of human bladder cancer. Phytomedicine (2023) 118:154943. doi: 10.1016/j.phymed.2023.154943
85. Orecchioni M, Ghosheh Y, Pramod AB, Ley K. Macrophage Polarization: Different Gene Signatures in M1(LPS+) vs. Classically and M2(LPS-) vs. Alternatively Activated Macrophages. Front Immunol (2019) 10:1084. doi: 10.3389/fimmu.2019.01084
86. Liu PS, Wang H, Li X, Chao T, Teav T, Christen S, et al. alpha-ketoglutarate orchestrates macrophage activation through metabolic and epigenetic reprogramming. Nat Immunol (2017) 18(9):985–94. doi: 10.1038/ni.3796
87. Zhou B, Liu J, Zeng L, Zhu S, Wang H, Billiar TR, et al. Extracellular SQSTM1 mediates bacterial septic death in mice through insulin receptor signalling. Nat Microbiol (2020) 5(12):1576–87. doi: 10.1038/s41564-020-00795-7
88. Sierra-Filardi E, Nieto C, Domínguez-Soto A, Barroso R, Sánchez-Mateos P, Puig-Kroger A, et al. CCL2 shapes macrophage polarization by GM-CSF and M-CSF: identification of CCL2/CCR2-dependent gene expression profile. J Immunol (2014) 192(8):3858–67. doi: 10.4049/jimmunol.1302821
89. Roca H, Varsos ZS, Sud S, Craig MJ, Ying C, Pienta KJ. CCL2 and interleukin-6 promote survival of human CD11b+ peripheral blood mononuclear cells and induce M2-type macrophage polarization. J Biol Chem (2009) 284(49):34342–54. doi: 10.1074/jbc.M109.042671
90. Xie Y, Chen Z, Zhong Q, Zheng Z, Chen Y, Shangguan W, et al. M2 macrophages secrete CXCL13 to promote renal cell carcinoma migration, invasion, and EMT. Cancer Cell Int (2021) 21(1):677. doi: 10.1186/s12935-021-02381-1
91. Wheeler KC, Jena MK, Pradhan BS, Nayak N, Das S, Hsu CD, et al. VEGF may contribute to macrophage recruitment and M2 polarization in the decidua. PLoS One (2018) 13(1):e0191040. doi: 10.1371/journal.pone.0191040
92. Wang X, Chen T, Leng L, Fan J, Cao K, Duan Z, et al. MIF produced by bone marrow-derived macrophages contributes to teratoma progression after embryonic stem cell transplantation. Cancer Res (2012) 72(11):2867–78. doi: 10.1158/0008-5472.CAN-11-3247
93. Zheng Y, Hao S, Xiang C, Han Y, Shang Y, Zhen Q, et al. The correlation between SPP1 and immune escape of EGFR mutant lung adenocarcinoma was explored by bioinformatics analysis. Front Oncol (2021) 11:592854. doi: 10.3389/fonc.2021.592854
94. Pang X, Shang H, Deng B, Wen F, Zhang Y. The interaction of adrenomedullin and macrophages induces ovarian cancer cell migration via activation of RhoA signaling pathway. Int J Mol Sci (2013) 14(2):2774–87. doi: 10.3390/ijms14022774
95. Walters DK, Mercher T, Gu TL, O'Hare T, Tyner JW, Loriaux M, et al. Activating alleles of JAK3 in acute megakaryoblastic leukemia. Cancer Cell (2006) 10(1):65–75. doi: 10.1016/j.ccr.2006.06.002
96. Koo GC, Tan SY, Tang T, Poon SL, Allen GE, Tan L, et al. Janus kinase 3-activating mutations identified in natural killer/T-cell lymphoma. Cancer Discovery (2012) 2(7):591–7. doi: 10.1158/2159-8290.CD-12-0028
97. McKinney M, Moffitt AB, Gaulard P, Travert M, De Leval L, Nicolae A, et al. The genetic basis of hepatosplenic T-cell lymphoma. Cancer Discovery (2017) 7(4):369–79. doi: 10.1158/2159-8290.CD-16-0330
98. Nairismagi ML, Tan J, Lim JQ, Nagarajan S, Ng CC, Rajasegaran V, et al. JAK-STAT and G-protein-coupled receptor signaling pathways are frequently altered in epitheliotropic intestinal T-cell lymphoma. Leukemia (2016) 30(6):1311–9. doi: 10.1038/leu.2016.13
99. Bains T, Heinrich MC, Loriaux MM, Beadling C, Nelson D, Warrick A, et al. Newly described activating JAK3 mutations in T-cell acute lymphoblastic leukemia. Leukemia (2012) 26(9):2144–6. doi: 10.1038/leu.2012.74
100. Lee CH, Wu CL, Shiau AL. Toll-like receptor 4 signaling promotes tumor growth. J Immunother (2010) 33(1):73–82. doi: 10.1097/CJI.0b013e3181b7a0a4
101. Karin M, Greten FR. NF-kappaB: linking inflammation and immunity to cancer development and progression. Nat Rev Immunol (2005) 5(10):749–59. doi: 10.1038/nri1703
102. Bernard V, Semaan A, Huang J, San Lucas FA, Mulu FC, Stephens BM, et al. Single-cell transcriptomics of pancreatic cancer precursors demonstrates epithelial and microenvironmental heterogeneity as an early event in neoplastic progression. Clin Cancer Res (2019) 25(7):2194–205. doi: 10.1158/1078-0432.CCR-18-1955
103. Ding X, Yang W, Shi X, Du P, Su L, Qin Z, et al. TNF receptor 1 mediates dendritic cell maturation and CD8 T cell response through two distinct mechanisms. J Immunol (2011) 187(3):1184–91. doi: 10.4049/jimmunol.1002902
104. Liu C, Wang L, Jiang Q, Zhang J, Zhu L, Lin L, et al. Hepatoma-derived growth factor and DDX5 promote carcinogenesis and progression of endometrial cancer by activating beta-catenin. Front Oncol (2019) 9:211. doi: 10.3389/fonc.2019.00211
105. Zhang C, Chang X, Chen D, Yang F, Li Z, Li D, et al. Downregulation of HDGF inhibits the tumorigenesis of bladder cancer cells by inactivating the PI3K-AKT signaling pathway. Cancer Manag Res (2019) 11:7909–23. doi: 10.2147/CMAR.S215341
106. Zhou Y, Zhou N, Fang W, Huo J. Overexpressed HDGF as an independent prognostic factor is involved in poor prognosis in Chinese patients with liver cancer. Diagn Pathol (2010) 5:58. doi: 10.1186/1746-1596-5-58
107. Abe R, Peng T, Sailors J, Bucala R, Metz CN. Regulation of the CTL response by macrophage migration inhibitory factor. J Immunol (2001) 166(2):747–53. doi: 10.4049/jimmunol.166.2.747
108. Tillmann S, Bernhagen J, Noels H. Arrest functions of the MIF ligand/receptor axes in atherogenesis. Front Immunol (2013) 4:115. doi: 10.3389/fimmu.2013.00115
109. Vargas-Caraveo A, Sayd A, Robledo-Montaña J, Caso JR, Madrigal JLM, García-Bueno B, et al. Toll-like receptor 4 agonist and antagonist lipopolysaccharides modify innate immune response in rat brain circumventricular organs. Neuroinflamm (2020) 17(1):6. doi: 10.1186s12974-019-1690-2
110. Ledeboer A, Liu T, Shumilla JA, Mahoney JH, Vijay S, Gross MI, et al. The glial modulatory drug AV411 attenuates mechanical allodynia in rat models of neuropathic pain. Neuron Glia Biol (2006) 2(4):279–91. doi: 10.1017/S1740925X0700035X
Keywords: scRNA-seq, PDAC, TIMs, necroptosis, TME
Citation: Dong W, Zhao H, Xiao S, Zheng L, Fan T, Wang L, Zhang H, Hu Y, Yang J, Wang T and Xiao W (2023) Single-cell RNA-seq analyses inform necroptosis-associated myeloid lineages influence the immune landscape of pancreas cancer. Front. Immunol. 14:1263633. doi: 10.3389/fimmu.2023.1263633
Received: 20 July 2023; Accepted: 22 November 2023;
Published: 12 December 2023.
Edited by:
Takaji Matsutani, Maruho, JapanReviewed by:
Qiming Gong, Affiliated Hospital of Youjiang Medical University for Nationalities, ChinaMing Zheng, Academy of Military Medical Sciences, China
Copyright © 2023 Dong, Zhao, Xiao, Zheng, Fan, Wang, Zhang, Hu, Yang, Wang and Xiao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wenhua Xiao, d19oeGlhb0Bob3RtYWlsLmNvbQ==; Tao Wang, d2FuZ3RhbzExMTdAZ21haWwuY29t
†These authors have contributed equally to this work