Quantitative single cell analysis of cell population dynamics during submandibular salivary gland development and differentiation

Decoding mTOR signalling heterogeneity in the tumour microenvironment using multiplexed imaging and graph convolutional networks

Evaluating the contribution of the tumour microenvironment (TME) in tumour progression has proven a complex challenge due to the intricate interactions within the TME. Multiplexed imaging is an emerging technology that allows concurrent assessment of multiple of these components simultaneously. Here we utilise a highly multiplexed dataset of 61 markers across 746 colorectal tumours to investigate how complex mTOR signalling in different tissue compartments influences patient prognosis. We found that the signalling of mTOR pathway can have heterogeneous activation patterns in tumour and immune compartments which correlate with patient prognosis. Using graph neural networks, we determined the most predictive features of mTOR activity in immune cells and identified relevant cellular subpopulations. We validated our observations using spatial transcriptomics data analysis in an independent patient cohort. Our work provides a framework for studying complex cell signalling and reveals important insights for developing mTOR-based therapies.

Submandibular gland epithelial development and the importance of junctions

Salivary glands consist of highly specialized epithelial cells that secrete the fluid, saliva, and/or transport saliva into the oral cavity. Saliva is essential to lubricate the oral cavity for food consumption and to maintain the hygiene of the oral cavity. In this review, we will focus on the formation of the epithelial cell lineage and the cell junctions that are essential for formation of saliva and maintenance of the epithelial barrier between the ducts that transport saliva and the extracellular environment.

Identifying fibrogenic cells following salivary gland obstructive injury

Fibrosis results from excess extracellular matrix accumulation, which alters normal tissue architecture and impedes function. In the salivary gland, fibrosis can be induced by irradiation treatment for cancer therapy, Sjögren's Disease, and other causes; however, it is unclear which stromal cells and signals participate in injury responses and disease progression. As hedgehog signaling has been implicated in fibrosis of the salivary gland and other organs, we examined contributions of the hedgehog effector, Gli1, to fibrotic responses in salivary glands. To experimentally induce a fibrotic response in female murine submandibular salivary glands, we performed ductal ligation surgery. We detected a progressive fibrotic response where both extracellular matrix accumulation and actively remodeled collagen significantly increased at 14 days post-ligation. Macrophages, which participate in extracellular matrix remodeling, and Gli1+ and PDGFRα+ stromal cells, which may deposit extracellular matrix, both increased with injury. Using single-cell RNA-sequencing, Gli1 + cells were not found in discrete clusters at embryonic day 16 but were found in clusters expressing the stromal genes Pdgfra and/or Pdgfrb. In adult mice, Gli1+ cells were similarly heterogeneous but more cells co-expressed PDGFRα and PDGFRβ. Using Gli1-CreERT2; ROSA26tdTomato lineage-tracing mice, we found that Gli1-derived cells expand with ductal ligation injury. Although some of the Gli1 lineage-traced tdTomato+ cells expressed vimentin and PDGFRβ following injury, there was no increase in the classic myofibroblast marker, smooth muscle alpha-actin. Additionally, there was little change in extracellular matrix area, remodeled collagen area, PDGFRα, PDGFRβ, endothelial cells, neurons, or macrophages in Gli1 null salivary glands following injury when compared with controls, suggesting that Gli1 signaling and Gli1+ cells have only a minor contribution to mechanical injury-induced fibrotic changes in the salivary gland. We used scRNA-seq to examine cell populations that expand with ligation and/or showed increased expression of matrisome genes. Some Pdgfra + /Pdgfrb + stromal cell subpopulations expanded in response to ligation, with two stromal cell subpopulations showing increased expression of Col1a1 and a greater diversity of matrisome genes, consistent with these cells being fibrogenic. However, only a few cells in these subpopulations expressed Gli1, consistent with a minor contribution of these cells to extracellular matrix production. Defining the signaling pathways driving fibrotic responses in stromal cell sub-types could reveal future therapeutic targets.

Autologous mesenchymal stem cells offer a new paradigm for salivary gland regeneration

Salivary gland (SG) dysfunction, due to radiotherapy, disease, or aging, is a clinical manifestation that has the potential to cause severe oral and/or systemic diseases and compromise quality of life. Currently, the standard-of-care for this condition remains palliative. A variety of approaches have been employed to restore saliva production, but they have largely failed due to damage to both secretory cells and the extracellular matrix (niche). Transplantation of allogeneic cells from healthy donors has been suggested as a potential solution, but no definitive population of SG stem cells, capable of regenerating the gland, has been identified. Alternatively, mesenchymal stem cells (MSCs) are abundant, well characterized, and during SG development/homeostasis engage in signaling crosstalk with the SG epithelium. Further, the trans-differentiation potential of these cells and their ability to regenerate SG tissues have been demonstrated. However, recent findings suggest that the "immuno-privileged" status of allogeneic adult MSCs may not reflect their status post-transplantation. In contrast, autologous MSCs can be recovered from healthy tissues and do not present a challenge to the recipient's immune system. With recent advances in our ability to expand MSCs in vitro on tissue-specific matrices, autologous MSCs may offer a new therapeutic paradigm for restoration of SG function.

Neurotrophin signaling is a central mechanism of salivary dysfunction after irradiation that disrupts myoepithelial cells

The mechanisms that prevent regeneration of irradiated (IR) salivary glands remain elusive. Bulk RNAseq of IR versus non-IR human salivary glands showed that neurotrophin signaling is highly disrupted post-radiation. Neurotrophin receptors (NTRs) were significantly upregulated in myoepithelial cells (MECs) post-IR, and single cell RNAseq revealed that MECs pericytes, and duct cells are the main sources of neurotrophin ligands. Using two ex vivo models, we show that nerve growth factor (NGF) induces expression of MEC genes during development, and upregulation of NTRs in adult MECs is associated with stress-induced plasticity and morphological abnormalities in IR human glands. As MECs are epithelial progenitors after gland damage and are required for proper acinar cell contraction and secretion, we propose that MEC-specific upregulation of NTRs post-IR disrupts MEC differentiation and potentially impedes the ability of the gland to regenerate.

Identifying Fibrogenic Cells Following Salivary Gland Obstructive Injury

Fibrosis results from excess extracellular matrix accumulation, which alters normal tissue architecture and impedes function. In the salivary gland, fibrosis can be induced by irradiation treatment for cancer therapy, Sjögren's Disease, and other causes; however, it is unclear which stromal cells and signals participate in injury responses and disease progression. As hedgehog signaling has been implicated in fibrosis of the salivary gland and other organs, we examined contributions of the hedgehog effector, Gli1, to fibrotic responses in salivary glands. To experimentally induce a fibrotic response in female murine submandibular salivary glands, we performed ductal ligation surgery. We detected a progressive fibrotic response where both extracellular matrix accumulation and actively remodeled collagen trended upwards at 7 days and significantly increased at 14 days post- ligation. Macrophages, which participate in extracellular matrix remodeling, Gli1 + and PDGFRα + stromal cells, which may deposit extracellular matrix, both increased with injury. Using single-cell RNA-sequencing, we found that a majority of Gli1 + cells at embryonic day 16 also express Pdgfra and/or Pdgfrb. However, in adult mice, only a small subset of Gli1 + cells express PDGFRα and/or PDGFRβ at the protein level. Using lineage-tracing mice, we found that Gli1-derived cells expand with ductal ligation injury. Although some of the Gli1 lineage-traced tdTomato + cells expressed vimentin and PDGFRβ following injury, there was no increase in the classic myofibroblast marker, smooth muscle alpha-actin. Additionally, there was little change in extracellular matrix area, remodeled collagen area, PDGFRα, PDGFRβ, endothelial cells, neurons, or macrophages in Gli1 null salivary glands following injury when compared with controls, suggesting that Gli1 signaling and Gli1 + cells have only a minor contribution to mechanical injury-induced fibrotic changes in the salivary gland. We used scRNA-seq to examine cell populations that expand with ligation and/or showed increased expression of matrisome genes. Pdgfra + /Pdgfrb + stromal cell subpopulations both expanded in response to ligation, showed increased expression and a greater diversity of matrisome genes expressed, consistent with these cells being fibrogenic. Defining the signaling pathways driving fibrotic responses in stromal cell sub-types could reveal future therapeutic targets.

