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

Specific 3-O-sulfated heparan sulfate domains regulate salivary gland basement membrane metabolism and epithelial differentiation

Heparan sulfate (HS) regulation of FGFR function, which is essential for salivary gland (SG) development, is determined by the immense structural diversity of sulfated HS domains. 3-O-sulfotransferases generate highly 3-O-sulfated HS domains (3-O-HS), and Hs3st3a1 and Hs3st3b1 are enriched in myoepithelial cells (MECs) that produce basement membrane (BM) and are a growth factor signaling hub. Hs3st3a1;Hs3st3b1 double-knockout (DKO) mice generated to investigate 3-O-HS regulation of MEC function and growth factor signaling show loss of specific highly 3-O-HS and increased FGF/FGFR complex binding to HS. During development, this increases FGFR-, BM- and MEC-related gene expression, while in adult, it reduces MECs, increases BM and disrupts acinar polarity, resulting in salivary hypofunction. Defined 3-O-HS added to FGFR pulldown assays and primary organ cultures modulates FGFR signaling to regulate MEC BM synthesis, which is critical for secretory unit homeostasis and acinar function. Understanding how sulfated HS regulates development will inform the use of HS mimetics in organ regeneration.

Salivary gland organoid transplantation as a therapeutic option for radiation-induced xerostomia

BACKGROUND

Xerostomia is a pathological condition characterized by decreased salivation due to salivary gland dysfunction and is frequently attributed to irreversible damage as a side effect of radiation therapy. Stem cell-derived organoid therapy has garnered attention as a promising avenue for resolving this issue. However, Matrigel, a hydrogel commonly used in organoid culture, is considered inappropriate for clinical use due to its undefined composition and immunogenicity. In this study, we aimed to develop a method for culturing collagen-based human salivary gland organoids (hSGOs) suitable for clinical applications and evaluated their therapeutic effectiveness.

METHODS

Human salivary gland stem cells were isolated from the salivary gland tissues and cultured in both Matrigel and collagen. We compared the gene and protein expression patterns of salivary gland-specific markers and measured amylase activity in the two types of hSGOs. To evaluate the therapeutic effects, we performed xenogeneic and allogeneic transplantation using human and mouse salivary gland organoids (hSGOs and mSGOs), respectively, in a mouse model of radiation-induced xerostomia.

RESULTS

hSGOs cultured in Matrigel exhibited self-renewal capacity and differentiated into acinar and ductal cell lineages. In collagen, they maintained a comparable self-renewal ability and more closely replicated the characteristics of salivary gland tissue following differentiation. Upon xenotransplantation of collagen-based hSGOs, we observed engraftment, which was verified by detecting human-specific nucleoli and E-cadherin expression. The expression of mucins, especially MUC5B, within the transplanted hSGOs suggested a potential improvement in the salivary composition. Moreover, the allograft procedure using mSGOs led to increased salivation, validating the efficacy of our approach.

CONCLUSIONS

This study showed that collagen-based hSGOs can be used appropriately in clinical settings and demonstrated the effectiveness of an allograft procedure. Our research has laid the groundwork for the future application of collagen-based hSGOs in allogeneic clinical trials.

Representing core gene expression activity relationships using the latent structure implicit in Bayesian networks

MOTIVATION

Many types of networks, such as co-expression or ChIP-seq-based gene-regulatory networks, provide useful information for biomedical studies. However, they are often too full of connections and difficult to interpret, forming "indecipherable hairballs."

RESULTS

To address this issue, we propose that a Bayesian network can summarize the core relationships between gene expression activities. This network, which we call the LatentDAG, is substantially simpler than conventional co-expression network and ChIP-seq networks (by two orders of magnitude). It provides clearer clusters, without extraneous cross-cluster connections, and clear separators between modules. Moreover, one can find a number of clear examples showing how it bridges the connection between steps in the transcriptional regulatory network and other networks (e.g. RNA-binding protein). In conjunction with a graph neural network, the LatentDAG works better than other biological networks in a variety of tasks, including prediction of gene conservation and clustering genes.

AVAILABILITY AND IMPLEMENTATION

Code is available at https://github.com/gersteinlab/LatentDAG.

Primary Sjögren's syndrome: new perspectives on salivary gland epithelial cells

Primary Sjögren's syndrome (pSS) is a chronic autoimmune disease primarily affecting exocrine glands such as the salivary glands, leading to impaired secretion and sicca symptoms. As the mainstay of salivation, salivary gland epithelial cells (SGECs) have an important role in the pathology of pSS. Emerging evidence suggests that the interplay between immunological factors and SGECs may not be the initial trigger or the sole mechanism responsible for xerostomia in pSS, challenging conventional perceptions. To deepen our understanding, current research regarding SGECs in pSS was reviewed. Among the extensive aberrations in cellular architecture and function, this review highlighted certain alterations of SGECs that were identified to occur independently of or in absence of lymphocytic infiltration. In particular, some of these alterations may serve as upstream factors of immuno-inflammatory responses. These findings underscore the significance of introspecting the pathogenesis of pSS and developing interventions targeting SGECs in the early stages of the disease.

Salivary Gland Tissue Recombination Can Modify Cell Fate

Although mesenchyme is essential for inducing the epithelium of ectodermal organs, its precise role in organ-specific epithelial fate determination remains poorly understood. To elucidate the roles of tissue interactions in cellular differentiation, we performed single-cell RNA sequencing and imaging analyses on recombined tissues, where mesenchyme and epithelium were switched ex vivo between two types of embryonic mouse salivary glands: the parotid gland (a serous gland) and the submandibular gland (a predominantly mucous gland). We found partial induction of molecules that define gland-specific acinar and myoepithelial cells in recombined salivary epithelium. The parotid epithelium recombined with submandibular mesenchyme began to express mucous acinar genes not intrinsic to the parotid gland. While myoepithelial cells do not normally line parotid acini, newly induced myoepithelial cells densely populated recombined parotid acini. However, mucous acinar and myoepithelial markers continued to be expressed in submandibular epithelial cells recombined with parotid mesenchyme. Consequently, some epithelial cells appeared to be plastic, such that their fate could still be modified in response to mesenchymal signaling, whereas other epithelial cells appeared to be already committed to a specific fate. We also discovered evidence for bidirectional induction: transcriptional changes were observed not only in the epithelium but also in the mesenchyme after heterotypic tissue recombination. For example, parotid epithelium induced the expression of muscle-related genes in submandibular fibroblasts that began to mimic parotid fibroblast gene expression. These studies provide the first comprehensive unbiased molecular characterization of tissue recombination approaches exploring the regulation of cell fate.

Analysis of the heterogeneity and complexity of murine extraorbital lacrimal gland via single-cell RNA sequencing

PURPOSE

The lacrimal gland is essential for maintaining ocular surface health and avoiding external damage by secreting an aqueous layer of the tear film. However, a healthy lacrimal gland's inventory of cell types and heterogeneity remains understudied.

METHODS

Here, 10X Genome-based single-cell RNA sequencing was used to generate an unbiased classification of cellular diversity in the extraorbital lacrimal gland (ELG) of C57BL/6J mice. From 43,850 high-quality cells, we produced an atlas of cell heterogeneity and defined cell types using classic marker genes. The possible functions of these cells were analyzed through bioinformatics analysis. Additionally, the CellChat was employed for a preliminary analysis of the cell-cell communication network in the ELG.

