Transcriptional profiling reveals gland-specific differential expression in the three major salivary glands of the adult mouse

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.

Nkx2.3 transcription factor is a key regulator of mucous cell identity in salivary glands

Saliva is vital to oral health, fulfilling multiple functions in the oral cavity. Three pairs of major salivary glands and hundreds of minor salivary glands contribute to saliva production. The secretory acinar cells within these glands include two distinct populations. Serous acinar cells secrete a watery saliva containing enzymes, while mucous acinar cells secrete a more viscous fluid containing highly glycosylated mucins. Despite their shared developmental origins, the parotid gland (PG) is comprised of only serous acinar cells, while the sublingual gland (SLG) contains predominantly mucous acinar cells. The instructive signals that govern the identity of serous versus mucous acinar cell phenotypes are not yet known. The homeobox transcription factor Nkx2.3 is uniquely expressed in the SLG. Disruption of the Nkx2.3 gene was reported to delay the maturation of SLG mucous acinar cells. To examine whether Nkx2.3 plays a role in directing the mucous cell phenotype, we analyzed SLG from Nkx2.3-/- mice using RNAseq, immunostaining and proteomic analysis of saliva. Our results indicate that Nkx2.3, most likely in concert with other transcription factors uniquely expressed in the SLG, is a key regulator of the molecular program that specifies the identity of mucous acinar cells.

Emergence of saliva protein genes in the secretory calcium-binding phosphoprotein (SCPP) locus and accelerated evolution in primates

Genes within the secretory calcium-binding phosphoprotein (SCPP) family evolved in conjunction with major evolutionary milestones: the formation of a calcified skeleton in vertebrates, the emergence of tooth enamel in fish, and the introduction of lactation in mammals. The SCPP gene family also contains genes expressed primarily and abundantly in human saliva. Here, we explored the evolution of the saliva-related SCPP genes by harnessing currently available genomic and transcriptomic resources. Our findings provide insights into the expansion and diversification of SCPP genes, notably identifying previously undocumented convergent gene duplications. In primate genomes, we found additional duplication and diversification events that affected genes coding for proteins secreted in saliva. These saliva-related SCPP genes exhibit signatures of positive selection in the primate lineage while the other genes in the same locus remain conserved. We found that regulatory shifts and gene turnover events facilitated the accelerated gain of salivary expression. Collectively, our results position the SCPP gene family as a hotbed of evolutionary innovation, suggesting the potential role of dietary and pathogenic pressures in the adaptive diversification of the saliva composition in primates, including humans.

Tlx1 regulates acinar and duct development in mouse salivary glands

Tlx1 encodes a transcription factor expressed in several craniofacial structures of developing mice. The role of Tlx1 in salivary gland development was examined using morphological and immunohistochemical analyses of Tlx1 null mice. Tlx1 is expressed in submandibular and sublingual glands but not parotid glands of neonatal and adult male and female C57Bl/6J (Tlx1+/+ ) mice. TLX1 protein was localized to the nuclei of terminal tubule cells, developing duct cells and mesenchymal cells in neonatal submandibular and sublingual glands, and to nuclei of duct cells and connective tissue cells in adult glands. Occasionally, TLX1 was observed in nuclei of epithelial cells in or adjacent to the acini. Submandibular glands were smaller and sublingual glands were larger in size in mutant mice (Tlx1-/- ) compared to wild-type mice. Differentiation of terminal tubule and proacinar cells of neonatal Tlx1-/- submandibular glands was abnormal; expression of their characteristic products, submandibular gland protein C and parotid secretory protein, respectively, was reduced. At 3 weeks postnatally, terminal tubule cells at the acinar-intercalated duct junction were poorly developed or absent in Tlx1-/- mice. Granular convoluted ducts in adult mutant mice were decreased, and epidermal growth factor and nerve growth factor expression were reduced. Along with normal acinar cell proteins, adult acinar cells of Tlx1-/- mice continued to express neonatal proteins and expressed parotid proteins not normally present in submandibular glands. Sublingual gland mucous acinar and serous demilune cell differentiation were altered. Tlx1 is necessary for proper differentiation of submandibular and sublingual gland acinar cells, and granular convoluted ducts. The mechanism(s) underlying Tlx1 regulation of salivary gland development and differentiation remains unknown.

Claudin-10 Decrease in the Submandibular Gland Contributes to Xerostomia