A cellular hierarchy of Notch and Kras signaling controls cell fate specification in the developing mouse salivary gland

The development of the mouse salivary gland involves a tip-driven process of branching morphogenesis that takes place in concert with differentiation into acinar, myoepithelial, and ductal (basal and luminal) sub-lineages. By combining clonal lineage tracing with a three-dimensional (3D) reconstruction of the branched epithelial network and single-cell RNA-seq analysis, we show that in tips, a heterogeneous population of renewing progenitors transition from a Krt14+ multipotent state to unipotent states via two transcriptionally distinct bipotent states, one restricted to the Krt14+ basal and myoepithelial lineage and the other to the Krt8+ acinar and luminal lineage. Using genetic perturbations, we show how the differential expression of Notch signaling correlates with spatial segregation, exits from multipotency, and promotes the Krt8+ lineage, whereas Kras activation promotes proacinar fate. These findings provide a mechanistic basis for how positional cues within growing tips regulate the process of lineage segregation and ductal patterning.

Bulk RNA-seq and scRNA-seq analysis reveal an activation of immune response and compromise of secretory function in major salivary glands of obese mice

Obesity affects the function of multiple organs/tissues including the exocrine organ salivary glands. However, the effects of obesity on transcriptomes and cell compositions in the salivary glands have yet been studied by bulk RNA-sequencing and single-cell RNA-sequencing. Besides, the cell types in the sublingual gland, one of the three major salivary glands, have yet been characterized by the approach of single-cell RNA-sequencing. In this report, we find that the histological structure of the three major salivary glands are not obviously affected in the obese mice. Bulk RNA-sequencing analysis shows that the most prominent changes observed in the three major salivary glands of the obese mice are the mobilization of transcriptomes related to the immune response and down-regulation of genes related to the secretory function of the salivary glands. Based on single-cell RNA-sequencing analysis, we identify and annotate 17 cell clusters in the sublingual gland for the first time, and find that obesity alters the relative compositions of immune cells and secretory cells in the major glands of obese mice. Integrative analysis of the bulk RNA-sequencing and single-cell RNA-sequencing data confirms the activation of immune response genes and compromise of secretory function in the three major salivary glands of obese mice. Consequently, the secretion of extracellular matrix proteins is significantly reduced in the three major salivary glands of obese mice. These results provide new molecular insights into understanding the effect of obesity on salivary glands.

The Semaphorin 3A-AKT axis-mediated cell proliferation in salivary gland morphogenesis and adenoid cystic carcinoma pathogenesis

We recently demonstrated that Semaphorin 3 A (Sema3A), the expression of which is negatively regulated by Wnt/β-catenin signaling, promotes odontogenic epithelial cell proliferation, suggesting the involvement of Sema3A in tooth germ development. Salivary glands have a similar developmental process to tooth germ development, in which reciprocal interactions between the oral epithelium and adjacent mesenchyme proceeds via stimulation with several growth factors; however, the role of Sema3A in the development of salivary glands is unknown. There may thus be a common mechanism between epithelial morphogenesis and pathogenesis; however, the role of Sema3A in salivary gland tumors is also unclear. The current study investigated the involvement of Sema3A in submandibular gland (SMG) development and its expression in adenoid cystic carcinoma (ACC) specimens. Quantitative RT-PCR and immunohistochemical analyses revealed that Sema3A was expressed both in epithelium and in mesenchyme in the initial developmental stages of SMG and their expressions were decreased during the developmental processes. Loss-of-function experiments using an inhibitor revealed that Sema3A was required for AKT activation-mediated cellular growth and formation of cleft and bud in SMG rudiment culture. In addition, Wnt/β-catenin signaling decreased the Sema3A expression in the rudiment culture. ACC arising from salivary glands frequently exhibits malignant potential. Immunohistochemical analyses of tissue specimens obtained from 10 ACC patients showed that Sema3A was hardly observed in non-tumor regions but was strongly expressed in tumor lesions, especially in myoepithelial neoplastic cells, at high frequencies where phosphorylated AKT expression was frequently detected. These results suggest that the Sema3A-AKT axis promotes cell growth, thereby contributing to morphogenesis and pathogenesis, at least in ACC, of salivary glands.

Thermoresponsive fiber-based microwells capable of formation and retrieval of salivary gland stem cell spheroids for the regeneration of irradiation-damaged salivary glands

Three-dimensional spheroid culture enhances cell-to-cell interactions among stem cells and promotes the expression of stem cell properties; however, subsequent retrieval and delivery of these cells remain a challenge. We fabricated a thermoresponsive fiber-based microwell scaffold by combining electrospinning and hydrogel micropatterning. The resultant scaffold appeared to facilitate the formation of cellular spheroids of uniform size and enabled the expression of more stem cell-secreting growth factor genes ( EGF, IGF-1, FGF1, FGF2, and HGF), pluripotent stem cell-related genes ( SOX2 and NANOG), and adult epithelial stem cell-related genes ( LGR4, LGR5, and LGR6) than salivary gland stem cells in a monolayer culture (SGSCmonolayer). The spheroids could be retrieved efficiently by decreasing temperature. SGSC-derived spheroid (SGSCspheroid) cells were then implanted into the submandibular glands of mice at 2 weeks after fractionated X-ray irradiation at a dose of 7.5 Gy/day. At 16 weeks post-irradiation, restoration of salivary function was detected only in SGSCspheroid-implanted mice. The production of submandibular acini specific mucin increased in SGSCspheroid-implanted mice, compared with PBS control. More MIST1+ mature acinar cells were preserved in the SGSCspheroid-implanted group than in the PBS control group. Intriguingly, SGSCspheroid-implanted mice exhibited greater amelioration of tissue damage and preservation of KRT7+ terminally differentiated luminal ductal cells than SGSCmonolayer-implanted mice. The SGSCspheroid-implanted mice also showed less DNA damage and apoptotic cell death than the SGSCmonolayer-implanted mice at 2 weeks post-implantation. Additionally, a significant increase in Ki67+AQP5+ proliferative acinar cells was noted only in SGSCspheroid-implanted mice. Our results suggest that a thermoresponsive fiber-based scaffold could be of use to facilitate the production of function-enhanced SGSCspheroid cells and their subsequent retrieval and delivery to damaged salivary glands to alleviate radiation-induced apoptotic cell death and promote salivary gland regeneration.

Salivary Gland Function, Development and Regeneration