RESULTS

Over 37 subclasses of cells were identified, including seven types of glandular epithelial cells, three types of fibroblasts, ten types of myeloid-derived immune cells, at least eleven types of lymphoid-derived immune cells, and five types of vascular-associated cell subsets. The cell-cell communication network analysis revealed that fibroblasts and immune cells play a pivotal role in the dense intercellular communication network within the mouse ELG.

CONCLUSIONS

This study provides a comprehensive transcriptome atlas and related database of the mouse ELG.

Salivary gland developmental mechanics

The salivary gland undergoes branching morphogenesis to elaborate into a tree-like structure with numerous saliva-secreting acinar units, all joined by a hierarchical ductal system. The expansive epithelial surface generated by branching morphogenesis serves as the structural basis for the efficient production and delivery of saliva. Here, we elucidate the process of salivary gland morphogenesis, emphasizing the role of mechanics. Structurally, the developing salivary gland is characterized by a stratified epithelium tightly encased by the basement membrane, which is in turn surrounded by a mesenchyme consisting of a dense network of interstitial matrix and mesenchymal cells. Diverse cell types and extracellular matrices bestow this developing organ with organized, yet spatially varied mechanical properties. For instance, the surface epithelial sheet of the bud is highly fluidic due to its high cell motility and weak cell-cell adhesion, rendering it highly pliable. In contrast, the inner core of the bud is more rigid, characterized by reduced cell motility and strong cell-cell adhesion, which likely provide structural support for the tissue. The interactions between the surface epithelial sheet and the inner core give rise to budding morphogenesis. Furthermore, the basement membrane and the mesenchyme offer mechanical constraints that could play a pivotal role in determining the higher-order architecture of a fully mature salivary gland.

Bioprinting salivary gland models and their regenerative applications

OBJECTIVE

Salivary gland (SG) hypofunction is a common clinical condition arising from radiotherapy to suppress head and neck cancers. The radiation often destroys the SG secretory acini, and glands are left with limited regenerative potential. Due to the complex architecture of SG acini and ducts, three-dimensional (3D) bioprinting platforms have emerged to spatially define these in vitro epithelial units and develop mini-organs or organoids for regeneration. Due to the limited body of evidence, this comprehensive review highlights the advantages and challenges of bioprinting platforms for SG regeneration.

METHODS

SG microtissue engineering strategies such as magnetic 3D bioassembly of cells and microfluidic coaxial 3D bioprinting of cell-laden microfibers and microtubes have been proposed to replace the damaged acinar units, avoid the use of xenogeneic matrices (like Matrigel), and restore salivary flow.

RESULTS

Replacing the SG damaged organ is challenging due to its complex architecture, which combines a ductal network with acinar epithelial units to facilitate a unidirectional flow of saliva. Our research group was the first to develop 3D bioassembly SG epithelial functional organoids with innervation to respond to both cholinergic and adrenergic stimulation. More recently, microtissue engineering using coaxial 3D bioprinting of hydrogel microfibers and microtubes could also supported the formation of viable epithelial units. Both bioprinting approaches could overcome the need for Matrigel by facilitating the assembly of adult stem cells, such as human dental pulp stem cells, and primary SG cells into micro-sized 3D constructs able to produce their own matrix and self-organize into micro-modular tissue clusters with lumenized areas. Furthermore, extracellular vesicle (EV) therapies from organoid-derived secretome were also designed and validated ex vivo for SG regeneration after radiation damage.

CONCLUSION

Magnetic 3D bioassembly and microfluidic coaxial bioprinting platforms have the potential to create SG mini-organs for regenerative applications via organoid transplantation or organoid-derived EV therapies.

ΔNp63 regulates Sfrp1 expression to direct salivary gland branching morphogenesis

Branching morphogenesis is a complex process shared by many organs including the lungs, kidney, prostate, as well as several exocrine organs including the salivary, mammary and lacrimal glands. This critical developmental program ensures the expansion of an organ's surface area thereby maximizing processes of cellular secretion or absorption. It is guided by reciprocal signaling from the epithelial and mesenchymal cells. While signaling pathways driving salivary gland branching morphogenesis have been relatively well-studied, our understanding of the underlying transcriptional regulatory mechanisms directing this program, is limited. Here, we performed in vivo and ex vivo studies of the embryonic mouse submandibular gland to determine the function of the transcription factor ΔNp63, in directing branching morphogenesis. Our studies show that loss of ΔNp63 results in alterations in the differentiation program of the ductal cells which is accompanied by a dramatic reduction in branching morphogenesis that is mediated by dysregulation of WNT signaling. We show that ΔNp63 modulates WNT signaling to promote branching morphogenesis by directly regulating Sfrp1 expression. Collectively, our findings have revealed a novel role for ΔNp63 in the regulation of this critical process and offers a better understanding of the transcriptional networks involved in branching morphogenesis.

Immunomodulation of salivary gland function due to cancer therapy

Functional salivary glands (SG) are essential for maintaining oral health, and salivary dysfunction is a persistent major clinical challenge. Several cancer therapies also have off-target effects leading to SG dysfunction. Recent advances highlight the role of SG immune populations in homeostasis, dysfunction and gland regeneration. Here, we review what is known about SG immune populations during development and postnatal homeostasis. We summarize recent findings of immune cell involvement in SG dysfunction following cancer treatments such as irradiation (IR) for head and neck cancers, immune transplant leading to graft-versus-host-disease (GVHD) and immune checkpoint inhibitor (ICI) treatment. The role of immune cells in SG in both homeostasis and disease, is an emerging field of research that may provide important clues to organ dysfunction and lead to novel therapeutic targets.

Single-cell analysis reveals the transcriptional alterations in the submandibular glands of aged mice

OBJECTIVES

Aging-related salivary gland changes, such as lymphocyte infiltration and acinar cell loss decrease saliva secretion, thereby affecting quality of life. The precise molecular mechanisms underlying these changes remain unclear.

METHODS

We here performed single-cell RNA sequencing to clarify gene expression changes in each cell type comprising the submandibular glands (SMGs) of adult and aged mice.

RESULTS

The proportion of acinar cells decreased in various epithelial clusters annotated with cell type-specific marker genes. Expression levels of the cellular senescence markers, Cdkn2a/p16 and Cdkn1a/p21, were increased in the basal and striated ducts of aged SMGs relative to their levels in those of adult SMGs. In contrast, senescence-associated secretory phenotype-related genes, except transforming growth factor-β, exhibited little change in expression in aged SMGs relative to adult SMGs.

CONCLUSIONS

Gene Ontology analysis revealed increased expression levels of genes encoding major histocompatibility complex (MHC) class I components in the ductal component cells of aged SMGs. MHC class I expression may thus be associated with salivary gland aging.

CFTR expression in human salivary gland acinar cells

