Transcriptomic and Single-Cell Analysis Reveals Regulatory Networks and Cellular Heterogeneity in Mouse Primary Sjögren's Syndrome Salivary Glands

Single-Cell Sequencing: High-Resolution Analysis of Cellular Heterogeneity in Autoimmune Diseases

Autoimmune diseases (AIDs) are complex in etiology and diverse in classification but clinically show similar symptoms such as joint pain and skin problems. As a result, the diagnosis is challenging, and usually, only broad treatments can be available. Consequently, the clinical responses in patients with different types of AIDs are unsatisfactory. Therefore, it is necessary to conduct more research to figure out the pathogenesis and therapeutic targets of AIDs. This requires research technologies with strong extraction and prediction capabilities. Single-cell sequencing technology analyses the genomic, epigenomic, or transcriptomic information at the single-cell level. It can define different cell types and states in greater detail, further revealing the molecular mechanisms that drive disease progression. These advantages enable cell biology research to achieve an unprecedented resolution and scale, bringing a whole new vision to life science research. In recent years, single-cell technology especially single-cell RNA sequencing (scRNA-seq) has been widely used in various disease research. In this paper, we present the innovations and applications of single-cell sequencing in the medical field and focus on the application contributing to the differential diagnosis and precise treatment of AIDs. Despite some limitations, single-cell sequencing has a wide range of applications in AIDs. We finally present a prospect for the development of single-cell sequencing. These ideas may provide some inspiration for subsequent research.

Advances in cellular and molecular pathways of salivary gland damage in Sjögren's syndrome

Sjögren's Syndrome (SS) is an autoimmune disorder characterized by dysfunction of exocrine glands. Primarily affected are the salivary glands, which exhibit the most frequent pathological changes. The pathogenesis involves susceptibility genes, non-genetic factors such as infections, immune cells-including T and B cells, macrophage, dendritic cells, and salivary gland epithelial cells. Inflammatory mediators such as autoantibodies, cytokines, and chemokines also play a critical role. Key signaling pathways activated include IFN, TLR, BAFF/BAFF-R, PI3K/Akt/mTOR, among others. Comprehensive understanding of these mechanisms is crucial for developing targeted therapeutic interventions. Thus, this study explores the cellular and molecular mechanisms underlying SS-related salivary gland damage, aiming to propose novel targeted therapeutic approaches.

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.

CXCL9 may serve as a potential biomarker for primary Sjögren's syndrome with extra-glandular manifestations

BACKGROUND

Primary Sjögren's syndrome (pSS) is an autoimmune condition that causes harm to exocrine glands and also has extra-glandular manifestations (EGM). pSS patients with EGM have a worse prognosis than those with only sicca symptoms. Previous studies have shown that the minor salivary glands (MSG) of pSS patients exhibit a unique profile of cytokines and chemokines compared to healthy controls. However, there is a lack of research comparing pSS with EGM (pSS-EGM) and pSS without EGM (pSS-non-EGM). This study aims to explore potential biomarkers associated with pSS, particularly pSS with EGM.

METHODS

By utilizing RNA sequencing, we conducted an analysis on the gene expression profiles of MSG in 63 patients diagnosed with pSS, as well as 12 non-pSS individuals. Furthermore, we also investigated the MSG of pSS patients, both with and without EGM. Through bioinformatics analysis, we identified genes with differential expression (DEGs) and determined the core hub genes using PPI network. We then analyzed the top 20 DEGs and their correlation with the patients' clinical characteristics, and validated our findings using peripheral blood plasma.

RESULTS

A total of 725 differentially expressed genes (DEGs) were identified in the comparison between pSS and non-pSS groups, and 727 DEGs were observed between pSS-EGM and pSS-non-EGM. It is noteworthy that the expression levels of CXCL9 were higher in both pSS patients and pSS-EGM when compared to the control group. Taking into consideration the significance of the top 20 DEGs in relation to clinical parameters and the central hub genes, we ultimately chose CXCL9. In comparison to the non-pSS group, pSS patients exhibited notably greater expression of the CXCL9 gene in the MSG, as well as higher levels of CXCL9 protein in their plasma (p < 0.001). Furthermore, the expression of the CXCL9 gene and levels of CXCL9 protein were notably higher in pSS patients accompanied by EGM and those with SSA antibodies. Additionally, a correlation was found between the expression of the CXCL9 gene and the EULAR Sjogren's Syndrome Disease Activity Index (ESSDAI), as well as with immunoglobulin G (IgG) levels and erythrocyte sedimentation rate (ESR). Meanwhile, the protein levels of CXCL9 were found to be correlated with IgG levels and ESSDAI.

