Abstract
Pancreatic duct adenocarcinoma (PDAC) is an aggressive and lethal malignancy. Pancreatic cancer stem cells (PCSCs) are assumed to contribute to the initiation and invasion of PDAC. In this study, we performed single-cell RNA sequencing (scRNA-seq) analysis of PDAC tumor samples from patients and control pancreas tissues to reveal the transformation process of cancer stem cell (CSC)-like ductal cells into ductal cells with invasive potential and we screened out CSC-related genes (CRGs). Subsequently, we applied LASSO and Cox regression models to identify five CRGs with potential prognostic values and constructed a risk prognostic model using the Cancer Genome Atlas datasets. The risk models were verified using Gene Expression Omnibus datasets. Patients in the high-risk group had a significantly poor overall survival (Pvalue<0.0001), as illustrated by the Kaplan-Meier survival curve, and the area under the curve confirmed the accuracy of predictions by our risk model. Tumor mutation burden variations were used to further explore the differences between the two risk cohorts. In addition, the Human Protein Atlas was used to investigate the protein expression of five hub CRGs. In brief, we utilized scRNA-seq to reveal the invasive trajectory of ductal cells and identified crucial CRGs in PDAC, which may help predict patient survival and provide potential clinical therapeutic targets against CSCs.