The key role of CFTR in secretory epithelia has been extensively documented. Additionally, CFTR plays a significant role in ion absorption in exocrine glands, including salivary and sweat glands. Most of the knowledge about CFTR expression comes from animal models such as the mouse or the rat, but there is limited information about CFTR expression in human tissues. In the present study, we assessed the expression of CFTR in human submandibular and parotid glands. Consistent with findings in rodent salivary glands, our immunolocalization studies show that CFTR is expressed in duct cells. However, CFTR expression in human salivary glands differs from that in rodents, as immunolocalization and single-cell RNA sequencing analysis from a previous study performed in the human parotid gland revealed the presence of CFTR protein and transcripts within a distinct cell cluster. Based on cell marker expression, this cluster corresponds to acinar cells. To obtain functional evidence supporting CFTR expression, we isolated human parotid acinar cells through collagenase digestion. Acinar cells displayed an anion conductance that was activated in response to cAMP-increasing agents and was effectively blocked by CFTRInh172, a known CFTR blocker. This study provides novel evidence of CFTR expression within acinar cells of human salivary glands. This finding challenges the established model positioning CFTR exclusively in duct cells from exocrine glands.NEW & NOTEWORTHY This study addresses the uncertainty about the impact of CFTR on human salivary gland function. We found CFTR transcripts in a subset of duct cells known as ionocytes, as well as in acinar cells. Isolated human parotid acinar cells exhibited Cl- conductance consistent with CFTR activity. This marks the first documented evidence of functional CFTR expression in human salivary gland acinar cells.

The Mayo Clinic Salivary Tissue-Organoid Biobanking: A Resource for Salivary Regeneration Research

The salivary gland (SG) is an essential organ that secretes saliva, which supports versatile oral function throughout life, and is maintained by elusive epithelial stem and progenitor cells (SGSPC). Unfortunately, aging, drugs, autoimmune disorders, and cancer treatments can lead to salivary dysfunction and associated health consequences. Despite many ongoing therapeutic efforts to mediate those conditions, investigating human SGSPC is challenging due to lack of standardized tissue collection, limited tissue access, and inadequate purification methods. Herein, we established a diverse and clinically annotated salivary regenerative biobanking at the Mayo Clinic, optimizing viable salivary cell isolation and clonal assays in both 2D and 3D-matrigel growth environments. Our analysis identified ductal epithelial cells in vitro enriched with SGSPC expressing the CD24/EpCAM/CD49f+ and PSMA- phenotype. We identified PSMA expression as a reliable SGSPC differentiation marker. Moreover, we identified progenitor cell types with shared phenotypes exhibiting three distinct clonal patterns of salivary differentiation in a 2D environment. Leveraging innovative label-free unbiased LC-MS/MS-based single-cell proteomics, we identified 819 proteins across 71 single cell proteome datasets from purified progenitor-enriched parotid gland (PG) and sub-mandibular gland (SMG) cultures. We identified distinctive co-expression of proteins, such as KRT1/5/13/14/15/17/23/76 and 79, exclusively observed in rare, scattered salivary ductal basal cells, indicating the potential de novo source of SGSPC. We also identified an entire class of peroxiredoxin peroxidases, enriched in PG than SMG, and attendant H2O2-dependent cell proliferation in vitro suggesting a potential role for PRDX-dependent floodgate oxidative signaling in salivary homeostasis. The distinctive clinical resources and research insights presented here offer a foundation for exploring personalized regenerative medicine.

Mesenchymal stem cell-derived HGF attenuates radiation-induced senescence in salivary glands via compensatory proliferation

BACKGROUND & AIM

Irradiation of the salivary glands during head and neck cancer treatment induces cellular senescence in response to DNA damage and contributes to radiation-induced hyposalivation by affecting the salivary gland stem/progenitor cell (SGSC) niche. Cellular senescence, such as that induced by radiation, is a state of cell-cycle arrest, accompanied by an altered pro-inflammatory secretome known as the senescence-associated secretory phenotype (SASP) with potential detrimental effects on the surrounding microenvironment. We hypothesized that the pro-regenerative properties of mesenchymal stem cells (MSCs) may attenuate cellular senescence post-irradiation. Therefore, here we evaluated the effects of adipose-derived MSCs (ADSCs) on the radiation-induced response of salivary gland organoids (SGOs).

METHODS

Proteomic analyses to identify soluble mediators released by ADSCs co-cultured with SGOS revealed secretion of hepatocyte growth factor (HGF) in ADSCs, suggesting a possible role in the stem cell crosstalk. Next, the effect of recombinant HGF in the culture media of ex vivo grown salivary gland cells was tested in 2D monolayers and 3D organoid models.

RESULTS

Treatment with HGF robustly increased salivary gland cell proliferation. Importantly, HGF supplementation post-irradiation enhanced proliferation at lower doses of radiation (0, 3, 7 Gy), but not at higher doses (10, 14 Gy) where most cells stained positive for senescence-associated beta-galactosidase. Furthermore, HGF had no effect on the senescence-associated secretory phenotype (SASP) of irradiated SGOs, suggesting there may be compensatory proliferation by cell-division competent cells instead of a reversal of cellular senescence after irradiation.

CONCLUSION

ADSCs may positively influence radiation recovery through HGF secretion and can promote the ex vivo expansion of salivary gland stem/progenitor cells to enhance the effects of co-transplanted SGSC.

A closer look into the cellular and molecular biology of myoepithelial cells across various exocrine glands

Myoepithelial cells (MECs) are a unique subset of epithelial cells that possess several smooth muscle cell characteristics, such as a high number of actin-myosin filaments and the ability to contract. These cells are primarily located around the secretory cells of exocrine glands, including the salivary, mammary, lacrimal, and sweat glands. Their primary functions involve the construction of the basement membrane and help with secretion of gland products through contraction. So far, no comparative analysis of MECs in different exocrine glands had ever evaluated their differences. In this review, we took advantage of the various publicly available scRNAseq data from mouse exocrine glands to identify their shared and unique characteristics. The aim of this review is to compare the role of MECs in maintaining healthy glandular function, their involvement in disease states, and their regenerative capacity, with a particular emphasis on the latest research findings in these areas.

Salivary Gland Bioengineering

Salivary gland dysfunction affects millions globally, and tissue engineering may provide a promising therapeutic avenue. This review delves into the current state of salivary gland tissue engineering research, starting with a study of normal salivary gland development and function. It discusses the impact of fibrosis and cellular senescence on salivary gland pathologies. A diverse range of cells suitable for tissue engineering including cell lines, primary salivary gland cells, and stem cells are examined. Moreover, the paper explores various supportive biomaterials and scaffold fabrication methodologies that enhance salivary gland cell survival, differentiation, and engraftment. Innovative engineering strategies for the improvement of vascularization, innervation, and engraftment of engineered salivary gland tissue, including bioprinting, microfluidic hydrogels, mesh electronics, and nanoparticles, are also evaluated. This review underscores the promising potential of this research field for the treatment of salivary gland dysfunction and suggests directions for future exploration.

CSF1R-dependent macrophages in the salivary gland are essential for epithelial regeneration after radiation-induced injury

