Tumor microenvironment characterization and Immune infiltration in head and neck squamous cell carcinoma

Objection: Head and neck squamous cell carcinoma (HNSCC), a basic malignant tumor of the head and neck distinct. As a famous heterogeneous disease, the jobs of immune cells inside the tumor-related as yet missing for HNSCC, particularly in current immunotherapy. Method: We explored the TME, TMB and evaluate the 22 TIICs subsets of immune response based on GEO and TCGA database of HNSCC to explore its relationship with atomic subpopulation, survival, function and expression difference and reveal potential targets and biomarkers for immunotherapy. Results: Observing the download of GSE6631 database contained 22 HNSCC samples and 22 normal samples and TCGA database contained 111 HNSCC and 12 normal tissues. The results suggested that the expression of macrophages M0 and T cells CD4 memory resting was significant difference and may plays an important role in regulate cancer progression (P<0.05). The result of tumor mutational burden revealed that the most common somatic mutations variant classification was missense mutation, the most common DNA sequence polymorphism type was SNP, the most common single nucleotide variants (SNV) class was C>T, the variants per sample median was 78 in HNSCC patients. Top 10 mutated genes that related to TMB was TP53, TTN, FAT1, MUC16, CDKN2A, CSMD3, SYNE1, LRP1B, NOTCH1 and PIK3CA. We portrayed the immune scene in detail, uncovering the awesome immune infiltration styles of various subtypes in HNSCC. Conclusion: The intricate connection between TIIC, TMB and genomic alterations was additionally set up. Our paintings advance the information of immune response and offers significant assets for research to enhance immunotherapy. performed utilizing GO (cellular components, biological and molecular functions) enrichment and KEGG pathway analysis via R package. Finally, we evaluated the related to tumor-infiltrating immune cells and TMB. restores diminished levels of glutathione on account of chemotherapy and radiotherapy and it encourages recovery processes of the intestine epithelium and immunological system . Our results also suggested that the expression of T cells regulatory (Tregs) was critical distinction between high TMB group and low TMB group indicating Tregs is essential in the tumor immune response. Meanwhile, our study has limitations, one is that our results not verified in clinical samples and cannot provides precise clinical data as well as the relatively few patients, another is that due to the histological types and multiple anatomical sites and of HNSCC the tumor-infiltrating immune cells may vary widely. We will do a lot of related research in the follow-up and eventually translated into the development of new strategies for precision cancer medicine. We portrayed the immune scene in detail, uncovering the awesome immune infiltration styles of various subtypes in HNSCC. The intricate connection between TIIC, TMB and genomic alterations was additionally set up. Our paintings advance the information of immune response and offers significant assets for research to enhance immunotherapy.

3 when detected of HNSCC is usually made at advanced stages and the 5-year survival rate is still under 50% now, while due to local recurrence and metastasis, which reduces survival rate to 35%. 4 The occurrences and progression of HNSCC is a complicated process involving multiple molecules. While, the malignant phenotypes of cancers are characterized by the intrinsic activities of tumor cells as well as by the immune cells recruited to and activated 5 . As a heterogeneous disease, the jobs of immune cells in the tumor-related as yet missing for HNSCC, particularly in current immunotherapy.
The tumor microenvironment (TME) comprises not only cell types (immune cells, fibroblasts and endothelial cells etc.) but also extracellular parts (hormones, cytokines, extracellular matrix, growth factors, etc.) that are encompassing tumor cells and fed by a vascular network 6 . The TME plays an indispensable role across tumor initiation, progression and even metastasis and effects therapeutic efficacy 7 .Tumor-infiltrating immune cells (TIIC) have been accounted for to be viably targeted by drugs and to connect with the clinical result by structure an ecosystem in the tumor microenvironment to manager cancer progression and have appeared potential prognostic value 8 .
There is a complicated interaction which considerable prognostic relevance as the immune system has a double job by means of assisting both host barrier and tumor progression 9 . Upregulation of inhibitory immune checkpoints with the aid of TIICs by tumor cells prompts tumor avoidance from host immunosurveillance 10 . Chromatin modifiers lead to CD8 + T-cell prohibition and antigen presentation limit of HNSCC 11 . Moreover, no past investigation has revealed insight into the prognostic estimation of these TIICs subpopulation. Past examinations assessed the estimation of these TIIC by flow cytometry and immunohistochemistry (IHC) 12 and are limited by the quantity of fluorescent channels accessible as well as few immune cell types can be assessed immediately 13 . A metagene instrument CIBERSORT utilizes deconvolution of mass gene expression information and a sophisticated algorithm for quantification the cellular element of immune reaction and cell kinds in heterogeneous samples has significantly extended the capability of the genomic database 14 − 4 15 .Tumor mutational burden (TMB) has been depicted as indicator of immunological reaction and tumor behavior 16 − 17 , mutation formation advances carcinogenesis by means of activation or inactivation of genes and related pathways, in this way creating novel peptide sequences which can animate immune reaction 18 .High TMB may additionally constitute a high underlying wide variety of drivers and imply a higher-risk tumor. As is a developing biomarker of affectability to immune checkpoint inhibitors and has been appeared to be more significantly connected with reaction to PD-1 and PD-L1 blockade immunotherapy than PD-1 or PD-L1 expression 19 . In this study, we carried out CIBERSORT to explored the TME, TMB and quantify the 22 TIICs subsets of immune reaction based on GEO and TCGA database of HNSCC to explore its association with molecular subpopulation, survival, function and expression difference and reveal potential targets and biomarkers for immunotherapy.

