Impact of IL-10 Promoter Polymorphisms and IL-10 Serum Levels on Advanced Laryngeal Squamous Cell Carcinoma and Survival Rate

Background/Aim: Prognosis of advanced stages of laryngeal squamous cell carcinoma (LSCC) remains poor. To clarify therapeutic targets and improve survival rate, identification of new specific and prognostic biomarkers of LSCC is required. The study aimed to evaluate the impact of IL-10:rs1800871, rs1800872, rs1800896 single nucleotide polymorphisms (SNPs), and IL-10 serum levels on LSCC development and determine associations of selected SNPs with patient survival rate. Patients and Methods: A total of 300 LSCC patients and 533 controls were included in the study. Genotyping was carried out using RT-PCR; IL-10 serum levels were analyzed by ELISA. Results: Significant associations were identified between IL-10 rs1800871 variants and advanced stage of LSCC patient group in the codominant, recessive and additive models (OR=0.473, p=0.027; OR=0.510, p=0.040; and OR=0.733; p=0.037). Significant variants of IL-10 rs1800872 were determined in the codominant, recessive and additive models (OR=0.473, p=0.027; OR=0.510, p=0.040; and OR=0.733, p=0.037). The distribution of IL-10 SNPs genotypes did not impact LSCC patient survival rate (respectively, p=0.952; p=0.952; p=0.991). Conclusion: IL-10:rs1800871 and rs1800872 SNPs are associated with advanced stage of LSCC. The genotypic distribution of IL-10 SNPs does not influence the survival rate of LSCC patients.

Interleukin-10 (IL-10) is a pleiotropic cytokine that is originally called cytokine synthesis inhibitory factor (CSIF) and has structural similarities to interferon-gamma (INFγ) (12,13). In humans, IL-10 is encoded by the IL-10 gene and located on chromosome 1 (1q31-1q32) (14). Various cells can produce IL-10, including macrophages, mast cells, T cells regulatory subsets, dendritic cells, B cells, a cluster of differentiation (CD) 4+ and CD8+ T cells, and natural killer (NK) cells (12,15,16). The receptor of IL-10 is a tworeceptor complex composed of two copies each of IL-10 specific chain -IL10R1 and IL10R2 proteins (17). Activated IL-10R complex induces signal transducer and activator of transcription 3 (STAT3) signaling via phosphorylation of the cytoplasmic tails of IL-10R1 and IL10R2 by Janus Kinase 1 (JAK1) and Tyrosine kinase 2 (Tyk2) (12,17). In unstimulated tissues, the expression of IL-10 is deficient, suggesting the importance of IL-10 for anti-inflammatory functions and indicating that provoking factors as massive inflammation are required for observation of elevated IL-10 levels (18).
IL-10 is known as an immunosuppressive cytokine that stimulates proliferation and metastasis of tumor cells (19). However, the role of IL-10 in carcinogenesis remains controversial and considered to be paradoxical (12,20,21). Various studies have demonstrated a positive correlation between IL-10 levels in serum and tumor tissues and a poor prognosis in melanoma, lung cancer, T/NK cells lymphomas, and HNSCC (22)(23)(24)(25)(26). Previous studies also suggested that tumor cells can produce IL-10 by themselves, and IL-10 might be a potential biomarker for IL-10-targeted therapy leading to the conclusion that antagonists of IL-10 might be used as a novel treatment of cancer patients (22-25, 27, 28). In contrast, the IL-10 anti-tumor effect is also demonstrated in murine tumor models where rapid tumor rejection that increased secretion of IL-10 was identified (12). Moreover, it is suggested that IL-10 plays a critical role in memory development and the function of appropriate CD8 T cells (29). Regarding the IL-10 role in anti-tumor immune responses, it can be presumed that the transcription activity of IL-10 or IL-10 protein structure could be modified by genetic polymorphisms, which could provide susceptibility to different oncological disorders involving HNSCC (30,31).
Several studies indicating the role of IL-10: rs1800872, rs1800871, rs1800896 polymorphisms in HNSCC have already been published (32)(33)(34)(35). However, there is still a lack of data analyzing the associations between these polymorphisms and the pure group of LSCC patients and the impact on patients' survival rate.
The objective of this study was to evaluate the impact of IL-10: rs1800871, rs1800871, and rs1800896 SNPs and IL-10 serum levels on LSCC development and to determine possible associations of selected SNPs with the survival rate of LSCC patients.

