Single-cell RNA sequencing reveals PDGFRα+ stromal cell subpopulations that promote proacinar cell differentiation in embryonic salivary gland organoids.

Stromal cells can direct the differentiation of epithelial progenitor cells during organ development. Fibroblast growth factor (FGF) signaling is essential for submandibular salivary gland development. Through stromal fibroblast cells, FGF2 can indirectly regulate proacinar cell differentiation in organoids, but the mechanisms are not understood. We performed single-cell RNA-sequencing and identified multiple stromal cell subsets, including Pdgfra+ stromal subsets expressing both Fgf2 and Fgf10. When combined with epithelial progenitor cells in organoids, magnetic-activated cell-sorted PDGFRα+ cells promoted proacinar cell differentiation similarly to total stroma. Gene expression analysis revealed that FGF2 increased the expression of multiple stromal genes, including Bmp2 and Bmp7. Both BMP2 and BMP7 synergized with FGF2, stimulating proacinar cell differentiation but not branching. However, stromal cells grown without FGF2 did not support proacinar organoid differentiation and instead differentiated into myofibroblasts. In organoids, TGFß1 treatment stimulated myofibroblast differentiation and inhibited the proacinar cell differentiation of epithelial progenitor cells. Conversely, FGF2 reversed the effects of TGFß1. We also demonstrated that adult salivary stromal cells were FGF2 responsive and could promote proacinar cell differentiation. These FGF2 signaling pathways may have applications in future regenerative therapies.

Raman microspectroscopy fingerprinting of organoid differentiation state.

BACKGROUND: Organoids, which are organs grown in a dish from stem or progenitor cells, model the structure and function of organs and can be used to define molecular events during organ formation, model human disease, assess drug responses, and perform grafting in vivo for regenerative medicine approaches. For therapeutic applications, there is a need for nondestructive methods to identify the differentiation state of unlabeled organoids in response to treatment with growth factors or pharmacologicals. METHODS: Using complex 3D submandibular salivary gland organoids developed from embryonic progenitor cells, which respond to EGF by proliferating and FGF2 by undergoing branching morphogenesis and proacinar differentiation, we developed Raman confocal microspectroscopy methods to define Raman signatures for each of these organoid states using both fixed and live organoids. RESULTS: Three separate quantitative comparisons, Raman spectral features, multivariate analysis, and machine learning, classified distinct organoid differentiation signatures and revealed that the Raman spectral signatures were predictive of organoid phenotype. CONCLUSIONS: As the organoids were unlabeled, intact, and hydrated at the time of imaging, Raman spectral fingerprints can be used to noninvasively distinguish between different organoid phenotypes for future applications in disease modeling, drug screening, and regenerative medicine.

RDH10-mediated retinol metabolism and RARα-mediated retinoic acid signaling are required for submandibular salivary gland initiation.

In mammals, the epithelial tissues of major salivary glands generate saliva and drain it into the oral cavity. For submandibular salivary glands (SMGs), the epithelial tissues arise during embryogenesis from naïve oral ectoderm adjacent to the base of the tongue, which begins to thicken, express SOX9 and invaginate into underlying mesenchyme. The developmental mechanisms initiating salivary gland development remain unexplored. In this study, we show that retinoic acid (RA) signaling activity at the site of gland initiation is colocalized with expression of retinol metabolic genes Rdh10 and Aldh1a2 in the underlying SMG mesenchyme. Utilizing a novel ex vivo assay for SMG initiation developed for this study, we show that RDH10 and RA are required for salivary gland initiation. Moreover, we show that the requirement for RA in gland initiation involves canonical signaling through retinoic acid receptors (RAR). Finally, we show that RA signaling essential for gland initiation is transduced specifically through RARα, with no contribution from other RAR isoforms. This is the first study to identify a molecular signal regulating mammalian salivary gland initiation.

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

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

Endothelial cell regulation of salivary gland epithelial patterning.

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

Mesenchymal Cells Affect Salivary Epithelial Cell Morphology on PGS/PLGA Core/Shell Nanofibers.