The salivary glands often become damaged in individuals receiving radiotherapy for head and neck cancer, resulting in chronic dry mouth. This leads to detrimental effects on their health and quality of life, for which there is no regenerative therapy. Macrophages are the predominant immune cell in the salivary glands and are attractive therapeutic targets due to their unrivaled capacity to drive tissue repair. Yet, the nature and role of macrophages in salivary gland homeostasis and how they may contribute to tissue repair after injury are not well understood. Here, we show that at least two phenotypically and transcriptionally distinct CX3CR1+ macrophage populations are present in the adult salivary gland, which occupy anatomically distinct niches. CD11c+CD206-CD163- macrophages typically associate with gland epithelium, whereas CD11c-CD206+CD163+ macrophages associate with blood vessels and nerves. Using a suite of complementary fate mapping systems, we show that there are highly dynamic changes in the ontogeny and composition of salivary gland macrophages with age. Using an in vivo model of radiation-induced salivary gland injury combined with genetic or antibody-mediated depletion of macrophages, we demonstrate an essential role for macrophages in clearance of cells with DNA damage. Furthermore, we show that epithelial-associated macrophages are indispensable for effective tissue repair and gland function after radiation-induced injury, with their depletion resulting in reduced saliva production. Our data, therefore, provide a strong case for exploring the therapeutic potential of manipulating macrophages to promote tissue repair and thus minimize salivary gland dysfunction after radiotherapy.

FGFR2 is essential for salivary gland duct homeostasis and MAPK-dependent seromucous acinar cell differentiation

Exocrine acinar cells in salivary glands (SG) are critical for oral health and loss of functional acinar cells is a major clinical challenge. Fibroblast growth factor receptors (FGFR) are essential for early development of multiple organs, including SG. However, the role of FGFR signaling in specific populations later in development and during acinar differentiation are unknown. Here, we use scRNAseq and conditional deletion of murine FGFRs in vivo to identify essential roles for FGFRs in craniofacial, early SG development and progenitor function during duct homeostasis. Importantly, we also discover that FGFR2 via MAPK signaling is critical for seromucous acinar differentiation and secretory gene expression, while FGFR1 is dispensable. We show that FGF7, expressed by myoepithelial cells (MEC), activates the FGFR2-dependent seromucous transcriptional program. Here, we propose a model where MEC-derived FGF7 drives seromucous acinar differentiation, providing a rationale for targeting FGFR2 signaling in regenerative therapies to restore acinar function.

Craniofacial developmental biology in the single-cell era

The evolution of a unique craniofacial complex in vertebrates made possible new ways of breathing, eating, communicating and sensing the environment. The head and face develop through interactions of all three germ layers, the endoderm, ectoderm and mesoderm, as well as the so-called fourth germ layer, the cranial neural crest. Over a century of experimental embryology and genetics have revealed an incredible diversity of cell types derived from each germ layer, signaling pathways and genes that coordinate craniofacial development, and how changes to these underlie human disease and vertebrate evolution. Yet for many diseases and congenital anomalies, we have an incomplete picture of the causative genomic changes, in particular how alterations to the non-coding genome might affect craniofacial gene expression. Emerging genomics and single-cell technologies provide an opportunity to obtain a more holistic view of the genes and gene regulatory elements orchestrating craniofacial development across vertebrates. These single-cell studies generate novel hypotheses that can be experimentally validated in vivo. In this Review, we highlight recent advances in single-cell studies of diverse craniofacial structures, as well as potential pitfalls and the need for extensive in vivo validation. We discuss how these studies inform the developmental sources and regulation of head structures, bringing new insights into the etiology of structural birth anomalies that affect the vertebrate head.

Single-cell Transcriptomic Analysis of Salivary Gland Endothelial Cells

Vascular endothelial cells have important functions in fibrosis via direct and indirect methods and in regeneration through secretion of tissue-specific, paracrineacting angiocrine factors. In the salivary gland, endothelial cells are required for proper development, but their roles within adult glands are largely unknown. The goal of this work was to identify ligand-receptor interactions between endothelial cells and other cell types that are important during homeostasis, fibrosis, and regeneration. To model salivary gland fibrosis and regeneration, we utilized a reversible ductal ligation. To induce injury, a clip was applied to the primary ducts for 14 days, and to induce a regenerative response, the clip was subsequently removed for 5 days. To identify endothelial cell-produced factors, we used single-cell RNA-sequencing of stromal-enriched cells from adult submandibular and sublingual salivary glands. Transcriptional profiles of homeostatic salivary gland endothelial cells were compared to endothelial cells of other organs. Salivary gland endothelial cells were found to express unique genes and displayed the highest overlap in gene expression with other fenestrated endothelial cells from the colon, small intestine, and kidney. Comparison of the 14-day ligated, mock ligated, and 5-day deligated stromal-enriched transcripts and lineage tracing were used to identify evidence for a partial endoMT phenotype, which was observed in a small number of endothelial cell subsets with ligation. CellChat was used to predict changes in ligand-receptor interactions in response to ligation and deligation. CellChat predicted that after ligation, endothelial cells are sources of protein tyrosine phosphatase receptor type m, tumor necrosis factor ligand superfamily member 13, and myelin protein zero signaling and targets for tumor necrosis factor signaling. Following deligation, CellChat predicted that endothelial cells are sources of chemokine (C-X-C motif) and EPH signaling to promote regenerative responses. These studies will inform future endothelial cell-based regenerative therapies.

Lipopolysaccharide increases bitter taste sensitivity via epigenetic changes in Tas2r gene clusters

T2R bitter receptors, encoded by Tas2r genes, are not only critical for bitter taste signal transduction but also important for defense against bacteria and parasites. However, little is known about whether and how Tas2r gene expression are regulated. Here, we show that in an inflammation model mimicking bacterial infection using lipopolysaccharide, the expression of many Tas2rs was significantly upregulated and mice displayed markedly increased neural and behavioral responses to bitter compounds. Using single-cell assays for transposase-accessible chromatin with sequencing (scATAC-seq), we found that the chromatin accessibility of Tas2rs was highly celltype specific and lipopolysaccharide increased the accessibility of many Tas2rs. scATAC-seq also revealed substantial chromatin remodeling in immune response genes in taste tissue stem cells, suggesting potential long-lasting effects. Together, our results suggest an epigenetic mechanism connecting inflammation, Tas2r gene regulation, and altered bitter taste, which may explain heightened bitter taste that can occur with infections and cancer treatments.

The First Transcriptomic Atlas of the Adult Lacrimal Gland Reveals Epithelial Complexity and Identifies Novel Progenitor Cells in Mice

The lacrimal gland (LG) secretes aqueous tears. Previous studies have provided insights into the cell lineage relationships during tissue morphogenesis. However, little is known about the cell types composing the adult LG and their progenitors. Using scRNAseq, we established the first comprehensive cell atlas of the adult mouse LG to investigate the cell hierarchy, its secretory repertoire, and the sex differences. Our analysis uncovered the complexity of the stromal landscape. Epithelium subclustering revealed myoepithelial cells, acinar subsets, and two novel acinar subpopulations: Tfrchi and Car6hi cells. The ductal compartment contained Wfdc2+ multilayered ducts and an Ltf+ cluster formed by luminal and intercalated duct cells. Kit+ progenitors were identified as: Krt14+ basal ductal cells, Aldh1a1+ cells of Ltf+ ducts, and Sox10+ cells of the Car6hi acinar and Ltf+ epithelial clusters. Lineage tracing experiments revealed that the Sox10+ adult populations contribute to the myoepithelial, acinar, and ductal lineages. Using scRNAseq data, we found that the postnatally developing LG epithelium harbored key features of putative adult progenitors. Finally, we showed that acinar cells produce most of the sex-biased lipocalins and secretoglobins detected in mouse tears. Our study provides a wealth of new data on LG maintenance and identifies the cellular origin of sex-biased tear components.

