RNA-seq based transcriptomic map reveals new insights into mouse salivary gland development and maturation.

BACKGROUND: Mouse 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. RESULTS: To 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. CONCLUSIONS: Our 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.

A global analysis of the complex landscape of isoforms and regulatory networks of p63 in human cells and tissues.

BACKGROUND: The transcription factor p63 belongs to the p53/p63/p73 family and plays key functional roles during normal epithelial development and differentiation and in pathological states such as squamous cell carcinomas. The human TP63 gene, located on chromosome 3q28 is driven by two promoters that generate the full-length transactivating (TA) and N-terminal truncated (ΔN) isoforms. Furthermore alternative splicing at the C-terminus gives rise to additional α, ß, γ and likely several other minor variants. Teasing out the expression and biological function of each p63 variant has been both the focus of, and a cause for contention in the p63 field. RESULTS: Here we have taken advantage of a burgeoning RNA-Seq based genomic data-sets to examine the global expression profiles of p63 isoforms across commonly utilized human cell-lines and major tissues and organs. Consistent with earlier studies, we find ΔNp63 transcripts, primarily that of the ΔNp63α isoforms, to be expressed in most cells of epithelial origin such as those of skin and oral tissues, mammary glands and squamous cell carcinomas. In contrast, TAp63 is not expressed in the majority of normal cell-types and tissues; rather it is selectively expressed at moderate to high levels in a subset of Burkitt's and diffuse large B-cell lymphoma cell lines. We verify this differential expression pattern of p63 isoforms by Western blot analysis, using newly developed ΔN and TA specific antibodies. Furthermore using unsupervised clustering of human cell lines, tissues and organs, we show that ΔNp63 and TAp63 driven transcriptional networks involve very distinct sets of molecular players, which may underlie their different biological functions. CONCLUSIONS: In this study we report comprehensive and global expression profiles of p63 isoforms and their relationship to p53/p73 and other potential transcriptional co-regulators. We curate publicly available data generated in part by consortiums such as ENCODE, FANTOM and Human Protein Atlas to delineate the vastly different transcriptomic landscapes of ΔNp63 and TAp63. Our studies help not only in dispelling prevailing myths and controversies on p63 expression in commonly used human cell lines but also augur new isoform- and cell type-specific activities of p63.

ΔNp63 knockout mice reveal its indispensable role as a master regulator of epithelial development and differentiation.

The transcription factor p63 is important in the development of the skin as p63-null mice exhibit striking defects in embryonic epidermal morphogenesis. Understanding the mechanisms that underlie this phenotype is complicated by the existence of multiple p63 isoforms, including TAp63 and ΔNp63. To investigate the role of ΔNp63 in epidermal morphogenesis we generated ΔNp63 knock-in mice in which the ΔNp63-specific exon is replaced by GFP. Homozygous ΔNp63(gfp/gfp) animals exhibit severe developmental anomalies including truncated forelimbs and the absence of hind limbs, largely phenocopying existing knockouts in which all p63 isoforms are deleted. ΔNp63-null animals show a poorly developed stratified epidermis comprising isolated clusters of disorganized epithelial cells. Despite the failure to develop a mature stratified epidermis, the patches of ΔNp63-null keratinocytes are able to stratify and undergo a program of terminal differentiation. However, we observe premature expression of markers associated with terminal differentiation, which is unique to ΔNp63-null animals and not evident in the skin of mice lacking all p63 isoforms. We posit that the dysregulated and accelerated keratinocyte differentiation phenotype is driven by significant alterations in the expression of key components of the Notch signaling pathway, some of which are direct transcriptional targets of ΔNp63 as demonstrated by ChIP experiments. The analysis of ΔNp63(gfp/gfp) knockout mice reaffirms the indispensable role of the ΔN isoform of p63 in epithelial biology and confirms that ΔNp63-null keratinocytes are capable of committing to an epidermal cell lineage, but are likely to suffer from diminished renewal capacity and an altered differentiation fate.

Abnormal hair follicle development and altered cell fate of follicular keratinocytes in transgenic mice expressing DeltaNp63alpha.