Salivary glands produce and secrete saliva, which is essential for maintaining oral health and overall health. Understanding both the unique structure and physiological function of salivary glands, as well as how they are affected by disease and injury will direct the development of therapy to repair and regenerate them. Significant recent advances, particularly in the OMICS field, increase our understanding of how salivary glands develop at the cellular, molecular and genetic levels; the signaling pathways involved, the dynamics of progenitor cell lineages in development, homeostasis and regeneration and the role of the extracellular matrix microenvironment. These provide a template for cell and gene therapies as well as bioengineering approaches to repair or regenerate salivary function.

In situ functional cell phenotyping reveals microdomain networks in colorectal cancer recurrence

Tumors are dynamic ecosystems comprising localized niches (microdomains), possessing distinct compositions and spatial configurations of cancer and non-cancer cell populations. Microdomain-specific network signaling coevolves with a continuum of cell states and functional plasticity associated with disease progression and therapeutic responses. We present LEAPH, an unsupervised machine learning algorithm for identifying cell phenotypes, which applies recursive steps of probabilistic clustering and spatial regularization to derive functional phenotypes (FPs) along a continuum. Combining LEAPH with pointwise mutual information and network biology analyses enables the discovery of outcome-associated microdomains visualized as distinct spatial configurations of heterogeneous FPs. Utilization of an immunofluorescence-based (51 biomarkers) image dataset of colorectal carcinoma primary tumors (n = 213) revealed microdomain-specific network dysregulation supporting cancer stem cell maintenance and immunosuppression that associated selectively with the recurrence phenotype. LEAPH enables an explainable artificial intelligence platform providing insights into pathophysiological mechanisms and novel drug targets to inform personalized therapeutic strategies.

Organoid Models for Salivary Gland Biology and Regenerative Medicine

The salivary gland is composed of an elegant epithelial network that secrets saliva and maintains oral homeostasis. While cell lines and animal models furthered our understanding of salivary gland biology, they cannot replicate key aspects of the human salivary gland tissue, particularly the complex architecture and microenvironmental features that dictate salivary gland function. Organoid cultures provide an alternative system to recapitulate salivary gland tissue in vitro, and salivary gland organoids have been generated from pluripotent stem cells and adult stem/progenitor cells. In this review, we describe salivary gland organoids, the advances and limitations, and the promising potential for regenerative medicine.

Comparing development and regeneration in the submandibular gland highlights distinct mechanisms

A common question in organ regeneration is the extent to which regeneration recapitulates embryonic development. To investigate this concept, we compared the expression of two highly interlinked and essential genes for salivary gland development, Sox9 and Fgf10, during submandibular gland development, homeostasis and regeneration. Salivary gland duct ligation/deligation model was used as a regenerative model. Fgf10 and Sox9 expression changed during regeneration compared to homeostasis, suggesting that these key developmental genes play important roles during regeneration, however, significantly both displayed different patterns of expression in the regenerating gland compared to the developing gland. Regenerating glands, which during homeostasis had very few weakly expressing Sox9-positive cells in the striated/granular ducts, displayed elevated expression of Sox9 within these ducts. This pattern is in contrast to embryonic development, where Sox9 expression was absent in the proximally developing ducts. However, similar to the elevated expression at the distal tip of the epithelium in developing salivary glands, regenerating glands displayed elevated expression in a subpopulation of acinar cells, which during homeostasis expressed Sox9 at lower levels. A shift in expression of Fgf10 was observed from a widespread mesenchymal pattern during organogenesis to a more limited and predominantly epithelial pattern during homeostasis in the adult. This restricted expression in epithelial cells was maintained during regeneration, with no clear upregulation in the surrounding mesenchyme, as might be expected if regeneration recapitulated development. As both Fgf10 and Sox9 were upregulated in proximal ducts during regeneration, this suggests that the positive regulation of Sox9 by Fgf10, essential during development, is partially reawakened during regeneration using this model. Together these data suggest that developmentally important genes play a key role in salivary gland regeneration but do not precisely mimic the roles observed during development.

Development of the venom delivery system in Elapidae snake species: Naja siamensis and Oxyuranus microlepidotus

Many advanced snakes possess a unique venom delivery system which they utilise to subdue prey and for defence. Despite extensive efforts, the evolutionary differences in this key system between advanced snake families remains enigmatic. The current study has investigated the development of the venom delivery system using two oviparous Elapidae models, Naja siamensis and Oxyuranus microlepidotus. The development stages of the embryos in both models were detailed using previously standardised characterisation. Variations in the days post-oviposition between these stages was observed, despite a continuous development trajectory. These differences also translated to the development of the venom delivery system.

Generation of a Single-Cell RNAseq Atlas of Murine Salivary Gland Development

Understanding the dynamic transcriptional landscape throughout organ development will provide a template for regenerative therapies. Here, we generated a single-cell RNA sequencing atlas of murine submandibular glands identifying transcriptional profiles that revealed cellular heterogeneity during landmark developmental events: end bud formation, branching morphogenesis, cytodifferentiation, maturation, and homeostasis. Trajectory inference analysis suggests plasticity among acinar and duct populations. We identify transcription factors correlated with acinar differentiation including Spdef, Etv1, and Xbp1, and loss of Ybx1, Eno1, Sox11, and Atf4. Furthermore, we characterize two intercalated duct populations defined by either Gfra3 and Kit, or Gstt1. This atlas can be used to investigate specific cell functions and comparative studies predicting common mechanisms involved in development of branching organs.

Cell Sheets Restore Secretory Function in Wounded Mouse Submandibular Glands

Thermoresponsive cell culture plates release cells as confluent living sheets in response to small changes in temperature, with recovered cell sheets retaining functional extracellular matrix proteins and tight junctions, both of which indicate formation of intact and functional tissue. Our recent studies demonstrated that cell sheets are highly effective in promoting mouse submandibular gland (SMG) cell differentiation and recovering tissue integrity. However, these studies were performed only at early time points and extension of the observation period is needed to investigate duration of the cell sheets. Thus, the goal of this study was to demonstrate that treatment of wounded mouse SMG with cell sheets is capable of increasing salivary epithelial integrity over extended time periods. The results indicate that cell sheets promote tissue organization as early as eight days after transplantation and that these effects endure through Day 20. Furthermore, cell sheet transplantation in wounded SMG induces a significant time-dependent enhancement of cell polarization, differentiation and ion transporter expression. Finally, this treatment restored saliva quantity to pre-wounding levels at both eight and twenty days post-surgery and significantly improved saliva quality at twenty days post-surgery. These data indicate that cell sheets engineered with thermoresponsive cell culture plates are useful for salivary gland regeneration and provide evidence for the long-term stability of cell sheets, thereby offering a potential new therapeutic strategy for treating hyposalivation.

Developmental role of hyaluronic acid and its application in salivary gland tissue engineering

Dry mouth, or xerostomia, caused by salivary gland dysfunction significantly impacts oral/systemic health and quality of life. Although in vitro-generated artificial salivary glands have been considered as the fundamental solution, its structural complexity is difficult to reproduce using current biomaterials. Therefore, understanding and recapitulating the roles of biomacromolecules in salivary gland organogenesis is needed to solve these problems. Hyaluronic acid (HA) is a macromolecule abundant during salivary gland organogenesis, but its role remains unknown. Here, we verify the effects of HA on salivary gland organogenesis and artificial organ germ formation in solubilized and substrate-immobilized forms. In embryonic submandibular glands (eSMG), we found dense HA layers encapsulating proliferative c-Kit+ progenitor cells that were expressing CD44, an HA receptor. The blockage of HA synthesis, or degradation of HA, impaired eSMG growth by ablating the c-Kit+ progenitor cell population. We also found that high-molecular-weight (HMW) HA has a significant role in eSMG growth. Based on these findings, we discovered that HA is also crucial for in vitro formation of salivary gland organ germs, one of the most promising candidates for salivary gland tissue regeneration. We significantly enhanced salivary gland organ germ formation by supplementing HMW HA in solution; this effect was further increased when the HMW HA was immobilized on the substrate by polydopamine/HA co-immobilization. Our study suggests that the current use of HA in salivary gland tissue engineering can be further optimized.