Patients and Methods
The present case-control study was carried out at the Department of Otorhinolaryngology, Lithuanian University of Health Sciences (LUHS), Kaunas, Lithuania, and the Laboratory of Ophthalmology, Neuroscience Institute, LUHS, Kaunas, Lithuania.
Ethics Statement. The study protocol was approved by the Kaunas Regional Ethics Committee for Biomedical Research, LUHS (authorization number BE-2-37). An Informed Consent Form, in accordance with the Declaration of Helsinki, was obtained from all subjects involved in the study.
Study population. A total of 833 subjects were selected for the present study: 300 patients with LSCC and 533 healthy controls as a reference group. The demographic characteristics of LSCC patients and healthy control subjects are presented in Table I.  (36). Patients diagnosed with another type and localization of oncological disease, acute or chronic infectious diseases, as well as individuals using psychomotor suppressants and anti-epileptic drugs and persons younger than 18 years old, were rejected from this study.
Healthy controls. Patients who had an otorhinolaryngologist's consultation at the Out-patient Office of the Department of Otorhinolaryngology, LUHS, and were selected for surgical treatment (nasal bones reposition, septoplasty, rhinoseptoplasty, tympanoplasty, ossiculoplasty, tympanostomy, uvulopalatopharyngoplasty, or radiofrequency thermoablation) were involved in the present study as healthy control group subjects. The collection of peripheral venous blood samples was accomplished from the same catheter to induce general anaesthesia. Patients who attended the family doctor's consultation for the general check-up and had a collected peripheral blood sample were also included in this study. All the patients with the diagnosed oncologic disease, acute or chronic infectious diseases, individuals using psychomotor suppressants and anti-epileptic drugs, and persons younger than 18 years old were rejected from this study.
Subjects in the healthy control group were adjusted by age (p=0.054) and sex (p=0.813) according to demographic characteristics of the LSCC group (Table I).
SNP selection. The pathogenesis of LSCC on the molecular level remains unrevealed. According to published data, environmental carcinogens (37)(38)(39)(40)(41)(42) can induce damage to DNA and cause genomic irregularity. Therefore, further studies based on genetic alterations are required.
Recently, the controversial role of IL-10 suggested that IL-10 might be a significant biomarker for IL-10 targeted treatment for oncological patients (28). In accordance with these findings, for the present study, we selected IL-10: rs1800872, rs1800871, rs1800896 SNPs to determine possible associations with LSCC development. Selected IL-10 SNPs are the most frequently analyzed in the scientific literature and are strongly related to malignancies (43)(44)(45).
A few studies analyzing IL-10: rs1800872, rs1800871, rs1800896 associations with HNSCC have been published (32)(33)(34)(35)60). However, there is no data in the literature analyzing the association of these SNPs with the development of LSCC in a pure group of LSCC patients and the impact of these SNPs on the survival rate of LSCC patients.
DNA extraction, genotyping, and analysis of serum IL-10 levels.
The DNA extraction and analysis of IL-10: rs1800872, rs1800871, rs1800896 gene polymorphisms were performed at the Laboratory of Ophthalmology, Neuroscience Institute, LUHS.
The DNA was extracted from peripheral venous blood samples (leucocytes) collected in 200 μl test-tubes utilizing the silica-based membrane technology using a genomic DNA extraction kit (GeneJET Genomic DNA Purification Kit, Thermo Fisher Scientific, Vilnius, Lithuania), according to manufacturer's recommendations. The genotyping of IL-10: rs1800872, rs1800871, rs1800896 was carried out using the real-time polymerase chain reaction (PCR) method. Identification of all single-nucleotide polymorphisms (SNPs) was performed using TaqMan ® Genotyping assays (Thermo Fisher Scientific, Inc, Pleasanton, USA). The genotyping was performed using a "StepOnePlus" real-time PCR quantification system (Thermo Fisher Scientific, Singapore). Data of individual genotypes were received using the Allelic Discrimination program during the real-time PCR.
IL-10 serum levels were evaluated in 20 control subjects and 20 LSCC patients. The assay was performed using the Invitrogen ELISA Kit (Cat. No. KHC0101) for human IL-10, standard curve sensibility range=0-500 pg/ml, sensitivity <1 pg/ml, following the manufacturer's instructions, and analyzed on the Multiskan FC Microplate Photometer (Thermo Scientific, Waltham, MA) at 450 nm.
Quality control of genotyping. The repetitive analysis of 5% randomly chosen initial samples was accomplished for selected SNPs to confirm the same rate of genotypes from initial and repetitive genotyping.
Survival rate. The LSCC group data about the mortality rate, including the survival period after diagnosis of LSCC, were collected from the Lithuanian State Register of Death Cases and Their Causes, Institute of Hygiene, Vilnius.
Statistical analysis. The data about the study participants' demographic characteristics are presented as absolute numbers with percentages in brackets and compared between the control group subjects and LSCC group using the Chi-square test. The frequencies of all SNPs genotypes and alleles are presented in percentages. The Analysis of Hardy-Weinberg using a chi-square test was performed in the control group to compare the analyzed and expected frequencies of IL-10: rs1800872, rs1800871, rs1800896 SNPs. Chi-square test was used to compare the distribution of IL-10: rs1800872, rs1800871, rs1800896 SNPs in the LSCC, and control groups. Binomial logistic regression analysis with an adjusted odds ratio (OR) and its 95% confidence interval (95% CI) was utilized to evaluate the influence of IL-10: rs1800872, rs1800871, rs1800896 genotypes and alleles on LSCC development and the risk prediction for LSCC patients with these polymorphisms. The binomial logistic regression analysis results are  represented as genetic models: codominant, dominant, recessive, overdominant, and additive. Akaike Information Criterion (AIC) was used for the best genetic model selection. Only statistically significant variables are presented in the tables. Statistical analysis was performed using the SPSS/W 20.0 software (Statistical Package for the Social Sciences for Windows, Inc., Chicago, Illinois, USA). The findings were considered statistically significant when p<0.05.