The transcription factor p63 plays an essential role in epidermal morphogenesis. Animals lacking p63 fail to form many ectodermal organs, including the skin and hair follicles. Although the indispensable role of p63 in stratified epithelial skin development is well established, relatively little is known about this transcriptional regulator in directing hair follicle morphogenesis. Here, using specific antibodies, we have established the expression pattern of DeltaNp63 in hair follicle development and cycling. DeltaNp63 is expressed in the developing hair placode, whereas in mature hair its expression is restricted to the outer root sheath (ORS), matrix cells and to the stem cells of the hair follicle bulge. To investigate the role of DeltaNp63 in hair follicle morphogenesis and cycling, we have utilized a Tet-inducible mouse model system with targeted expression of this isoform to the ORS of the hair follicle. DeltaNp63 transgenic animals display dramatic defects in hair follicle development and cycling, eventually leading to severe hair loss. Strikingly, expression of DeltaNp63 leads to a switch in cell fate of hair follicle keratinocytes, causing them to adopt an interfollicular epidermal (IFE) cell identity. Moreover, DeltaNp63 transgenic animals exhibit a depleted hair follicle stem-cell niche, which further contributes to the overall cycling defects observed in the mutant animals. Finally, global transcriptome analysis of transgenic skin identified altered expression levels of crucial mediators of hair morphogenesis, including key members of the Wnt/beta-catenin signaling pathway, which, in part, account for these effects. Our data provide evidence supporting a role for DeltaNp63alpha in actively suppressing hair follicle differentiation and directing IFE cell lineage commitment.

ΔNp63 maintains the fidelity of the myoepithelial cell lineage and directs cell differentiation programs in the murine salivary gland.

Salivary glands consist of several epithelial cell types of distinct lineages and functional characteristics that are established by directed differentiation programs of resident stem and progenitor cells. We have shown that ΔNp63, a crucial transcriptional regulator of stem/progenitor cells, is enriched in both the basal and myoepithelial cell (MEC) populations and that ΔNp63 positive cells maintain all the descendent epithelial cell lineages of the adult mouse salivary glands (mSGs). Although this pivotal role of ΔNp63 in driving the broader epithelial cell fate and identity in the mSG has been demonstrated, how ΔNp63 functions specifically in the commitment and differentiation of the MEC population is less understood. Using multiple genetic mouse models that allow for cell tracing, we show that ΔNp63 is critical in maintaining and renewing MECs, in part through the transcriptional regulation of Acta2 gene expression, a defining marker of this cell population. We demonstrate that during adult mSG homeostasis, ΔNp63 enriched MECs function as bipotent progenitor cells that maintain not only the MEC population, but also the distinctly different ductal cell lineages. The fidelity of this process is dependent on ΔNp63 expression, since MEC-specific ablation of ΔNp63 results in altered MEC differentiation and affects cellular plasticity resulting in aberrant differentiation of the intercalated ducts and acinar cells. In contrast, we find that the contribution of MECs to ductal and acinar cell regeneration following severe injury is independent of ΔNp63. Our observations offer new insights into cellular mechanisms driving MEC fate choices and differentiation programs in the context of salivary gland homeostasis and in response to injury and regeneration. Long term, these findings have implications for better treatment of salivary gland dysfunction through stem cell-based approaches.

Genome-wide search identifies Ccnd2 as a direct transcriptional target of Elf5 in mouse mammary gland.

BACKGROUND: The ETS transcription factor Elf5 (also known as ESE-2) is highly expressed in the mammary gland and plays an important role in its development and differentiation. Indeed studies in mice have illustrated an essential role for Elf5 in directing alveologenesis during pregnancy. Although the molecular mechanisms that underlie the developmental block in Elf5 null mammary glands are beginning to be unraveled, this investigation has been hampered by limited information about the identity of Elf5-target genes. To address this shortcoming, in this study we have performed ChIP-cloning experiments to identify the specific genomic segments that are occupied by Elf5 in pregnant mouse mammary glands. RESULTS: Sequencing and genomic localization of cis-regulatory regions bound by Elf5 in vivo has identified several potential target genes covering broad functional categories. A subset of these target genes demonstrates higher expression levels in Elf5-null mammary glands suggesting a repressive functional role for this transcription factor. Here we focus on one putative target of Elf5, the Ccnd2 gene that appeared in our screen. We identify a novel Elf5-binding segment upstream of the Ccnd2 gene and demonstrate that Elf5 can transcriptionally repress Ccnd2 by directly binding to the proximal promoter region. Finally, using Elf5-null mammary epithelial cells and mammary glands, we show that loss of Elf5 in vivo leads to up regulation of Ccnd2 and an altered expression pattern in luminal cells. CONCLUSIONS: Identification of Elf5-targets is an essential first step in elucidating the transcriptional landscape that is shaped by this important regulator. Our studies offer new toolbox in examining the biological role of Elf5 in mammary gland development and differentiation.

