Identification of the immune cell infiltration landscape in pancreatic cancer to assist immunotherapy.

Background: A malignant tumor's immune environment, including infiltrating immune cell status, can be critical to patient outcomes. Recent studies have shown that immune cell infiltration (ICI) in pancreatic cancer (PC) is highly correlated with the response to immunotherapy and patient prognosis. Therefore, we aimed to create an ICI score that accurately predicts patient outcomes and immunotherapeutic efficacy. Methods: The ICI statuses of patients with PC were estimated from the publicly available The Cancer Genome Atlas (TCGA) pancreatic ductal adenocarcinoma and GSE57495 gene expression datasets using two computational algorithms (CIBERSORT and ESTIMATE). ICI and transcriptome subsets were defined using a clustering algorithm, and survival analysis was also performed. Principal component analysis was used to calculate the novel ICI score, and gene set enrichment analysis was performed to identify the pathways underlying the defined clusters. The tumor mutational burden (TMB) was further explored in TCGA cohort, and survival analysis was used to assess the capability of the ICI and TMB scores to predict overall survival. Additionally, common driver gene mutations and their differential expression in the different ICI score group were investigated. Results: The ICI landscapes of 240 patients were generated using the devised algorithm, revealing three ICI and three gene clusters whose use improved the prediction of overall survival (p = 0.019 and p < 0.001, respectively). Crucial immune checkpoint genes were differentially expressed among these subtypes; the RIG-I-LIKE and NOD-LIKE receptor signaling pathways were enriched in samples with low ICI scores (p < 0.05). We also found that the TMB scores could predict survival outcomes, whereas the ICI scores also could predict prognoses independent of TMB. Notably, ICI scores could effectively predict responses to immunotherapy. KRAS, TP53, CDKN2A, SMAD4 and TTN remained the most commonly mutated genes in PC; moreover, KRAS and TP53 mutation rates were significantly different between the two ICI score groups. Conclusions: We developed a novel ICI score that could independently predict the response to immunotherapy and survival of patients with PC. Evaluation of the ICI landscape in a larger cohort could clarify the interactions between these infiltrating cells, the tumor microenvironment and response to immunotherapy.