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.

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.

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.

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.