Evaluating the transcriptional landscape and cell-cell communication networks in chronically irradiated parotid glands

Understanding the transcriptional landscape that results in chronic salivary hypofunction after irradiation will help identify injury mechanisms and develop regenerative therapies. We present scRNA-seq analysis from control and irradiated murine parotid glands collected 10 months after irradiation. We identify a population of secretory cells defined by specific expression of Etv1, which may be an acinar cell precursor. Acinar and Etv1+ secretory express Ntrk2 and Erbb3, respectively while the ligands for these receptors are expressed in myoepithelial and stromal cells. Furthermore, our data suggests that secretory cells and CD4+CD8+T-cells are the most transcriptionally affected during chronic injury with radiation, suggesting active immune involvement. Lastly, evaluation of cell-cell communication networks predicts that neurotrophin, neuregulin, ECM, and immune signaling are dysregulated after irradiation, and thus may play a role in the lack of repair. This resource will be helpful to understand cell-specific pathways that may be targeted to repair chronic damage in irradiated glands.

Major depression-related factor NEGR1 controls salivary secretion in mouse submandibular glands

Salivary gland cells, which secrete water in response to neuronal stimulation, are closely connected to other neurons. Transcriptomic studies show that salivary glands also express some proteins responsible for neuronal function. However, the physiological functions of these common neuro-exocrine factors in salivary glands are largely unknown. Here, we studied the function of Neuronal growth regulator 1 (NEGR1) in the salivary gland cells. NEGR1 was also expressed in mouse and human salivary glands. The structure of salivary glands of Negr1 knockout (KO) mice was normal. Negr1 KO mice showed tempered carbachol- or thapsigargin-induced intracellular Ca²⁺ increases and store-operated Ca²⁺ entry. Of interest, the activity of the large-conductance Ca²⁺-activated K⁺ channel (BK channel) was increased, whereas Ca²⁺-activated Cl^- channel ANO1 channel activity was not altered in Negr1 KO mice. Pilocarpine- and carbachol-induced salivation was decreased in Negr1 KO mice. These results suggest that NEGR1 influence salivary secretion though the muscarinic Ca²⁺ signaling.

High-Resolution Transcriptomic Landscape of the Human Submandibular Gland

Saliva-secreting and transporting cells are part of the complex cellular milieu of the human salivary gland, where they play important roles in normal glandular physiology and diseased states. However, comprehensive molecular characterization, particularly at single-cell resolution, is still incomplete, in part due to difficulty in procuring normal human tissues. Here, we perform an in-depth analysis of male and female adult human submandibular gland (SMG) samples by bulk RNA sequencing (RNA-seq) and examine the molecular underpinnings of the heterogeneous cell populations by single-cell (sc) RNA-seq. Our results from scRNA-seq highlight the remarkable diversity of clusters of epithelial and nonepithelial cells that reside in the SMG that is also faithfully recapitulated by deconvolution of the bulk-RNA data sets. Our analyses reveal complex transcriptomic heterogeneity within both the ductal and acinar subpopulations and identify atypical SMG cell types, such as mucoacinar cells that are unique to humans and ionocytes that have been recently described in the mouse. We use CellChat to explore ligand-receptor interactome predictions that likely mediate crucial cell-cell communications between the various cell clusters. Finally, we apply a trajectory inference method to investigate specific cellular branching points and topology that offers insights into the dynamic and complex differentiation process of the adult SMG. The data sets and the analyses herein comprise an extensive wealth of high-resolution information and a valuable resource for a deeper mechanistic understanding of human SMG biology and pathophysiology.

Nuclear receptor Nr5a2 promotes diverse connective tissue fates in the jaw

Organ development involves the sustained production of diverse cell types with spatiotemporal precision. In the vertebrate jaw, neural-crest-derived progenitors produce not only skeletal tissues but also later-forming tendons and salivary glands. Here we identify the pluripotency factor Nr5a2 as essential for cell-fate decisions in the jaw. In zebrafish and mice, we observe transient expression of Nr5a2 in a subset of mandibular postmigratory neural-crest-derived cells. In zebrafish nr5a2 mutants, nr5a2-expressing cells that would normally form tendons generate excess jaw cartilage. In mice, neural-crest-specific Nr5a2 loss results in analogous skeletal and tendon defects in the jaw and middle ear, as well as salivary gland loss. Single-cell profiling shows that Nr5a2, distinct from its roles in pluripotency, promotes jaw-specific chromatin accessibility and gene expression that is essential for tendon and gland fates. Thus, repurposing of Nr5a2 promotes connective tissue fates to generate the full repertoire of derivatives required for jaw and middle ear 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.

Transcriptomic profile comparison reveals conservation of ionocytes across multiple organs

Single-cell RNA sequencing has recently led to the identification of a flurry of rare, new cell types, such as the CFTR-high ionocytes in the airway epithelium. Ionocytes appear to be specifically responsible for fluid osmolarity and pH regulation. Similar cells exist in multiple other organs and have received various names, including intercalated cell in the kidney, mitochondria-rich cell in the inner ear, clear cell in the epididymis, and ionocyte in the salivary gland. Here, we compare the previously published transcriptomic profile of cells expressing FOXI1, the signature transcription factor expressed in airway ionocytes. Such FOXI1+ cells were found in datasets representing human and/or murine kidney, airway, epididymis, thymus, skin, inner ear, salivary gland, and prostate. This allowed us to assess the similarities between these cells and identify the core transcriptomic signature of this ionocyte 'family'. Our results demonstrate that, across all these organs, ionocytes maintain the expression of a characteristic set of genes, including FOXI1, KRT7, and ATP6V1B1. We conclude that the ionocyte signature defines a class of closely related cell types across multiple mammalian organs.

Genetic Study of Elf5 and Ehf in the Mouse Salivary Gland

Salivary gland (SG) development, maturation, and homeostasis require coordinated roles of transcription factors (TFs) that dictate specific cell identities and fate. The ETS family of proteins are important transcriptional drivers of diverse cell lineages, tissue development, and differentiation programs and hence are also likely to play an important role in the SG. Here we have leveraged genomic and epigenomic data of the SG to examine the expression profile of ETS genes and identified 2 closely related paralogs, Elf5 and Ehf, that are highly expressed in distinct epithelial subpopulations. By using a well-defined mouse knockout model of Elf5, we show that Elf5, despite its enriched expression in the acinar cells, is functionally dispensable for maintaining the homeostatic state of the adult SG epithelium. The lack of a discernible phenotype of the Elf5-null SG might be due to possible functional redundancy with Ehf or other ETS factors. To probe this possibility and to examine the specific consequences of Ehf loss in the SG, we used CRISPR-Cas9 to generate mice in which the DNA-binding ETS domain of Ehf is disrupted due to an insertion mutation. We demonstrate that the Ehf mutant (EhfMut) mice exhibit a distinct cellular phenotype with decreased granular convoluted tubules that are accompanied by an increased accumulation of the intercalated Sox9-positive ductal cell population. Interestingly, the ductal phenotype of the EhfMut animals is highly pronounced in males, reaffirming the established sexual dimorphism of the SG that exists in rodents. Our results show that unlike Elf5, Ehf plays a nonredundant role in directing ductal cell differentiation of the SG and highlights the phenotypic subtlety in mutant mice of closely related TFs and the importance of careful consideration of cell type-specific studies.

