Disease Markers
Volume 2019 (2019), Article ID 6523837, 10 pages
https://doi.org/10.1155/2019/6523837

Review Article

Genetic Association between NFKBIA and NFKB1 Gene Polymorphisms and the Susceptibility to Head and Neck Cancer: A Meta-Analysis

Lin Li and Zhong-Ti Zhang

VIP Department, School of Stomatology, China Medical University, Shenyang 110002, China

Correspondence should be addressed to Zhong-Ti Zhang; zhangzhongti@163.com

Received 2 May 2019; Revised 21 August 2019; Accepted 1 September 2019; Published 12 September 2019

Academic Editor: Małgorzata Knaś

Copyright © 2019 Lin Li and Zhong-Ti Zhang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Background. The role of the NFKB1 gene rs28362491 polymorphism and NFKBIA gene rs2233406 polymorphism in the development of head and neck cancer (HNC) remains controversial. This meta-analysis was performed to assess the relationship between the gene polymorphisms and HNC quantitatively. Methods. PubMed, Embase, Web of Science, WanFang Data, and China National Knowledge databases were used to search for eligible articles. The relationship was evaluated by STATA 11.0. Results. Eight eligible articles were included in our study. Nine case-control studies from the eight included articles were correlated with rs28362491 polymorphism. Four articles were related to rs2233406 polymorphism. Overall, a significant correlation was observed between the rs28362491 polymorphism and a decreased risk of HNCs (, for DD vs. II; , for DD vs. DI+II). In subgroup analyses, the rs28362491 polymorphism was associated with the risk of nasopharyngeal carcinoma (NC), but not with oral cancer (OC). In addition, no statistical correlation was found between the polymorphism of rs2233406 and HNCs. Conclusion. rs28362491 polymorphism was significantly associated with the risk of HNCs, especially with NC. Additionally, our results showed that no association was discovered between rs2233406 polymorphism and HNCs.

1. Introduction

Head and neck cancer (HNC) is the sixth most common cancer [1], which arises in the oral cavity, pharynx, nasal cavity, and larynx. Appropriately 650,000 new cases of HNCs and 350,000 deaths occur each year worldwide according to the statistics [2]. As a cancer with a high morbidity rate and low survival rate, it reduces patient quality of life [3–5]. Although the primary etiology of HNC is related to human papillomavirus (HPV), smoking, alcohol, and Epstein-Barr virus (EBV) infection [6–10], only a small number of subjects exposed to these risk factors will suffer from a HNC, which suggests that individual genetic susceptibility might also be strongly correlated with the development of HNCs [3, 11, 12].

Nuclear factor-κB (NF-κB), which was identified initially in 1986, acts as an essential transcription factor associated with cell survival, cell growth, cell replication, and cell apoptosis [13–15]. NF-κB also plays an important role in the development of carcinogenesis and tumor cell’s resistance to radiotherapy and chemotherapy [16]. The NF-κB family consists of five members: p65/Rel A, p50/p105, p52/p100, c-Rel, and Rel B. The heterodimer of the p65/Rel A and p50/p105 subunits, encoded by the NFKB2 and NFKB1 genes, respectively, is the most common form of NF-κB. Moreover, NF-κB activity would be influenced by NF-κB inhibitor alpha (NFKBIA) [17–19]. Inappropriate expression of NF-κB has been reported in certain human tumors, such as colorectal cancer, breast cancer, colon cancer, multiple myeloma, melanoma, prostate cancer, and Hodgkin’s lymphoma [20–24].

The NFKB1 and NFKBIA genes, located on chromosome 4q24 and 14q13, respectively, are associated with the development of cancers, including oral cancer (OC) [25], gastric cancer [26], colorectal cancer [27], and Hodgkin’s lymphoma [28]. In recent years, it has been reported that some polymorphisms of the NFKB1 gene (rs28362491) and NFKBIA gene (rs2233406) might be associated with HNCs [25, 29–35]. However, the results are conflicting. It might be because of small sample sizes, publication bias, and different ethnic backgrounds. Therefore, it is necessary to perform a comprehensive meta-analysis with a large sample size and high statistical power to assess the association between the polymorphisms of the NFKB1 and NFKBIA genes and HNCs. Additionally, trial sequential analysis (TSA) was also conducted to minimize random errors and strengthen the reliability of our conclusions.

2. Materials and Methods

2.1. Selection of Relevant Studies

The meta-analysis was designed and implemented according to the recommendations of the Preferred Reporting Items for Systematic Reviews and Meta-Analyses statement [36].

PubMed, Embase, Web of Science, WanFang Data, and China National Knowledge databases were searched comprehensively prior to May 1, 2019 with the following terms: “Nuclear factor-κB”, or “Nuclear factor kappa B”, or “NF-κB”, or “NFKB1”, or “rs28362491”, or “NFKBIA”, or “Nuclear factor kappa B inhibitor”, or “rs2233406” AND “carcinoma”, or “cancer”, or “neoplasm”, or “tumor”, or “carcinogenesis” AND “polymorphism”, or “variants”, or “variant”, or “mutation”, or “mutations”, or “polymorphisms” AND “head and neck”, or “HNC”, or “oral”, or “oral cavity”, or “pharyngeal”, or “laryngeal”, or “nasopharyngeal”, or “oropharyngeal”, or “laryngopharyngeal”, or “hypopharyngeal”. Finally, references cited by all the included studies were scanned to identify relevant articles.

