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.

Primary Sjögren’s syndrome (pSS) is a systemic autoimmune disease characterized primarily by immune-mediated destruction of exocrine tissues, such as those of the salivary and lacrimal glands, resulting in the loss of saliva and tear production, respectively. This disease predominantly affects middle-aged women, often in an insidious manner with the accumulation of subtle changes in glandular function occurring over many years. Patients commonly suffer from pSS symptoms for years before receiving a diagnosis. Currently, there is no effective cure for pSS and treatment options and targeted therapy approaches are limited due to a lack of our overall understanding of the disease etiology and its underlying pathology. To better elucidate the underlying molecular nature of this disease, we have performed RNA-sequencing to generate a comprehensive global gene expression profile of minor salivary glands from an ethnically diverse cohort of patients with pSS. Gene expression analysis has identified a number of pathways and networks that are relevant in pSS pathogenesis. Moreover, our detailed integrative analysis has revealed a primary Sjögren’s syndrome molecular signature that may represent important players acting as potential drivers of this disease. Finally, we have established that the global transcriptomic changes in pSS are likely to be attributed not only to various immune cell types within the salivary gland but also epithelial cells which are likely playing a contributing role. Overall, our comprehensive studies provide a database-enriched framework and resource for the identification and examination of key pathways, mediators, and new biomarkers important in the pathogenesis of this disease with the long-term goals of facilitating earlier diagnosis of pSS and to mitigate or abrogate the progression of this debilitating disease.

Sjögren’s Syndrome (SS) is a chronic autoimmune disease of unknown etiology which primarily affects the salivary and lacrimal glands resulting in the loss of secretory function. Treatment options for SS have been hampered due to the lack of a better understanding of the underlying gene regulatory circuitry and the interplay between the myriad pathological cellular states that contribute to salivary gland dysfunction. To better elucidate the molecular nature of SS, we have performed RNA-sequencing analysis of the submandibular glands (SMG) of a well-established primary Sjögren’s Syndrome (pSS) mouse model. Our comprehensive examination of global gene expression and comparative analyses with additional SS mouse models and human datasets, have identified a number of important pathways and regulatory networks that are relevant in SS pathobiology. To complement these studies, we have performed single-cell RNA sequencing to examine and identify the molecular and cellular heterogeneity of the diseased cell populations of the mouse SMG. Interrogation of the single-cell transcriptomes has shed light on the diversity of immune cells that are dysregulated in SS and importantly, revealed an activated state of the salivary gland epithelial cells that contribute to the global immune mediated responses. Overall, our broad studies have not only revealed key pathways, mediators and new biomarkers, but have also uncovered the complex nature of the cellular populations in the SMG that are likely to drive the progression of SS. These newly discovered insights into the underlying molecular mechanisms and cellular states of SS will better inform targeted therapeutic discoveries.

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.

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.