Engineering salivary glands is of interest due to the damaging effects of radiation therapy and the autoimmune disease Sjögren's syndrome on salivary gland function. One of the current problems in tissue engineering is that in vitro studies often fail to predict in vivo regeneration due to failure of cells to interact with scaffolds and of the single cell types that are typically used for these studies. Although poly (lactic co glycolic acid) (PLGA) nanofiber scaffolds have been used for in vitro growth of epithelial cells, PLGA has low compliance and cells do not penetrate the scaffolds. Using a core-shell electrospinning technique, we incorporated poly (glycerol sebacate) (PGS) into PLGA scaffolds to increase the compliance and decrease hydrophobicity. PGS/PLGA scaffolds promoted epithelial cell penetration into the scaffold and apical localization of tight junction proteins, which is necessary for epithelial cell function. Additionally, co-culture of the salivary epithelial cells with NIH3T3 mesenchymal cells on PGS/PLGA scaffolds facilitated epithelial tissue reorganization and apical localization of tight junction proteins significantly more than in the absence of the mesenchyme. These data demonstrate the applicability of PGS/PLGA nanofibers for epithelial cell self-organization and facilitation of co-culture cell interactions that promote tissue self-organization in vitro.

Heterotypic control of basement membrane dynamics during branching morphogenesis.

Many mammalian organs undergo branching morphogenesis to create highly arborized structures with maximized surface area for specialized organ function. Cooperative cell-cell and cell-matrix adhesions that sculpt the emerging tissue architecture are guided by dynamic basement membranes. Properties of the basement membrane are reciprocally controlled by the interacting epithelial and mesenchymal cell populations. Here we discuss how basement membrane remodeling is required for branching morphogenesis to regulate cell-matrix and cell-cell adhesions that are required for cell patterning during morphogenesis and how basement membrane impacts morphogenesis by stimulation of cell patterning, force generation, and mechanotransduction. We suggest that in addition to creating mature epithelial architecture, remodeling of the epithelial basement membrane during branching morphogenesis is also essential to promote maturation of the stromal mesenchyme to create mature organ structure. Recapitulation of developmental cell-matrix and cell-cell interactions are of critical importance in tissue engineering and regeneration strategies that seek to restore organ function.

Integrins promote cytokinesis through the RSK signaling axis.

Cytokinesis is the final stage in cell division. Although integrins can regulate cytokinesis, the mechanisms involved are not fully understood. In this study, we demonstrate that integrin-regulated ERK (extracellular signal-related kinase) and RSK (p90 ribosomal S6 kinase) signaling promotes successful cytokinesis. Inhibiting the activation of ERK and RSK in CHO cells by a mutation in the integrin ß1 cytoplasmic tail or with pharmacological inhibitors results in the accumulation of cells with midbodies and the formation of binucleated cells. Activation of ERK and RSK signaling by the expression of constitutively active RAF1 suppresses the mutant phenotype in a RSK-dependent manner. Constitutively active RSK2 also restores cytokinesis inhibited by the mutant integrin. Importantly, the regulatory role of the RSK pathway is not specific to CHO cells. MCF-10A human mammary epithelial cells and HPNE human pancreatic ductal epithelial cells exhibit a similar dependence on RSK for successful cytokinesis. In addition, depriving mitotic MCF10A cells of integrin-mediated adhesion by incubating them in suspension suppressed ERK and RSK activation and resulted in a failure of cytokinesis. Furthermore, inhibition of RSK or integrins within the 3D context of a developing salivary gland organ explant also leads to an accumulation of epithelial cells with midbodies, suggesting a similar defect in cytokinesis. Interestingly, neither ERK nor RSK regulates cytokinesis in human fibroblasts, suggesting cell-type specificity. Taken together, our results identify the integrin-RSK signaling axis as an important regulator of cytokinesis in epithelial cells. We propose that the proper interaction of cells with their microenvironment through integrins contributes to the maintenance of genomic stability by promoting the successful completion of cytokinesis.

Highly multiplexed single-cell analysis of formalin-fixed, paraffin-embedded cancer tissue.