2.2. Inclusion and Exclusion Criteria

No limitations were imposed on language or publication date. Articles were enrolled if they met all the following criteria: (1) the paper was case-controlled; (2) the paper assessed the risk of HNC and the rs28362491 polymorphism of the NFKB1 gene and/or the rs2233406 polymorphism of the NFKBIA gene; and (3) information provided by the paper was sufficient to estimate the correlation. Articles which satisfied one of the following criteria were excluded: (1) duplicated papers, reviews, editorials, case reports, commentaries, and nonhuman studies; (2) full-text was not found; and (3) no sufficient data were reported.

2.3. Data Extraction and Quality Assessment

Two writers collected the information of each enrolled article according to the inclusion and exclusion criteria independently. First author, year of publication, country, ethnicity, genotyping method, source of control, distributions of genotype, and cancer type were fetched. Quality assessment of the enrolled papers in the analysis was conducted by the Newcastle-Ottawa Scale (NOS) criteria [37].

2.4. Trial Sequential Analysis (TSA)

Considering that meta-analysis might result in false positive (type-1 error) or false negative (type-2 error) results due to the sparse data [38], TSA was applied to minimize random errors and increase the reliability of our conclusions. In the present study, TSA was conducted with an overall 5% type-1 error risk and 20% type-2 error risk [39]. According to the required information size (RIS) and risk for type-1 error and type-2 error, TSA monitoring boundaries were constructed. If the -curve crossed with the TSA monitoring boundary before the RIS was reached or total sample size included in the study reached the RIS, a robust conclusion might have been affirmed and further study would be unnecessary.

2.5. False-Positive Report Probability (FPRP) Analysis

We calculated the false-positive report probability (FPRP) [40] to assess the significant results. We set 0.5 as the FPRP threshold and assigned a prior probability of 0.1 to detect an odds ratio (OR) of 1.50 for an association with genotypes under investigation. Only a significant finding with a was considered a noteworthy result.

2.6. Statistical Method

Statistical analysis was performed using STATA 11.0 (StataCorp LP, College Station, TX, United States). The strength of correlation between the gene polymorphisms and HNCs was evaluated through pooled odds ratios (ORs) along with the corresponding 95% confidence intervals (CIs). The ORs and CIs of each gene polymorphism were calculated for the allelic, homozygote, heterozygote, dominant, and recessive genetic models of each article, respectively. The combined ORs were calculated by -test, and a value of was regarded as statistically significant.

The between-article heterogeneity was estimated by Cochran’s -test and statistics. If the value in the and , the Mantel-Haenszel fixed-effects model was used to estimate the pooled ORs. Otherwise, the DerSimonian and Laird random-effects model was performed to evaluate the correlation between the polymorphisms of the NFKB1 gene (rs28362491), NFKBIA gene (rs2233406), and HNCs. Hardy-Weinberg equilibrium (HWE) in the controls was calculated by the chi-square test [41]. In order to explore the possible heterogeneity among the articles, subgroup analyses according to the genotyping method and tumor type were performed. Sensitivity analysis was conducted by deleting a single article every time in order to assess the reliability of the overall results. Publication bias was detected by Begg’s test and Egger’s linear regression.

3. Results

3.1. Article Characteristics

Eight articles [25, 29–35] involving 4434 cases and 4913 controls were finally enrolled in the meta-analysis. The process of study selection is demonstrated in Figure 1.

Figure 1: Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) flow diagram for a study selection process.

Nine case-control studies from eight included articles [25, 29–35] were correlated with rs28362491 polymorphism. Among the eight studies, four studies [25, 30, 34, 35] were related to rs2233406 polymorphism. Overall, the meta-analysis included three OC articles, three nasopharyngeal carcinoma (NC) articles, and two HNC articles. All studies showed that the genotype distribution in the controls was conformed to HWE. The detailed information of these enrolled papers is demonstrated in Table 1.

Table 1: Characteristics of articles enrolled in the meta-analysis.

FA (year)
NOS
Ethnicity
Genotyping method
SC
Cases
Controls
Cancer type
HWE

rs28362491

ins/ins
ins/del
del/del
ins/ins
ins/del
del/del

Lin (2006)
6
Asian
PCR-PAGE
HB
59
103
50
43
100
58
OC
0.993

Lehnerdt (2008)
7
Caucasian
Germany
PCR-RFLP
HB
132
179
53
118
141
48
HNC
0.586

Zhou (2009)
8
Asian
PCR-PAGE
HB
74
67
22
71
90
42
NC
0.177

Lin (2012)
8
Asian
TaqMan
HB
116
246
100
81
271
168
OC
0.099

Liu (2012)
7
Asian
PCR
PB
269
467
170
280
433
193
NC
0.289

Liu (2015)
7
Asian
TaqMan
HB
316
438
152
336
512
224
NC
0.262

