A single cell atlas of the murine pancreatic ductal tree identifies novel cell populations with potential implications in pancreas regeneration and exocrine pathogenesis.

Pancreatic ducts form an intricate network of tubules that secrete bicarbonate and drive acinar secretions into the duodenum. This network is formed by centroacinar cells, terminal, intercalated, intracalated ducts, and the main pancreatic duct. Ductal heterogeneity at the single-cell level has been poorly characterized. Here, we used scRNA-seq to comprehensively characterize mouse ductal heterogeneity at single-cell resolution of the entire ductal epithelium from centroacinar cells to the main duct. Moreover, we used organoid cultures, injury models and pancreatic tumor samples to interrogate the role of novel ductal populations in pancreas regeneration and exocrine pathogenesis. In our study, we have identified the coexistence of 15 ductal populations within the healthy pancreas and characterized their organoid formation capacity and endocrine differentiation potential. Cluster isolation and subsequent culturing let us identify ductal cell populations with high organoid formation capacity and endocrine and exocrine differentiation potential in vitro, including Wnt-responsive-population, ciliated-population and FLRT3+ cells. Moreover, we have characterized the location of these novel ductal populations in healthy pancreas, chronic pancreatitis and tumor samples, hightlihgting a putative role of WNT-responsive, IFN-responsive and EMT-populations in pancreatic exocrine pathogenesis as their expression inceases in chronic pancreatitis and PanIN lesions. In light of our discovery of previously unidentified ductal populations, we unmask the potential roles of specific ductal populations in pancreas regeneration and exocrine pathogenesis. Overall design: Tg(Sox9-EGFP) EB209 Gsat obtained from J.Ferrer laboratory (CRG) were utilized in this study. Sox9:eGFP positive cells were sorted from small, medoim and large ducts via FACS. Ten thousand sorted cells from each ductal fraction were pooled from 3 male mice in two independent experiments and used for scRNA-seq.