ΔNp63 is upregulated during salivary gland regeneration following duct ligation and irradiation in mice

The transcription factor p63, a component of the p53 family, has important functions in development, homeostasis, and regeneration of epithelial tissues. However, the role of p63 in the regeneration of exocrine glands, including the salivary glands (SGs), has not been fully investigated. We investigated p63 expression in SG regeneration induced by duct ligation and irradiation. The expression of ΔNp63, a p63 isoform, increased and was colocalized with keratin 5 positive cells were myoepithelial cells. Furthermore, ΔNp63 expression was regulated by FGF7 stimulation via p38 MAPK phosphorylation and affected SG morphogenesis. These results suggest that ΔNp63 is essential for SG regeneration and may be a new target for regenerative treatment.

VDR in salivary gland homeostasis and cancer

The vitamin D receptor (VDR) and its ligand 1,25(OH)₂D₃ (1,25D) impact differentiation and exert anti-tumor effects in many tissues, but its role in salivary gland has yet to be defined. Using immunohistochemistry (IHC), we have detected strong VDR expression in murine and human salivary gland ducts. Compared to normal gland, VDR protein expression was retained in differentiated human pleomorphic adenoma (PA) but was undetectable in undifferentiated PA and in carcinomas, suggesting deregulation of VDR during salivary cancer progression. To gain insight into the potential role of VDR in salivary cancer, we assessed the effects of vitamin D in vivo and in vitro. Despite the presence of VDR in salivary gland, chronic dietary vitamin D restriction did not alter morphology of the salivary epithelium in C57/Bl6 mice. The localization of VDR in ductal epithelium prompted us to examine the effects of 1,25D in an established cell line (mSGc) derived from normal murine submandibular gland (SMG). This previously characterized cell line consists of multiple stem, progenitor and differentiated cell types as determined by mutually exclusive cellular expression of basal, ductal and myoepithelial markers. We demonstrated VDR expression and regulation of VDR target genes Vdr and Postn by 1,25D in mSGc, indicating functional ligand-mediated transcriptional activity. The effect of VDR signaling on epithelial differentiation markers was assessed by qPCR and IHC in mSGc cells treated with 1,25D. We found that 1,25D reduced mRNA expression of the basal cell progenitor marker keratin 5 (K5) and increased expression of the differentiated ductal cell marker keratin 7 (K7). Further, we found that 1,25D significantly decreased the number of proliferating cells, including proliferating K5⁺ cells. Characterization of cell cycle by Muse cytometry indicated 1,25D treatment decreased cells in S, G₂, and M phase. The inhibition of K5⁺ cell proliferation by 1,25D is of particular interest because K5⁺ basal cells contribute to a wide variety of salivary tumor types. Our studies suggest that 1,25D alters cancer-relevant progenitor and differentiation markers in the salivary gland.

Tracing tumorigenesis in a solid tumor model at single-cell resolution

Characterizing the complex composition of solid tumors is fundamental for understanding tumor initiation, progression and metastasis. While patient-derived samples provide valuable insight, they are heterogeneous on multiple molecular levels, and often originate from advanced tumor stages. Here, we use single-cell transcriptome and epitope profiling together with pathway and lineage analyses to study tumorigenesis from a developmental perspective in a mouse model of salivary gland squamous cell carcinoma. We provide a comprehensive cell atlas and characterize tumor-specific cells. We find that these cells are connected along a reproducible developmental trajectory: initiated in basal cells exhibiting an epithelial-to-mesenchymal transition signature, tumorigenesis proceeds through Wnt-differential cancer stem cell-like subpopulations before differentiating into luminal-like cells. Our work provides unbiased insights into tumor-specific cellular identities in a whole tissue environment, and emphasizes the power of using defined genetic model systems.

Understanding tumour development at a granular level is a challenge in solid tumours. Here, the authors provide a cell atlas across tumour development in a genetic model of salivary gland squamous cell carcinoma using single-cell transcriptome and epitope profiling.

Regional Differences following Partial Salivary Gland Resection

Regenerative medicine aims to repair, replace, or restore function to tissues damaged by aging, disease, or injury. Partial organ resection is not only a common clinical approach in cancer therapy but also an experimental injury model used to examine mechanisms of regeneration and repair in organs. We performed a partial resection, or partial sialoadenectomy, in the female murine submandibular salivary gland (SMG) to establish a model for investigation of repair mechanisms in salivary glands (SGs). After partial sialoadenectomy, we performed whole-gland measurements over a period of 56 d and found that the gland increased slightly in size. We used microarray analysis and immunohistochemistry (IHC) to examine messenger RNA and protein changes in glands over time. Microarray analysis identified dynamic changes in the transcriptome 3 d after injury that were largely resolved by day 14. At the 3-d time point, we detected gene signatures for cell cycle regulation, inflammatory/repair response, and extracellular matrix (ECM) remodeling in the partially resected glands. Using quantitative IHC, we identified a transient proliferative response throughout the gland. Both secretory epithelial and stromal cells expressed Ki67 that was detectable at day 3 and largely resolved by day 14. IHC also revealed that while most of the gland underwent a wound-healing response that resolved by day 14, a small region of the gland showed an aberrant sustained fibrotic response characterized by increased levels of ECM deposition, sustained Ki67 levels in stromal cells, and a persistent M2 macrophage response through day 56. The partial submandibular salivary gland resection model provides an opportunity to examine a normal healing response and an aberrant fibrotic response within the same gland to uncover mechanisms that prevent wound healing and regeneration in mammals. Understanding regional differences in the wound-healing responses may ultimately affect regenerative therapies for patients.

ROCK inhibitor increases proacinar cells in adult salivary gland organoids

Salisphere-derived adult epithelial cells have been used to improve saliva production of irradiated mouse salivary glands. Importantly, optimization of the cellular composition of salispheres could improve their regenerative capabilities. The Rho Kinase (ROCK) inhibitor, Y27632, has been used to increase the proliferation and reduce apoptosis of progenitor cells grown in vitro. In this study, we investigated whether Y27632 could be used to improve expansion of adult submandibular salivary epithelial progenitor cells or to affect their differentiation potential in different media contexts. Application of Y27632 in medium used previously to grow salispheres promoted expansion of Kit⁺ and Mist1⁺ cells, while in simple serum-containing medium Y27632 increased the number of cells that expressed the K5 basal progenitor marker. Salispheres derived from Mist1^CreERT2; R26^TdTomato mice grown in salisphere media with Y27632 included Mist1-derived cells. When these salispheres were incorporated into 3D organoids, inclusion of Y27632 in the salisphere stage increased the contribution of Mist1-derived cells expressing the proacinar/acinar marker, Aquaporin 5 (AQP5), in response to FGF2-dependent mesenchymal signals. Optimization of the cellular composition of salispheres and organoids can be used to improve the application of adult salivary progenitor cells in regenerative medicine strategies.

Establishment of a Murine Pro-acinar Cell Line to Characterize Roles for FGF2 and α3β1 Integrins in Regulating Pro-acinar Characteristics

