Single-cell transcriptome profiling of primary tumors and paired organoids of pancreatobiliary cancer

Single-cell RNA-seq (scRNA-seq) and cancer organoid model have shown promise in investigating tumor microenvironment heterogeneity and facilitating chemotherapeutic drug testing to inform treatment selection. It is still unknown whether the scRNA-seq results based on organoid can faithfully reflect the heterogeneity of primary pancreatobiliary cancer. To reveal the similarities and differences between primary tumors and their matched organoids at transcriptome level, we conducted scRNA-seq for paired primary tumors and organoids from one cholangiocarcinoma (CCA) and two pancreatic ductal adenocarcinoma (PDAC) patients. We identified inter-patient and intra-tumor heterogeneity and found that the organoids retained copy number variation (CNV) patterns of primary tumors. There was no significant difference in cancer stem cell (CSC) properties between the primary tumors and the organoids, whereas organoid from one PDAC case had increased mesenchymal-score and decreased epithelial-score compared with the primary tumors. All organoids showed a transition tendency from the classical subtype to the basal-like subtype in the transcriptional level. Organoids and primary tumors differed in metabolic and unfolded protein response (UPR) signatures. In addition, we revealed the heterogeneity of cancer associated fibroblasts (CAFs) and T cells, and explored the developmental trajectory of T cells. Our findings facilitate further understanding of organoid model and confirm its application prospects in pancreatobiliary cancer.