Stem and progenitor cells of the submandibular salivary gland (SMG) give rise to, maintain, and regenerate the multiple lineages of mature epithelial cells including those belonging to the ductal, acinar, basal and myoepithelial subtypes. Here we have exploited single cell RNA-sequencing and in vivo genetic lineage tracing technologies to generate a detailed map of the cell fate trajectories and branch points of the basal and myoepithelial cell populations of the mouse SMG during embryonic development and in adults. Our studies show that the transcription factor p63 and alpha-smooth muscle actin (SMA) serve as faithful markers of the basal and myoepithelial cell lineages, respectively and that both cell types are endowed with progenitor cell properties. However, p63+ basal and SMA+ myoepithelial cells exhibit distinct cell fates by virtue of maintaining different cellular lineages during morphogenesis and in adults. Collectively, our results reveal the dynamic and complex nature of the diverse SMG cell populations and highlight the distinct differentiation potential of the p63 and SMA expressing subtypes in the stem and progenitor cell hierarchy. Long term these findings have profound implications towards a better understanding of the molecular mechanisms that dictate lineage commitment and differentiation programs during development and adult gland maintenance.

BackgroundMouse models have served a valuable role in deciphering various facets of Salivary Gland (SG) biology, from normal developmental programs to diseased states. To facilitate such studies, gene expression profiling maps have been generated for various stages of SG organogenesis. However these prior studies fall short of capturing the transcriptional complexity due to the limited scope of gene-centric microarray-based technology. Compared to microarray, RNA-sequencing (RNA-seq) offers unbiased detection of novel transcripts, broader dynamic range and high specificity and sensitivity for detection of genes, transcripts, and differential gene expression. Although RNA-seq data, particularly under the auspices of the ENCODE project, have covered a large number of biological specimens, studies on the SG have been lacking.ResultsTo better appreciate the wide spectrum of gene expression profiles, we isolated RNA from mouse submandibular salivary glands at different embryonic and adult stages. In parallel, we processed RNA-seq data for 24 organs and tissues obtained from the mouse ENCODE consortium and calculated the average gene expression values. To identify molecular players and pathways likely to be relevant for SG biology, we performed functional gene enrichment analysis, network construction and hierarchal clustering of the RNA-seq datasets obtained from different stages of SG development and maturation, and other mouse organs and tissues. Our bioinformatics-based data analysis not only reaffirmed known modulators of SG morphogenesis but revealed novel transcription factors and signaling pathways unique to mouse SG biology and function. Finally we demonstrated that the unique SG gene signature obtained from our mouse studies is also well conserved and can demarcate features of the human SG transcriptome that is different from other tissues.ConclusionsOur RNA-seq based Atlas has revealed a high-resolution cartographic view of the dynamic transcriptomic landscape of the mouse SG at various stages. These RNA-seq datasets will complement pre-existing microarray based datasets, including the Salivary Gland Molecular Anatomy Project by offering a broader systems-biology based perspective rather than the classical gene-centric view. Ultimately such resources will be valuable in providing a useful toolkit to better understand how the diverse cell population of the SG are organized and controlled during development and differentiation.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3228-7) contains supplementary material, which is available to authorized users.

Atopic dermatitis (AD) is the most common inflammatory skin disease with no well-delineated cause or effective cure. Here we show that the p53 family member p63, specifically the ΔNp63, isoform has a key role in driving keratinocyte activation in AD. We find that overexpression of ΔNp63 in transgenic mouse epidermis results in a severe skin phenotype that shares many of the key clinical, histological and molecular features associated with human AD. This includes pruritus, epidermal hyperplasia, aberrant keratinocyte differentiation, enhanced expression of selected cytokines and chemokines and the infiltration of large numbers of inflammatory cells including type 2  T-helper cells – features that are highly representative of AD dermatopathology. We further demonstrate several of these mediators to be direct transcriptional targets of ΔNp63 in keratinocytes. Of particular significance are two p63 target genes, IL-31 and IL-33, both of which are key players in the signaling pathways implicated in AD. Importantly, we find these observations to be in good agreement with elevated levels of ΔNp63 in skin lesions of human patients with AD. Our studies reveal an important role for ΔNp63 in the pathogenesis of AD and offer new insights into its etiology and possible therapeutic targets.

Summary Multipotent ΔNp63-positive cells maintain all epithelial cell lineages of the embryonic and adult salivary gland (SG). However, the molecular mechanisms by which ΔNp63 regulates stem/progenitor (SP) cell populations in the SG remains elusive. To understand the role of ΔNp63 in directing cell fate choices in this gland, we have generated ΔNp63-deleted adult mice and primary salivary cell cultures to probe alterations in SP cell differentiation and function. In parallel, we have leveraged RNA-seq and ChIP-seq-based characterization of the ΔNp63-driven cistrome and scRNA-seq analysis to molecularly interrogate altered SG cellular identities and differentiation states dependent on ΔNp63. Our studies reveal that ablation of ΔNp63 results in a loss of the SP cell population and skewed differentiation that is mediated by Follistatin-dependent dysregulated TGF-β/Activin signaling. These findings offer new revelations into the SP cell gene regulatory networks that are likely to be relevant for normal or diseased SG states.

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.