Abstract
The tumor microenvironment (TME) is comprised of tumor cells, infiltrating immune cells, and stroma. Multiple reports suggest that the immune cell infiltration (ICI) in TME is strongly associated with responsiveness to immunotherapy and prognosis of certain cancers. Thus far, the ICI profile of pancreatic carcinoma (PC) remains unclear. Here, we employed two algorithms to characterize the ICI profile of PC patients. Based on our results, we identified 2 ICI patterns and calculated the ICI score by using principal component analysis. Furthermore, we revealed that patients with low ICI scores had a better prognosis, compared to high ICI scores. Moreover, we discovered that a low tumor mutation burden (TMB) offered better overall survival (OS), relative to high TMB. In this study, a high ICI score referred to elevated PD-L1/TGF-β levels, increased activation of cell cycle pathway and DNA repair pathway, as well as reduced expression of immune-activation-related genes. We also demonstrated that three metabolic pathways were suppressed in the low ICI score group. These data may explain why a high ICI score equates to a poor prognosis. Based on our analysis, the ICI score can be used as an effective predictor of PC prognosis. Hence, establishing an ICI profile, based on a large patient population, will not only enhance our knowledge of TME but also aid in the development of immunotherapies specific to PC.
Introduction
Pancreatic carcinoma (PC) is an aggressive cancer that provides a 5-year overall survival (OS) rate to <9% of affected patients (–). Despite advancements in the field, PC prognosis remains poor, due to late diagnosis and restrictive treatment strategies (, ). Immunotherapy is an approach that enhances host immunity, which, in turn, targets and eliminates tumor cells (–). It has ushered in new advancements in anti-tumor treatment, including PC therapy, carrying enormous survival benefits (, , ). However, only a few PC patients are responsive to immunotherapy (–). Thus, exploring new predictive biomarkers for PC patients’ prognosis and developing novel therapeutic strategies are of urgent need.
The lack of patient response to immunotherapy likely stems from our limited understanding of the tumor microenvironment (TME). TME is comprised of tumor-infiltrating immune cells (TIICs) and stroma, and it is specific to individual cancers (). Hence, PC-specific TME has a unique immune cell infiltration (ICI) and characteristically desmoplastic stroma (, ). Several studies on TME showed a strong correlation between ICI and tumor growth, metastasis, and sensitivity to immunotherapy (). One example is tumor-associated macrophages (TAMs), which drive tumorigenesis, via the release of immunosuppressive cytokines, like TGF-β, thereby promoting poor prognosis (, ). Alternately, elevated levels of CD4+ T cells and CD8+ T cell infiltration are associated with better prognosis and responsiveness to immunotherapy (). However, TIICs alone do not regulate sensitivity to immunotherapy. In fact, in some studies, patients with high TIICs were shown to be immunotherapy-resistant, likely due to TIICs dysfunction brought on by TAM-regulated immunosuppressive cytokines (, ). In addition, a massive infiltration of stromal cells can further block TIICs infiltration into tumor tissue (). These data suggest that immunotherapy responsiveness is strongly modulated by intrinsic events in TME and not by the action of an individual cell population. Thus far, there are few reports on an extensive ICI profile in TME of PC patients. In the last few decades, advancements in next-generation sequencing (NGS) technology, particularly the NGS algorithm technology, have enabled the extraction of large amounts of biological data on PC formation and metastasis ().
Here, we employed two algorithms, CIBERSORT and ESTIMATE, to assess the gene-expression profiles of a large number of tumor samples and establish an extensive overview of the PC immune landscape (, ). In addition, we categorized PC patients into two distinct subgroups, based on the ICI patterns. Finally, we established an ICI score to characterize the immune cell infiltration profile of PC, which could accurately predict patient prognosis and be of benefit to the selection of subsequent immunotherapeutic strategies.
Materials and Methods
An outline of the study, along with the data collection procedure, is summarized in Figure 1.
Figure 1
PC Datasets and Samples
The datasets analyzed in the current study are available at The Cancer Genome Atlas (TCGA) at https://portal.gdc.cancer.gov, reference number TCGA-PAAD, and Gene Expression Omnibus (GEO) at https://www.ncbi.nlm.nih.gov/geo/, under the accession number GSE57495.
Consensus Clustering for TIICs
The “CIBERSORT” R package, carrying an LM22 signature and 1000 permutations, was employed for the quantification of individual immune cell infiltrating levels in PC (). ESTIMATE was used to establish the immune and stromal scores in PC samples (). The hierarchical agglomerative clustering of PC was completed per ICI pattern for individual PC samples. Subsequently, the unsupervised clustering “Pam” technique, closely following the Euclidean and Ward’s linkage, was employed, using the “ConsensuClusterPlus” R package (), and replicated 1,000X to verify classification stability.
Differentially Expressed Genes (DEGs) Related to the ICI Phenotype
Patients were grouped into ICI clusters, according to their ICI profile, to identify genes associated with the ICI patterns. DEGs, among ICI subtypes, were identified by using the “limma” R package, and with the following criteria: P < 0.05 and absolute fold-change > 1.4.
Dimension Reduction and ICI Score Computation
Next, using unsupervised clustering, the patients were stratified, according to their DEG values. The positively- and negatively-regulated DEG values were further classified as ICI signature gene (ICISG) A (ICISG-A) or ICISG-B, respectively. Moreover, the Boruta algorithm was employed for the dimension reduction of ICISG-A and ICISG-B (), and the principal component 1 was retrieved as the signature score, using PCA. Finally, the ICI scores of individual patients were obtained using a technique similar to the gene expression grade index (), as described below:
Somatic Data Collection
First, we downloaded the TCGA cohort’s genetic information from the TCGA data portal. Next, we calculated the sum of non-synonymous mutations in PC. The somatic changes, in PC driver genes, were next assessed against high and low ICI scores. The PC driver genes were recognized using the “maftool” R package (). Lastly, the leading 20 driver DEGs were chosen for further analysis.
Statistical Analyses
GraphPad Prism version7.0 or SPSS version 21.0 (IBM Corporation, Armonk, NY, USA) was employed for all statistical analyses. Comparison between groups was performed with the Wilcoxon test. The X-tile software, which iteratively selects potential threshold to maximize rank statistic, was employed for PC patient classification into subgroups and to decrease computational batch outcome (). The Kaplan-Meier plotter produced the OS curves for each examined category. The log-rank test analyzed significance. The chi-square test examined the association between the ICI scores and somatic mutation frequency. The Spearman analysis calculated the correlation coefficient. Two-tailed P< 0.05 was the statistical significance threshold.
Results
The ICI Profile in PC TME
245 tumor samples were obtained from the TCGA-PAAD and GEO-GSE57495 cohorts. Infiltrating immune cells were quantified in PC tissues, using both CIBERSORT and ESTIMATE algorithms. Next, these patients were categorized into discrete subgroups, via unsupervised clustering, using the “ConsesusClusterPlus” package of the R software.
Based on our analysis, two isolated ICI subtypes were found (Figures 2A, B). To delineate the intrinsic changes that lead to the distinct phenotypes, we further analyzed the ICI pattern of PC TME. As depicted in Figure 2C, ICI cluster A had a large population of regulatory T cells, M0 macrophages, and activated mast cells infiltration, whereas ICI cluster B had more naive B cells, CD8+ T cells, resting memory CD4+ T cells, activated memory CD4+ T cells, resting NK cells, monocytes, M1 macrophages, resting dendritic cells, and resting mast cells. Additionally, the ICI cluster B phenotype also exhibited elevated immune and stromal scores. We, next, compiled the correlation coefficient heatmap to illustrate the ICI profile in PC TME (Figure 2D).
Figure 2
Identification of the Immune Genetics Subtype
To elucidate the genetic features of discrete immune-related phenotypes, we conducted genetic differential analyses, using the “limma” program of R software. Using unsupervised clustering of 114 DEGs (obtained from prior differential analyses), we categorized the TCGA and GEO cohorts into two genomic clusters, namely gene clusters A and B (Figure 3A). Among the 114 DEGs, the 22 genes that were positively correlated with the gene cluster in PC patients were placed in the ICISG-A category and the remaining DEGs were placed in the ICISG-B category. Simultaneously, to eliminate noise and the presence of repetitive genes, we employed the Boruta algorithm to conduct dimension reduction in ICISG-A and ICISG-B. We, also, generated a heatmap illustrating the gene expression profile of 114 DEGs within ICISG-A and ICISG-B (Figure 3B), with the “clusterProfiler” R program. As illustrated in Figures 3C, D, we summarized significantly enriched biological processes in both ICISG-A and ICISG-B, using the Gene Ontology (GO) enrichment analysis. Furthermore, gene cluster A contained a large population of CD8+ T cells, monocytes, M1 macrophages, resting dendritic cells, and resting mast cells infiltration. Alternately, gene cluster B had massive infiltration of M0 macrophages, activated dendritic cells, and activated mast cells (Figure 3E).
Figure 3
ICI Score Computation and Association Between ICI Score and Somatic Variation
To establish quantitative indicators for the ICI profile, we employed the principal-component analysis (PCA) to calculate 2 cluster scores: (1) the ICI score A from ICISG-A and (2) the ICI score B from ICISG-B. The scores of ICI score A and ICI score B were calculated for individual patients, based on a compilation of related scores. Subsequently, we obtained the ICI score, which was used as an estimator of PC prognosis. Patients in the TCGA and GEO cohorts were then classified into 2 categories, according to an optimal ICI score threshold, produced by the X-tile software. To assess the immunological activities and drug sensitivity of each group, we chose for evaluation CD274 (PD-L1) and TGFB1 (TGF-β) as the immune checkpoint-related genes, and CCL19, CCR7, CD3D, CD3E, CD79B, IL33, CD8A, CXCL9, and IDO1 as immune-activation-related genes. Accordingly, we demonstrated that immune-activation-related genes were markedly upregulated in the low ICI score group, relative to the high ICI score group (Figure 4A). Additionally, the high ICI group exhibited increased levels of TGF-β and PD-L1 relative to the low ICI score group, as evidenced by the Wilcoxon test (Figure 4A). In addition, the gene set enrichment analysis (GSEA) revealed that the cell cycle and DNA repair pathways, particularly, mismatch repair and nucleotide excision repair, were significantly enriched and activated in the high ICI score group, as opposed to the low ICI score group. Conversely, three metabolic pathways (glycerolipid metabolism; fatty acid metabolism; glycine, serine and threonine metabolism) were significantly enriched and inhibited in the low ICI score group (Figure 4B).
Figure 4
Previous evidence revealed that tumors with a high tumor burden mutation (TMB, non-synonymous variants) contributed to an elevated number of infiltrating CD8+ T cells that target and destroy tumors. Based on this, TMB may be an influencing factor in patient prognosis and responsiveness to cancer immunotherapy (, ). The KEYNOTE 012 clinical trial showed that a high TMB was correlated with increased sensitivity to PD-1 blockades and a favorable prognosis in PC patients (). Here, we evaluated somatic PC driver gene variants’ distribution in low and high ICI score subgroups, using “maftools”. The top 20 DEGs were then selected for further analysis (Figures 4C, D). We demonstrated that the levels of KRAS, TP53, CDKN2A, and SMAD4 were markedly altered in the low and high ICI score subgroups (chi-square test, P<0.05). Given the clinical significance of TMB, we, next, examined the intrinsic relationship between the TMB and ICI score to assess genetic imprints of individual ICI score subgroups. As illustrated in Figure 4E, our correlation analysis validated that the ICI score was negatively associated with TMB (Spearman coefficient: R = -0.25, P = 0.0025).
The Significance of ICI Score in PC Prognosis
We, next, examined the prognostic capabilities of the ICI score. Using the Kaplan-Meier plot, we demonstrated that the low ICI score group had markedly improved OS rate, compared to the high ICI score group (median survival time 738 vs. 460days, log-rank test, P < 0.0001). Additionally, the ICI score prognosis accuracy was further confirmed in TCGA-PAAD and GEO-GSE57495, respectively (Figures 5A–C).
Figure 5
We also assigned patients to distinct groups, according to the TMB immune threshold, as described previously. Based on our analysis, patients with low TMB experienced improved OS, relative to high TMB (log-rank test, P = 0.005, Figure 5D). Given that TMB may have an impact on the prognostic capabilities of the ICI score, we next examined the collaborative impact of ICI score and TMB in predicting prognosis. The stratified survival analysis demonstrated that TMB failed to regulate ICI score-based prognosis estimation. In fact, the ICI score exhibited marked OS differences in both high and low TMB groups (log-rank test, Low TMB & Low ICI score (LL) versus Low TMB & High ICI score (LH); High TMB & Low ICI score (HL) versus High TMB & High ICI score (HH), Figure 5E). Collectively, these data are indicative of the predictive potential of the ICI score in PC prognosis, independent of TMB.
Discussion
Immunotherapy has remarkable efficacy in multiple malignancies. However, only a small population of PC patients respond positively to immunotherapy, owing to the complicated TME (). Therefore, it is urgent and necessary to explore the TME in PC patients in order to predict prognosis, design personalized therapy, and aid in the development of new PC therapies. In a previous study, researchers examined the ICI profile of PC, stratified PC patients according to infiltrating T-cell activity in TME and correlated cytolytic immune activity with mutational, structural, and neoepitope features of the tumor (). The previous study provided insight into the intrinsic activities of PC TME. However, we aimed to explore the PC TME from a new perspective. Here, we developed a methodology for quantifying a comprehensive PC TME. Based on our analysis, we recommend the ICI score to be a promising prognostic biomarker for PC.
Different types of tumors vary widely in their TME contexture, especially in ICI (). The TME heterogeneity, which impacts tumor progression and prognosis, has been identified in both human and murine PC (, –). PC is classified into different subtypes based on the TME heterogeneity. Moreover, an abundance of intratumoral CD8+ T cells is closely correlated with the response to immunotherapy in PC (–). However, multiple clinical and genomic studies have suggested that PC, unlike other types of cancer, is correlated with low ICI, especially CD8+ T cells (, , ). This may explain why only a minority of PC patients are sensitive to immunotherapy, relative to other types of tumors with high CD8+ T cell infiltration. In this study, we evaluated the ICI in 245 PC samples and assigned the samples to two discrete immune categories and gene clusters. We demonstrated that the levels of CD4+ T cells, CD8+ T cells, M1 macrophages, and the high immune score didn’t predict PC prognosis, unlike previous reports (, ). The PC TME is complex and may interfere with the intercellular interactions between infiltrating immune cells, thereby affecting immunity tolerance and activity (, ). Given these circumstances, the immune phenotypes and gene clusters may not serve as an appropriate biomarker to aid in predict PC prognosis and sensitivity to immunotherapy (Figures S1A, B).
Given the diverse nature of the immune milieu in PC patients, it was critical to characterize the ICI patterns of individual PC patients. The individual-based model, derived from tumor subtype-specific biomarkers, has been widely used in predicting prognosis in breast and colorectal cancers (, ). Here, using the Boruta algorithm, we identified PC biomarkers and generated an individualized ICI score to assess ICI patterns. With the help of GSEA, we also discovered three metabolic pathways (glycerolipid metabolism; fatty acid metabolism; glycine, serine and threonine metabolism) were enriched and these metabolism‐related genes were lowly expressed in the low ICI score group, relative to the high ICI score group. Emerging evidence has revealed that oncogenes can induce metabolic alterations in tumor cells and TIICs, that can restrict immune responses to cancer therapy and promote poor prognosis (, ). Previous studies have suggested that, in TME, competition for nutrients between cancer cells and T cells contributes to immunosuppression (, ). Moreover, a recent study evaluated metabolic features of tumor cell types in vivo and revealed that individual cell populations had distinct programs of nutrient uptake that might serve an important role in the development of future cancer therapies by altering TME (). In addition, we, also, observed markedly elevated expression of cell cycle- and DNA repair- (mismatch repair and nucleotide excision repair) related genes in the high ICI score group. Recent studies indicated that the cell cycle activity of both cancer and immune cells in TME could mediate antitumor immunity (–). Meanwhile, cell cycle inhibition can result in both cell autonomous (e.g., dysregulation of the antigen presentation machinery), and non-cell autonomous (e.g., release of SASP-associated factors and T cell recruiting chemokines) mechanisms of antitumor immunity (). This finding has strong implications for PC therapy, as it suggests that existing drugs that modulate the cell cycle, may potentially have an added (and untapped) benefit of sensitizing tumors to immunotherapy. A complex DNA repair machinery has evolved to protect genomic integrity in the face of numerous sources of DNA damage. When DNA repair fails, this damage can lead to carcinogenesis and tumor genomic instability (). Although the coordinated activities of DNA repair pathways can quickly correct most DNA damages, delayed or improper repairs can lead to changes in the tumor genome, thereby reconstructing the TME (, ). A large amount of evidence shows that DNA repair plays an important role in driving sensitivity and response to immunotherapy (, ).
High PD-L1 levels are common in tumor cells, but not normal cells (, ). Several cancer types can exhibit immunosuppressed TME, along with elevated levels of PD-L1, which inhibits T-cell-mediated cytotoxicity of tumor cells (57). This may contribute to the correlation of elevated PD-L1 levels to poor prognosis. In a meta-analysis study involving PC patients, 19-62.5% of patients exhibited elevated PD-L1 levels and corresponding a poor prognosis, as opposed to those with low PD-L1 levels (58). In accordance with other studies, we also demonstrated that a high ICI score, with elevated PD-L1 levels, had a poor prognosis, relative to a low ICI score. Moreover, immunotherapies like PD-1/PD-L1 antibodies were found to be highly efficient in treating other forms of cancers but exhibited little success in PC therapy (, 59). The specific insensitivity to immunotherapy, observed with PC patients, may be due to the PC-specific TME (, 59). In PC, abundant extracellular matrix (ECM) surrounds tumor cells and creates a physical barrier that blocks entry of drugs and cytotoxic T cells (60). As such, PC therapy must first involve manipulation of the PC stroma to allow for the infiltration of T cells and progression of antitumor immunity. Alternately, TGF-β is involved in the pancreatic stellate cell-mediated secretion of the stiff fibrillary ECM, that sustains tumor survival (61, 62). Moreover, recently, it was suggested that TGF-β also protected tumor cells by restricting infiltrating T cells (63). Based on our analysis, TGF-β was highly expressed in the high ICI score group, which may have contributed to subsequent poor PC prognosis and immunotherapy failure. Single-agent immunotherapies are therefore ineffective in PC. Several reports have suggested that in tumors with low T cell infiltration, the combined blockade of the PD-L1 checkpoint and TGF-β signaling pathway can enhance CD8+ T cell infiltration in TME and stimulate strong anti-tumor immunity (64, 65). As such, the synergistic suppression of TGF-β and PD-L1 pathways may be an effective therapy to curb PC.
In addition, in most tumors, a high mutation load is equivalent to an increased amount of tumor antigens, which can regulate survival benefits by reshaping the TME, and can also be used as a biomarker for immunotherapy responsiveness and prognosis (66, 67). In terms of mechanism, high TMB provides more opportunities for the generation of “non-self” neoantigens, which can activate the enrichment of immune cells (68). However, these theories have only been confirmed in some immunotherapeutic hot tumors. However, in immunotherapeutic cold tumors, such as PC, these rules may not apply. PC is known to have a low tumor mutation load, which may contribute to the unresponsiveness of PC immunotherapy (69). Based on our analysis, the PC patients with a low mutation load had a better prognosis. This may be attributed to the unique population of dysfunctional antigen-presenting dendritic cells in PC, which failed to induce anti-cancer immunity. Instead, the Th17 cells released IL-17 to modulate TME in order to promote tumor growth and metastasis (70). Based on our analysis, the association between the ICI score and TMB was -0.25. Using stratified analysis, we demonstrated that the prognosis, using the ICI score, was not regulated by TMB in PC. Furthermore, the absence of an ICI score-TMB correlation is indicative of these factors modulating individual aspects of tumor immunobiology. In our analysis, the ICI score was found to be predictive of PC patient prognosis, independent of TMB.
There were several limitations to our study. First, the available public dataset for PC is limited. So, in this study 245 tumor samples were obtained from the TCGA-PAAD and GEO-GSE57495 cohorts. Hence, the sample size might be relatively small. A larger sample population will aid in a better understanding of the PC ICI profile. Second, we performed bioinformatic analysis of PC samples and developed an index, the ICI score, to stratify patients. Specifically, the ICI score may serve as a tool for the classification of PC patients and may contribute to the clinical treatment of patients. Some expressions of immune-related genes and enrichments of molecular pathways were different between the ICI score subgroups. The subgroup-specific differences could also serve as potential therapeutic targets. Further in vitro and in vivo validations will be helpful to discover these potential therapeutic targets.
Conclusion
In conclusion, we extensively evaluated the ICI profile in PC patients, enabling a clearer understanding of the anti-/pro-tumor immune response in PC. According to our analysis, the ICI score can serve as a valid prognostic biomarker of PC. Therefore, a systematic assessment of tumor ICI patterns in PC patients is essential to individualized PC therapy.
Funding
This work was supported by the Natural Science Youth Fund of Southwest Medical University [Grant Number: 2019ZQN072].
Statements
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: TCGA at https://portal.gdc.cancer.gov, reference number TCGA-PAAD, and GEO at https://www.ncbi.nlm.nih.gov/geo/, under the accession number GSE57495.
Author contributions
SS and DL conceived and designed the whole project. DL and FK analyzed the data and wrote the manuscript. JD and MZ carried out data interpretations. FP, YG, and CF provided specialized expertise and collaboration in data analysis. BL and XY revised the manuscript. All authors contributed to the article and approved the submitted version.
Conflict of interest
The 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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.677609/full#supplementary-material
Supplementary Figure 1(A) Kaplan-Meier curves for the discrete ICI clusters of the TCGA and GEO cohorts, Log-rank test, P=0.995. (B) Kaplan-Meier curves for the discrete gene clusters of the TCGA and GEO cohorts, Log-rank test, P=0.256.
References
1
NeoptolemosJPKleeffJMichlPCostelloEGreenhalfWPalmerDH. Therapeutic Developments in Pancreatic Cancer: Current and Future Perspectives. Nat Rev Gastroenterol Hepatol (2018) 15(6):333–48. doi: 10.1038/s41575-018-0005-x
2
Ilic MIlicI. Epidemiology of Pancreatic Cancer. World J Gastroenterol (2016) 22(44):9694–705. doi: 10.3748/wjg.v22.i44.9694
3
McGuiganAKellyPTurkingtonRCJonesCColemanHGMcCainRS. Pancreatic Cancer: A Review of Clinical Diagnosis, Epidemiology, Treatment and Outcomes. World J Gastroenterol (2018) 24(43):4846–61. doi: 10.3748/wjg.v24.i43.4846
4
MooreADonahueT. Pancreatic Cancer. Jama (2019) 322(14):1426. doi: 10.1001/jama.2019.14699
5
RileyRSJuneCHLangerRMitchellMJ. Delivery Technologies for Cancer Immunotherapy. Nat Rev Drug Discov (2019) 18(3):175–96. doi: 10.1038/s41573-018-0006-z
6
KrugerSIlmerMKoboldSCadilhaBLEndresSOrmannsSet al. Advances in Cancer Immunotherapy 2019 - Latest Trends. J Exp Clin Cancer Res (2019) 38(1):268. doi: 10.1186/s13046-019-1266-0
7
SanmamedMFChenL. A Paradigm Shift in Cancer Immunotherapy: From Enhancement to Normalization. Cell (2018) 175(2):313–26. doi: 10.1016/j.cell.2018.09.035
8
MorrisonAHByrneKTVonderheideRH. Immunotherapy and Prevention of Pancreatic Cancer. Trends Cancer (2018) 4(6):418–28. doi: 10.1016/j.trecan.2018.04.001
9
BanerjeeKKumarSRossKAGautamSPoelaertBNasserMWet al. Emerging Trends in the Immunotherapy of Pancreatic Cancer. Cancer Lett (2018) 417:35–46. doi: 10.1016/j.canlet.2017.12.012
10
FengMXiongGCaoZYangGZhengSSongXet al. PD-1/PD-L1 and Immunotherapy for Pancreatic Cancer. Cancer Lett (2017) 407:57–65. doi: 10.1016/j.canlet.2017.08.006
11
RenBCuiMYangGWangHFengMYouLet al. Tumor Microenvironment Participates in Metastasis of Pancreatic Cancer. Mol Cancer (2018) 17(1):108. doi: 10.1186/s12943-018-0858-1
12
LigorioMSilSMalagon-LopezJNiemanLTMisaleSDi PilatoMet al. Stromal Microenvironment Shapes the Intratumoral Architecture of Pancreatic Cancer. Cell (2019) 178(1):160–75.e27. doi: 10.1016/j.cell.2019.05.012
13
HinshawDCShevdeLA. The Tumor Microenvironment Innately Modulates Cancer Progression. Cancer Res (2019) 79(18):4557–66. doi: 10.1158/0008-5472.CAN-18-3962
14
MantovaniAMarchesiFMalesciALaghiLAllavenaP. Tumour-Associated Macrophages as Treatment Targets in Oncology. Nat Rev Clin Oncol (2017) 14(7):399–416. doi: 10.1038/nrclinonc.2016.217
15
De PalmaMLewisCE. Macrophage Regulation of Tumor Responses to Anticancer Therapies. Cancer Cell (2013) 23(3):277–86. doi: 10.1016/j.ccr.2013.02.013
16
VassilakopoulouMAvgerisMVelchetiVKotoulaVRampiasTChatzopoulosKet al. Evaluation of PD-L1 Expression and Associated Tumor-Infiltrating Lymphocytes in Laryngeal Squamous Cell Carcinoma. Clin Cancer Res (2016) 22(3):704–13. doi: 10.1158/1078-0432.CCR-15-1543
17
AyersMLuncefordJNebozhynMMurphyELobodaAKaufmanDRet al. IFN-γ-Related Mrna Profile Predicts Clinical Response to PD-1 Blockade. J Clin Invest (2017) 127(8):2930–40. doi: 10.1172/JCI91190
18
ŞenbabaoğluYGejmanRSWinerAGLiuMVan AllenEMde VelascoGet al. Tumor Immune Microenvironment Characterization in Clear Cell Renal Cell Carcinoma Identifies Prognostic and Immunotherapeutically Relevant Messenger RNA Signatures. Genome Biol (2016) 17(1):231. doi: 10.1186/s13059-016-1092-z
19
ChenDSMellmanI. Elements of Cancer Immunity and the Cancer-Immune Set Point. Nature (2017) 541(7637):321–30. doi: 10.1038/nature21349
20
MüllerSRaulefsSBrunsPAfonso-GrunzFPlötnerAThermannRet al. Next-Generation Sequencing Reveals Novel Differentially Regulated Mrnas, Lncrnas, Mirnas, Sdrnas and a Pirna in Pancreatic Cancer. Mol Cancer (2015) 14:94. doi: 10.1186/s12943-015-0401-6
21
NewmanAMLiuCLGreenMRGentlesAJFengWXuYet al. Robust Enumeration of Cell Subsets From Tissue Expression Profiles. Nat Methods (2015) 12(5):453–7. doi: 10.1038/nmeth.3337
22
KursaMBRudnickiWR. Feature Selection With the Boruta Package. J Stat Software (2010) 36(11):1–13. doi: 10.18637/jss.v036.i11
23
YoshiharaKShahmoradgoliMMartínezEVegesnaRKimHTorres-GarciaWet al. Inferring Tumour Purity and Stromal and Immune Cell Admixture From Expression Data. Nat Commun (2013) 4:2612. doi: 10.1038/ncomms3612
24
YuGWangLGHanYHeQY. Clusterprofiler: An R Package for Comparing Biological Themes Among Gene Clusters. Omics (2012) 16(5):284–7. doi: 10.1089/omi.2011.0118
25
SotiriouCWirapatiPLoiSHarrisAFoxSSmedsJet al. Gene Expression Profiling in Breast Cancer: Understanding the Molecular Basis of Histologic Grade to Improve Prognosis. J Natl Cancer Inst (2006) 98(4):262–72. doi: 10.1093/jnci/djj052
26
MayakondaALinDCAssenovYPlassCKoefflerHP. Maftools: Efficient and Comprehensive Analysis of Somatic Variants in Cancer. Genome Res (2018) 28(11):1747–56. doi: 10.1101/gr.239244.118
27
CampRLDolled-FilhartMRimmDL. X-Tile: A New Bio-Informatics Tool for Biomarker Assessment and Outcome-Based Cut-Point Optimization. Clin Cancer Res (2004) 10(21):7252–9. doi: 10.1158/1078-0432.CCR-04-0713
28
RizviNAHellmannMDSnyderAKvistborgPMakarovVHavelJJet al. Cancer Immunology. Mutational Landscape Determines Sensitivity to PD-1 Blockade in Non-Small Cell Lung Cancer. Science (2015) 348(6230):124–8. doi: 10.1126/science.aaa1348
29
McGranahanNFurnessAJRosenthalRRamskovSLyngaaRSainiSKet al. Clonal Neoantigens Elicit T Cell Immunoreactivity and Sensitivity to Immune Checkpoint Blockade. Science (2016) 351(6280):1463–9. doi: 10.1126/science.aaf1490
30
SeiwertTYBurtnessBMehraRWeissJBergerREderJPet al. Safety and Clinical Activity of Pembrolizumab for Treatment of Recurrent or Metastatic Squamous Cell Carcinoma of the Head and Neck (KEYNOTE-012): An Open-Label, Multicentre, Phase 1b Trial. Lancet Oncol (2016) 17(7):956–65. doi: 10.1016/S1470-2045(16)30066-3
31
BalliDRechAJStangerBZVonderheideRH. Immune Cytolytic Activity Stratifies Molecular Subsets of Human Pancreatic Cancer. Clin Cancer Res (2017) 23(12):3129–38. doi: 10.1158/1078-0432.CCR-16-2128
32
QuailDFJoyceJA. Microenvironmental Regulation of Tumor Progression and Metastasis. Nat Med (2013) 19(11):1423–37. doi: 10.1038/nm.3394
33
TsujikawaTKumarSBorkarRNAzimiVThibaultGChangYHet al. Quantitative Multiplex Immunohistochemistry Reveals Myeloid-Inflamed Tumor-Immune Complexity Associated With Poor Prognosis. Cell Rep (2017) 19(1):203–17. doi: 10.1016/j.celrep.2017.03.037
34
LiJByrneKTYanFYamazoeTChenZBaslanTet al. Tumor Cell-Intrinsic Factors Underlie Heterogeneity of Immune Cell Infiltration and Response to Immunotherapy. Immunity (2018) 49(1):178–93.e7. doi: 10.1016/j.immuni.2018.06.006
35
BaileyPChangDKNonesKJohnsALPatchAMGingrasMCet al. Genomic Analyses Identify Molecular Subtypes of Pancreatic Cancer. Nature (2016) 531(7592):47–52. doi: 10.1038/nature16965
36
BlandoJSharmaAHigaMGZhaoHVenceLYadavSSet al. Comparison of Immune Infiltrates in Melanoma and Pancreatic Cancer Highlights VISTA as a Potential Target in Pancreatic Cancer. Proc Natl Acad Sci USA (2019) 116(5):1692–7. doi: 10.1073/pnas.1811067116
37
MacherlaSLaksSNaqashARBulumulleAZervosEMuzaffarM. Emerging Role of Immune Checkpoint Blockade in Pancreatic Cancer. Int J Mol Sci (2018) 19(11):3505. doi: 10.3390/ijms19113505
38
YangJJiangYHeRLiuWYangMTaoLet al. DKK2 Impairs Tumor Immunity Infiltration and Correlates With Poor Prognosis in Pancreatic Ductal Adenocarcinoma. J Immunol Res (2019) 2019:8656282. doi: 10.1155/2019/8656282
39
ChangJHJiangYPillarisettyVG. Role of Immune Cells in Pancreatic Cancer From Bench to Clinical Application: An Updated Review. Med (Baltimore) (2016) 95(49):e5541. doi: 10.1097/MD.0000000000005541
40
WangSLiYXingCDingCZhangHChenLet al. Tumor Microenvironment in Chemoresistance, Metastasis and Immunotherapy of Pancreatic Cancer. Am J Cancer Res (2020) 10(7):1937–53.
41
CallariMCappellettiVD’AiutoFMusellaVLemboAPetelFet al. Subtype-Specific Metagene-Based Prediction of Outcome After Neoadjuvant and Adjuvant Treatment in Breast Cancer. Clin Cancer Res (2016) 22(2):337–45. doi: 10.1158/1078-0432.CCR-15-0757
42
BramsenJBRasmussenMHOngenHMattesenTBØrntoftMWÁrnadóttirSSet al. Molecular-Subtype-Specific Biomarkers Improve Prediction of Prognosis in Colorectal Cancer. Cell Rep (2017) 19(6):1268–80. doi: 10.1016/j.celrep.2017.04.045
43
BaderJEVossKRathmellJC. Targeting Metabolism to Improve the Tumor Microenvironment for Cancer Immunotherapy. Mol Cell (2020) 78(6):1019–33. doi: 10.1016/j.molcel.2020.05.034
44
RingelAEDrijversJMBakerGJCatozziAGarcía-CañaverasJCGassawayBMet al. Obesity Shapes Metabolism in the Tumor Microenvironment to Suppress Anti-Tumor Immunity. Cell (2020) 183(7):1848–66.e26. doi: 10.1016/j.cell.2020.11.009
45
YinZBaiLLiWZengTTianHCuiJ. Targeting T Cell Metabolism in the Tumor Microenvironment: An Anti-Cancer Therapeutic Strategy. J Exp Clin Cancer Res (2019) 38(1):403. doi: 10.1186/s13046-019-1409-3
46
JiangSYanW. T-Cell Immunometabolism Against Cancer. Cancer Lett (2016) 382(2):255–8. doi: 10.1016/j.canlet.2016.09.003
47
ReinfeldBIMaddenMZWolfMMChytilABaderJEPattersonARet al. Cell-Programmed Nutrient Partitioning in the Tumour Microenvironment. Nature (2021). doi: 10.1038/s41586-021-03442-1
48
JinXDingDYanYLiHWangBMaLet al. Phosphorylated RB Promotes Cancer Immunity by Inhibiting NF-κb Activation and PD-L1 Expression. Mol Cell (2019) 73(1):22–35.e6. doi: 10.1016/j.molcel.2018.10.034
49
ZhangJBuXWangHZhuYGengYNihiraNTet al. Cyclin D-CDK4 Kinase Destabilizes PD-L1 Via Cullin 3-SPOP to Control Cancer Immune Surveillance. Nature (2018) 553(7686):91–5. doi: 10.1038/nature25015
50
LiJStangerBZ. Cell Cycle Regulation Meets Tumor Immunosuppression. Trends Immunol (2020) 41(10):859–63. doi: 10.1016/j.it.2020.07.010
51
HelledayTPetermannELundinCHodgsonBSharmaRA. DNA Repair Pathways as Targets for Cancer Therapy. Nat Rev Cancer (2008) 8(3):193–204. doi: 10.1038/nrc2342
52
MengXDuanCPangHChenQHanBZhaCet al. DNA Damage Repair Alterations Modulate M2 Polarization of Microglia to Remodel the Tumor Microenvironment Via the P53-Mediated MDK Expression in Glioma. EBioMedicine (2019) 41:185–99. doi: 10.1016/j.ebiom.2019.01.067
53
BeverKMLeDT. DNA Repair Defects and Implications for Immunotherapy. J Clin Invest (2018) 128(10):4236–42. doi: 10.1172/JCI122010
54
MouwKWD’AndreaAD. DNA Repair Deficiency and Immunotherapy Response. J Clin Oncol (2018) 36(17):1710–3. doi: 10.1200/JCO.2018.78.2425
55
WangXTengFKongLYuJ. PD-L1 Expression in Human Cancers and Its Association With Clinical Outcomes. Onco Targets Ther (2016) 9:5023–39. doi: 10.2147/OTT.S105862
56
LarsenTVHussmannDNielsenAL. PD-L1 and PD-L2 Expression Correlated Genes in Non-Small-Cell Lung Cancer. Cancer Commun (Lond) (2019) 39(1):30. doi: 10.1186/s40880-019-0376-6
57
MoRJHanZDLiangYKYeJHWuSLLinSXet al. Expression of PD-L1 in Tumor-Associated Nerves Correlates With Reduced CD8(+) Tumor-Associated Lymphocytes and Poor Prognosis in Prostate Cancer. Int J Cancer (2019) 144(12):3099–110. doi: 10.1002/ijc.32061
58
GaoHLLiuLQiZHXuHXWangWQWuCTet al. The Clinicopathological and Prognostic Significance of PD-L1 Expression in Pancreatic Cancer: A Meta-Analysis. Hepatobiliary Pancreat Dis Int (2018) 17(2):95–100. doi: 10.1016/j.hbpd.2018.03.007
59
HenriksenADyhl-PolkAChenINielsenD. Checkpoint Inhibitors in Pancreatic Cancer. Cancer Treat Rev (2019) 78:17–30. doi: 10.1016/j.ctrv.2019.06.005
60
HoWJJaffeeEMZhengL. The Tumour Microenvironment in Pancreatic Cancer - Clinical Challenges and Opportunities. Nat Rev Clin Oncol (2020) 17(9):527–40. doi: 10.1038/s41571-020-0363-5
61
KangYRoifeDLeeYLvHSuzukiRLingJet al. Transforming Growth Factor-β Limits Secretion of Lumican by Activated Stellate Cells Within Primary Pancreatic Adenocarcinoma Tumors. Clin Cancer Res (2016) 22(19):4934–46. doi: 10.1158/1078-0432.CCR-15-2780
62
PapageorgisPStylianopoulosT. Role of TGFβ in Regulation of the Tumor Microenvironment and Drug Delivery (Review). Int J Oncol (2015) 46(3):933–43. doi: 10.3892/ijo.2015.2816
63
MariathasanSTurleySJNicklesDCastiglioniAYuenKWangYet al. TGFβ Attenuates Tumour Response to PD-L1 Blockade by Contributing to Exclusion of T Cells. Nature (2018) 554(7693):544–8.doi: 10.1038/nature25501
64
WangYGaoZDuXChenSZhangWWangJet al. Co-Inhibition of the TGF-β Pathway and the PD-L1 Checkpoint by Ph-Responsive Clustered Nanoparticles for Pancreatic Cancer Microenvironment Regulation and Anti-Tumor Immunotherapy. Biomater Sci (2020) 8(18):5121–32. doi: 10.1039/D0BM00916D
65
LindHGameiroSRJochemsCDonahueRNStraussJGulleyJMet al. Dual Targeting of TGF-β and PD-L1 Via a Bifunctional Anti-PD-L1/TGF-βrii Agent: Status of Preclinical and Clinical Advances. J Immunother Cancer (2020) 8(1):e000433. doi: 10.1136/jitc-2019-000433
66
ChanTAYarchoanMJaffeeESwantonCQuezadaSAStenzingerAet al. Development of Tumor Mutation Burden as an Immunotherapy Biomarker: Utility for the Oncology Clinic. Ann Oncol (2019) 30(1):44–56. doi: 10.1093/annonc/mdy495
67
YarchoanMHopkinsAJaffeeEM. Tumor Mutational Burden and Response Rate to PD-1 Inhibition. N Engl J Med (2017) 377(25):2500–1. doi: 10.1056/NEJMc1713444
68
ChalmersZRConnellyCFFabrizioDGayLAliSMEnnisRet al. Analysis of 100,000 Human Cancer Genomes Reveals the Landscape of Tumor Mutational Burden. Genome Med (2017) 9(1):34. doi: 10.1186/s13073-017-0424-2
69
OttPABangYJPiha-PaulSARazakARABennounaJSoriaJCet al. T-Cell-Inflamed Gene-Expression Profile, Programmed Death Ligand 1 Expression, and Tumor Mutational Burden Predict Efficacy in Patients Treated With Pembrolizumab Across 20 Cancers: KEYNOTE-028. J Clin Oncol (2019) 37(4):318–27. doi: 10.1200/JCO.2018.78.2276
70
HegdeSKrisnawanVEHerzogBHZuoCBredenMAKnolhoffBLet 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
Summary
Keywords
immune cell infiltration, prognostic biomarker, pancreatic carcinoma, tumor microenvironment, immunotherapy
Citation
Luo D, Kuang F, Du J, Zhou M, Peng F, Gan Y, Fang C, Yang X, Li B and Su S (2021) Characterization of the Immune Cell Infiltration Profile in Pancreatic Carcinoma to Aid in Immunotherapy. Front. Oncol. 11:677609. doi: 10.3389/fonc.2021.677609
Received
08 March 2021
Accepted
19 April 2021
Published
13 May 2021
Volume
11 - 2021
Edited by
Liangrong Shi, Central South University, China
Reviewed by
Jinyang Li, The Rockefeller University, United States; Xiaodong Li, First People’s Hospital of Changzhou, China
Updates
Copyright
© 2021 Luo, Kuang, Du, Zhou, Peng, Gan, Fang, Yang, Li and Su.
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: Song Su, 13882778554@163.com
†These authors have contributed equally to this work
This article was submitted to Cancer Immunity and Immunotherapy, a section of the journal Frontiers in Oncology
Disclaimer
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.