Tight junction proteins play a crucial role in paracellular transport in salivary gland epithelia. It is clear that severe xerostomia in patients with HELIX syndrome is caused by mutations in the claudin-10 gene. However, little is known about the expression pattern and role of claudin-10 in saliva secretion in physical and disease conditions. In the present study, we found that only claudin-10b transcript was expressed in human and mouse submandibular gland (SMG) tissues, and claudin-10 protein was dominantly distributed at the apicolateral membranes of acini in human, rat, and mouse SMGs. Overexpression of claudin-10 significantly reduced transepithelial electrical resistance and increased paracellular transport of dextran and Na+ in SMG-C6 cells. In C57BL/6 mice, pilocarpine stimulation promoted secretion and cation concentration in saliva in a dose-dependent increase. Assembly of claudin-10 to the most apicolateral portions in acini of SMGs was observed in the lower pilocarpine (1 mg/kg)-treated group, and this phenomenon was much obvious in the higher pilocarpine (10 mg/kg)-treated group. Furthermore, 7-, 14-, and 21-wk-old nonobese diabetic (NOD) and BALB/c mice were used to mimic the progression of hyposalivation in Sjögren syndrome. Intensity of claudin-10 protein was obviously lower in SMGs of 14- and 21-wk-old NOD mice compared with that of age-matched BALB/c mice. In the cultured mouse SMG tissues, interferon-γ (IFN-γ) downregulated claudin-10 expression. In claudin-10-overexpressed SMG-C6 cells, paracellular permeability was decreased. Furthermore, IFN-γ stimulation increased p-STAT1 level, whereas pretreatment with JAK/STAT1 antagonist significantly alleviated the IFN-γ-induced claudin-10 downregulation. These results indicate that claudin-10 functions as a pore-forming component in acinar epithelia of SMGs, assembly of claudin-10 is required for saliva secretion, and downregulation of claudin-10 induces hyposecretion. These findings may provide new clues to novel therapeutic targets on hyposalivation.

The role and mechanism of tight junctions in the regulation of salivary gland secretion

Tight junctions (TJs) are cell-cell interactions that localize at the most apical portion of epithelial/endothelial cells. One of the predominant functions of TJs is to regulate material transport through paracellular pathway, which serves as a selective barrier. In recent years, the expression and function of TJs in salivary glands has attracted great interest. The characteristics of multiple salivary gland TJ proteins have been identified. During salivation, the activation of muscarinic acetylcholine receptor and transient receptor potential vanilloid subtype 1, as well as other stimuli, promote the opening of acinar TJs by inducing internalization of TJs, thereby contributing to increased paracellular permeability. Besides, endothelial TJs are also redistributed with leakage of blood vessels in cholinergic-stimulated submandibular glands. Furthermore, under pathological conditions, such as Sjögren's syndrome, diabetes mellitus, immunoglobulin G4-related sialadenitis, and autotransplantation, the integrity and barrier function of TJ complex are impaired and may contribute to hyposalivation. Moreover, in submandibular glands of Sjögren's syndrome mouse model and patients, the endothelial barrier is disrupted and involved in hyposecretion and lymphocytic infiltration. These findings enrich our understanding of the secretory mechanisms that link the importance of epithelial and endothelial TJ functions to salivation under both physiological and pathophysiological conditions.

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.

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.

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.

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.

A mechanism of gene evolution generating mucin function

How novel gene functions evolve is a fundamental question in biology. Mucin proteins, a functionally but not evolutionarily defined group of proteins, allow the study of convergent evolution of gene function. By analyzing the genomic variation of mucins across a wide range of mammalian genomes, we propose that exonic repeats and their copy number variation contribute substantially to the de novo evolution of new gene functions. By integrating bioinformatic, phylogenetic, proteomic, and immunohistochemical approaches, we identified 15 undescribed instances of evolutionary convergence, where novel mucins originated by gaining densely O-glycosylated exonic repeat domains. Our results suggest that secreted proteins rich in proline are natural precursors for acquiring mucin function. Our findings have broad implications for understanding the role of exonic repeats in the parallel evolution of new gene functions, especially those involving protein glycosylation.

Salivary gland organoid culture maintains distinct glandular properties of murine and human major salivary glands

Salivary glands that produce and secrete saliva, which is essential for lubrication, digestion, immunity, and oral homeostasis, consist of diverse cells. The long-term maintenance of diverse salivary gland cells in organoids remains problematic. Here, we establish long-term murine and human salivary gland organoid cultures. Murine and human salivary gland organoids express gland-specific genes and proteins of acinar, myoepithelial, and duct cells, and exhibit gland functions when stimulated with neurotransmitters. Furthermore, human salivary gland organoids are established from isolated basal or luminal cells, retaining their characteristics. Single-cell RNA sequencing also indicates that human salivary gland organoids contain heterogeneous cell types and replicate glandular diversity. Our protocol also enables the generation of tumoroid cultures from benign and malignant salivary gland tumor types, in which tumor-specific gene signatures are well-conserved. In this study, we provide an experimental platform for the exploration of precision medicine in the era of tissue regeneration and anticancer treatment.