Liu (2015)
7
Asian
TaqMan
HB
236
331
117
274
438
195
NC
0.169

Chen (2018)
8
Asian
MassARRAY
HB
124
197
100
163
230
90
OC
0.577

Gupta (2017)
7
Asian
PCR-RFLP
HB
132
162
18
144
146
22
HNC
0.064

rs2233406

CC
CT
TT
CC
CT
TT

Lin (2012)
8
Asian
TaqMan
HB
351
101
10
438
78
4
OC
0.797

Liu (2012)
7
Asian
PCR
PB
706
185
15
694
201
11
NC
0.402

Liu (2015)
7
Asian
TaqMan
HB
701
188
17
813
244
15
NC
0.492

Chen (2018)
8
Asian
MassARRAY
HB
308
108
9
365
110
8
OC
0.931

FA: first author; SC: source of control; OC: oral cancer; NC: nasopharyngeal carcinoma; HNC: head and neck cancer; HB: hospital-based study; PB: population-based study; HWE: Hardy-Weinberg equilibrium; PCR-RFLP: polymerase chain reaction-restriction fragment length polymorphism; PCR: polymerase chain reaction; PCR-PAGE: polymerase chain reaction-polyacrylamide gel electrophoresis.

3.2. Meta-Analysis Results

The results of the overall analyses and subgroup analyses for gene polymorphisms (rs28362491 and rs2233406) with HNCs are shown in Table 2.

Table 2: Results of overall and subgroup analyses for gene polymorphisms (rs28362491 and rs2233406).

rs28362491
No
D versus I
DD versus II
DI versus II
DD+DI versus II
DD versus DI+II

OR
95% CI

OR
(95% CI)

OR
(95% CI)

OR
(95% CI)

OR
(95% CI)

Overall
9
0.89
0.79-1.00
0.059
0.76
0.60-0.97
0.026
0.95
0.82-1.08
0.419
0.89
0.76-1.05
0.165
0.80
0.68-0.95
0.010

PCR-RFLP
2
1.04
0.89-1.22
0.630
0.96
0.65-1.40
0.813
1.17
0.93-1.48
0.179
1.13
0.91-1.41
0.274
0.88
0.62-1.26
0.495

PCR-PAGE
2
0.74
0.61-0.91
0.004
0.57
0.38-0.85
0.007
0.73
0.53-1.02
0.063
0.67
0.49-0.92
0.012
0.69
0.49-0.98
0.039

TaqMan
3
0.79
0.69-0.91
0.001
0.61
0.45-0.83
0.002
0.83
0.69-1.00
0.047
0.75
0.60-0.94
0.014
0.70
0.60-0.83
<0.001

OC
3
0.86
0.58-1.27
0.457
0.73
0.32-1.65
0.446
0.82
0.56-1.20
0.314
0.79
0.46-1.33
0.369
0.84
0.49-1.46
0.543

HNC
2
1.04
0.88-1.22
0.630
0.96
0.65-1.40
0.813
1.17
0.93-1.48
0.179
1.13
0.91-1.41
0.274
0.88
0.62-1.26
0.495

NC
4
0.87
0.78-0.96
0.007
0.75
0.62-0.90
0.002
0.94
0.81-1.09
0.400
0.87
0.74-1.02
0.096
0.78
0.68-0.89
<0.001

rs2233406
No
T versus C
TT versus CC
TC versus CC
TT+TC versus CC
TT versus CC+TC

OR
95% CI

OR
(95% CI)

OR
(95% CI)

OR
(95% CI)

OR
(95% CI)

Overall
4
1.13
0.90-1.42
0.296
1.51
0.99-2.30
0.058
1.08
0.85-1.39
0.523
1.12
0.86-1.44
0.399
1.50
0.98-2.29
0.061

TaqMan
2
1.25
0.72-2.16
0.431
1.68
0.93-3.04
0.084
1.19
0.66-2.12
0.565
1.23
0.68-2.23
0.497
1.67
0.92-3.02
0.089

OC
2
1.39
0.97-1.98
0.074
1.92
0.93-3.98
0.080
1.36
0.99-1.88
0.059
1.40
0.98-2.00
0.062
1.81
0.87-3.75
0.111

NC
2
0.96
0.83-1.10
0.520
1.33
0.79-2.24
0.290
0.90
0.77-1.05
0.180
0.92
0.79-1.07
0.296
1.36
0.81-2.29
0.251

3.2.1. Analysis for rs28362491 Polymorphism

For the existence of significant heterogeneity in the -test (value of in or ), the random-effects model was performed to analyze the correlation between rs28362491 polymorphism and HNCs under all genetic models. Overall, rs28362491 polymorphism was significantly correlated with a decreased risk of HNCs under homozygote and recessive genetic models (, for DD vs. II, Figure 2; , for DD vs. DI+II). In subgroup analyses, a significant association was discovered for the polymerase chain reaction-polyacrylamide gel electrophoresis (PCR-PAGE) genotyping method subgroup (, for D vs. I; , for DD vs. II; , for DD+DI vs. II;, for DD vs. DI+II), TaqMan genotyping method subgroup (, for D vs. I; , for DD vs. II; , for DI vs. II; , for DD+DI vs. II; , for DD vs. DI+II), and NC (, for D vs. I; , for DD vs. II; , for DD vs. DI+II).