p63 and Its Target Follistatin Maintain Salivary Gland Stem/Progenitor Cell Function through TGF-ß/Activin Signaling.

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.

Functional characterization and genomic studies of a novel murine submandibular gland epithelial cell line.

A better understanding of the normal and diseased biology of salivary glands (SG) has been hampered, in part, due to difficulties in cultivating and maintaining salivary epithelial cells. Towards this end, we have generated a mouse salivary gland epithelial cell (mSGc) culture system that is well-suited for the molecular characterization of SG cells and their differentiation program. We demonstrate that mSGc can be maintained for multiple passages without a loss of proliferation potential, readily form 3D-spheroids and importantly express a panel of well-established salivary gland epithelial cell markers. Moreover, mSGc 3D-spheroids also exhibit functional maturation as evident by robust agonist-induced intracellular calcium signaling. Finally, transcriptomic characterization of mSGc by RNA-seq and hierarchical clustering analysis with adult organ RNA-seq datasets reveal that mSGc retain most of the molecular attributes of adult mouse salivary gland. This well-characterized mouse salivary gland cell line will fill a critical void in the field by offering a valuable resource to examine various mechanistic aspects of mouse salivary gland biology.

Genetic and scRNA-seq Analysis Reveals Distinct Cell Populations that Contribute to Salivary Gland Development and Maintenance.

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.

Elf5 inhibits the epithelial-mesenchymal transition in mammary gland development and breast cancer metastasis by transcriptionally repressing Snail2.

The epithelial-mesenchymal transition (EMT) is a complex process that occurs during organogenesis and in cancer metastasis. Despite recent progress, the molecular pathways connecting the physiological and pathological functions of EMT need to be better defined. Here we show that the transcription factor Elf5, a key regulator of mammary gland alveologenesis, controls EMT in both mammary gland development and metastasis. We uncovered this role for Elf5 through analyses of Elf5 conditional knockout animals, various in vitro and in vivo models of EMT and metastasis, an MMTV-neu transgenic model of mammary tumour progression and clinical breast cancer samples. Furthermore, we demonstrate that Elf5 suppresses EMT by directly repressing the transcription of Snail2, a master regulator of mammary stem cells and a known inducer of EMT. These findings establish Elf5 not only as a key cell lineage regulator during normal mammary gland development, but also as a suppressor of EMT and metastasis in breast cancer.

An active role of the DeltaN isoform of p63 in regulating basal keratin genes K5 and K14 and directing epidermal cell fate.

BACKGROUND: One major defining characteristic of the basal keratinocytes of the stratified epithelium is the expression of the keratin genes K5 and K14. The temporal and spatial expression of these two genes is usually tightly and coordinately regulated at the transcriptional level. This ensures the obligate pairing of K5 and K14 proteins to generate an intermediate filament (IF) network that is essential for the structure and function of the proliferative keratinocytes. Our previous studies have shown that the basal-keratinocyte restricted transcription factor p63 is a direct regulator of K14 gene. METHODOLOGY/PRINCIPAL FINDINGS: Here we provide evidence that p63, specifically the DeltaN isoform also regulates the expression of the K5 gene by binding to a conserved enhancer within the 5' upstream region. By using specific antibodies against DeltaNp63, we show a concordance in the expression between basal keratins and DeltaNp63 proteins but not the TAp63 isoforms during early embryonic skin development. We demonstrate, that contrary to a previous report, transgenic mice expressing DeltaNp63 in lung epithelium exhibit squamous metaplasia with de novo induction of K5 and K14 as well as transdifferentiation to the epidermal cell lineage. Interestingly, the in vivo epidermal inductive properties of DeltaNp63 do not require the C-terminal SAM domain. Finally, we show that DeltaNp63 alone can restore the expression of the basal keratins and reinitiate the failed epidermal differentiation program in the skin of p63 null animals. SIGNIFICANCE: DeltaNp63 is a critical mediator of keratinocyte stratification program and directly regulates the basal keratin genes.