CONCLUSION

CXCL9 proves to be a valuable biomarker in pSS, specifically due to its strong ability to differentiate between pSS patients with EGM and those without EGM. There is a significant correlation between CXCL9 and various clinical parameters both at the gene and protein level. Therefore, CXCL9 could be a potential target for future treatment of pSS.

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.

Single-Cell RNA-Seq Analysis Identifies Angiotensinogen and Galanin as Unique Molecular Markers of Acinar Cells in Murine Salivary Glands

The submandibular gland (SMG) and sublingual gland (SLG) are two of three major salivary glands in mammals and comprise serous and mucous acinar cells. The two glands share some functional properties, which are largely dependent on the types of acinar cells. In recent years, while ScRNA-seq (single-cell sequencing) with a 10 × platform has been used to explore molecular markers in salivary glands, few studies have examined the acinar heterogeneity and unique molecular markers between SMG and SLG. This study aimed to identify the molecular markers of acinar cells in the SLG and SMG. We performed ScRNA-seq analyses in 4-week-old mice and verified the screened molecular markers using reverse transcription-quantitative real-time PCR, immunohistochemistry, and immunofluorescence. Our results showed prominently heterogeneous acinar cells, although there was great similarity in the cluster composition between the two glands at 4 weeks. Furthermore, we demonstrated that Agt is a specific marker of SMG serous acinar cells, whereas Gal is a specific marker of SLG mucous acinar cells. Trajectory inference revealed that Agt and Gal represent two types of differential acinar cell clusters during late development in adults. Thus, we reveal previously unknown specific markers for salivary acinar cell diversity, which has extensive implications for their further functional research.

Application of omics in Sjögren's syndrome

OBJECTIVE

The pathogenesis, diagnosis, and treatment of Sjögren's syndrome (SS) face many challenges, and there is an urgent need to develop new technologies to improve our understanding of SS.

METHODS

By searching the literature published domestically and internationally in the past 20 years, this artical reviewed the research of various omics techniques in SS.

RESULTS

Omics technology provided valuable insights into the pathogenesis, early diagnosis, condition and efficacy evaluation of SS. It is helpful to reveal the pathogenesis of the disease and explore new treatment schemes, which will open a new era for the study of SS.

CONCLUSION

At present, omics research has made some gratifying achievements, but there are still many uncertainties. Therefore, in the future, we should improve research techniques, standardize the collection of samples, and adopt a combination of multi-omics techniques to jointly study the pathogenesis of SS and provide new schemes for its treatment.

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.

Review of genes potentially related to hyposecretion in male non-obese diabetic (NOD) mice, a Sjögren's syndrome model

BACKGROUND

Sjögren's syndrome (SS) is known to cause dry eyes and mouth due to inflammation of the lacrimal and salivary glands. However, some reports imply that other factors trigger dry eyes and mouth. We previously investigated various factors using RNA-sequencing analysis of lacrimal glands from male non-obese diabetic (NOD) mice, an SS model. In this review, we described (1) the exocrine features of male and female NOD mice, (2) the up- and down-regulated genes in the lacrimal glands of male NOD mice as revealed by our RNA-sequencing data, and (3) comparisons between these genes and data in the Salivary Gland Gene Expression Atlas.

HIGHLIGHTS

Male NOD mice exhibit a steady worsening of lacrimal hyposecretion and dacryoadenitis, whereas females exhibit a complex pathophysiological condition that includes diabetic disease, salivary hyposecretion, and sialadenitis. Ctss, an up-regulated gene, is a potential inducer of lacrimal hyposecretion and is also expressed in salivary glands. Two other up-regulated genes, Ccl5 and Cxcl13, may worsen the inflammation of SS in both the lacrimal and salivary glands. The genes Esp23, Obp1a, and Spc25 were detected as down-regulated, but judging the relationship between these genes and hyposecretion is difficult as only limited information is available. Another down-regulated gene, Arg1, is involved in lacrimal hyposecretion, and it also has the potential to cause salivary hyposecretion in NOD mice.

CONCLUSION

In NOD mice, males may be better than females at evaluating the pathophysiology of SS. Some regulated genes revealed by our RNA-sequencing data might be potential therapeutic targets for SS.

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.

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.

Adipose-derived mesenchymal stem cells promote salivary duct regeneration via a paracrine effect