Limitations on the number of unique protein and DNA molecules that can be characterized microscopically in a single tissue specimen impede advances in understanding the biological basis of health and disease. Here we present a multiplexed fluorescence microscopy method (MxIF) for quantitative, single-cell, and subcellular characterization of multiple analytes in formalin-fixed paraffin-embedded tissue. Chemical inactivation of fluorescent dyes after each image acquisition round allows reuse of common dyes in iterative staining and imaging cycles. The mild inactivation chemistry is compatible with total and phosphoprotein detection, as well as DNA FISH. Accurate computational registration of sequential images is achieved by aligning nuclear counterstain-derived fiducial points. Individual cells, plasma membrane, cytoplasm, nucleus, tumor, and stromal regions are segmented to achieve cellular and subcellular quantification of multiplexed targets. In a comparison of pathologist scoring of diaminobenzidine staining of serial sections and automated MxIF scoring of a single section, human epidermal growth factor receptor 2, estrogen receptor, p53, and androgen receptor staining by diaminobenzidine and MxIF methods yielded similar results. Single-cell staining patterns of 61 protein antigens by MxIF in 747 colorectal cancer subjects reveals extensive tumor heterogeneity, and cluster analysis of divergent signaling through ERK1/2, S6 kinase 1, and 4E binding protein 1 provides insights into the spatial organization of mechanistic target of rapamycin and MAPK signal transduction. Our results suggest MxIF should be broadly applicable to problems in the fields of basic biological research, drug discovery and development, and clinical diagnostics.

ROCK1-directed basement membrane positioning coordinates epithelial tissue polarity.

The basement membrane is crucial for epithelial tissue organization and function. However, the mechanisms by which basement membrane is restricted to the basal periphery of epithelial tissues and the basement membrane-mediated signals that regulate coordinated tissue organization are not well defined. Here, we report that Rho kinase (ROCK) controls coordinated tissue organization by restricting basement membrane to the epithelial basal periphery in developing mouse submandibular salivary glands, and that ROCK inhibition results in accumulation of ectopic basement membrane throughout the epithelial compartment. ROCK-regulated restriction of PAR-1b (MARK2) localization in the outer basal epithelial cell layer is required for basement membrane positioning at the tissue periphery. PAR-1b is specifically required for basement membrane deposition, as inhibition of PAR-1b kinase activity prevents basement membrane deposition and disrupts overall tissue organization, and suppression of PAR-1b together with ROCK inhibition prevents interior accumulations of basement membrane. Conversely, ectopic overexpression of wild-type PAR-1b results in ectopic interior basement membrane deposition. Significantly, culture of salivary epithelial cells on exogenous basement membrane rescues epithelial organization in the presence of ROCK1 or PAR-1b inhibition, and this basement membrane-mediated rescue requires functional integrin ß1 to maintain epithelial cell-cell adhesions. Taken together, these studies indicate that ROCK1/PAR-1b-dependent regulation of basement membrane placement is required for the coordination of tissue polarity and the elaboration of tissue structure in the developing submandibular salivary gland.

Cell-based multi-parametric model of cleft progression during submandibular salivary gland branching morphogenesis.

Cleft formation during submandibular salivary gland branching morphogenesis is the critical step initiating the growth and development of the complex adult organ. Previous experimental studies indicated requirements for several epithelial cellular processes, such as proliferation, migration, cell-cell adhesion, cell-extracellular matrix (matrix) adhesion, and cellular contraction in cleft formation; however, the relative contribution of each of these processes is not fully understood since it is not possible to experimentally manipulate each factor independently. We present here a comprehensive analysis of several cellular parameters regulating cleft progression during branching morphogenesis in the epithelial tissue of an early embryonic salivary gland at a local scale using an on lattice Monte-Carlo simulation model, the Glazier-Graner-Hogeweg model. We utilized measurements from time-lapse images of mouse submandibular gland organ explants to construct a temporally and spatially relevant cell-based 2D model. Our model simulates the effect of cellular proliferation, actomyosin contractility, cell-cell and cell-matrix adhesions on cleft progression, and it was used to test specific hypotheses regarding the function of these parameters in branching morphogenesis. We use innovative features capturing several aspects of cleft morphology and quantitatively analyze clefts formed during functional modification of the cellular parameters. Our simulations predict that a low epithelial mitosis rate and moderate level of actomyosin contractility in the cleft cells promote cleft progression. Raising or lowering levels of contractility and mitosis rate resulted in non-progressive clefts. We also show that lowered cell-cell adhesion in the cleft region and increased cleft cell-matrix adhesions are required for cleft progression. Using a classifier-based analysis, the relative importance of these four contributing cellular factors for effective cleft progression was determined as follows: cleft cell contractility, cleft region cell-cell adhesion strength, epithelial cell mitosis rate, and cell-matrix adhesion strength.

