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.

Purpose

Human papillomavirus (HPV)-related head and neck squamous cell carcinoma (HNSCC) is associated with daily marijuana use and is also increasing in parallel with increased marijuana use in the United States. Our study is designed to define the interaction between cannabinoids and HPV-positive HNSCC.

Experimental design

The expression of cannabinoid receptors CNR1 and CNR2 was analyzed using The Cancer Genome Atlas (TCGA) HNSCC data. We used agonists, antagonists, siRNAs, or shRNA-based models to explore the roles of CNR1 and CNR2 in HPV-positive HNSCC cell lines and animal models. Cannabinoid downstream pathways involved were determined by Western blotting and analyzed in a primary HPV HNSCC cohort with single-sample gene set enrichment analysis (ssGSEA) and the OncoGenome Positioning System (Onco-GPS).

Results

In TCGA cohort, the expression of CNR1 and CNR2 was elevated in HPV-positive HNSCC compared with HPV-negative HNSCC, and knockdown of CNR1/CNR2 expression inhibited proliferation in HPV-positive HNSCC cell lines. Specific CNR1 and CNR2 activation as well as nonselective cannabinoid receptor activation in cell lines and animal models promoted cell growth, migration, and inhibited apoptosis through p38 MAPK pathway activation. CNR1/CNR2 antagonists suppressed cell proliferation and migration and induced apoptosis. Using whole-genome expression analysis in a primary HPV HNSCC cohort, we identified specific p38 MAPK pathway activation signature in tumors from HPV HNSCC patients with objective measurement of concurrent cannabinoid exposure.

Conclusions

Cannabinoids can promote progression of HPV-positive HNSCC through p38 MAPK pathway activation.