Results
Hardy-Weinberg equilibrium (HWE) test was carried out between the control group subjects. The results demonstrated that genotypes of SNPs in the control group did not deviate from HWE (p>0.05) (Table II).
Association between IL-10 SNPs and LSCC. In the present study, we performed genotyping of IL-10: rs1800872, rs1800871, and rs1800896 SNPs in 833 subjects (300 LSCC patients and 533 control group subjects) and studied possible associations between selected SNPs and LSCC development. Table III lists the frequencies of selected SNPs genotype profiles between the control group and LSCC patients. However, we did not find any statistically significant differences between the control group and LSCC patients (Table III).
Interestingly, we noticed that genotypes at position -819 (rs1800871) were changed synchronously with -592 (rs1800872) position. To evaluate the impact of IL-10 rs1800871, rs1800872, and rs1800896 on the LSCC development, a binomial logistic regression was applied. However, our results did not show any statistically significant results.
Genotypic and allelic distribution according to the development of LSCC by stages. To analyze the possible genotypic differences to the progress of LSCC, we divided the LSCC group into two subgroups: the early-stage LSCC group (composed of I and II stages LSCC patients) and the advanced-stage LSCC group (formed of III and IV stages LSCC patients) and compared them with control group subjects and between each other.
We found statistically significant differences in the allele distribution of IL-10 rs1800871 and rs1800872 in the advanced-stage LSCC subgroup compared to the control group: the G allele at rs1800871 was less frequent in the LSCC group than in controls, while the A allele was statistically significantly more frequent in the advanced stage of LSCC subgroup than in the control group (accordingly: 70.6% vs. 81.6% and 29.4% vs. 18.4%; p=0.030). Respectively, the G allele at 1800872 was less common in the advanced-stage LSCC subgroup than in controls, although the T allele at rs1800872 was statistically significantly more frequent in the advanced stage of LSCC subgroup than in the control group (accordingly: 29.4% vs. 18 SNPs analysis between the early and advanced stages of LSCC revealed statistically significant differences in the genotypic distribution of IL-10 rs1800871 and rs1800872. We found out that the GG genotype at rs1800871 was more frequent in the early stage compared to the advanced stage of the LSCC subgroups (63.9% vs. 52.7%; p=0.049), while the AA genotype was less frequent in the early stage compared to the advanced stage of the LSCC subgroups (3.6% vs. 11.5%; p=0.008). Also, the G allele at rs1800871 was more frequent in the early stage LSCC subgroup compared to the advanced stage of LSCC patients (80.2% vs. 70.6%; p=0.007). The rs1800872 GG genotype was more frequent in the early stage compared to the advanced stage of the LSCC subgroup (63.9% vs. 52.7%; p=0.049), while the TT genotype was less frequent in the early stage compared to the advanced stage of the LSCC subgroup (3.6% vs. 11.5%, p=0.008). Also, the G allele at rs1800872 was more frequent in the early stage of LSCC comparing to the advanced stage of LSCC patients (80.2% vs. 70.6%, p=0.007) (Table IV).
Binomial logistic regression analysis was performed to evaluate the impact of selected IL-10 SNPs on LSCC development in the disease's advanced stage. It revealed that IL-10 rs1800871 AA genotype is associated with about 0.5fold decreased odds of LSCC development in the codominant and recessive models (OR=0.473; 95% CI=0.244-0.918; p=0.027 and OR=0.510; 95% CI=0.268-0.971; p=0.040, respectively). Also, each allele A is associated with 0.7-fold decreased odds of LSCC development in the additive model (OR=0.733; 95% CI=0.548-0.982; p=0.037) (Table V).
Genotypic and allelic distribution according to the metastasis to the neck lymph nodes. To better understand the development of LSCC, we divided LSCC patients group into two subgroups: patients with no metastasis to the neck lymph nodes (NMNLN) subgroup and patients with metastasis to the neck lymph nodes (WMNLN) subgroup. The subgroups were compared with the control group separately and between each other. After statistical analysis, it was revealed that the GG genotype at rs1800871 was less frequent in the WMNLN subgroup than in control subjects (44.8% vs. 60.2%; p=0.024), while the AA genotype was more frequent in the WMNLN subgroup than in control subjects (17.2% vs. 6.2%; p=0.002). The GG genotype at rs1800872 was less frequent in the WMNLN subgroup than in control subjects (44.8% vs. 60.2%; p=0.024), while the TT genotype was more frequent in the WMNLN subgroup than in control  (Table VI).
Comparing the NMNLN and WMNLN subgroups between each other, we identified that the genotypes and alleles at rs1800871, and respectively at rs1800872, were more frequent in the NMNLN subgroup than in the WMNLN subgroup (p=0.001).
The GG genotype at rs1800871 was less frequent in the WMNLN subgroup than in the NMNLN subgroup ( (Table VI).
After binomial logistic regression analysis, the statistically significant differences in the impact of IL-10 rs1800871 and rs1800872 on the WMNLN patients' subgroup were revealed. It was clarified that genotype AA at IL-10 rs1800871 was associated with 3.    Association between SNPs and five-year survival rate. The five-year survival rate of selected 300 LSCC patients was 65%. Analyzing LSCC patients' five-year survival rate and the genotype distribution, we did not find any statistically significant differences between the genotypic distribution of selected SNPs and the survival rate of the LSCC patients (Figures 2 and 3).

Discussion
Recent studies suggested that polymorphisms in pro-and antiinflammatory cytokines promote different diseases, including carcinogenesis (34,61). It has been indicated that immune cytokines such as IL-10 and inflammation play an essential role during carcinogenesis by inducing angiogenesis, immune suppression, and cell metastasis (62,63). According to that, the molecular effect of IL-10 gene promoter SNPs is not precisely understood, and the results remain conflicting (17,34). For this study, we selected three IL-10: rs1800871, rs1800872, and rs1800896 SNPs to identify their possible role in LSCC development and determine associations of genetic distribution of selected SNPs and patients' survival rate; IL-10 serum levels were measured as well. To the best of for knowledge, this is the first study analyzing the impact of before mentioned SNPs on the survival rate of LSCC patients. The immunosuppressive characteristics of IL-10, corresponding with the high levels of IL-10 identified in tumor tissues are well determined in the literature (20). In tumor cells, IL-10 inhibits the expression of major histocompatibility complex (MHC) class II in the antigenpresenting cell, MHC class I expression, and natural killer group 2D (NKG2D) ligand expression (12,64). Conditions as carcinogenesis trigger upregulation of NKG2D, inducing the recognition of the NKG2D receptor (20). This results in the activation of signaling pathways, leading to stimulation of antitumor functions of NK cells (65). Moreover, blockage of NKG2D pathways presents increased susceptibility of induced malignant transformations of mice (66). Hence, IL-10 induced inhibition of NKG2D ligand expression influences tumor immune surveillance (64). However, during the down-regulation of MHC class I expression, it stimulates human leukocyte antigen-G (HLA-G) molecules expression that is usually observed in pathological mechanisms, as carcinogenesis (20). HLA-G is associated with tumor immune escape, metastasis, and poor prognosis (67). Referring to the critical role of IL-10 in immune response regulation in malignant proceses, a huge interest has been focused on the identification of naturally occurring gene polymorphisms in IL-10 and flanking regions (17). It is suggested that alterations in IL-10 promoter activity might be one of the mechanisms provoking IL-10 production modifications (68). However, the molecular effects of IL-10 rs1800871, rs1800872, and rs1800896 polymorphisms genotypic distribution in malignancy pathogenesis remain not completely understood (17). According to previous studies and focusing on histological aspects that LSCC belongs to squamous cell carcinomas, significant associations between IL-10 rs1800871, rs1800872, and rs1800896 polymorphisms with lung, thyroid, gastric, cervical, and prostate cancer are identified (55,(57)(58)(59). Several studies announced that IL-10 has pro-tumorigenic and anti-tumorigenic effects (12,(22)(23)(24)(25)(26). For this reason, numerous studies evaluating the association between polymorphisms of IL-10 as potential biomarkers and different types of cancer have been reported. From all IL-10 promoter polymorphisms, the rs1800896 (-1082A/G) is the one most widely studied. A meta-analysis by Khan et al. analyzed twenty-three case-control relevant studies of IL-10 promoter-1082 A>G (rs1800896) polymorphism significance in the increased risk of carcinogenesis involving 4753 cases (576 cases of HNSCC) and 6,086 control subjects. It concluded that IL-10 promoter-1082 polymorphism is not associated under a heterozygous (GG vs. AG) model with cancer developing risk (32). However, in a meta-analysis conducted by Niu et al., a systematic review of 9 relevant studies including 2,258 cases of only HNSCC and 2,887 control subjects (623 cases of nasopharyngeal carcinoma, 1,357 -oral SCC, and 278 cases of all HNSCC without specific location) was performed and suggested that IL-10 -1082A>G (rs1800896) polymorphism might be associated with increased risk of HNSCC development (34). In the present study, we identified argumentative results: all genotypes of IL-10 rs1800896 (TT, TC, and CC) were very similar in the LSCC and control groups (23.3 vs. 27.8%, 54.3 vs. 50.5%, and 22.3 vs. 21,8%; p=0.364) demonstrating that there was no statistically significant difference between IL-10 rs1800896 in both investigated groups.
The impact of IL-10 promoter polymorphisms on vocal fold leukoplakia and LSCC was investigated by Zhou et al. (35). The study involved 146 patients with LSCC (97.3% male and 2.7% female) and 119 control subjects (95.8% male and 4.2% female) and revealed that the -592A/C (rs1800872), -819T/C (rs1800871) and -1082A/G (rs1800896) polymorphisms and elevated IL-10 plasma levels are associated with an increased risk of LSCC and vocal fold leukoplakia development. Moreover, higher IL-10 plasma levels and higher OR values of selected SNPs were identified in the advanced stages of LSCC patients and the patients with metastasis to the neck lymph nodes. In the present study, our LSCC group was two times larger compared to Zhou et al. study and was composed of 300 patients. Furthermore, our control group was more than fourtimes larger than in the previously mentioned study and was formed of 533 subjects adjusted by age and sex. Nevertheless, we did not identify any statistically significant differences according to genotype distribution after the binomial logistic regression analysis was accomplished. However, results of our study indicated associations of IL-  suggesting that increasing the study population could lead to more statistically reliable results (35). Our study results suggest that IL-10 rs1800871 and rs1800872 polymorphisms play a significant role only in the advanced stages of LSCC development compared to healthy controls and comparing them to the early stage of the LSCC group. Hence, it leads us to the presumption that LSCC might develop slower in patients who carry the AA genotype at rs1800871 and the TT at rs1800872 with about 0.5-fold decreased odds of LSCC development under the codominant and recessive (p=0.027, p=0.040 and p=0.027, p=0.040, respectively) models. Also, each allele A at rs1800871 and as well each T allele at rs1800872 is associated with 0.7-fold decreased odds of LSCC development under the additive model (p=0.037 and p=0.037, respectively).
Makni et al. investigated possible IL-10 promoter gene variants in 194 HNSCC patients -137 nasopharyngeal cancer (NPC), 57 LSCC and in 263 control subjects and found out that all three SNPs were significant in NPC development, but only IL-10 rs1800871 and AT haplotype may be used for identification of higher risk of LSCC (33). In this study, the synchronous changing of genotypes at position -819 (rs1800871) with -592 (rs1800872) position was not identified. However, this was noticed both in our study and in the study designed by Zhou et al. (35). Also, in the study performed by Makni et al. the dominant sex in the control group was female: 171 (65.8%) female from 263 control subjects; respectively 6 (10.5%) females from 57 LSCC patients. In the present study, we selected and adjusted   both investigated groups for age and sex (LSCC group: 95.7% males and 4.3% females; control group: 95.3% males and 4.7% females; p=0.813), with a reference that reported predominant sex of LSCC is male (3,33).
In the literature, huge interest has been given on elevated IL-10 plasma levels as a significant marker and adverse prognostic factor in HNSCC (35,69,70). Increased IL-10 serum levels of HNSCC patients have been observed at advanced stages of the disease (19.3±1.1.2 pg/ml, n=29 vs. 9.1±0.9 pg/ml, n=31, p<0.0005) and in patients with nodal metastasis (19.1±1.6 pg/ml, n=29 vs. 9.4±1.4 pg/ml, n=31) (26). However, some researchers suggest that serum cytokine levels cannot be used as potential biomarker in LSCC (71). In the present study, the IL-10 plasma level was slightly lower in LSCC patients than in control subjects; however, no statistically significant differences were observed (10.3±15.99 pg/ml, n=20 vs. 14.6±2.61 pg/ml, n=20; p=0.579). The median IL-10 plasma level in HNSCC patients was also lower comparing to the control group in the study carried out by Makni  Regarding the poor prognosis of HNSCC and urgently needed prognostic biomarkers to diagnose these disorders earlier, the present study was the first to investigate possible associations of IL-10: rs1800871, rs1800872, rs1800896 SNPs with LSCC patients' survival rate. As a result, we identified that the distribution of selected SNPs genotypes had no significant impact on LSCC patients' survival rate (respectively, p=0.952; p=0.952, and p=0.991).
This study's strength was the inclusion of a large study population (833 subjects in total), pure LSCC patients group, and precise selection of investigated groups -both groups were adjusted for age and sex.
Generally, the results of the present study are in concordance with the data of the literature. Some discrepancies between our results and literature data can be possibly explained by the significant diversity of malignances traditionally united under the umbrella of the HNSCC term because this term includes malignant tumors of different localizations (oral, pharyngeal, laryngeal regions) as well as of different etiology, biological and clinical behavior. Therefore, pooling all cancer types together may mask the possible significant associations of IL-10: rs1800871, rs1800872, and rs1800896 with individual cancer types. In the present study, a pure cohort of LSCC patients was investigated. Results revealed no statistically significant differences in the distribution of IL-10: rs1800871, rs1800872, rs1800896 variants between the LSCC patients and control groups. Presumably, these findings reflect the biological and clinical peculiarities of the LSCC tumor.
To the best of our knowledge, this is the first report that associates the role of IL-10: rs1800871, rs1800872, rs1800896 SNPs, and development of LSCC in large, pure and homogenous LSCC patients cohort and age-and sex-matched control subjects. This particuliarity allowed us to accomplish a precise analysis of associations between selected SNPs and LSCC development, with a concrete tumor in one anatomical region of head and neck. It is well known from clinical practice that LSCC (notably, the glottis cancer) is characterized as less aggressive in accordance with local spreading and low metastatic rate compared to other HNSCC (72,73). Thus, not identified discrepancies of IL-10: rs1800871, rs1800872, rs1800896 variants distribution between this low aggressive LSCC patients' and sex-matched control subjects is comprehensible. However, a statistically significant association between IL-10 rs1800871 and rs1800872 gene polymorphisms and LSCC in patients with the advanced stage of LSCC and patients with metastasis to the neck lymph nodes suggests the importance of these SNPs in the development of this peculiar tumor.
Certain limitations of the present study, as involvement of environmental factors should be deliberated; smoking and alcohol consumption habits were not investigated. However, this is a targeted task for future investigation.

Conclusion
The present study indicated a significant association between IL-10 rs1800871 and rs1800872 gene polymorphisms with an increased risk for LSCC development in advanced stages of LSCC and patients with metastasis to the neck lymph nodes. These results might be useful for selecting candidate prognostic factors in identifying and elaborating novel targeted therapies and immunotherapy strategies and could improve LSCC patients' survival rate. However, further studies involving risk factors such as smoking and alcohol consumption are necessary to define these SNPs' role in LSCC carcinogenesis.

Conflicts of Interest
The Authors declare no potential conflicts of interest associated with this study.