Radiation therapy for head and neck cancers results in permanent damage to the saliva producing acinar compartment of the salivary gland. To date, a pure pro-acinar cell line to study underlying mechanisms of acinar cell differentiation in culture has not been described. Here, we report the establishment of a pro-acinar (mSG-PAC1) and ductal (mSG-DUC1) cell line, from the murine submandibular salivary gland (SMG), which recapitulate developmental milestones in differentiation. mSG-DUC1 cells express the ductal markers, keratin-7 and keratin-19, and form lumenized spheroids. mSG-PAC1 cells express the pro-acinar markers SOX10 and aquaporin-5. Using the mSG-PAC1 cell line, we demonstrate that FGF2 regulates specific steps during acinar cell maturation. FGF2 up-regulates aquaporin-5 and the expression of the α3 and α6 subunits of the α3β1 and α6β1 integrins that are known to promote SMG morphogenesis and differentiation. mSG-DUC1 and mSG-PAC1 cells were derived from genetically modified mice, homozygous for floxed alleles of the integrin α3 subunit. Similar to SMGs from α3-null mice, deletion of α3 alleles in mSG-PAC1 cells results in the up-regulation of E-cadherin and the down-regulation of CDC42. Our data indicate that mSG-DUC1 and mSG-PAC1 cells will serve as important tools to gain mechanistic insight into salivary gland morphogenesis and differentiation.

Concise Review: A Critical Evaluation of Criteria Used to Define Salivary Gland Stem Cells

In the effort to develop cell-based therapies to treat salivary gland dysfunction, many different populations of cells in the adult salivary glands have been proposed as stem cells. These cell populations vary, depending on the assay used, and are often nonoverlapping, leading to the conclusion that salivary glands harbor multiple stem cells. The goal of this review is to critically appraise the assays and properties used to identify stem cells in the adult salivary gland, and to consider the caveats of each. Re-evaluation of the defining criteria may help to reconcile the many potential stem cell populations described in the salivary gland, in order to increase comparability between studies and build consensus in the field. Stem Cells 2019;37:1144-1150.

Distribution of aquaporin-5, transforming growth factor-β1 and laminin during regeneration of atrophic rat submandibular glands after duct ligation

The aim of this study was to determine the localization of aquaporin-5 (AQP5), transforming growth factor-β₁ (TGF-β₁) and laminin during regeneration of the rat submandibular gland. After duct ligation for 7 days, the regenerating glands were collected on days 0, 1, 3, 7, and 14 after ligation release to study the process of regeneration. Immunohistochemical staining revealed apical expression of AQP5 in many acinar cells, strong expression in intercalated ducts (ICDs) of the normal submandibular gland at Day 14, and strong expression in duct-like structures (DLSs) during regeneration from Day 0 to 7. However, a few AQP5-negative acinar cells were detected during regeneration. At Day 0, immunopositivity for TGF-β₁ was detected in connective tissue. At Days 3 and 7 during regeneration, TGF-β₁ immunostaining was observed in DLSs, which were surrounded by α-smooth muscle actin-positive thickened myoepithelial cells. Laminin staining was strong in the thickened basement membrane of DLSs at Day 3 during regeneration, but weak around acinar cells at Day 14. These findings suggest that TGF-β₁ is involved in the environment around DLSs, myoepithelial cells and laminin, that DLSs have the same functional properties as ICDs, and that AQP5-negative acinar cells may be mucous cells.

Generating Embryonic Salivary Gland Organoids

Organoids are important research tools for studying organ morphogenesis and differentiation because they recapitulate ex vivo the native 3D organization of cells that is essential for proper cell and organ function. The composition of organoids can be manipulated to incorporate specific cell types to facilitate molecular interrogation of cell-cell interactions during organoid formation. A method for generating organoids derived from both embryonic salivary gland epithelial progenitor cells and mesenchymal support cells is described. Methods for isolating enriched populations of the epithelial cells as clusters and the mesenchyme cells as single cells from mouse embryonic submandibular salivary glands are also provided. Separating the epithelial and mesenchymal cell populations allows for independent molecular manipulation of each cell type. In addition, methods for lentiviral transduction of the mesenchyme cells and quantitative image analysis of organoids are provided. The methods described here are useful for exploring mechanisms driving organ formation. © 2018 by John Wiley & Sons, Inc.

Postnatal changes in the development of rat submandibular glands in offspring of diabetic mothers: Biochemical, histological and ultrastructural study

Development and maturation of submandibular salivary glands are influenced by intrauterine diabetic environment. Several studies investigated the effects of diabetes on the salivary glands. However, the effects of maternal diabetes on the submandibular glands of the offspring was not properly examined. Therefore, the present study was designed to describe the changes in the development of the submandibular glands of the offspring of diabetic mothers. The submandibular glands of the offspring of Streptozotocin (STZ)-induced diabetic female rats were examined at two and four weeks after birth. Detection of mRNA demonstrated that maternal diabetes affects the level of different indicators. The reduction of expression of epidermal growth factor (EGF); a protein mitogen, cytokeratin 5 (CK5); an epithelial cell progenitor, CK7 and aquaporin 5 (AQP5); differentiation markers and B cell lymphoma 2 (Bcl2); an antiapoptotic marker were found. Increase in Bcl2-associated X protein (Bax); an apoptotic marker was detected. These changes indicate their effects on saliva secretion, glands tumorigenesis, growth of normal oral flora and oral microbes, with decreased protein synthesis and production of xerostomia and dental caries. Loss of normal glandular architecture, significant increase in fibrosis, by the detection of collagen fibers, and stagnation of secretory granules were found with atrophic changes in the acinar cells. Marked defect of polysaccharides in the acinar cells, denoting functional changes, was manifested by significant reduction of the intensity of periodic acid-Schiff (PAS) reaction. The positive immunoreactivity of caspase-3, denoting cellular apoptosis, and minimal reaction of alpha-smooth muscle actin (α SMA) and proliferating cell nuclear antigen (PCNA) were evident in the offspring of diabetic mothers. We conclude that maternal diabetes produces degenerative effects in the structure and function of the submandibular salivary glands of the offspring, reflecting possible influences on their secretory activity affecting oral and digestive health.

c-Kit+ Cells in Adult Salivary Glands do not Function as Tissue Stem Cells

A rare population of salivary gland cells isolated based on c-Kit immunoreactivity are thought to represent tissue stem cells since they exhibit the most robust proliferative and differentiation capacity ex vivo. Despite their high promise for cell-based therapies aimed at restoring salivary function, the precise location and in vivo function of c-Kit⁺ stem cells remain unclear. Here, by combining immunostaining with c-KitCre^ERT2-based genetic labeling and lineage tracing in the adult mouse salivary glands, we show that c-Kit is expressed in a relatively large and heterogeneous cell population that consists mostly of differentiated cells. Moreover, c-Kit does not mark ductal stem cells that are known to express cytokeratin K14. Tracking the fate of in vivo-labeled c-Kit⁺ or that of K14⁺ cells in spheroid cultures reveals a limited proliferative potential for c-Kit⁺ cells and identifies K14⁺ cells as the major source of salispheres in these cultures. Long-term in vivo lineage tracing studies indicate that although c-Kit marks at least two discrete ductal cell lineages, c-Kit⁺ cells do not contribute to the normal maintenance of any other cell lineages. Our results indicate that c-Kit is not a reliable marker for salivary gland stem cells, which has important implications for salivary gland regenerative therapies.

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.

Opportunities and Challenges in Implementation of Multiparameter Single Cell Analysis Platforms for Clinical Translation