BPIFA2 as a Novel Early Biomarker to Identify Fatal Radiation Injury After Radiation Exposure

Background

Current dosimeters cannot cope with the two tasks of medical rescue in the early stage of nuclear accident, the accurate determination of radiation exposure and the identification of patients with fatal radiation injury. As radiation can cause alterations in serum components, it is feasible to develop biomarkers for radiation injury from serum. This study aims to investigate whether serum BPIFA2 could be used as a potential biomarker of predicting fatal radiation injury in the early stage after nuclear accident.

Methods

A rabbit anti-mouse BPIFA2 polyclonal antibody was prepared to detect the expression of BPIFA2. C57BL/6J female mice were exposed to total body radiation (TBI) at different dose and Partial body radiation (PBI) at lethal dose to detect the dynamic changes of BPIFA2 in serum at different time points after irradiation by Western blot assay.

Results

BPIFA2 in mice serum were significantly increased at 1–12 h post-irradiation at .5–10 Gy, and increased again significantly at 3 d after 10 Gy irradiation with associated with mortality closely. It also increased rapidly after PBI and was closely related to injury degree, regardless whether the salivary glands were irradiated.

Conclusions

The increase of serum BPIFA2 is a novel early biomarker not only for identifying radiation exposure, but also for fatal radiation injury playing a vital role in rational use of medical resources, and greater efficiency of medical treatment to minimize casualties.

Altered Rheological Properties of Saliva with Aging in Mouse Sublingual Gland

Mucin in saliva plays a critical role in the hydration and lubrication of the oral mucosa by retaining water molecules, and its impaired function may be associated with hyposalivation-independent xerostomia. Age-dependent effects on salivary gland function and rheological properties of secreted saliva are not fully understood as aging is a complex and multifactorial process. We aimed to evaluate age-related changes in the rheological properties of saliva and elucidate the underlying mechanism. We performed ex vivo submandibular gland (SMG) and sublingual gland (SLG) perfusion experiments to collect saliva from isolated glands of young (12 wk old) and aged (27 mo old) female C57BL/6J mice and investigate the rheological properties by determining the spinnbarkeit (viscoelasticity). While fluid secretion was comparable in SMG and SLG of both mice, spinnbarkeit showed a significant decrease in SLG saliva of aged mice than that of young mice. There were no significant differences in GalNAc concentration between young and aged SLG saliva. Liquid chromatography/tandem mass spectrometry analysis of SLG saliva revealed that (Hex)₁ (HexNAc)₁ (NeuAc)₁ at m/z 793.31 was the most abundant O-glycan structure in SLG saliva commonly detected in both mice. Lectin staining of salivary gland tissue showed that SLG stained strongly with Maackia amurensis lectin II (MAL II) while Sambucus nigra agglutinin (SNA) stained little, if any, SLG. The messenger RNA expression of St3gal1 that encodes an α-2,3 sialic acid sialyltransferase SIAT4-A showed a decrease in SLG of aged mice, confirmed by a Western blot analysis. Lectin blot analysis in SLG saliva revealed that the relative signal intensity detected by MAL II was significantly lower in aged SLG. Our results suggest that spinnbarkeit decreases in SLG of aging mice due to downregulation of sialic acid linked to α-2,3 sialic acid sialyltransferase expression.

Differential regulation of transposable elements (TEs) during the murine submandibular gland development

The submandibular gland (SG) is a relatively simple organ formed by three cell types: acinar, myoepithelial, and an intricate network of duct-forming epithelial cells, that together fulfills several physiological functions from assisting food digestion to acting as an immune barrier against pathogens. Successful SG organogenesis is the product of highly controlled and orchestrated genetic and transcriptional programs. Mounting evidence links Transposable Elements (TEs), originally thought to be selfish genetic elements, to different aspects of gene regulation in mammalian development and disease. To our knowledge, the role of TEs during murine SG organogenesis has not been studied. Using novel bioinformatic tools and publicly available RNA-Seq datasets, our results indicate that a significant number of genic and intergenic TEs are differentially expressed during the SG development. Furthermore, changes in expression of specific TEs correlated with that of genes involved in cellular division and differentiation, critical aspects for SG maturation. Altogether, we propose that TEs modulate gene networks that operate during SG development.

A Global Vista of the Epigenomic State of the Mouse Submandibular Gland

