Single-cell transcriptome reveals the heterogeneity of malignant ductal cells and the prognostic value of REG4 and SPINK1 in primary pancreatic ductal adenocarcinoma

Abstract

Background Pancreatic ductal adenocarcinoma (PDAC) is one of the leading causes of cancer-related deaths, with very limited therapeutic options available. This study aims to comprehensively depict the heterogeneity and identify prognostic targets for PDAC with single-cell RNA sequencing (scRNA-seq) analysis. Methods ScRNA-seq analysis was performed on 16 primary PDAC and three adjacent lesions. A series of analytical methods were applied for analysis in cell clustering, gene profiling, lineage trajectory analysis and cell-to-cell interactions. In vitro experiments including colony formation, wound healing and sphere formation assay were performed to assess the role of makers. Results A total of 32,480 cells were clustered into six major populations, among which the ductal cell cluster expressing high copy number variants (CNVs) was defined as malignant cells. Malignant cells were further subtyped into five subgroups which exhibited specific features in immunologic and metabolic activities. Pseudotime trajectory analysis indicated that components of various oncogenic pathways were differentially expressed along tumor progression. Furthermore, intensive substantial crosstalk between ductal cells and stromal cells was identified. Finally, genes (REG4 and SPINK1) screened out of differentially expressed genes (DEGs) were upregulated in PDAC cell lines. Silencing either of them significantly impaired proliferation, invasion, migration and stemness of PDAC cells. Conclusions Our findings offer a valuable resource for deciphering the heterogeneity of malignant ductal cells in PDAC. REG4 and SPINK1 are expected to be promising targets for PDAC therapy.