The high-content interrogation of single cells with platforms optimized for the multiparameter characterization of cells in liquid and solid biopsy samples can enable characterization of heterogeneous populations of cells ex vivo. Doing so will advance the diagnosis, prognosis, and treatment of cancer and other diseases. However, it is important to understand the unique issues in resolving heterogeneity and variability at the single cell level before navigating the validation and regulatory requirements in order for these technologies to impact patient care. Since 2013, leading experts representing industry, academia, and government have been brought together as part of the Foundation for the National Institutes of Health (FNIH) Biomarkers Consortium to foster the potential of high-content data integration for clinical translation.

A novel role for cilia-dependent sonic hedgehog signaling during submandibular gland development

BACKGROUND

Submandibular glands (SMGs) are specialized epithelial structures which generate saliva necessary for mastication and digestion. Loss of SMGs can lead to inflammation, oral lesions, fungal infections, problems with chewing/swallowing, and tooth decay. Understanding the development of the SMG is important for developing therapeutic options for patients with impaired SMG function. Recent studies have suggested Sonic hedgehog (Shh) signaling in the epithelium plays an integral role in SMG development; however, the mechanism by which Shh influences gland development remains nebulous.

RESULTS

Using the Kif3a^f/f ;Wnt1-Cre ciliopathic mouse model to prevent Shh signal transduction by means of the loss of primary cilia in neural crest cells, we report that mesenchymal Shh activity is necessary for gland development. Furthermore, using a variety of murine transgenic lines with aberrant mesenchymal Shh signal transduction, we determine that loss of Shh activity, by means of loss of the Gli activator, rather than gain of Gli repressor, is sufficient to cause the SMG aplasia. Finally, we determine that loss of the SMG correlates with reduced Neuregulin1 (Nrg1) expression and lack of innervation of the SMG epithelium.

CONCLUSIONS

Together, these data suggest a novel mechanistic role for mesenchymal Shh signaling during SMG development. Developmental Dynamics 247:818-831, 2018. © 2018 Wiley Periodicals, Inc.

Salivary glands regenerate after radiation injury through SOX2-mediated secretory cell replacement

Salivary gland acinar cells are routinely destroyed during radiation treatment for head and neck cancer that results in a lifetime of hyposalivation and co-morbidities. A potential regenerative strategy for replacing injured tissue is the reactivation of endogenous stem cells by targeted therapeutics. However, the identity of these cells, whether they are capable of regenerating the tissue, and the mechanisms by which they are regulated are unknown. Using in vivo and ex vivo models, in combination with genetic lineage tracing and human tissue, we discover a SOX2+ stem cell population essential to acinar cell maintenance that is capable of replenishing acini after radiation. Furthermore, we show that acinar cell replacement is nerve dependent and that addition of a muscarinic mimetic is sufficient to drive regeneration. Moreover, we show that SOX2 is diminished in irradiated human salivary gland, along with parasympathetic nerves, suggesting that tissue degeneration is due to loss of progenitors and their regulators. Thus, we establish a new paradigm that salivary glands can regenerate after genotoxic shock and do so through a SOX2 nerve-dependent mechanism.

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.

FGF2-dependent mesenchyme and laminin-111 are niche factors in salivary gland organoids

Epithelial progenitor cells are dependent upon a complex 3D niche to promote their proliferation and differentiation during development, which can be recapitulated in organoids. The specific requirements of the niche remain unclear for many cell types, including the proacinar cells that give rise to secretory acinar epithelial cells that produce saliva. Here, using ex vivo cultures of E16 primary mouse submandibular salivary gland epithelial cell clusters, we investigated the requirement for mesenchymal cells and other factors in producing salivary organoids in culture. Native E16 salivary mesenchyme, but not NIH3T3 cells or mesenchymal cell conditioned medium, supported robust protein expression of the progenitor marker Kit and the acinar/proacinar marker AQP5, with a requirement for FGF2 expression by the mesenchyme. Enriched salivary epithelial clusters that were grown in laminin-enriched basement membrane extract or laminin-111 together with exogenous FGF2, but not with EGF, underwent morphogenesis to form organoids that displayed robust expression of AQP5 in terminal buds. Knockdown of FGF2 in the mesenchyme or depletion of mesenchyme cells from the organoids significantly reduced AQP5 levels even in the presence of FGF2, suggesting a requirement for autocrine FGF2 signaling in the mesenchyme cells for AQP5 expression. We conclude that basement membrane proteins and mesenchyme cells function as niche factors in salivary organoids.

Immunolocalization patterns of cytokeratins during salivary acinar cell development in mice

Embryonic development of the mouse salivary glands begins with epithelial thickening and continues with sequential changes from the pre-bud to terminal bud stages. After birth, morphogenesis proceeds, and the glands develop into a highly branched epithelial structure that terminates with saliva-producing acinar cells at the adult stage. Acinar cells derived from the epithelium are differentiated into serous, mucous, and seromucous types. During differentiation, cytokeratins, intermediate filaments found in most epithelial cells, play vital roles. Although the localization patterns and developmental roles of cytokeratins in different epithelial organs, including the mammary glands, circumvallate papilla, and sweat glands, have been well studied, their stage-specific localization and morphogenetic roles during salivary gland development have yet to be elucidated. Therefore, the aim of this study was to determine the stage and acinar cell type-specific localization pattern of cytokeratins 4, 5, 7, 8, 13, 14, 18, and 19 in the major salivary glands (submandibular, sublingual, and parotid glands) of the mouse at the E15.5, PN0, PN10, and adult stages. In addition, cell physiology, including cell proliferation, was examined during development via immunostaining for Ki67 to understand the cellular mechanisms that govern acinar cell differentiation during salivary gland morphogenesis. The distinct localization patterns of cytokeratins in conjunction with cell physiology will reveal the roles of epithelial cells in salivary gland formation during the differentiation of serous, mucous or seromucous salivary glands.

RARα and RARγ reciprocally control K5+ progenitor cell expansion in developing salivary glands

Understanding the mechanisms of controlled expansion and differentiation of basal progenitor cell populations during organogenesis is essential for developing targeted regenerative therapies. Since the cytokeratin 5-positive (K5+) basal epithelial cell population in the salivary gland is regulated by retinoic acid signaling, we interrogated how isoform-specific retinoic acid receptor (RAR) signaling impacts the K5+ cell population during salivary gland organogenesis to identify RAR isoform-specific mechanisms that could be exploited in future regenerative therapies. In this study, we utilized RAR isoform-specific inhibitors and agonists with murine submandibular salivary gland organ explants. We determined that RARα and RARγ have opposing effects on K5+ cell cycle progression and cell distribution. RARα negatively regulates K5+ cells in both whole organ explants and in isolated epithelial rudiments. In contrast, RARγ is necessary but not sufficient to positively maintain K5+ cells, as agonism of RARγ alone failed to significantly expand the population. Although retinoids are known to stimulate differentiation, K5 levels were not inversely correlated with differentiated ductal cytokeratins. Instead, RARα agonism and RARγ inhibition, corresponding with reduced K5, resulted in premature lumenization, as marked by prominin-1. With lineage tracing, we demonstrated that K5+ cells have the capacity to become prominin-1+ cells. We conclude that RARα and RARγ reciprocally control K5+ progenitor cells endogenously in the developing submandibular salivary epithelium, in a cell cycle-dependent manner, controlling lumenization independently of keratinizing differentiation. Based on these data, isoform-specific targeting RARα may be more effective than pan-RAR inhibitors for regenerative therapies that seek to expand the K5+ progenitor cell pool.

