Branching organs, including the salivary and mammary glands, lung, and kidney, arise as epithelial buds that are morphologically very similar. However, the mesenchyme is known to guide epithelial morphogenesis and to help govern cell fate and eventual organ specificity. We performed single-cell transcriptome analyses of 14,441 cells from embryonic day 12 submandibular and parotid salivary glands to characterize their molecular identities during bud initiation. The mesenchymal cells were considerably more heterogeneous by clustering analysis than the epithelial cells. Nonetheless, distinct clusters were evident among even the epithelial cells, where unique molecular markers separated presumptive bud and duct cells. Mesenchymal cells formed separate, well-defined clusters specific to each gland. Neuronal and muscle cells of the 2 glands in particular showed different markers and localization patterns. Several gland-specific genes were characteristic of different rhombomeres. A muscle cluster was prominent in the parotid, which was not myoepithelial or vascular smooth muscle. Instead, the muscle cluster expressed genes that mediate skeletal muscle differentiation and function. Striated muscle was indeed found later in development surrounding the parotid gland. Distinct spatial localization patterns of neuronal and muscle cells in embryonic stages appear to foreshadow later differences in adult organ function. These findings demonstrate that the establishment of transcriptional identities emerges early in development, primarily in the mesenchyme of developing salivary glands. We present the first comprehensive description of molecular signatures that define specific cellular landmarks for the bud initiation stage, when the neural crest–derived ectomesenchyme predominates in the salivary mesenchyme that immediately surrounds the budding epithelium. We also provide the first transcriptome data for the largely understudied embryonic parotid gland as compared with the submandibular gland, focusing on the mesenchymal cell populations.
mTORC1 activation in lung mesenchyme drives sex- and age-dependent pulmonary structure and function decline
