Molecular characterization of the salivary adenoid cystic carcinoma immune landscape by anatomic subsites.

Adenoid cystic carcinoma (AdCC) is a slow-growing salivary gland malignancy that relapses frequently. AdCCs of the submandibular gland exhibit unique differences in prognosis and treatment response to adjuvant radiotherapy compared to other sites, yet the role of tumor anatomic subsite on gene expression and tumor immune microenvironment (TIME) composition remains unclear. We used 87 samples, including 48 samples (27 AdCC and 21 normal salivary gland tissue samples) from 4 publicly available AdCC RNA sequencing datasets, a validation set of 33 minor gland AdCCs, and 39 samples from an in-house cohort (30 AdCC and 9 normal salivary gland samples). RNA sequencing data were used for single sample gene set enrichment analysis and TIME deconvolution. Quantitative PCR and multiplex immunofluorescence were performed on the in-house cohort. Wilcoxon rank-sum, nonparametric equality-of-medians tests and linear regression models were used to evaluate tumor subsite differences. AdCCs of different anatomic subsites including parotid, submandibular, sublingual, and minor salivary glands differed with respect to expression of several key tumorigenic pathways. Among the three major salivary glands, the reactive oxygen species (ROS)/nuclear factor erythroid 2-related factor 2 (NRF2) pathway signature was significantly underexpressed in AdCC of submandibular compared to parotid and sublingual glands while this association was not observed among normal glands. Additionally, the NRF2 pathway, whose expression was associated with favorable overall survival, was overexpressed in AdCCs of parotid gland compared to minor and submandibular glands. The TIME deconvolution identified differences in CD4