FGFR2b is essential for salivary gland duct homeostasis and MAPK-dependent seromucous acinar cell differentiation

Exocrine secretory acinar cells in salivary glands (SG) are critical for oral health and loss of functional acinar cells is a major clinical challenge. Fibroblast growth factor receptors (FGFR) are essential for early development of multiple organs, including SG. However, the role of FGFR signaling in specific epithelial SG populations later in development and during acinar differentiation are unknown. Here, we predicted FGFR dependence in specific populations using scRNAseq data and conditional mouse models to delete FGFRs in vivo. We identifed essential roles for FGFRs in craniofacial and early SG development, as well as progenitor function during duct homeostasis. Importantly, we discovered that FGFR2b was critical for seromucous and serous acinar cell differentiation and secretory gene expression (Bpifa2 and Lpo) via MAPK signaling, while FGFR1b was dispensable. We show that FGF7, expressed by myoepithelial cells (MEC), activated the FGFR2b-dependent seromucous transcriptional program. We propose a model where MEC-derived FGF7 drives seromucous acinar differentiaton, providing a rationale for targeting FGFR2b signaling in regenerative therapies to restore acinar function.

Inflammation induces bitter taste oversensitization via epigenetic changes in Tas2r gene clusters

T2R bitter receptors, encoded by Tas2r genes, are not only critical for bitter taste signal transduction but also important for defense against bacteria and parasites. However, little is known about whether and how Tas2r gene expression are regulated. Here we show that, in an inflammation model mimicking bacterial infection, the expression of many Tas2rs are significantly up-regulated and mice displayed markedly increased neural and behavioral responses to bitter compounds. Using single-cell assays for transposase-accessible chromatin with sequencing (scATAC-seq), we found that the chromatin accessibility of Tas2rs was highly cell type specific and inflammation increased the accessibility of many Tas2rs . scATAC-seq also revealed substantial chromatin remodeling in immune response genes in taste tissue stem cells, suggesting potential long-term effects. Together, our results suggest an epigenetic mechanism connecting inflammation, Tas2r gene regulation, and altered bitter taste, which may explain heightened bitter taste that can occur with infections and cancer treatments.

Δ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.

Derivation of Human Salivary Epithelial Progenitors from Pluripotent Stem Cells via Activation of RA and Wnt Signaling

Derivation of salivary gland epithelial progenitors (SGEPs) from human pluripotent stem cells (hPSCs) has great potential in developmental biology and regenerative medicine. At present, no efficient method is available to generate salivary gland cells from hPSCs. Here, we described for the first time a robust protocol for direct differentiation of hPSCs into SGEPs by mimicking retinoic acid and Wnt signaling. These hPSC-derived SGEPs expressed SOX9, KRT5, and KRT19, important progenitor markers of developing salivary glands. CD24 and α-SMA positive cells, capable of restoring the functions of injured salivary glands, were also present in SGEP cultures. Importantly, RNA-sequencing revealed that the SGEPs resembled the transcript profiles of human fetal submandibular glands. Therefore, we provided an efficient protocol to induce hPSCs differentiation into SGEPs. Our study provides a foundation for generating functional hPSCs derived salivary gland acinar cells and three-dimensional organoids, potentially serving as new models for basic study and future translational research.

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.

The submandibular and sublingual glands maintain oral microbial homeostasis through multiple antimicrobial proteins

Introduction

Oral microbial homeostasis is a key factor affecting oral health, and saliva plays a significant role in maintaining oral microbial homeostasis. The submandibular gland (SMG) and sublingual gland (SLG) together produce the most saliva at rest. Organic ingredients, including antimicrobial proteins, are rich and distinctive and depend on the type of acinar cells in the SMG and SLG. However, the functions of the SMG and SLG in maintaining oral microbial homeostasis have been difficult to identify and distinguish, given their unique anatomical structures.

Methods

In this study, we independently removed either the SMG or SLG from mouse models. SMGs were aseptically removed in three mice in the SMG-removal group, and SLGs were aseptically removed in three mice in the SLG-removal group. Three mice from the sham-operated group were only anesthetized and incised the skin. After one month, we analyzed their oral microbiome through 16S rRNA sequencing. And then, we analyzed each gland using proteomics and single-cell RNA sequencing.

Results

Our study revealed that the microbiome balance was significantly disturbed, with decreased bacterial richness, diversity, and uniformity in the groups with the SMG or SLG removed compared with the sham-operated group. We identified eight secreted proteins in the SMG and two in the SLG that could be involved in maintaining oral microbial homeostasis. Finally, we identified multiple types of cells in the SMG and SLG (including serous acinar, mucinous acinar, ductal epithelial, mesenchymal, and immune cells) that express potential microbiota homeostasis regulatory proteins. Our results suggest that both the SMG and SLG play crucial roles in maintaining oral microbial homeostasis via excretion. Furthermore, the contribution of the SMG in maintaining oral microbial homeostasis appears to be superior to that of the SLG. These findings also revealed the possible antimicrobial function of gland secreta.

Discussion

Our results suggest that control of oral microbial dysbiosis is necessary when the secretory function of the SMG or SLG is impaired. Our study could be the basis for further research on the prevention of oral diseases caused by microbial dysbiosis.

[Developments in Research on Salivary Biomarkers in the Diagnosis of Systemic Diseases]

Saliva, an important biological fluid secreted by oral glands, serves multiple functions. It performs cleaning and protective functions for oral tissues, safeguarding against biological, mechanical and chemical stimuli, while allowing for the sensory perception of taste and temperature. It is also responsible for the preliminary digestion of food. These functions and properties of saliva are attributed to the presence of electrolytes, buffers, proteins, glycoproteins, and lipids in saliva. Recent studies have found that saliva contains biomarkers that are closely connected with the pathophysiological status of the human body, suggesting that saliva makes an ideal biological fluid for drug monitoring and biomarker screening. Therefore, salivary biomarkers can be used as an instrument for physical monitoring and localization of the occurrence of diseases, thereby accomplishing early diagnosis of diseases and assessment of the overall health status of patients. However, the actual application of salivary biomarkers in the diagnosis and treatment of systemic diseases is still not widely available, and the establishment of evaluation criteria and the exploration of its mechanism are not sufficiently investigated. Herein, we reviewed the latest research findings on applying the salivary biomarkers in the diagnosis of systemic diseases.

Neuronal-epithelial cross-talk drives acinar specification via NRG1-ERBB3-mTORC2 signaling