OBJECTIVES

The self-regeneration of exocrine tissues, including salivary glands, is limited and their regeneration mechanism has not yet been fully elucidated. Here we identify the role of adipose-derived mesenchymal stem cells (AMSCs) in salivary gland regeneration.

METHODS

AMSCs expressing mesenchymal stem cell markers were applied to a submandibular gland injury model and the mechanism of salivary gland repair and regeneration was analyzed.

RESULTS

Transplanted green fluorescent protein (GFP)-labeled AMSCs grew tightly together and promoted ductal regeneration in the regenerative nodule, with slight infiltration of nonspecific immune cells. A comprehensive gene analysis through RNA-sequencing revealed increased expression of bone morphogenetic protein (BMP), transforming growth factor (TGF), and Wnt in AMSC-transplanted regenerative nodules. The factors released from AMSCs scavenge hydrogen peroxidase-induced reactive oxygen species (ROS) through Wnt promoter activity in vitro. Furthermore, AMSC-conditioned medium recovered the growth of the hydrogen peroxidase-damaged primordium of the submandibular gland culture ex vivo.

CONCLUSIONS

These results suggest that AMSC-released factors scavenge ROS and maintain salivary gland repair and regeneration via paracrine effects. Thus, AMSCs could be a practical and applicable tool for use in salivary gland regeneration.

Single-Cell Genomics for Investigating Pathogenesis of Inflammatory Diseases

Recent technical advances have enabled unbiased transcriptomic and epigenetic analysis of each cell, known as "single-cell analysis". Single-cell analysis has a variety of technical approaches to investigate the state of each cell, including mRNA levels (transcriptome), the immune repertoire (immune repertoire analysis), cell surface proteins (surface proteome analysis), chromatin accessibility (epigenome), and accordance with genome variants (eQTLs; expression quantitative trait loci). As an effective tool for investigating robust immune responses in coronavirus disease 2019 (COVID-19), many researchers performed single-cell analysis to capture the diverse, unbiased immune cell activation and differentiation. Despite challenges elucidating the complicated immune microenvironments of chronic inflammatory diseases using existing experimental methods, it is now possible to capture the simultaneous immune features of different cell types across inflamed tissues using various single-cell tools. In this review, we introduce patient-based and experimental mouse model research utilizing single-cell analyses in the field of chronic inflammatory diseases, as well as multi-organ atlas targeting immune cells.

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.

Exploiting the role of T cells in the pathogenesis of Sjögren's syndrome for therapeutic treatment

Sjögrens syndrome (SS) is caused by autoantibodies that attack proprioceptive salivary and lacrimal gland tissues. Damage to the glands leads to dry mouth and eyes and affects multiple systems and organs. In severe cases, SS is life-threatening because it can lead to interstitial lung disease, renal insufficiency, and lymphoma. Histological examination of the labial minor salivary glands of patients with SS reveals focal lymphocyte aggregation of T and B cells. More studies have been conducted on the role of B cells in the pathogenesis of SS, whereas the role of T cells has only recently attracted the attention of researchers. This review focusses on the role of various populations of T cells in the pathogenesis of SS and the progress made in research to therapeutically targeting T cells for the treatment of patients with SS.

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.

Spatial transcriptomics of the lacrimal gland features macrophage activity and epithelium metabolism as key alterations during chronic inflammation

The lacrimal gland (LG) is an exocrine gland that produces the watery part of the tear film that lubricates the ocular surface. Chronic inflammation, such as Sjögren's syndrome (SS), is one of the leading causes of aqueous-deficiency dry eye (ADDE) disease worldwide. In this study we analyzed the chronic inflammation in the LGs of the NOD.B10Sn-H2b/J (NOD.H-2b) mice, a mouse model of SS, utilizing bulk RNAseq and Visium spatial gene expression. With Seurat we performed unsupervised clustering and analyzed the spatial cell distribution and gene expression changes in all cell clusters within the LG sections. Moreover, for the first time, we analyzed and validated specific pathways defined by bulk RNAseq using Visium technology to determine activation of these pathways within the LG sections. This analysis suggests that altered metabolism and the hallmarks of inflammatory responses from both epithelial and immune cells drive inflammation. The most significant pathway enriched in upregulated DEGs was the "TYROBP Causal Network", that has not been described previously in SS. We also noted a significant decrease in lipid metabolism in the LG of the NOD.H-2b mice. Our data suggests that modulation of these pathways can provide a therapeutic strategy to treat ADDE.