Figure 2: Forest plot for the association of rs28362491 polymorphism and HNC risk under a homozygote genetic model.

3.2.2. Analysis for rs2233406 Polymorphism

With respect to rs2233406 polymorphism, the random-effects model was conducted under all genetic models except the recessive model because the statistical heterogeneity between the articles was substantial (value of in or ). We discovered that there was no relationship between the rs2233406 polymorphism and the risk of HNCs under all genetic models (, for T vs. C; , for TT vs. CC; , for TC vs. CC; , for TT+TC vs. CC; , for TT vs. CC+TC). In the stratified analyses by genotyping method and tumor type, no association was observed in all genetic models ().

3.3. Sensitivity Analyses and Publication Bias

In the sensitivity analyses, no substantive change was discovered in the combined ORs after excluding one paper at a time. Egger linear regression tests and Begg’s funnel plots were conducted to evaluate publication bias. In all the genetic models, no remarkable publication bias was found by the value in the Egger test (D vs. I: ; DD vs. II: ; DI vs. II: ; DD+DI vs. II: ; DD vs. DI+II: ; T vs. C: ; TT vs. CC: ; TC vs. CC: ; TT+TC vs. CC: ; TT vs. TC+CC: ) and Begg’s funnel plot for rs28362491 and rs2233406 polymorphisms.

3.4. TSA and FPRP Analysis Results

Nine trials were used to evaluate the association between the rs28362491 gene polymorphism and the HNC susceptibility. The results of TSA analysis (taking the homozygous model data for example) showed that the cumulative -curve had not crossed the trial monitoring boundary before the RIS was reached (Figure 3), indicating that the cumulative evidence might be insufficient and further studies are needed to strengthen the conclusion. However, the total sample size of another four genetic models included in the study all reached the RIS, suggesting that the cumulative evidence might be sufficient and further studies would be unnecessary (Figure 4 for the recessive model in rs28362491 gene polymorphism). As for rs2233406, the results of TSA analysis demonstrated that the cumulative -curve had not crossed the trial monitoring boundary before the RIS was reached under all the genetic models, indicating that further trials are necessary to evaluate the correlation between rs2233406 polymorphism and HNCs (figures were not shown).

Figure 3: The required information size to demonstrate the association of rs28362491 polymorphism and HNC susceptibility under a homozygote genetic model.

Figure 4: The required information size to demonstrate the association of rs28362491 polymorphism and HNC susceptibility under a recessive genetic model.

The FPRP values for all the discovered significant results are demonstrated in Table 3. For a prior probability of 0.1, the FPRP values were all less than 0.50 in the significant findings, indicating that these significant correlations were noteworthy.

Table 3: False-positive report probability values for associations between the rs28362491 polymorphism and the HNC risk.

Variables
OR (95% CI)
a
Powerb
Prior probability

0.25
0.1
0.01
0.001
0.0001

D versus I

PCR-PAGE
0.74 (0.61-0.91)
0.004
0.839
0.015
0.044
0.338
0.837
0.981

TaqMan
0.79 (0.69-0.91)
0.001
0.991
0.003
0.010
0.098
0.523
0.916

NC
0.87 (0.78-0.96)
0.006
1.000
0.016
0.048
0.355
0.847
0.982

DD versus II

Overall
0.76 (0.60-0.97)
0.027
0.854
0.088
0.225
0.761
0.970
0.997

PCR-PAGE
0.57 (0.38-0.85)
0.006
0.221
0.073
0.192
0.723
0.963
0.996

TaqMan
0.61 (0.45-0.83)
0.002
0.286
0.017
0.050
0.364
0.853
0.983

NC
0.75 (0.62-0.90)
0.002
0.897
0.007
0.020
0.180
0.688
0.957

DI versus II

TaqMan
0.83 (0.69-1.00)
0.050
0.989
0.132
0.313
0.833
0.981
0.998

DD+DI versus II

PCR-PAGE
0.67 (0.49-0.92)
0.013
0.512
0.072
0.190
0.720
0.963
0.996

TaqMan
0.75 (0.60-0.94)
0.013
0.847
0.042
0.117
0.594
0.937
0.993

DD versus DI+II

Overall
0.80 (0.68-0.95)
0.011
0.981
0.032
0.091
0.524
0.918
0.991

PCR-PAGE
0.69 (0.49-0.98)
0.038
0.576
0.166
0.374
0.868
0.985
0.998

TaqMan
0.70 (0.60-0.83)
0.00004
0.713
<0.001
0.001
0.006
0.054
0.363

NC
0.78 (0.68-0.89)
0.0002
0.990
0.001
0.002
0.022
0.184
0.693

aChi-square test was adopted to calculate the genotype frequency distributions. bStatistical power was calculated using the number of observations in the subgroup and the OR and values in this table.

4. Discussion