The parotid, submandibular, and sublingual glands represent a trio of oral secretory glands whose primary function is to produce saliva, facilitate digestion of food, provide protection against microbes, and maintain oral health. While recent studies have begun to shed light on the global gene expression patterns and profiles of salivary glands, particularly those of mice, relatively little is known about the location and identity of transcriptional control elements. Here we have established the epigenomic landscape of the mouse submandibular salivary gland (SMG) by performing chromatin immunoprecipitation sequencing experiments for 4 key histone marks. Our analysis of the comprehensive SMG data sets and comparisons with those from other adult organs have identified critical enhancers and super-enhancers of the mouse SMG. By further integrating these findings with complementary RNA-sequencing based gene expression data, we have unearthed a number of molecular regulators such as members of the Fox family of transcription factors that are enriched and likely to be functionally relevant for SMG biology. Overall, our studies provide a powerful atlas of cis-regulatory elements that can be leveraged for better understanding the transcriptional control mechanisms of the mouse SMG, discovery of novel genetic switches, and modulating tissue-specific gene expression in a targeted fashion.

An ancient, conserved gene regulatory network led to the rise of oral venom systems

Although oral venom systems are ecologically important characters, how they originated is still unclear. In this study, we show that oral venom systems likely originated from a gene regulatory network conserved across amniotes. This network, which we term the “metavenom network,” comprises over 3,000 housekeeping genes coexpressed with venom and play a role in protein folding and modification. Comparative transcriptomics revealed that the network is conserved between venom glands of snakes and salivary glands of mammals. This suggests that while these tissues have evolved different functions, they share a common regulatory core, that persisted since their common ancestor. We propose several evolutionary mechanisms that can utilize this common regulatory core to give rise to venomous animals from their nonvenomous ancestors.

Oral venom systems evolved multiple times in numerous vertebrates enabling the exploitation of unique predatory niches. Yet how and when they evolved remains poorly understood. Up to now, most research on venom evolution has focused strictly on the toxins. However, using toxins present in modern day animals to trace the origin of the venom system is difficult, since they tend to evolve rapidly, show complex patterns of expression, and were incorporated into the venom arsenal relatively recently. Here we focus on gene regulatory networks associated with the production of toxins in snakes, rather than the toxins themselves. We found that overall venom gland gene expression was surprisingly well conserved when compared to salivary glands of other amniotes. We characterized the “metavenom network,” a network of ∼3,000 nonsecreted housekeeping genes that are strongly coexpressed with the toxins, and are primarily involved in protein folding and modification. Conserved across amniotes, this network was coopted for venom evolution by exaptation of existing members and the recruitment of new toxin genes. For instance, starting from this common molecular foundation, Heloderma lizards, shrews, and solenodon, evolved venoms in parallel by overexpression of kallikreins, which were common in ancestral saliva and induce vasodilation when injected, causing circulatory shock. Derived venoms, such as those of snakes, incorporated novel toxins, though still rely on hypotension for prey immobilization. These similarities suggest repeated cooption of shared molecular machinery for the evolution of oral venom in mammals and reptiles, blurring the line between truly venomous animals and their ancestors.

Development of a functional salivary gland tissue chip with potential for high-content drug screening

Radiation therapy for head and neck cancers causes salivary gland dysfunction leading to permanent xerostomia. Limited progress in the discovery of new therapeutic strategies is attributed to the lack of in vitro models that mimic salivary gland function and allow high-throughput drug screening. We address this limitation by combining engineered extracellular matrices with microbubble (MB) array technology to develop functional tissue mimetics for mouse and human salivary glands. We demonstrate that mouse and human salivary tissues encapsulated within matrix metalloproteinase-degradable poly(ethylene glycol) hydrogels formed in MB arrays are viable, express key salivary gland markers, and exhibit polarized localization of functional proteins. The salivary gland mimetics (SGm) respond to calcium signaling agonists and secrete salivary proteins. SGm were then used to evaluate radiosensitivity and mitigation of radiation damage using a radioprotective compound. Altogether, SGm exhibit phenotypic and functional parameters of salivary glands, and provide an enabling technology for high-content/throughput drug testing.

The salivary protein BPIFA2 differentially regulates sodium preference and blood pressure in male and female mice

BPIFA2 (PSP, SPLUNC2, C20orf70) is a major salivary protein of uncertain physiological function. BPIFA2 is downregulated in salivary glands of spontaneously hypertensive rats, pointing to a role in blood pressure regulation. This study used a novel Bpifa2 knockout mouse model to test the role of BPIFA2 in sodium preference and blood pressure. Blood pressure did not differ between wild-type male and female mice but was significantly lower in male knockout mice compared to male wild-type mice. In contrast, blood pressure was increased in female knockout mice compared to female wild-type mice. Female wild-type mice showed a significant preference for 0.9% saline compared to male mice. This difference was reduced in the knockout mice. BPIFA2 is an LPS-binding protein but it remains to be determined if the reported effects are mediated by the LPS-binding activity of BPIFA2.

