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.

A Global Vista of the Epigenomic State of the Mouse Submandibular Gland.

The parotid, submandibular, and sublingual glands represent a trio of oral secretory glands whose primary function is to produce saliva, facilitate digestion of food, provide protection against microbes, and maintain oral health. While recent studies have begun to shed light on the global gene expression patterns and profiles of salivary glands, particularly those of mice, relatively little is known about the location and identity of transcriptional control elements. Here we have established the epigenomic landscape of the mouse submandibular salivary gland (SMG) by performing chromatin immunoprecipitation sequencing experiments for 4 key histone marks. Our analysis of the comprehensive SMG data sets and comparisons with those from other adult organs have identified critical enhancers and super-enhancers of the mouse SMG. By further integrating these findings with complementary RNA-sequencing based gene expression data, we have unearthed a number of molecular regulators such as members of the Fox family of transcription factors that are enriched and likely to be functionally relevant for SMG biology. Overall, our studies provide a powerful atlas of cis-regulatory elements that can be leveraged for better understanding the transcriptional control mechanisms of the mouse SMG, discovery of novel genetic switches, and modulating tissue-specific gene expression in a targeted fashion.

Transcriptomic and Single-Cell Analysis of the Murine Parotid Gland.

The salivary complex of mammals consists of 3 major pairs of glands: the parotid, submandibular, and sublingual glands. While the 3 glands share similar functional properties, such as saliva secretion, their differences are largely based on the types of secretions they produce. While recent studies have begun to shed light on the underlying molecular differences among the glands, few have examined the global transcriptional repertoire over various stages of gland maturation. To better elucidate the molecular nature of the parotid gland, we have performed RNA sequencing to generate comprehensive and global gene expression profiles of this gland at different stages of maturation. Our transcriptomic characterization and hierarchical clustering analysis with adult organ RNA sequencing data sets has identified a number of molecular players and pathways that are relevant for parotid gland biology. Moreover, our detailed analysis has revealed a unique parotid gland-specific gene signature that may represent important players that could impart parotid gland-specific biological properties. To complement our transcriptomic studies, we have performed single-cell RNA sequencing to map the transcriptomes of parotid epithelial cells. Interrogation of the single-cell transcriptomes revealed the degree of molecular and cellular heterogeneity of the various epithelial cell types within the parotid gland. Moreover, we uncovered a mixed-lineage population of cells that may reflect molecular priming of differentiation potentials. Overall our comprehensive studies provide a powerful tool for the discovery of novel molecular players important in parotid gland biology.