NF-κB plays an essential role in immune response, cell apoptosis, cell proliferation, and the development of cancers [42–44]. Moreover, many inflammatory cytokines, including IL-6, IL-8, and TNF-α would be influenced by NF-κB in regulating their biological effects. Therefore, an abnormality in NF-κB function would disturb these various biological behaviors and eventually lead to tumorigenesis [45].

The -94 ins/del ATTG (rs28362491) promoter polymorphism of the NFKB1 gene encoding p50/p105 subunit of the NF-κB family could modulate the transcription, production, and function of the p50/p105 subunit [45, 46]. NF-κB inhibitor I kappa B (IκB) would inhibit the function of NF-κB by binding to it. NF-κB inhibitor alpha (NFKBIA), which is the most common protein of the IκB family, is encoded by the NFKBIA gene. As a functional polymorphism of the NFKBIA gene, the -826 C/T (rs2233406) polymorphism might play an important role in influencing the function of NF-κB. Recently, various articles have explored the relationship between rs28362491 and rs2233406 polymorphisms and HNCs [25, 29–35]. However, the results were inconclusive.

To the best of our knowledge, this is the first meta-analysis to evaluate the association between rs28362491 and rs2233406 polymorphisms and HNCs. In the present study, nine case-control studies from eight eligible papers were enrolled to assess the correlation between -94 ins/del ATTG polymorphism and HNCs. Overall, the results showed that -94 ins/del ATTG (rs28362491) polymorphism was significantly correlated with a decreased risk of HNCs under homozygote and recessive genetic models (, for DD vs. II, Figure 2; , for DD vs. DI+II). Although the mechanism remains unknown, considering the essential role of NF-κB in the development of cancer, we hypothesize that the possible mechanisms are as follows. NF-κB plays a critical role in various biological pathways, which brings us to the idea that ATTG deletion results in the loss of binding to nuclear proteins, which causes lower promoter activity [47]. The promoter sequence with -94 del allele leads to decreased transcriptional activity and thereby results in decreased p50/p105 expression levels, which plays an essential role in transactivating antiapoptosis genes and restraining cell apoptosis, thereby promoting cellular proliferation. Besides, compared to the -94 ins allele carriers, the p50 in -94 del carriers could form lesser heterodimers with p65 to regulate the inflammatory pathway. Therefore, the -94 del allele might act as a protective role in HNC risk. However, no association between the rs2233406 polymorphism and the risk of HNCs was observed.

In the subgroup analyses of rs28362491 polymorphism based on the genotyping method, a statistically significant correlation was discovered in the PCR-PAGE genotyping method subgroup under allelic, homozygote, dominant, and recessive genetic models and in the TaqMan genotyping method subgroup under all the genetic models, but not in any other genotyping method subgroups. This might be because various genotyping methods would influence the relationship, indicating that it is necessary to confirm a genotyping method with a high specificity and sensitivity to increase the reliability of the results.

Since HNCs include OC, pharyngeal cancer, NC, and laryngeal neoplasm, further stratification analysis by tumor type was performed. When stratified by tumor type, we found that the rs28362491 polymorphism was significantly correlated with a decreased risk of NC under allelic, homozygote, and recessive genetic models. No association was discovered between the gene polymorphism and OC under all genetic models. This might be caused by various microenvironments, because the same gene could play an individual role in different tumor sites [48]. Surprisingly, our results were different from the results presented by four previous meta-analyses [49–52], which suggested that a significant correlation was seen between the rs28362491 polymorphism and OC. These contradictory results might be because all previous analyses only included two studies associated with OC. Therefore, the results suggested by the previous meta-analyses with limited sample sizes were unreliable.

In the subgroup analyses of rs2233406 polymorphism according to the genotyping method and tumor type, there was no relationship between the rs2233406 polymorphism and HNCs.

However, our study has some inevitable limitations. Firstly, some potential articles that have not been published were not enrolled in the present study, so a publication bias might exist. Secondly, our meta-analysis had a relatively small sample size in each subgroup, so the results of the subgroup analyses might not have enough power to identify the association. Thirdly, the environmental factors, such as smoking and alcohol, also play an essential role in the development of HNCs. Unfortunately, subgroup analyses according to smoking or alcohol consumption could not be conducted since there were no sufficient relevant data from most of the enrolled studies. Finally, the results of TSA analysis (taking the homozygous model data for example) showed that the cumulative evidence might be insufficient and further studies would be needed to strengthen the conclusion of rs2233406 polymorphism and the homozygous model in rs28362491 polymorphism. Therefore, further studies with a larger sample size and more information are required to observe the role of rs28362491 and rs2233406 polymorphisms in HNCs.