The parotid secretory protein BPIFA2 is a salivary surfactant that affects lipopolysaccharide action

NEW FINDINGS

What is the central question of this study? The salivary protein BPIFA2 binds lipopolysaccharide, but its physiological function is not known. This study uses a new knockout mouse model to explore the physiological role of BPIFA2 in the oral cavity and systemic physiology. What is the main finding and its importance? BPIFA2 is a crucial surfactant in mouse saliva. In its absence, saliva exhibits the surface tension of water. Depletion of BPIFA2 affects salivary and ingested lipopolysaccharide and leads to systemic sequelae that include increased insulin secretion and metabolomic changes. These results suggest that the lipopolysaccharide-binding activity of BPIFA2 affects the activity of ingested lipopolysaccharide in the intestine and that BPIFA2 depletion causes mild metabolic endotoxaemia.

ABSTRACT

Saliva plays important roles in the mastication, swallowing and digestion of food, speech and lubrication of the oral mucosa, antimicrobial and anti-inflammatory activities, and the control of body temperature in grooming animals. The salivary protein BPIFA [BPI fold containing family A member 2; former names: parotid secretory protein (PSP), SPLUN2 and C20orf70] is related to lipid-binding and lipopolysaccharide (LPS)-binding proteins expressed in the mucosa. Indeed, BPIFA2 binds LPS, but the physiological role of BPIFA2 remains to be determined. To address this question, Bpifa2 knockout (Bpifa2^{tm1(KOMP)Vlcg} ) (KO) mice were phenotyped, with emphasis on the saliva and salivary glands. Stimulated whole saliva collected from KO mice was less able to spread on a hydrophobic surface than wild-type saliva, and the surface tension of KO saliva was close to that of water. These data suggest that BPIFA2 is a salivary surfactant that is mainly responsible for the low surface tension of mouse saliva. The reduced surfactant activity of KO saliva did not affect consumption of dry food or grooming, but saliva from KO mice contained less LPS than wild-type saliva. Indeed, mice lacking BPIFA2 responded to ingested LPS with an increased stool frequency, suggesting that BPIFA2 plays a role in the solubilization and activity of ingested LPS. Consistent with these findings, BPIFA2-depleted mice also showed increased insulin secretion and metabolomic changes that were consistent with a mild endotoxaemia. These results support the distal physiological function of a salivary protein and reinforce the connection between oral biology and systemic disease.

Ano6 disruption impairs acinar cell regulatory volume decrease and protein secretion in murine submandibular salivary glands

The widely expressed Anoctamin 6 (Ano6) supports different Ca²⁺ -dependent functions, but little is known about its role in salivary glands. Mouse submandibular gland (SMG) acinar cells exhibited a robust regulatory volume decrease (RVD) following cell swelling that was reduced approximately 70% in Ano6^-/- mice. Ca²⁺ -free conditions nearly eliminated the RVD response suggesting that Ano6 is an obligatory component of the cell volume-activated, Ca²⁺ -dependent RVD pathway in salivary gland acinar cells. Ex vivo agonist-stimulated secretion of water and ions was unaffected by Ano6 disruption under both isotonic and hypotonic conditions suggesting that Ano6 does not play a major role in fluid and electrolyte secretion. In contrast, the total amount of β-adrenergic-dependent protein secretion by the SMG was significantly reduced in Ano6^-/- mice. Closer inspection of these latter results revealed that protein secretion was affected only in the female SMG by Ano6 disruption. These results indicate that Ano6 modulates the RVD response and protein secretion by salivary gland acinar cells.

Transcriptomic and Single-Cell Analysis of the Murine Parotid Gland

The salivary complex of mammals consists of 3 major pairs of glands: the parotid, submandibular, and sublingual glands. While the 3 glands share similar functional properties, such as saliva secretion, their differences are largely based on the types of secretions they produce. While recent studies have begun to shed light on the underlying molecular differences among the glands, few have examined the global transcriptional repertoire over various stages of gland maturation. To better elucidate the molecular nature of the parotid gland, we have performed RNA sequencing to generate comprehensive and global gene expression profiles of this gland at different stages of maturation. Our transcriptomic characterization and hierarchical clustering analysis with adult organ RNA sequencing data sets has identified a number of molecular players and pathways that are relevant for parotid gland biology. Moreover, our detailed analysis has revealed a unique parotid gland-specific gene signature that may represent important players that could impart parotid gland-specific biological properties. To complement our transcriptomic studies, we have performed single-cell RNA sequencing to map the transcriptomes of parotid epithelial cells. Interrogation of the single-cell transcriptomes revealed the degree of molecular and cellular heterogeneity of the various epithelial cell types within the parotid gland. Moreover, we uncovered a mixed-lineage population of cells that may reflect molecular priming of differentiation potentials. Overall our comprehensive studies provide a powerful tool for the discovery of novel molecular players important in parotid gland biology.