SUMMARY STATEMENT

RARα and RARγ reciprocally control K5+ progenitor cell proliferation and distribution in the developing submandibular salivary epithelium in a cell cycle-dependent manner while regulating lumenization independently of keratinizing differentiation.

Endothelial cell regulation of salivary gland epithelial patterning

Perfusion-independent regulation of epithelial pattern formation by the vasculature during organ development and regeneration is of considerable interest for application in restoring organ function. During murine submandibular salivary gland development, the vasculature co-develops with the epithelium during branching morphogenesis; however, it is not known whether the vasculature has instructive effects on the epithelium. Using pharmacological inhibitors and siRNA knockdown in embryonic organ explants, we determined that VEGFR2-dependent signaling is required for salivary gland epithelial patterning. To test directly for a requirement for endothelial cells in instructive epithelial patterning, we developed a novel ex vivo cell fractionation/reconstitution assay. Immuno-depletion of CD31⁺ endothelial cells in this assay confirmed a requirement for endothelial cells in epithelial patterning of the gland. Depletion of endothelial cells or inhibition of VEGFR2 signaling in organ explants caused an aberrant increase in cells expressing the ductal proteins K19 and K7, with a reduction in Kit⁺ progenitor cells in the endbuds of reconstituted glands. Addition of exogenous endothelial cells to reconstituted glands restored epithelial patterning, as did supplementation with the endothelial cell-regulated mesenchymal factors IGFBP2 and IGFBP3. Our results demonstrate that endothelial cells promote expansion of Kit⁺ progenitor cells and suppress premature ductal differentiation in early developing embryonic submandibular salivary gland buds.

SOX2 regulates acinar cell development in the salivary gland

Acinar cells play an essential role in the secretory function of exocrine organs. Despite this requirement, how acinar cells are generated during organogenesis is unclear. Using the acini-ductal network of the developing human and murine salivary gland, we demonstrate an unexpected role for SOX2 and parasympathetic nerves in generating the acinar lineage that has broad implications for epithelial morphogenesis. Despite SOX2 being expressed by progenitors that give rise to both acinar and duct cells, genetic ablation of SOX2 results in a failure to establish acini but not ducts. Furthermore, we show that SOX2 targets acinar-specific genes and is essential for the survival of acinar but not ductal cells. Finally, we illustrate an unexpected and novel role for peripheral nerves in the creation of acini throughout development via regulation of SOX2. Thus, SOX2 is a master regulator of the acinar cell lineage essential to the establishment of a functional organ.

DOI:http://dx.doi.org/10.7554/eLife.26620.001

The salivary glands produce fluid that contains enzymes to help us to digest our food. These glands contain a tree-like network of cells – known as acinar cells – that produce the fluid, and cells that form ducts to transport the fluid out of the glands. Both types of cells form from stem cells as animal embryos develop. Like all developing organs, the salivary glands receive many different signals that guide how they grow. However, the identity of the cues that instruct a stem cell to produce a new acinar cell or duct cell are not known.

Emmerson et al. studied how the salivary glands develop in mouse embryos. The experiments show that a protein called SOX2 – which is an essential regulator of stem cells in embryos – is required for acinar cells to form. Loss of SOX2 inhibited the production of acinar but not duct cells. Furthermore, nerves that surround the gland provide support to cells that produce SOX2 and promote the formation of acinar cells.

Further experiments suggest that the nerves also play the same role in humans. Adult organs often use developmental signals to repair or regenerate tissue. As such, understanding how an organ develops may lead to new therapies that can stimulate salivary glands and other organs to regenerate after they have been damaged in adults.

DOI:http://dx.doi.org/10.7554/eLife.26620.002

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3228-7) contains supplementary material, which is available to authorized users.

Par-1b is required for morphogenesis and differentiation of myoepithelial cells during salivary gland development

The salivary epithelium initiates as a solid mass of epithelial cells that are organized into a primary bud that undergoes morphogenesis and differentiation to yield bilayered acini consisting of interior secretory acinar cells that are surrounded by contractile myoepithelial cells in mature salivary glands. How the primary bud transitions into acini has not been previously documented. We document here that the outer epithelial cells subsequently undergo a vertical compression as they express smooth muscle α-actin and differentiate into myoepithelial cells. The outermost layer of polarized epithelial cells assemble and organize the basal deposition of basement membrane, which requires basal positioning of the polarity protein, Par-1b. Whether Par-1b is required for the vertical compression and differentiation of the myoepithelial cells is unknown. Following manipulation of Par-1b in salivary gland organ explants, Par-1b-inhibited explants showed both a reduced vertical compression of differentiating myoepithelial cells and reduced levels of smooth muscle α-actin. Rac1 knockdown and inhibition of Rac GTPase function also inhibited branching morphogenesis. Since Rac regulates cellular morphology, we investigated a contribution for Rac in myoepithelial cell differentiation. Inhibition of Rac GTPase activity showed a similar reduction in vertical compression and smooth muscle α-actin levels while decreasing the levels of Par-1b protein and altering its basal localization in the outer cells. Inhibition of ROCK, which is required for basal positioning of Par-1b, resulted in mislocalization of Par-1b and loss of vertical cellular compression, but did not significantly alter levels of smooth muscle α-actin in these cells. Overexpression of Par-1b in the presence of Rac inhibition restored basement membrane protein levels and localization. Our results indicate that the basal localization of Par-1b in the outer epithelial cells is required for myoepithelial cell compression, and Par-1b is required for myoepithelial differentiation, regardless of its localization.

The WNT/MYB pathway suppresses KIT expression to control the timing of salivary proacinar differentiation and duct formation

Growth factor signaling is involved in the development of various organs, but how signaling regulates organ morphogenesis and differentiation in a coordinated manner remains to be clarified. Here, we show how WNT signaling controls epithelial morphogenetic changes and differentiation using the salivary gland as a model. Experiments using genetically manipulated mice and organ cultures revealed that WNT signaling at an early stage (E12-E15) of submandibular salivary gland (SMG) development inhibits end bud morphogenesis and differentiation into proacini by suppressing Kit expression through the upregulation of the transcription factor MYB, and concomitantly increasing the expression of distal progenitor markers. In addition, WNT signaling at the early stage of SMG development promoted end bud cell proliferation, leading to duct formation. WNT signaling reduction at a late stage (E16-E18) of SMG development promoted end bud maturation and suppressed duct formation. Thus, WNT signaling controls the timing of SMG organogenesis by keeping end bud cells in an undifferentiated bipotent state.

The role of human fibronectin- or placenta basement membrane extract-based gels in favouring the formation of polarized salivary acinar-like structures

