Numerous studies have characterized the existence of cell subtypes, along with their corresponding transcriptional profiles, within the developing mouse pancreas. The upstream mechanisms that initiate and maintain gene expression programs across cell states, however, remain largely unknown. Here, we generate single-nucleus ATAC-Sequencing data of developing murine pancreas and perform an integrated, multi-omic analysis of both chromatin accessibility and RNA expression to describe the chromatin landscape of the developing pancreas at both E14.5 and E17.5 at single-cell resolution. We identify candidate transcription factors regulating cell fate and construct gene regulatory networks of active transcription factor binding to regulatory regions of downstream target genes. This work serves as a valuable resource for the field of pancreatic biology in general and contributes to our understanding of lineage plasticity among endocrine cell types. In addition, these data identify which epigenetic states should be represented in the differentiation of stem cells to the pancreatic beta cell fate to best recapitulate in vitro the gene regulatory networks that are critical for progression along the beta cell lineage in vivo.

Pancreatic development requires spatially and temporally controlled expression of growth factors derived from mesenchyme. Here, we report that in mice the secreted factor Fgf9 is expressed principally by mesenchyme and then mesothelium during early development, then subsequently by both mesothelium and rare epithelial cells by E12.5 and onwards. Global knockout of the Fgf9 gene resulted in the reduction of pancreas and stomach size, as well as complete asplenia. The number of early Pdx1+ pancreatic progenitors was reduced at E10.5, as was proliferation of mesenchyme at E11.5. Although loss of Fgf9 did not interfere with differentiation of later epithelial lineages, single-cell RNA-Sequencing identified transcriptional programs perturbed upon loss of Fgf9 during pancreatic development, including loss of the transcription factor Barx1. Lastly, we identified conserved expression patterns of FGF9 and receptors in human fetal pancreas, suggesting that FGF9 expressed by pancreatic mesenchyme may similarly affect the development of the human pancreas.