Despite these shortcomings, our meta-analysis has several highlights. To our knowledge, this present meta-analysis is the first one to conclude a relationship between rs28362491 polymorphism and HNCs and discover that rs28362491 polymorphism is associated with the development of HNCs, especially NC, although there were four meta-analyses [49–52] which focused on the association between the rs28362491 polymorphism and OC. Compared with these, our meta-analysis has the following advantages: Firstly, nine case-control studies from eight eligible papers were enrolled to assess the correlation between -94 ins/del ATTG polymorphism and HNCs; all the previous analyses only included two studies associated with OC. Secondly, the results presented by the four previous meta-analyses [49–52] were different from our results, which suggested that no significant correlation was seen between the rs28362491 polymorphism and OC. These contradictory results might be because all the previous analyses only included two studies associated with OC. Therefore, the results suggested by our meta-analyses with a higher sample size were more reliable than these previous meta-analyses. Thirdly, our study demonstrated that there was no association between the rs2233406 polymorphism and the risk of HNCs, which was not mentioned in previous articles. Finally, in order to make our conclusion more credible, publication bias analysis and sensitivity analysis were performed. Egger’s linear regression tests and Begg’s funnel plots demonstrated that no obvious publication bias was found. The sensitivity analysis suggested that our conclusions are reliable. Additionally, we performed FPRP analysis, and the results showed that the significant findings of our study are robust. Besides, the meta-analysis clearly shows that the rs28362491 gene polymorphism could be generally applied as a novel prognostic biomarker for HNCs, especially for NC, and it might play a protective role in populations with HNCs, which provides guidance for precise prediction of prognosis and/or therapeutic response.

5. Conclusions

In conclusion, this meta-analysis explored that the risk of HNCs was significantly correlated with rs28362491 polymorphism, but not with rs2233406 polymorphism. Moreover, the significant correlation between rs28362491 polymorphism and the susceptibility to NC was identified for the first time. However, more gene-environment and gene-gene interaction papers with larger sample sizes should be addressed to evaluate the relationship between rs28362491 and rs2233406 polymorphisms and HNCs.

Conflicts of Interest

The authors declare that they have no conflict of interest.

Acknowledgments

The present article was funded with the following grants: the Clinical Diagnosis and Treatment Capacity Building Project of Liaoning Province (grant no. LNCCC-D17-2015) and the Natural Science Foundation of Liaoning Province, China (grant no. 2015020558).

References

A. Jemal, F. Bray, M. M. Center, J. Ferlay, E. Ward, and D. Forman, “Global cancer statistics,” CA: A Cancer Journal for Clinicians, vol. 61, no. 2, pp. 69–90, 2011.

A. Argiris, M. V. Karamouzis, D. Raben, and R. L. Ferris, “Head and neck cancer,” The Lancet, vol. 371, no. 9625, pp. 1695–1709, 2008.

M. Brunotto, A. M. Zarate, A. Bono, J. L. Barra, and S. Berra, “Risk genes in head and neck cancer: a systematic review and meta-analysis of last 5 years,” Oral Oncology, vol. 50, no. 3, pp. 178–188, 2014.

R. C. Chai, D. Lambie, M. Verma, and C. Punyadeera, “Current trends in the etiology and diagnosis of HPV-related head and neck cancers,” Cancer Medicine, vol. 4, no. 4, pp. 596–607, 2015.

Global Burden of Disease Study 2013 Collaborators, “Global, regional, and national incidence, prevalence, and years lived with disability for 301 acute and chronic diseases and injuries in 188 countries, 1990-2013: a systematic analysis for the Global Burden of Disease Study 2013,” The Lancet, vol. 386, no. 9995, pp. 743–800, 2015.

M. Marron, P. Boffetta, Z. F. Zhang et al., “Cessation of alcohol drinking, tobacco smoking and the reversal of head and neck cancer risk,” International Journal of Epidemiology, vol. 39, no. 1, pp. 182–196, 2010.

C. R. Leemans, B. J. M. Braakhuis, and R. H. Brakenhoff, “The molecular biology of head and neck cancer,” Nature Reviews Cancer, vol. 11, no. 1, pp. 9–22, 2011.

L. Kim, T. King, and M. Agulnik, “Head and neck cancer: changing epidemiology and public health implications,” Oncology, vol. 24, no. 10, pp. 915–919, 2010, 924.

M. Hashibe, P. Brennan, S. Benhamou et al., “Alcohol drinking in never users of tobacco, cigarette smoking in never drinkers, and the risk of head and neck cancer: pooled analysis in the International Head and Neck Cancer Epidemiology consortium,” Journal of the National Cancer Institute, vol. 99, no. 10, pp. 777–789, 2007.

R. Sankaranarayanan, E. Masuyer, R. Swaminathan, J. Ferlay, and S. Whelan, “Head and neck cancer: a global perspective on epidemiology and prognosis,” Anticancer Research, vol. 18, no. 6b, pp. 4779–4786, 1998.

M. Lacko, B. J. M. Braakhuis, E. M. Sturgis et al., “Genetic susceptibility to head and neck squamous cell carcinoma,” International Journal of Radiation Oncology, Biology, Physics, vol. 89, no. 1, pp. 38–48, 2014.

T. Munshi, C. J. Heckman, and S. Darlow, “Association between tobacco waterpipe smoking and head and neck conditions: a systematic review,” Journal of the American Dental Association, vol. 146, no. 10, pp. 760–766, 2015.

S. C. Lin, S. Y. Lu, S. Y. Lee, C. Y. Lin, C. H. Chen, and K. W. Chang, “Areca (betel) nut extract activates mitogen-activated protein kinases and NF-κB in oral keratinocytes,” International Journal of Cancer, vol. 116, no. 4, pp. 526–535, 2005.