Slc4a11 disruption causes duct cell loss and impairs NaCl reabsorption in female mouse submandibular glands

Slc4a11, a member of the Slc4 HCO3- transporter family, has a wide tissue distribution. In mouse salivary glands, the expression of Slc4a11 mRNA was more than eightfold greater than the other nine members of the Slc4 gene family. The Slc4a11 protein displayed a diffuse subcellular distribution in both the acinar and duct cells of mouse submandibular glands (SMG). Slc4a11 disruption induced a significant increase in the Na+ and Cl- concentrations of stimulated SMG saliva, whereas it did not affect the fluid secretion rate in response to either β-adrenergic or cholinergic receptor stimulation. Heterologous expressed mouse Slc4a11 acted as a H+ /OH- transporter that was uncoupled of Na+ or Cl- movement, and this activity was blocked by ethyl-isopropyl amiloride (EIPA) but not 4,4'-Diisothiocyanato-2,2'-stilbenedisulfonic acid (DIDS). Slc4a11 disruption revealed that Slc4a11 does not play a major role in intracellular pH regulation in mouse salivary gland cells. In contrast, NaCl reabsorption was impaired in the SMG saliva of female compared to male Slc4a11 null mice, which correlated with the loss of duct cells and a decrease in expression of the duct-cell-specific transcription factor Ascl3. Together, our results suggest that Slc4a11 expression regulates the number of ducts cells in the mouse SMG and consequently NaCl reabsorption.

Slc26a6 is an apical membrane anion exchanger that drives HCO3--dependent fluid secretion in murine pancreatic acinar cells

The nonselective anion exchanger Slc26a6, also known as putative anion transporter 1 and chloride/formate exchanger, is thought to play a major role in HCO 3 - transport in exocrine glands. In this study, Slc26a6 null mice were used to explore the function of Slc26a6 in the exocrine pancreas. Slc26a6 primarily localized to the apical membrane of pancreatic exocrine acinar cells. The volume of stimulated juice secretion by the ex vivo pancreas was significantly reduced ~35% in Slc26a6-/- mice, but no changes occurred in the gross structure or gland weights of Slc26a6 null mice. The secretion of pancreatic juice by Slc26a6+/+ mice was dependent on HCO 3 - while, in contrast, fluid secretion by Slc26a6-/- mice was independent of HCO 3 - , suggesting that Slc26a6 mediates the HCO 3 - -dependent component of fluid secretion. Consistent with these observations, disruption of Slc26a6 also significantly reduced HCO 3 - secretion by the pancreas ~35%. Taken together, these results demonstrate that the apical Slc26a6 anion exchanger in acinar cells is involved in HCO 3 - -dependent fluid secretion but that another major HCO 3 - -independent pathway is the primary driver of the fluid secretion process in the mouse pancreas.

Single-Cell RNA-seq Identifies Cell Diversity in Embryonic Salivary Glands

Branching organs, including the salivary and mammary glands, lung, and kidney, arise as epithelial buds that are morphologically very similar. However, the mesenchyme is known to guide epithelial morphogenesis and to help govern cell fate and eventual organ specificity. We performed single-cell transcriptome analyses of 14,441 cells from embryonic day 12 submandibular and parotid salivary glands to characterize their molecular identities during bud initiation. The mesenchymal cells were considerably more heterogeneous by clustering analysis than the epithelial cells. Nonetheless, distinct clusters were evident among even the epithelial cells, where unique molecular markers separated presumptive bud and duct cells. Mesenchymal cells formed separate, well-defined clusters specific to each gland. Neuronal and muscle cells of the 2 glands in particular showed different markers and localization patterns. Several gland-specific genes were characteristic of different rhombomeres. A muscle cluster was prominent in the parotid, which was not myoepithelial or vascular smooth muscle. Instead, the muscle cluster expressed genes that mediate skeletal muscle differentiation and function. Striated muscle was indeed found later in development surrounding the parotid gland. Distinct spatial localization patterns of neuronal and muscle cells in embryonic stages appear to foreshadow later differences in adult organ function. These findings demonstrate that the establishment of transcriptional identities emerges early in development, primarily in the mesenchyme of developing salivary glands. We present the first comprehensive description of molecular signatures that define specific cellular landmarks for the bud initiation stage, when the neural crest-derived ectomesenchyme predominates in the salivary mesenchyme that immediately surrounds the budding epithelium. We also provide the first transcriptome data for the largely understudied embryonic parotid gland as compared with the submandibular gland, focusing on the mesenchymal cell populations.

