Summary
Understanding the dynamic transcriptional landscape throughout organ development will provide a template for regenerative therapies. Here, we generated a single-cell RNA sequencing atlas of murine submandibular glands identifying transcriptional profiles that revealed cellular heterogeneity during landmark developmental events: end bud formation, branching morphogenesis, cytodifferentiation, maturation, and homeostasis. Trajectory inference analysis suggests plasticity among acinar and duct populations. We identify transcription factors correlated with acinar differentiation including
Spdef, Etv1
, and
Xbp1
, and loss of
Ybx1
,
Eno1, Sox11,
and
Atf4
. Furthermore, we characterize two intercalated duct populations defined by either
Gfra3
and
Kit
, or
Gstt1
. This atlas can be used to investigate specific cell functions and comparative studies predicting common mechanisms involved in development of branching organs.