M. S. Hayden and S. Ghosh, “Shared principles in NF-κB signaling,” Cell, vol. 132, no. 3, pp. 344–362, 2008.

W. E. Naugler and M. Karin, “NF-κB and cancer — identifying targets and mechanisms,” Current Opinion in Genetics & Development, vol. 18, no. 1, pp. 19–26, 2008.

M. Koti, R. J. Gooding, P. Nuin et al., “Identification of the IGF1/PI3K/NF κB/ERK gene signalling networks associated with chemotherapy resistance and treatment response in high-grade serous epithelial ovarian cancer,” BMC Cancer, vol. 13, no. 1, p. 549, 2013.

S. Schreiber, S. Nikolaus, and J. Hampe, “Activation of nuclear factor κB inflammatory bowel disease,” Gut, vol. 42, no. 4, pp. 477–484, 1998.

G. Rogler, K. Brand, D. Vogl et al., “Nuclear factor κB is activated in macrophages and epithelial cells of inflamed intestinal mucosa,” Gastroenterology, vol. 115, no. 2, pp. 357–369, 1998.

M. F. Neurath, I. Fuss, G. Schurmann et al., “Cytokine gene transcription by NF-κB family members in patients with inflammatory bowel disease,” Annals of the New York Academy of Sciences, vol. 859, no. 1, pp. 149–159, 1998.

A. Jana, N. L. Krett, G. Guzman et al., “NFkB is essential for activin-induced colorectal cancer migration via upregulation of PI3K-MDM2 pathway,” Oncotarget, vol. 8, no. 23, pp. 37377–37393, 2017.

J. Domingo-Domenech, C. Oliva, A. Rovira et al., “Interleukin 6, a nuclear factor-κB target, predicts resistance to docetaxel in hormone-independent prostate cancer and nuclear factor-κB inhibition by PS-1145 enhances docetaxel antitumor activity,” Clinical Cancer Research, vol. 12, no. 18, pp. 5578–5586, 2006.

T. Mukhopadhyay, J. A. Roth, and S. A. Maxwell, “Altered expression of the p50 subunit of the NF-kappa B transcription factor complex in non-small cell lung carcinoma,” Oncogene, vol. 11, no. 5, pp. 999–1003, 1995.

R. C. Bargou, C. Leng, D. Krappmann et al., “High-level nuclear NF-kappa B and Oct-2 is a common feature of cultured Hodgkin/Reed-Sternberg cells,” Blood, vol. 87, no. 10, pp. 4340–4347, 1996.

V. Bours, E. Dejardin, F. Goujon-Letawe, M. P. Merville, and V. Castronovo, “The NF-κB transcription factor and cancer: high expression of NF-κB- and IκB-related proteins in tumor cell lines,” Biochemical Pharmacology, vol. 47, no. 1, pp. 145–149, 1994.

C. W. Lin, Y. S. Hsieh, C. H. Hsin et al., “Effects of NFKB1 and NFKBIA gene polymorphisms on susceptibility to environmental factors and the clinicopathologic development of oral cancer,” PLoS One, vol. 7, no. 4, article e35078, 2012.

S. S. Lo, J. H. Chen, C. W. Wu, and W. Y. Lui, “Functional polymorphism of NFKB1 promoter may correlate to the susceptibility of gastric cancer in aged patients,” Surgery, vol. 145, no. 3, pp. 280–285, 2009.

J. Gao, D. Pfeifer, L. J. He et al., “Association of NFKBIA polymorphism with colorectal cancer risk and prognosis in Swedish and Chinese populations,” Scandinavian Journal of Gastroenterology, vol. 42, no. 3, pp. 345–350, 2007.

E. T. Chang, B. M. Birmann, J. L. Kasperzyk et al., “Polymorphic variation in NFKB1 and other aspirin-related genes and risk of Hodgkin lymphoma,” Cancer Epidemiology, Biomarkers & Prevention, vol. 18, no. 3, pp. 976–986, 2009.

S. C. Lin, C. J. Liu, W. I. Yeh, M. T. Lui, K. W. Chang, and C. S. Chang, “Functional polymorphism in NFKB1 promoter is related to the risks of oral squamous cell carcinoma occurring on older male areca (betel) chewers,” Cancer Letters, vol. 243, no. 1, pp. 47–54, 2006.

F. Chen, F. Liu, L. Yan et al., “A functional haplotype of NFKB1 influence susceptibility to oral cancer: a population-based and in vitro study,” Cancer Medicine, vol. 7, no. 5, pp. 2211–2218, 2018.

B. Zhou, L. Rao, Y. Li et al., “A functional insertion/deletion polymorphism in the promoter region of NFKB1 gene increases susceptibility for nasopharyngeal carcinoma,” Cancer Letters, vol. 275, no. 1, pp. 72–76, 2009.

A. Gupta, V. Agnihotri, R. Kumar et al., “Effects of tobacco habits on the polymorphism of NFKB1 and NFKB1A gene of head and neck squamous cell carcinoma in Indian population,” Asian Pacific Journal of Cancer Prevention, vol. 18, no. 7, pp. 1855–1859, 2017.

