Identification of Candidate Therapeutic Target Genes and Profiling of Tumor-Infiltrating Immune Cells in Pancreatic Cancer via Integrated Transcriptomic Analysis

Abstract

Pancreatic cancer (PC) has a dismal prognosis despite advancing scientific and technological knowledge. The exploration of novel genes is critical to improving current therapeutic measures. This research is aimed at selecting hub genes that can act as candidate therapeutic target genes and as prognostic biomarkers in PC. Gene expression profiles of datasets GSE101448, GSE15471, and GSE62452 were extracted from the GEO database. The “limma” package was performed to select differentially expressed genes (DEGs) between PC and normal tissue samples in each dataset. Robust rank aggregation (RRA) algorithm was conducted to integrate multiple expression profiles and identify robust DEGs. GO analysis and KEGG analysis were conducted to identify the functional correlation of the DEGs. The CIBERSORT algorithm was conducted to estimate the immune cell composition of each tissue sample. STRING and Cytoscape were used to establish the protein-protein interaction (PPI) network. The cytoHubba plugin in Cytoscape was performed to identify hub genes. Survival analysis based on hub gene expression was performed with clinical information from TCGA database. 566 robust DEGs (338 upregulated genes and 226 downregulated genes) were identified. Tumor tissue had a higher infiltration of resting dendritic cells and tumor-associated macrophages (TAM), including M0, M1, and M2 macrophages, while infiltration levels of B memory cells, plasma cells, T cells CD8, T follicular helper cells, and NK cells in normal tissue were relatively higher. GO terms and KEGG pathway analysis results revealed enrichment in tumor-associated pathways, including the extracellular matrix organization, cell−substrate adhesion cytokine−cytokine receptor interaction, calcium signaling pathway, and glycine, serine, and threonine metabolism, to name a few. Finally, FN1, MSLN, PLAU, and VCAN were selected as hub genes. High expression of FN1, MSLN, PLAU, and VCAN in PC significantly correlated with poor prognosis. Integrated transcriptomic analysis was used to provide new insights into PC pathogenesis. FN1, MSLN, PLAU, and VCAN may be considered as novel biomarkers of PC.