Single-Cell Transcriptomic Analysis of Salivary Gland Endothelial Cells

Abstract

Vascular endothelial cells have important tissue-specific functions in fibrosis and regeneration. In the salivary gland, endothelial cells are required for proper development, but their roles within adult glands are largely unknown. To identify ligand–receptor interactions between endothelial cells and other cell types that may be important during fibrosis and regeneration, we used a reversible ductal ligation injury. To induce injury, a clip was applied to the primary ducts for 14 d, and to induce a regenerative response, the clip was subsequently removed for 5 d. To identify endothelial cell-produced factors, we used single-cell RNA sequencing of stromal-enriched cells from adult female 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 expressed many 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-d ligated, mock-ligated, and 5-d deligated stromal-enriched transcripts and lineage tracing revealed that endothelial cells retain their identity following ligation and recovery from injury. CellChat and NATMI were used to predict changes in ligand–receptor interactions from endothelial cells to other cells in response to ligation and deligation. CellChat and NATMI predicted that after ligation, interactions with fibroblasts, epithelial cells, and glial cells were increased, and following deligation, interactions with pericyte, glia, fibroblasts, and immune cells were increased. Some of the highest-ranked interactions predicted in ligated compared to mock endothelial cells were between glial cells via Col4a2-Cd93 and Jag2-Notch1, as well as epithelial cells via Pecam1-Cd38, while in deligated compared to ligated endothelial cells, the top interactions were between fibroblasts via Ntf3-Ntrk2, glial cells via Hspg2-Itgb1, and pericytes via Jam2-F11r. Understanding salivary gland endothelial cell signaling will inform future endothelial cell-based regenerative therapies.