AUTHOR=Meng Liangliang , He Xiaoxi , Hong Quan , Qiao Bo , Zhang Xiao , Wu Bin , Zhang Xiaobo , Wei Yingtian , Li Jing , Ye Zhaoxiang , Xiao Yueyong TITLE=CCR4, CCR8, and P2RY14 as Prognostic Factors in Head and Neck Squamous Cell Carcinoma Are Involved in the Remodeling of the Tumor Microenvironment JOURNAL=Frontiers in Oncology VOLUME=11 YEAR=2021 URL=https://www.frontiersin.org/journals/oncology/articles/10.3389/fonc.2021.618187 DOI=10.3389/fonc.2021.618187 ISSN=2234-943X ABSTRACT=

The tumor microenvironment (TME) plays a critical role in the initiation and progression of cancer. However, the specific mechanism of its regulation in head and neck squamous cell carcinoma (HNSCC) remains unclear. In this study, we first applied the ESTIMATE method to calculate the immune and stromal scores in patients’ tumor tissues from The Cancer Genome Atlas (TCGA) database. GSE41613, GSE30784, and GSE37991 data sets from the Gene Expression Omnibus (GEO) database were recruited for further validation. Differentially expressed genes (DEGs) were identified and then analyzed by Cox regression analysis and protein-protein interaction (PPI) network construction. DEGs significantly associated with prognosis and TME will be identified as hub genes. These genes were also validated at the protein level by immunohistochemical analysis of 10 pairs of primary tumor tissues and the adjacent normal tissues from our institution. The relationship between hub genes expression and immune cell fraction estimated by CIBERSORT software was also examined. 275 DEGs were significantly associated with TME. CCR4, CCR8, and P2RY14 have then identified as hub genes by intersection Cox and PPI analysis. Further investigation revealed that the expression of CCR4, CCR8, and P2RY14 was negatively correlated with clinicopathological characteristics (clinical stage, T stage) and positively associated with survival in HNSCC patients, especially in male patients. The expression of CCR8 and P2RY14 was lower in males than in females. CCR8 and P2RY14 were differentially expressed in tumor tissues than normal tissues, and the results were validated at the protein level by immunohistochemistry experiments. Gene set enrichment analysis (GSEA) showed that the high expression groups’ hub genes were mainly enriched for immune-related activities. In the low-expression groups, genes were primarily enriched in metabolic pathways. CIBERSORT results showed that the expression of these genes was all negatively correlated with the fraction of memory B cells and positively correlated with the fraction of the other four cells, including naive B cells, resting T cells CD4 memory, T cells follicular helper, and T cells regulatory (Tregs). The results suggest that CCR4, CCR8, and P2RY14 may be responsible for maintaining the immune dominance of TME, thus leading to a better prognosis.