The apical Na+ -HCO3 - cotransporter Slc4a7 (NBCn1) does not contribute to bicarbonate transport by mouse salivary gland ducts

The HCO₃ ^- secretion mechanism in salivary glands is unclear but is thought to rely on the co-ordinated activity of multiple ion transport proteins including members of the Slc4 family of bicarbonate transporters. Slc4a7 was immunolocalized to the apical membrane of mouse submandibular duct cells. In contrast, Slc4a7 was not detected in acinar cells, and correspondingly, Slc4a7 disruption did not affect fluid secretion in response to cholinergic or β-adrenergic stimulation in the submandibular gland (SMG). Much of the Na ⁺ -dependent intracellular pH (pH _i ) regulation in SMG duct cells was insensitive to 4,4'-diisothiocyano-2,2'-stilbenedisulfonic acid, S0859, and to the removal of extracellular HCO ₃ ^- . Consistent with these latter observations, the Slc4a7 null mutation had no impact on HCO ₃ ^- secretion nor on pH _i regulation in duct cells. Taken together, our results revealed that Slc4a7 targets to the apical membrane of mouse SMG duct cells where it contributes little if any to pH _i regulation or stimulated HCO ₃ ^- secretion.

Sexual dimorphisms in the transcriptomes of murine salivary glands

Transcriptional profiling identified 933 sexually dimorphic genes out of the 14 371 protein‐coding genes expressed in the three major murine salivary glands: parotid, sublingual, and submandibular. Most (89%) sex‐specific genes were enriched in a single gland, while only 0.5% of the sexually dimorphic genes were enriched in all glands. The sublingual gland displayed a strong male sex bias (94% of sex‐enriched genes), while a sex preference was not obvious in the parotid or submandibular glands. A subset of transcription factor genes was correlated with the expression of gland‐specific, sex‐enriched genes. Higher expression of Cftr chloride and Scnn1 sodium channels in the male submandibular correlated with greater NaCl reabsorption. In conclusion, adult salivary glands display sex‐ and gland‐specific differences in gene expression that reflect their unique functional properties.

Physiological cAMP-elevating secretagogues differentially regulate fluid and protein secretions in mouse submandibular and sublingual glands

The mechanisms underlying the functional differences in sympathetic and parasympathetic regulation of the major salivary glands have received little attention. The acute effects of parasympathetic muscarinic (carbachol)-dependent and combined parasympathetic-dependent plus cAMP-dependent pathways on fluid secretion rates, ion composition, and protein content were assessed using a newly developed ex vivo preparation that allows the simultaneous perfusion of the mouse submandibular (SMGs) and sublingual glands (SLGs). Our results confirm that the muscarinic-dependent pathway accounts for the bulk of salivation in SMGs and SLGs, whereas costimulation with a cAMP-increasing agent (forskolin, isoproterenol, or vasoactive intestinal peptide) did not increase the flow rate. Costimulation with carbachol plus the β-adrenergic agonist isoproterenol decreased the concentration of NaCl and produced a substantial increase in the protein and Ca²⁺ content of SMG but not SLG saliva, consistent with a sparse sympathetic innervation of the SLGs. On the other hand, forskolin, which bypasses receptors to increase intracellular cAMP by directly activating the enzyme adenylate cyclase, enhanced the secretion of protein and Ca²⁺ by both the SMGs and SLGs. In contrast, isoproterenol and vasoactive intestinal peptide specifically stimulated protein secretion in SMG and SLG salivas, respectively. In summary, cAMP-dependent signaling does not play a major role in the stimulation of fluid secretion in SMGs and SLGs, whereas each cAMP-increasing agonist behaves differently in a gland-specific manner suggesting differential expression of G protein-coupled receptors in the epithelial cells of SMGs and SLGs.

Integrated Epigenetic Mapping of Human and Mouse Salivary Gene Regulation