G. F. Lehnerdt, A. Bankfalvi, S. Grehl et al., “No association of the NF-κB1 -94ins/delATTG promoter polymorphism with relapse-free and overall survival in patients with squamous cell carcinomas of the head and neck region,” International Journal of Immunopathology and Pharmacology, vol. 21, no. 4, pp. 827–832, 2008.

Y. Liu, F. Qiu, L. Yang et al., “Polymorphisms of NFκB1 and IκBα and their synergistic effect on nasopharyngeal carcinoma susceptibility,” BioMed Research International, vol. 2015, Article ID 362542, 9 pages, 2015.

Y. H. Liu, The association between the genetic variants in NFKB1 and NFΚBIA and risk of Nasopharyngeal Carcinoma, Dissertation, Guangzhou Medical University, 2012.

A. Liberati, D. G. Altman, J. Tetzlaff et al., “The PRISMA statement for reporting systematic reviews and meta-analyses of studies that evaluate healthcare interventions: explanation and elaboration,” BMJ, vol. 339, no. 1, article b2700, 2009.

A. Stang, “Critical evaluation of the Newcastle-Ottawa scale for the assessment of the quality of nonrandomized studies in meta-analyses,” European Journal of Epidemiology, vol. 25, no. 9, pp. 603–605, 2010.

A. Afshari, J. Wetterslev, and A. F. Smith, “Can systematic reviews with sparse data be trusted?” Anaesthesia, vol. 72, no. 1, pp. 12–16, 2017.

J. Wetterslev, K. Thorlund, J. Brok, and C. Gluud, “Trial sequential analysis may establish when firm evidence is reached in cumulative meta-analysis,” Journal of Clinical Epidemiology, vol. 61, no. 1, pp. 64–75, 2008.

S. Wacholder, S. Chanock, M. Garcia-Closas, L. El Ghormli, and N. Rothman, “Assessing the probability that a positive report is false: an approach for molecular epidemiology studies,” Journal of the National Cancer Institute, vol. 96, no. 6, pp. 434–442, 2004.

M. Verhagen, A. van der Meij, P. A. M. van Deurzen et al., “Meta-analysis of the BDNF Val66Met polymorphism in major depressive disorder: effects of gender and ethnicity,” Molecular Psychiatry, vol. 15, no. 3, pp. 260–271, 2010.

Y. Yu, Y. Wan, and C. Huang, “The biological functions of NF-κB1 (p50) and its potential as an anti-cancer target,” Current Cancer Drug Targets, vol. 9, no. 4, pp. 566–571, 2009.

M. Karin, “Nuclear factor-κB in cancer development and progression,” Nature, vol. 441, no. 7092, pp. 431–436, 2006.

Q. Li and I. M. Verma, “NF-κB regulation in the immune system,” Nature Reviews Immunology, vol. 2, no. 10, pp. 725–734, 2002.

M. S. Mohd Suzairi, S. C. Tan, A. A. Ahmad Aizat et al., “The functional -94 insertion/deletion ATTG polymorphism in the promoter region of NFKB1 gene increases the risk of sporadic colorectal cancer,” Cancer Epidemiology, vol. 37, no. 5, pp. 634–638, 2013.

A. S. Karban, T. Okazaki, C. I. M. Panhuysen et al., “Functional annotation of a novel NFKB1 promoter polymorphism that increases risk for ulcerative colitis,” Human Molecular Genetics, vol. 13, no. 1, pp. 35–45, 2004.

Y. M. Oltulu, E. Coskunpinar, G. Ozkan et al., “Investigation of NF-κB1 and NF-κBIA gene polymorphism in non-small cell lung cancer,” BioMed Research International, vol. 2014, Article ID 530381, 6 pages, 2014.

Z. Saidak, R. Mentaverri, and E. M. Brown, “The role of the calcium-sensing receptor in the development and progression of cancer,” Endocrine Reviews, vol. 30, no. 2, pp. 178–195, 2009.

X. Yang, P. Li, J. Tao et al., “Association between NFKB1 -94ins/del ATTG promoter polymorphism and cancer susceptibility: an updated meta-analysis,” International Journal of Genomics, vol. 2014, Article ID 612972, 8 pages, 2014.

W. Duan, E. Wang, F. Zhang, T. Wang, X. You, and B. Qiao, “Association between the NFKB1-94ins/del ATTG polymorphism and cancer risk: an updated meta-analysis,” Cancer Investigation, vol. 32, no. 7, pp. 311–320, 2014.

Y. Q. Luo, D. Wang, T. Gong, and J. Zhu, “An updated meta-analysis of 37 case-control studies on the association between NFKB1 -94ins/del ATTG promoter polymorphism and cancer susceptibility,” Oncotarget, vol. 7, no. 36, pp. 58659–58670, 2016.

D. Wang, T. Xie, J. Xu et al., “Genetic association between NFKB1 -94 ins/del ATTG promoter polymorphism and cancer risk: a meta-analysis of 42 case-control studies,” Scientific Reports, vol. 6, no. 1, article 30220, 2016.