Acinar cells are the principal secretory units of multiple exocrine organs. A single-cell, layered, lumenized acinus forms from a large cohort of epithelial progenitors that must initiate and coordinate three cellular programs of acinar specification, namely, lineage progression, secretion, and polarization. Despite this well-known outcome, the mechanism(s) that regulate these complex programs are unknown. Here, we demonstrate that neuronal-epithelial cross-talk drives acinar specification through neuregulin (NRG1)-ERBB3-mTORC2 signaling. Using single-cell and global RNA sequencing of developing murine salivary glands, we identified NRG1-ERBB3 to precisely overlap with acinar specification during gland development. Genetic deletion of Erbb3 prevented cell lineage progression and the establishment of lumenized, secretory acini. Conversely, NRG1 treatment of isolated epithelia was sufficient to recapitulate the development of secretory acini. Mechanistically, we found that NRG1-ERBB3 regulates each developmental program through an mTORC2 signaling pathway. Thus, we reveal that a neuronal-epithelial (NRG1/ERBB3/mTORC2) mechanism orchestrates the creation of functional acini.

Single-cell transcriptomic analysis of the tumor ecosystem of adenoid cystic carcinoma

Adenoid cystic carcinoma (ACC) is a malignant tumor that originates from exocrine gland epithelial cells. We profiled the transcriptomes of 49,948 cells from paracarcinoma and carcinoma tissues of three patients using single-cell RNA sequencing. Three main types of the epithelial cells were identified into myoepithelial-like cells, intercalated duct-like cells, and duct-like cells by marker genes. And part of intercalated duct-like cells with special copy number variations which altered with MYB family gene and EN1 transcriptomes were identified as premalignant cells. Developmental pseudo-time analysis showed that the premalignant cells eventually transformed into malignant cells. Furthermore, MYB and MYBL1 were found to belong to two different gene modules and were expressed in a mutually exclusive manner. The two gene modules drove ACC progression into different directions. Our findings provide novel evidence to explain the high recurrence rate of ACC and its characteristic biological behavior.

Role of Snai2 and Notch signaling in salivary gland myoepithelial cell fate

Myoepithelial (ME) cells in exocrine glands exhibit both epithelial and mesenchymal features, contributing to fluid secretion through contraction. However, the regulation mechanism of behind this unique phenotype in salivary glands remains unclear. We established a flow cytometry-based purification method using cell surface molecules, epithelial cell adhesion molecule (EpCAM) and alpha 6 integrin (CD49f), to characterize ME cells. EpCAM⁺CD49f^high cells showed relatively high expression of ME cell-marker genes, such as alpha-smooth muscle actin (α-SMA). For lineage tracing and strict isolation, tdTomato⁺EpCAM⁺CD49f^high-ME cells were obtained from myosin heavy chain 11 (Myh11) -CreER^T2/tdTomato mice. Transcriptome analysis revealed that expression of genes involved in the epithelial-mesenchymal transition, including Snai2, were upregulated in the ME cell-enriched subset. Snai2 suppression in stable ME cells decreased α-SMA and increased Krt14 expression, suggesting that ME cell features may be controlled by the epithelial-mesenchymal balance regulated by Snai2. In contrast, ME cells showed reduced ME properties and expressed the ductal markers Krt18/19 under sphere culture conditions. Notch signaling was activated under sphere culture conditions; excessive activation of Notch signaling accelerated Krt18/19 expression, but reduced α-SMA and Snai2 expression, suggesting that the behavior of Snai2-expressing ME cells may be controlled by Notch signaling.

A Roadmap for the Human Oral and Craniofacial Cell Atlas

Oral and craniofacial tissues are uniquely adapted for continuous and intricate functioning, including breathing, feeding, and communication. To achieve these vital processes, this complex is supported by incredible tissue diversity, variously composed of epithelia, vessels, cartilage, bone, teeth, ligaments, and muscles, as well as mesenchymal, adipose, and peripheral nervous tissue. Recent single cell and spatial multiomics assays-specifically, genomics, epigenomics, transcriptomics, proteomics, and metabolomics-have annotated known and new cell types and cell states in human tissues and animal models, but these concepts remain limitedly explored in the human postnatal oral and craniofacial complex. Here, we highlight the collaborative and coordinated efforts of the newly established Oral and Craniofacial Bionetwork as part of the Human Cell Atlas, which aims to leverage single cell and spatial multiomics approaches to first understand the cellular and molecular makeup of human oral and craniofacial tissues in health and to then address common and rare diseases. These powerful assays have already revealed the cell types that support oral tissues, and they will unravel cell types and molecular networks utilized across development, maintenance, and aging as well as those affected in diseases of the craniofacial complex. This level of integration and cell annotation with partner laboratories across the globe will be critical for understanding how multiple variables, such as age, sex, race, and ancestry, influence these oral and craniofacial niches. Here, we 1) highlight these recent collaborative efforts to employ new single cell and spatial approaches to resolve our collective biology at a higher resolution in health and disease, 2) discuss the vision behind the Oral and Craniofacial Bionetwork, 3) outline the stakeholders who contribute to and will benefit from this network, and 4) outline directions for creating the first Human Oral and Craniofacial Cell Atlas.

Sci-Seq of Human Fetal Salivary Tissue Introduces Human Transcriptional Paradigms and a Novel Cell Population

Multiple pathologies and non-pathological factors can disrupt the function of the non-regenerative human salivary gland including cancer and cancer therapeutics, autoimmune diseases, infections, pharmaceutical side effects, and traumatic injury. Despite the wide range of pathologies, no therapeutic or regenerative approaches exist to address salivary gland loss, likely due to significant gaps in our understanding of salivary gland development. Moreover, identifying the tissue of origin when diagnosing salivary carcinomas requires an understanding of human fetal development. Using computational tools, we identify developmental branchpoints, a novel stem cell-like population, and key signaling pathways in the human developing salivary glands by analyzing our human fetal single-cell sequencing data. Trajectory and transcriptional analysis suggest that the earliest progenitors yield excretory duct and myoepithelial cells and a transitional population that will yield later ductal cell types. Importantly, this single-cell analysis revealed a previously undescribed population of stem cell-like cells that are derived from SD and expresses high levels of genes associated with stem cell-like function. We have observed these rare cells, not in a single niche location but dispersed within the developing duct at later developmental stages. Our studies introduce new human-specific developmental paradigms for the salivary gland and lay the groundwork for the development of translational human therapeutics.

Understanding the development of oral epithelial organs through single cell transcriptomic analysis

During craniofacial development, the oral epithelium begins as a morphologically homogeneous tissue that gives rise to locally complex structures, including the teeth, salivary glands and taste buds. How the epithelium is initially patterned and specified to generate diverse cell types remains largely unknown. To elucidate the genetic programs that direct the formation of distinct oral epithelial populations, we mapped the transcriptional landscape of embryonic day 12 mouse mandibular epithelia at single cell resolution. Our analysis identified key transcription factors and gene regulatory networks that define different epithelial cell types. By examining the spatiotemporal patterning process along the oral-aboral axis, our results propose a model in which the dental field is progressively confined to its position by the formation of the aboral epithelium anteriorly and the non-dental oral epithelium posteriorly. Using our data, we also identified Ntrk2 as a proliferation driver in the forming incisor, contributing to its invagination. Together, our results provide a detailed transcriptional atlas of the embryonic mandibular epithelium, and unveil new genetic markers and regulators that are present during the specification of various oral epithelial structures.

Single‐cell RNA sequencing in oral science: Current awareness and perspectives