Head and neck cancer patients treated with radiotherapy commonly experience hyposalivation and oral/tooth infections, leading to a reduced quality of life. Clinical management is currently unsatisfactory for dry mouth. Thus, there is a need for growing salivary fluid-secreting (acinar) cells for these patients. However, functionally-grown salivary acinar cells are cultured in Matrigel, a product that cannot be used clinically, owing to its source from a mouse sarcoma. Therefore, finding a gel suitable for clinical use and possessing properties similar to that of Matrigel would allow biopsied salivary cells to be expanded in vitro and transplanted into the mouths of xerostomic patients. This study tested gels made with human placenta basement membrane extract (BME) or fibronectin for the growth and differentiation of human salivary biopsies into acinar cells. We report here that, following expansion of primary human salivary gland epithelial cells (huSGs) in serum-free medium, using these gels (made from human proteins) allowed morphological and functional differentiation of salivary ductal cells into acinar-like cells. These (human) gels gave comparable results to Matrigel, such as differentiation into polarized acinar 3D units or monolayers with tight junction proteins (claudin-1, -2, -3) and exhibiting adequate transepithelial electrical resistance, acinar proteins (AQP5, α-amylase, mucin-1, NKCC1) and acinar adhesion-related cell markers (CD44, CD166). Ultrastructural, mRNA and protein analyses confirmed the formation of differentiated acinar polarized cells. The mitotic activity was highest with human placenta BME gel. This human culture model provided a reproducible approach to studying human salivary cell expansion and differentiation for tissue engineering. Copyright © 2016 John Wiley & Sons, Ltd.

Multiscale analysis of the murine intestine for modeling human diseases

When functioning properly, the intestine is one of the key interfaces between the human body and its environment. It is responsible for extracting nutrients from our food and excreting our waste products. It provides an environment for a host of healthful microbes and serves as a first defense against pathogenic ones. These processes require tight homeostatic controls, which are provided by the interactions of a complex mix of epithelial, stromal, neural and immune cells, as well as the resident microflora. This homeostasis can be disrupted by invasive microbes, genetic lesions, and carcinogens, resulting in diseases such Clostridium difficile infection, inflammatory bowel disease (IBD) and cancer. Enormous strides have been made in understanding how this important organ functions in health and disease using everything from cell culture systems to animal models to human tissue samples. This has resulted in better therapies for all of these diseases, but there is still significant room for improvement. In the United States alone, 14,000 people per year die of C. difficile, up to 1.6 million people suffer from IBD, and more than 50,000 people die every year from colon cancer. Because these and other intestinal diseases arise from complex interactions between the different components of the gut ecosystem, we propose that systems approaches that address this complexity in an integrative manner may eventually lead to improved therapeutics that deliver lasting cures. This review will discuss the use of systems biology for studying intestinal diseases in vivo with particular emphasis on mouse models. Additionally, it will focus on established experimental techniques that have been used to drive this systems-level analysis, and emerging techniques that will push this field forward in the future.

Identification of Stem Cells in the Secretory Complex of Salivary Glands

Salivary glands have an essential secretory function for maintaining oral and overall health. The epithelial compartment of the gland is composed of several highly specialized cell types that cooperate to secrete and deliver saliva to the oral cavity. The mouse submandibular gland has been used as a model for major salivary glands in human. The secretory complex in this model is composed of 2 secretory compartments, including acini and granular ducts connected by intercalated ducts. Contractile myoepithelial cells surround the secretory complex to facilitate salivary flow. Whether differentiated cells in the secretory complex are maintained by self-duplication or contribution from stem cells has remained an open question. Here, in analyzing the expression of basal cytokeratin (K) 14 in the secretory complex, we discovered a subset of K14(+) ductal cells in the intercalated ducts of the adult gland. These cells are distinct from the K14-expressing basal/myoepithelial cells, proliferate at a significantly higher rate than any other epithelial cell type in the gland, and reside in a spatially defined domain within the intercalated duct. Using inducible genetic lineage tracing, we show that K14(+) ductal cells represent a long-lived yet cycling population of stem cells that are established during development and contribute to the formation and maintenance of the granular ducts throughout life. Our data provide direct evidence for the existence of stem cells contributing to homeostasis of salivary glands, as well as new insights into glandular pathobiology.

Immunohistochemical localization of keratin 5 in the submandibular gland in adult and postnatal developing mice

Keratin 5 (K5) is a marker of basal progenitor cells in the epithelia of a number of organs. During prenatal development of the submandibular gland (SMG) in mice, K5(+) progenitor cells in the developing epithelia play important roles in its organogenesis. Although K5(+) cells are also present in the adult mouse SMG and may function in tissue regeneration, their histological localization has not yet investigated in detail. In the present study, we examined the immunohistochemical localization of K5 in the SMG in adult and postnatal developing mice. At birth, K5 immunoreactivity was detected in the entire duct system, in which it was localized in the basal cells of a double-layered epithelium, but was not detected in the terminal tubule or myoepithelial cells. At postnatal weeks 1-3, with the development of intercalated ducts (ID), striated ducts (SD), and excretory ducts (ED), K5-immunoreactive basal cells were gradually restricted to the ED and the proximal double-layered portions of the ID connecting to the SD. At the same time, K5 immunoreactivity appeared in myoepithelial cells, in which its positive ratio gradually increased. In adults, K5 immunoreactivity was localized to most myoepithelial cells, most basal cells in the ED, and a small number of ID cells at the boundary between the ID and SD in the female SMG or between the ID and granular convoluted tubules in the male SMG. These results suggest that K5 is a marker of differentiated myoepithelial cells and duct progenitor cells in the mouse SMG.

Changes in the Submandibular Salivary Gland Epithelial Cell Subpopulations During Progression of Sjögren's Syndrome-Like Disease in the NOD/ShiLtJ Mouse Model

Sjögren's syndrome (SS), an autoimmune exocrinopathy, is associated with dysfunction of the secretory salivary gland epithelium, leading to xerostomia. The etiology of SS disease progression is poorly understood as it is typically not diagnosed until late stage. Since mouse models allow the study of disease progression, we investigated the NOD/ShiLtJ mouse to explore temporal changes to the salivary epithelium. In the NOD/ShiLtJ model, SS presents secondary to autoimmune diabetes, and SS disease is reportedly fully established by 20 weeks. We compared epithelial morphology in the submandibular salivary glands (SMG) of NOD/ShiLtJ mice with SMGs from the parental strain at 12, 18, and 22 weeks of age and used immunofluorescence to detect epithelial proteins, including the acinar marker, aquaporin 5, ductal cell marker, cytokeratin 7, myoepithelial cell marker, smooth muscle α-actin, and the basal cell marker, cytokeratin 5, while confirming immune infiltrates with CD45R. We also compared these proteins in the labial salivary glands of human SS patients with control tissues. In the NOD/ShiLtJ SMG, regions of lymphocytic infiltrates were not associated with widespread epithelial tissue degradation; however, there was a decrease in the area of the gland occupied by secretory epithelial cells in favor of ductal epithelial cells. We observed an expansion of cells expressing cytokeratin 5 within the ducts and within the smooth muscle α-actin(+) basal myoepithelial population. The altered acinar/ductal ratio within the NOD/ShiLtJ SMG likely contributes to salivary hypofunction, while the expansion of cytokeratin 5 positive-basal cells may reflect loss of function or indicate a regenerative response.

Regenerating Salivary Glands in the Microenvironment of Induced Pluripotent Stem Cells

This report describes our initial attempt to regenerate salivary glands using induced pluripotent stem (iPS) cells in vivo and in vitro. Glandular tissues that were similar to the adult submandibular glands (SMGs) and sublingual glands could be partially produced by the transplantation of iPS cells into mouse salivary glands. However, the tumorigenicity of iPS cells has not been resolved yet. It is well known that stem cells affect their microenvironment, known as a stem cell niche. We focused on the niche and the interaction between iPS cells and salivary gland cells in our study on salivary gland regeneration. Coculture of embryonic SMG cells and iPS cells have better-developed epithelial structures and fewer undifferentiated specific markers than monoculture of embryonic SMG cells in vitro. These results suggest that iPS cells have a potential ability to accelerate differentiation for salivary gland development and regeneration.