High-Resolution Transcriptomic Landscape of the Human Submandibular Gland.

Saliva-secreting and transporting cells are part of the complex cellular milieu of the human salivary gland, where they play important roles in normal glandular physiology and diseased states. However, comprehensive molecular characterization, particularly at single-cell resolution, is still incomplete, in part due to difficulty in procuring normal human tissues. Here, we perform an in-depth analysis of male and female adult human submandibular gland (SMG) samples by bulk RNA sequencing (RNA-seq) and examine the molecular underpinnings of the heterogeneous cell populations by single-cell (sc) RNA-seq. Our results from scRNA-seq highlight the remarkable diversity of clusters of epithelial and nonepithelial cells that reside in the SMG that is also faithfully recapitulated by deconvolution of the bulk-RNA data sets. Our analyses reveal complex transcriptomic heterogeneity within both the ductal and acinar subpopulations and identify atypical SMG cell types, such as mucoacinar cells that are unique to humans and ionocytes that have been recently described in the mouse. We use CellChat to explore ligand-receptor interactome predictions that likely mediate crucial cell-cell communications between the various cell clusters. Finally, we apply a trajectory inference method to investigate specific cellular branching points and topology that offers insights into the dynamic and complex differentiation process of the adult SMG. The data sets and the analyses herein comprise an extensive wealth of high-resolution information and a valuable resource for a deeper mechanistic understanding of human SMG biology and pathophysiology.

Genetic Study of Elf5 and Ehf in the Mouse Salivary Gland.

Salivary gland (SG) development, maturation, and homeostasis require coordinated roles of transcription factors (TFs) that dictate specific cell identities and fate. The ETS family of proteins are important transcriptional drivers of diverse cell lineages, tissue development, and differentiation programs and hence are also likely to play an important role in the SG. Here we have leveraged genomic and epigenomic data of the SG to examine the expression profile of ETS genes and identified 2 closely related paralogs, Elf5 and Ehf, that are highly expressed in distinct epithelial subpopulations. By using a well-defined mouse knockout model of Elf5, we show that Elf5, despite its enriched expression in the acinar cells, is functionally dispensable for maintaining the homeostatic state of the adult SG epithelium. The lack of a discernible phenotype of the Elf5-null SG might be due to possible functional redundancy with Ehf or other ETS factors. To probe this possibility and to examine the specific consequences of Ehf loss in the SG, we used CRISPR-Cas9 to generate mice in which the DNA-binding ETS domain of Ehf is disrupted due to an insertion mutation. We demonstrate that the Ehf mutant (EhfMut) mice exhibit a distinct cellular phenotype with decreased granular convoluted tubules that are accompanied by an increased accumulation of the intercalated Sox9-positive ductal cell population. Interestingly, the ductal phenotype of the EhfMut animals is highly pronounced in males, reaffirming the established sexual dimorphism of the SG that exists in rodents. Our results show that unlike Elf5, Ehf plays a nonredundant role in directing ductal cell differentiation of the SG and highlights the phenotypic subtlety in mutant mice of closely related TFs and the importance of careful consideration of cell type-specific studies.