The emergence of single‐cell RNA sequencing enables simultaneous sequencing of thousands of cells, making the analysis of cell population heterogeneity more efficient. In recent years, single‐cell RNA sequencing has been used in the investigation of heterogeneous cell populations, cellular developmental trajectories, stochastic gene transcriptional kinetics, and gene regulatory networks, providing strong support in life science research. However, the application of single‐cell RNA sequencing in the field of oral science has not been reviewed comprehensively yet. Therefore, this paper reviews the development and application of single‐cell RNA sequencing in oral science, including fields of tissue development, teeth and jaws diseases, maxillofacial tumors, infections, etc., providing reference and prospects for using single‐cell RNA sequencing in studying the oral diseases, tissue development, and regeneration.

The Cellular and Molecular Landscape of Synchronous Pediatric Sialoblastoma and Hepatoblastoma

Sialoblastoma (SBL) is an infrequent embryonal malignant tumor originating from the salivary gland, resembling primitive salivary gland anlage, whereas hepatoblastoma (HB) is the most common pediatric liver malignancy. The simultaneous occurrence of both tumors is extremely rare. Here we reported a case of a 6-month-old infant diagnosed with synchronous SBL and HB. The patient received neoadjuvant chemotherapy followed by surgical resection. Fresh tissues of both tumors were collected before and after chemotherapy, which were further profiled by whole exome sequencing (WES) and single-cell RNA sequencing (scRNA-seq). WES analysis revealed potential somatic driver mutation PIK3CA p.Glu454Lys for SBL and canonical mutation CTNNB1 p.Ser45Pro for HB. No shared somatic variants or common copy number alterations were found between SBL and HB primary tumor samples. Though scRNA-seq, single-cell atlases were constructed for both tumors. SBL may recapitulate a pre-acinar stage in the development of salivary gland, including basaloid, duct-like, myoepithelial-like, and cycling phenotypes. In the meantime, HB was composed of tumor cells resembling different stages of the liver, including hepatocyte-like, hepatic progenitor-like, and hepatoblast-like cells. After chemotherapy, both tumors were induced into a more mature phenotype. In terms of transcriptional signatures, SBL and HB showed enhanced expression of epithelial markers KRT8, KRT18, and essential embryo development genes SDC1, MDK, indicating the disruption of normal embryo epithelium development. Finally, heterozygous deleterious germline mutation BLM and FANCI were identified which could predispose the patient to higher cancer risk. It partially explained the reason for the co-occurrence of SBL and HB. Taken together, we provided valuable resources for deciphering cellular heterogeneity and adaptive change of tumor cells after chemotherapy for synchronous SBL and HB, providing insights into the mechanisms leading to synchronous pediatric tumors.

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.

Optimizing Soluble Cues for Salivary Gland Tissue Mimetics Using a Design of Experiments (DoE) Approach

The development of therapies to prevent or treat salivary gland dysfunction has been limited by a lack of functional in vitro models. Specifically, critical markers of salivary gland secretory phenotype downregulate rapidly ex vivo. Here, we utilize a salivary gland tissue chip model to conduct a design of experiments (DoE) approach to test combinations of seven soluble cues that were previously shown to maintain or improve salivary gland cell function. This approach uses statistical techniques to improve efficiency and accuracy of combinations of factors. The DoE-designed culture conditions improve markers of salivary gland function. Data show that the EGFR inhibitor, EKI-785, maintains relative mRNA expression of Mist1, a key acinar cell transcription factor, while FGF10 and neurturin promote mRNA expression of Aqp5 and Tmem16a, channel proteins involved in secretion. Mist1 mRNA expression correlates with increased secretory function, including calcium signaling and mucin (PAS-AB) staining. Overall, this study demonstrates that media conditions can be efficiently optimized to support secretory function in vitro using a DoE approach.

Reversely immortalized mouse salivary gland cells presented a promising metabolic and fibrotic response upon BMP9/Gdf2 stimulation

The submandibular gland (SMG) and the sublingual gland (SLG) are two of the three major salivary glands in mammals. In mice, they are adjacent to each other and open into the oral cavity, producing saliva to lubricate the mouth and aid in food digestion. Though salivary gland dysfunction accompanied with fibrosis and metabolic disturbance is common in clinic, in-depth mechanistic research is lacking. Currently, research on how to rescue salivary function is challenging, as it must resort to using terminally differentiated acinar cells or precursor acinar cells with unknown differentiation. In this study, we established reversely immortalized mouse primary SMG cells (iSMGCs) and SLG cells (iSLGCs) on the first postnatal day (P0). The iSMGCs and iSLGCs grew well, exhibited many salivary gland characteristics, and retained the metabolism-related genes derived from the original tissue as demonstrated using transcriptome sequencing (RNA-seq) analysis. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of these two cell lines, which overlapped with those of the SMG and SLG, were enriched in cysteine and methionine metabolism. Furthermore, we investigated the role of bone morphogenetic protein 9 (BMP9), also known as growth differentiation factor 2(Gdf2), on metabolic and fibrotic functions in the SMG and SLG. We demonstrated that iSMGCs and iSLGCs presented promising adipogenic and fibrotic responses upon BMP9/Gdf2 stimulation. Thus, our findings indicate that iSMGCs and iSLGCs faithfully reproduce characteristics of SMG and SLG cells and present a promising prospect for use in future study of salivary gland metabolism and fibrosis upon BMP9/Gdf2 stimulation.

Bioengineering in salivary gland regeneration

Salivary gland (SG) dysfunction impairs the life quality of many patients, such as patients with radiation therapy for head and neck cancer and patients with Sjögren’s syndrome. Multiple SG engineering strategies have been considered for SG regeneration, repair, or whole organ replacement. An in-depth understanding of the development and differentiation of epithelial stem and progenitor cells niche during SG branching morphogenesis and signaling pathways involved in cell–cell communication constitute a prerequisite to the development of suitable bioengineering solutions. This review summarizes the essential bioengineering features to be considered to fabricate an engineered functional SG model using various cell types, biomaterials, active agents, and matrix fabrication methods. Furthermore, recent innovative and promising approaches to engineering SG models are described. Finally, this review discusses the different challenges and future perspectives in SG bioengineering.

Regulation of myoepithelial differentiation

The salivary gland can be permanently impaired by radiation treatment for head and neck cancers. Efforts at tissue regeneration have focused on saliva-producing acinar cells. However, myoepithelial cells are also critical to gland function, but mechanisms that regulate their differentiation are poorly defined. To study myoepithelial differentiation, we employed mSG-PAC1 murine salivary gland epithelial cells. We demonstrate that mSG-PAC1 spheroids exhibit phenotypic plasticity between pro-acinar and myoepithelial cell fates. Increased expression of pro-acinar/acinar or myoepithelial RNAs was identified from spheroids cultured under different media conditions by microarray followed by gene-set enrichment analysis. Spheroids cultured with different medium components expressed proteins typical of either acinar or myoepithelial cells, as detected by immunocytochemistry. We demonstrate that the pattern of TAZ expression in the epithelial compartment of the differentiating murine salivary gland correlates with the expression of the myoepithelial marker alpha-SMA, as is the case for TAZ expression in mSG-PAC1 spheroids. Our analysis also indicates that YAP/TAZ target genes are upregulated together with myoepithelial markers. Importantly, siRNA targeting of TAZ expression in mSG-PAC1 spheroids diminished the expression of myoepithelial markers. Our results in this in vitro cell model implicate TAZ signaling in myoepithelial differentiation.