Senescence and fibrosis in salivary gland aging and disease.

Salivary gland hypofunction is highly prevalent in aged and diseased individuals leading to significant discomfort and morbidity. One factor that contributes to salivary gland hypofunction is cellular aging, or senescence. Senescent cells can impair gland function by secreting paracrine-acting growth factors and cytokines, known as senescence-associated secretory phenotype (SASP) factors. These SASP factors stimulate inflammation, propagate the senescent phenotype through the bystander effect, and stimulate fibrosis. As senotherapeutics that target senescent cells have shown effectiveness in limiting disease manifestations in other conditions, there is interest in the use of these drugs to treat salivary gland hypofunction. In this review, we highlight the contribution of senescence and fibrosis to salivary gland pathologies. We also discuss therapeutic approaches to eliminate or modulate the senescent SASP phenotype for treating age-related salivary gland diseases and extending health span.

Raman hyperspectroscopy of saliva and machine learning for Sjögren's disease diagnostics.

Sjögren's disease is an autoimmune disorder affecting exocrine glands, causing dry eyes and mouth and other morbidities. Polypharmacy or a history of radiation to the head and neck can also lead to dry mouth. Sjogren's disease is often underdiagnosed due to its non-specific symptoms, limited awareness among healthcare professionals, and the complexity of diagnostic criteria, limiting the ability to provide therapy early. Current diagnostic methods suffer from limitations including the variation in individuals, the absence of a single diagnostic marker, and the low sensitivity and specificity, high cost, complexity, and invasiveness of current procedures. Here we utilized Raman hyperspectroscopy combined with machine learning to develop a novel screening test for Sjögren's disease. The method effectively distinguished Sjögren's disease patients from healthy controls and radiation patients. This technique shows potential for development of a single non-invasive, efficient, rapid, and inexpensive medical screening test for Sjögren's disease using a Raman hyper-spectral signature.

Evaluation of Alginate Hydrogel Microstrands for Stromal Cell Encapsulation and Maintenance.

Mesenchymal stromal cells (MSCs) have displayed potential in regenerating organ function due to their anti-fibrotic, anti-inflammatory, and regenerative properties. However, there is a need for delivery systems to enhance MSC retention while maintaining their anti-fibrotic characteristics. This study investigates the feasibility of using alginate hydrogel microstrands as a cell delivery vehicle to maintain MSC viability and phenotype. To accommodate cell implantation needs, we invented a Syringe-in-Syringe approach to reproducibly fabricate microstrands in small numbers with a diameter of around 200 µm and a porous structure, which would allow for transporting nutrients to cells by diffusion. Using murine NIH 3T3 fibroblasts and primary embryonic 16 (E16) salivary mesenchyme cells as primary stromal cell models, we assessed cell viability, growth, and expression of mesenchymal and fibrotic markers in microstrands. Cell viability remained higher than 90% for both cell types. To determine cell number within the microstrands prior to in vivo implantation, we have further optimized the alamarBlue assay to measure viable cell growth in microstrands. We have shown the effect of initial cell seeding density and culture period on cell viability and growth to accommodate future stromal cell delivery and implantation. Additionally, we confirmed homeostatic phenotype maintenance for E16 mesenchyme cells in microstrands.

The matrix reorganized: extracellular matrix remodeling and integrin signaling.

Via integrins, cells can sense dimensionality and other physical and biochemical properties of the extracellular matrix (ECM). Cells respond differently to two-dimensional substrates and three-dimensional environments, activating distinct signaling pathways for each. Direct integrin signaling and indirect integrin modulation of growth factor and other intracellular signaling pathways regulate ECM remodeling and control subsequent cell behavior and tissue organization. ECM remodeling is critical for many developmental processes, and remodeled ECM contributes to tumorigenesis. These recent advances in the field provide new insights and raise new questions about the mechanisms of ECM synthesis and proteolytic degradation, as well as the roles of integrins and tension in ECM remodeling.

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.

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.

Identifying fibrogenic cells following salivary gland obstructive injury.

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

Identifying Fibrogenic Cells Following Salivary Gland Obstructive Injury.

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