Materials And Methods
Data acquisition. TMB is the total number of mutations per megabase of tumor tissue. In general terms, the mutation density of a tumor gene, that is, the average number of mutations in the tumor genome, including the total number of gene coding errors, base substitutions, gene insertions, or deletion errors. The larger the TMB, the easier it is to be discovered by immune cells, and the easier it is to become a target for tumor immunity, so that the more likely it is to be effective for immunotherapy. We described the copy number and somatic mutations characteristic based on TCGA databased and then analyzed the connection between TMB and patients' clinical factors, including age, gender, grade and TNM stage and we also analyzed the effect of TMB on survival. TCGA workflow type based on VarScan2 Variant Aggregation and Masking. Differentially expressed genes (DEG) associated with TMB were resolved using the R package limma 21 , which implements an empirical Bayesian way to estimate geneexpression modifications the use of moderated t tests. DEGs among TMB were controlled by using significance standards as applied in the R package limma. The adjusted P value for multiple testing was calculated the usage of the Benjamini-Hochberg correction 22 .To examine the function of the identified DMGs, biological analyses were performed utilizing GO (cellular components, biological processes and molecular functions) enrichment and KEGG pathway analysis via R package. Finally, we evaluated the related to tumor-infiltrating immune cells and TMB.

Statistical analysis.
Statistical analyses have been conducted the use of R version 3.5.3 and Bioconductor. For examinations of two gatherings, statistical significance for normally distributed variables was evaluated by unpaired Student t tests, and nonnormally distributed variables were revealed by Wilcoxon tests. Each dataset was handled by a weighted average approach to contrast the differences in the composition of TIIC and using boxplot, heatmap, corHeatmap and vioplot to visualization the difference in normal and tumor samples. Overall survival (OS) was characterized as the time interval from the date of diagnosis to the date of death. The listwise deletion technique was utilized to deal with missing data, which excluded the entire sample from the investigation if any single value was absent. Wilcox test analysis was performed to evaluate the differences not only in the gene expression of immune checkpoint molecules but also TMB clinical information between tumor and normal tissues. To recognize differential genes in the GEGs evaluation, we applied the Benjamini-Hochberg way to transform the P values to FDRs. For every single statistical analysis, a Pvalue < 0.05 was viewed as significant.

Results
Clinical information.
Observing the download of GSE6631 database contained 22 HNSCC samples and 22 normal samples and after immune cell matrix filtering according to P value less than 0.5, the compliance with the inclusion criteria contained 12 HNSCC samples and 10 normal samples. TCGA database contained 111 HNSCC and 12 normal tissues and after immune cell matrix filtering according to P value less than 0.5, the compliance with the inclusion criteria contained 98 HNSCC and 1 normal tissue, while the clinical patients contained 117 patients and the clinicopathological attributes of these samples are appeared in Table 1.

Discussion
In this study, we carried out CIBERSORT to conducted a comprehensive and detailed assessment of the TME, TMB and quantify the 22 TIICs subsets of immune response based on GEO and TCGA database in HNSCC to explore its association with molecular subpopulation, survival, function and 9 expression difference. We also detailed that the constitute of TIIC differs substantially in numerous subtypes and clinical factors of HNSCC. Tumor mutational burden firmly add to the generation of novel tumor epitopes. Knowing whether a more highly-immunogenic tumor consists of a right away link to mutational burden help us understand the mechanistic clarification for observed clinical survival patterns. In our outcomes, TMB was firmly related with TMB among TCGA data, supporting TMB's position as a strong proxy for TME. This ought to be of terrific interest for perspective on the present improvement of immunomodulatory therapies and may be uncover potential biomarkers or targets for immunotherapy or determinants for prognosis. Resting memory CD4 + cells are subject to signals from touch with IL-7 and IL-15, yet not MHC class II, for their survival and intermittent homeostatic proliferation 30 . Our results suggested that the expression of macrophages M0 and T cells CD4 memory resting was significant difference and may plays an important role in regulate cancer progression. While the other immune cells infiltration was not obvious may be due to the too few samples size and needs more researches to conduct. We also analysis the interaction between different 22 TIICs and provided information for future immune cells co-expression study in HNSCC researches.
The result of tumor mutational burden revealed that the most common somatic mutations variant classification was missense mutation, the most common DNA sequence polymorphism type was SNP, the most common single nucleotide variants (SNV) class was C > T. Changes that modify amino acid sequences are known as missense mutations and these mutations may influence protein stability/structure and damage protein interactions with different biomolecules, translation proteins non-functional and advancing tumor progression 31  In order to analysis the relationship between TMB and differentially expressed genes, we divide TMB into high TMB and low TMB according to the median value and explored the functions of this DEGs.
Biological processes of DEGs suggested that these genes may play a significant role in tissue development and the response to acid chemical, molecular function were mainly enhanced in extracellular matrix structural constituent and mainly affect cellular structure. KEGG analysis revealed that the DMGs were for the most part enriched in glutathione metabolism. Glutamine (GLN) usually 11 known as a significant metabolite utilized for the development of malignancy cells. GLN consumption restores diminished levels of glutathione on account of chemotherapy and radiotherapy and it encourages recovery processes of the intestine epithelium and immunological system 41 . Our results also suggested that the expression of T cells regulatory (Tregs) was critical distinction between high TMB group and low TMB group indicating Tregs is essential in the tumor immune response.
Meanwhile, our study has limitations, one is that our results not verified in clinical samples and cannot provides precise clinical data as well as the relatively few patients, another is that due to the Declarations RNA-seq data and corresponding clinical data were acquired from the data portal for TCGA (https://portal. gdc.cancer.gov/)

Ethics approval and consent to participate
Not applicable.

Consent for publication
Not applicable.