Defining the relationship of salivary gland malignancies to novel cell subpopulations in human salivary glands using single nucleus RNA-sequencing.

Salivary glands have essential roles in maintaining oral health, mastication, taste and speech, by secreting saliva. Salivary glands are composed of several types of cells, and each cell type is predicted to be involved in the carcinogenesis of different types of cancers including adenoid cystic carcinoma (ACC), acinic cell carcinoma (AciCC), salivary duct carcinoma (SDC), myoepithelial carcinoma (MECA) and other histology. In our study, we performed single nucleus RNA-seq on three human salivary gland samples to clarify the gene expression profile of each complex cellular component of the salivary glands and related these expression patterns to expression found in salivary gland cancers (SGC) to infer cell of origin. By single nucleus RNA-seq, salivary gland cells were stratified into four clusters: acinar cells, ductal cells 1, ductal cells 2 and myoepithelial cells/stromal cells. The localization of each cell group was verified by IHC of each cluster marker gene, and one group of ductal cells was found to represent intercalated ductal cells labeled with HES1. Furthermore, in comparison with SGC RNA-seq data, acinar cell markers were upregulated in AciCC, but downregulated in ACC and ductal cell markers were upregulated in SDC but downregulated in MECA, suggesting that markers of origin are highly expressed in some SGC. Cell type expressions in specific SGC histology are similar to those found in normal salivary gland populations, indicating a potential etiologic relationship.