Significant effort has been applied to identify the genome-wide gene expression profiles associated with salivary gland development and pathophysiology. However, relatively little is known about the regulators that control salivary gland gene expression. We integrated data from DNase1 digital genomic footprinting, RNA-seq, and gene expression microarrays to comprehensively characterize the cis- and trans-regulatory components controlling gene expression of the healthy submandibular salivary gland. Analysis of 32 human tissues and 87 mouse tissues was performed to identify the highly expressed and tissue-enriched transcription factors driving salivary gland gene expression. Following RNA analysis, protein expression levels and subcellular localization of 39 salivary transcription factors were confirmed by immunohistochemistry. These expression analyses revealed that the salivary gland highly expresses transcription factors associated with endoplasmic reticulum stress, human T-cell lymphotrophic virus 1 expression, and Epstein-Barr virus reactivation. DNase1 digital genomic footprinting to a depth of 333,426,353 reads was performed and utilized to generate a salivary gland gene regulatory network describing the genome-wide chromatin accessibility and transcription factor binding of the salivary gland at a single-nucleotide resolution. Analysis of the DNase1 gene regulatory network identified dense interconnectivity among PLAG1, MYB, and 13 other transcription factors associated with balanced chromosomal translocations and salivary gland tumors. Collectively, these analyses provide a comprehensive atlas of the cis- and trans-regulators of the salivary gland and highlight known aberrantly regulated pathways of diseases affecting the salivary glands.

Defective NaCl Reabsorption in Salivary Glands of Eda-Null X-LHED Mice

Mutations in the ectodysplasin A gene ( EDA) cause X-LHED (X-linked hypohidrotic ectodermal dysplasia), the most common human form of ectodermal dysplasia. Defective EDA signaling is linked to hypoplastic development of epithelial tissues, resulting in hypotrichosis, hypodontia, hypohidrosis, and xerostomia. The primary objective of the present study was to better understand the salivary gland dysfunction associated with ectodermal dysplasia using the analogous murine disorder. The salivary flow rate and ion composition of the 3 major salivary glands were determined in adult Eda-deficient Tabby hemizygous male (Ta/Y) and heterozygous female (Ta/X) mice. Submandibular and sublingual glands of Eda-mutant mice were smaller than wild-type littermates, while parotid gland weight was not significantly altered. Fluid secretion by the 3 major salivary glands was essentially unchanged, but the decrease in submandibular gland size was associated with a dramatic loss of ducts in Ta/Y and Ta/X mice. Reabsorption of Na+ and Cl-, previously linked in salivary glands to Scnn1 Na+ channels and Cftr Cl- channels, respectively, was markedly reduced at high flow rates in the ex vivo submandibular glands of Ta/Y mice (~60%) and, to a lesser extent, Ta/X mice (Na+ by 14%). Consistent with decreased Na+ reabsorption in Ta/Y mice, quantitative polymerase chain reaction analysis detected decreased mRNA expression for Scnn1b and Scnn1g, genes encoding the β and γ subunits, respectively. Moreover, the Na+ channel blocker amiloride significantly inhibited Na+ and Cl- reabsorption by wild-type male submandibular glands to levels comparable to those observed in Ta/Y mice. In summary, fluid secretion was intact in the salivary glands of Eda-deficient mice but displayed marked Na+ and Cl- reabsorption defects that correlated with the loss of duct cells and decreased Scnn1 Na+ channel expression. These results provide a likely mechanism for the elevated NaCl concentration observed in the saliva of affected male and female patients with X-LHED.

The apical anion exchanger Slc26a6 promotes oxalate secretion by murine submandibular gland acinar cells

The solute carrier family 26 (SLC26) gene family encodes at least 10 different anion exchangers. SLC26 member 6 (SLC26A6 or CFEX/PAT-1) and the cystic fibrosis transmembrane conductance regulator (CFTR) co-localize to the apical membrane of pancreatic duct cells, where they act in concert to drive HCO3- and fluid secretion. In contrast, in the small intestine, SLC26A6 serves as the major pathway for oxalate secretion. However, little is known about the function of Slc26a6 in murine salivary glands. Here, RNA sequencing-based transcriptional profiling and Western blots revealed that Slc26a6 is highly expressed in mouse submandibular and sublingual salivary glands. Slc26a6 localized to the apical membrane of salivary gland acinar cells with no detectable immunostaining in the ducts. CHO-K1 cells transfected with mouse Slc26a6 exchanged Cl- for oxalate and HCO3-, whereas two other anion exchangers known to be expressed in salivary gland acinar cells, Slc4a4 and Slc4a9, mediated little, if any, Cl-/oxalate exchange. Of note, both Cl-/oxalate exchange and Cl-/HCO3- exchange were significantly reduced in acinar cells isolated from the submandibular glands of Slc26a6-/- mice. Oxalate secretion in submandibular saliva also decreased significantly in Slc26a6-/- mice, but HCO3- secretion was unaffected. Taken together, our findings indicate that Slc26a6 is located at the apical membrane of salivary gland acinar cells, where it mediates Cl-/oxalate exchange and plays a critical role in the secretion of oxalate into saliva.