2013 Volume 229 Issue 4 Pages 245-254
Tumor necrosis factor-alpha (TNF-α) has been proposed to contribute to the development of upper aerodigestive tract (UADT) cancer that is characterized by poor prognosis. The G-to-A nucleotide change at -308 of the TNF-α gene (-308G/A polymorphism) can increase the expression level of TNF-α and thus may affect the genetic susceptibility of UADT cancer. The association between the -308G/A polymorphism and UADT cancer has been widely studied, but the results published are quite controversial. To obtain a more precise conclusion, we performed a meta-analysis including 1,751 patients and 3,345 controls. The results indicated that the AA genotype of the -308G/A polymorphism had a 54%-increased risk of UADT cancer, compared with the G carriers (GG and GA genotypes) [odds ratio (OR) = 1.54, 95% confidence interval (CI): 1.07-2.21]. After stratified by ethnicity, the AA genotype was associated with increased risk of UADT cancers in South Asians (OR = 33.18 and 95% CI: 1.92-573.62 for AA vs. GA+GG) but not in Caucasians or East Asians. After stratified by tumor site, the -308G/A polymorphism was associated with increased risks of oropharynx cancer (OR = 2.68 and 95% CI: 1.34-5.35 for AA vs. GA+GG) but not associated with esophagus or larynx cancer. After stratified by histological type, the -308G/A polymorphism was associated with increased risks of squamous cell carcinoma (OR = 1.81 and 95% CI: 1.15-2.84 for AA vs. GA+GG) but not associated with adenocarcinoma. Our results indicate that the -308G/A polymorphism might contribute to an increased risk of UADT cancer susceptibility.
Cancers of the upper aerodigestive tract (UADT), including malignant neoplasms of the oral cavity, pharynx, larynx, and esophagus, account for approximately 5% of all new cancer cases worldwide (Ferlay et al. 2010). UADT cancers usually have poor prognosis, and they are often accompanied with severe functional defect after treatment. Many risk factors are involved in the development of UADT cancers. Among them, tobacco, areca quid and alcohol are supposed to be the most important ones (Ansary-Moghaddam et al. 2009). It has been proved that exposing to these factors would increase the risk of the disease; however, not all individuals exposed to these factors would develop UADT cancers. This indicates that genetic differences may also contribute to the pathogenesis of UADT cancers. Therefore, a large number of studies have been carried out to explore the association of gene polymorphisms with UADT cancers (Lazarus and Park 2000; Yu et al. 2012), and many susceptibility genes have been identified, such as hypoxia-inducible factor-1 alpha (Mera-Menendez et al. 2012) and microsomal epoxide hydrolase (Li et al. 2011).
Tumor necrosis factor-alpha (TNF-α) was initially identified as a serum factor which could induce necrosis of transplanted tumors in mice (Carswell et al. 1975), and was found to have antitumor effects in some clinical studies (Gutman et al. 1997). However, some studies also illustrated that TNF-α might promote the development and spreading of cancer (Ardizzoia et al. 1992). TNF-α expression is mostly regulated at transcriptional level. The TNF-α gene is located within the MHC region on chromosome 6p21.3 (Noguchi et al. 2002). Several single-nucleotide polymorphisms have been identified in this region, and the -308G/A polymorphism of the TNF-α gene is the most frequently studied one (Guo et al. 2005; Whiteman et al. 2010; Yang et al. 2011b). It has been illustrated that the G to A nucleotide change at -308 alters the transcription factor binding site, resulting in increased production of TNF-α (Abraham and Kroeger 1999; Elahi et al. 2009). Numerous studies have been carried out to investigate the association of this polymorphism with UADT cancer risk (Chiu et al. 2001; Chen et al. 2005; Guo et al. 2005; Liu et al. 2005; Gupta et al. 2008; Canova et al. 2009; Yapijakis et al. 2009; Oh et al. 2010; Whiteman et al. 2010; Yang et al. 2011a); however, the results of these studies remain controversial rather than conclusive. This might partly result from the small sample size of each study. In order to get a more precise conclusion of the relationship between the -308G/A polymorphism and UADT cancer susceptibility, we performed a comprehensive meta-analysis of all eligible case-control studies.
A comprehensive literature search in PubMed, EMBASE and China National Knowledge Infrastructure (CNKI) was carried out by two reviewers independently, to identify articles that had investigated the association of UADT cancer risk with TNF-α polymorphisms. The last updated search was performed on September 1, 2012. The search terms used were as follows: (oral or pharyngeal or pharynx or oropharyngeal or larynx or esophagus) and (cancer or carcinoma) and (Tumor Necrosis Factor or TNF or tumor necrosis factor-α or TNF-α) and (Polymorphism or mutation or variant). Manual search of the references of original articles or review articles was performed to identify additional articles on this topic. There was no language restriction. Inclusion criteria of this meta-analysis were defined as follows: (1) evaluation of the polymorphisms of -308G/A and UADT cancer risk; (2) use of a case-control design based on unrelated individuals; (3) sufficient genotype distributions for cases and controls so that an odds ratio (OR) with 95% confidence interval (CI) could be estimated; (4) genotype distribution of control population should be consistent with Hardy-Weinberg equilibrium (HWE). If more than one article were published using the same case series, studies with the largest sample size were included.
Data extractionAll the identified studies had been checked by two independent reviewers and a consensus on all items had been reached. The following items were extracted using a standardized data form: first author, publication year, ethnicity, definition of cases, source of controls, genotyping method, tumor site, tumor histology, total number of cases and controls, as well as genotype distributions in cases and controls.
Statistical AnalysisHardy-Weinberg equilibrium (HWE) in control group was firstly tested using Person’s χ 2 test. The strength of association of -308G/A polymorphism with UADT cancer risk was measured by OR and 95% CI. The statistical significance of summary OR was determined using Z-test. The pooled ORs were estimated for recessive model (AA vs. GA+GG) and dominant model (AA+GA vs. GG). The ORs of variant genotype AA compared with wild-type GG homozygote (AA vs. GG), allele A vs. allele G and GA vs. GG were also evaluated.
Heterogeneity was evaluated with a χ 2-based Q statistic, and statistical significance was assumed for P value less than 0.10. When P value was more than 0.10, OR was pooled using the fixed-effect model, otherwise the random effect model was used. The statistical significance of OR was analyzed using Z test, P < 0.05 was considered as statistically significant. To evaluate the ethnicity-specific, cancer-site specific and cancer-histology specific effects, subgroup analyses were performed according to ethnicity, tumor site and tumor histology. Sensitivity analyses were carried out by sequentially excluding individual studies to access the stability of the results. The possible publication bias was visually examined using the Begg’s funnel plot and the degree of asymmetry was tested by Egger’s test (Egger et al. 1997). All statistical tests were performed using STATA11.0.
The process of selecting studies was outlined in Fig. 1. After searching PubMed, Embase and CNKI electronic database, 276 articles were identified. After reading the titles and abstracts, 259 articles were excluded for no relevance to TNF-α gene and UADT cancer risk. After reading full texts, 6 articles were excluded for no relevance to the -308G/A polymorphism. One article was excluded for data overlapping. Three case-control studies were excluded for that the distribution of different genotypes in control group was not consistent with HWE. Finally, 12 case-control studies in 7 articles were identified (Chiu et al. 2001; Chen et al. 2005; Liu et al. 2005; Gupta et al. 2008; Oh et al. 2010; Whiteman et al. 2010; Yang et al. 2011a). In total, these studies included 1,751 cases and 3,345 controls. There were 5 studies in East Asian population, 1 study in South Asian population, 2 studies in Caucasian population, and 4 studies in mixed populations. Six of them investigated oropharynx cancer, 5 of them investigated esophagus cancer and 1 investigated larynx cancer. Nine studies investigated squamous cell carcinoma, 2 study investigated adenocarcinoma and 1 study investigated both cancer types. The characteristics of included case-controls were listed in Table 1. Genotype and allele distributions of each case-control study are listed in Table 2.
Flow diagram of included/excluded studies.
Studies | Country | Ethnicity | Cancer Site |
Cancer Histology |
Cases/ Controls |
Case gender (M/F) |
Case age | Exposure factors |
---|---|---|---|---|---|---|---|---|
Chiu CJ 2001 | China | East Asian | Oropharynx | SCC | 102/284 | 102/0 | 30-60 | Areca |
Liu CJ 2005 | China | East Asian | Oropharynx | SCC | 192/146 | 173/19 | 54.1 ± 9.8 | Unknown |
Chen WC 2005 | China | East Asian | Oropharynx | SCC | 137/102 | 125/12 | 53.3 ± 12.3 | Toabacco |
Gupta R 2008 | India | South Asian | Oropharynx | SCC | 94/133 | 76/18 | 28-72 | Toabacco |
Oh SS 2010 (E,M) | China | East Asian | Esophagus | Mixed | 194/340 | 129/65 | ≥ 20 | Tobacco Alcohol |
Oh SS 2010 (O,S) | USA | Mixed | Oropharynx | SCC | 214/849 | 183/31 | 18-65 | Tobacco Alcohol Marijuana |
Oh SS 2010 (L,S) | USA | Mixed | Larynx | SCC | 67/849 | 54/13 | 18-65 | Tobacco Alcohol Marijuana |
Oh SS 2010 (E,S) | USA | Mixed | Esophagus | SCC | 27/849 | 23/4 | 18-65 | Tobacco Alcohol Marijuana |
Oh SS 2010 (E,A) | USA | Mixed | Esophagus | ACC | 59/849 | 48/11 | 18-65 | Tobacco Alcohol Marijuana |
Whiteman DC 2010 (E,A) | Australia | Caucasian | Esophagus | ACC | 253/1,293 | 232/21 | 18-79 | Alcohol |
Whiteman DC 2010 (E,S) | Australia | Caucasian | Esophagus | SCC | 207/1,293 | 120/87 | 18-79 | Alcohol |
Yang CM 2011 | China | East Asian | Oropharynx | SCC | 205/198 | 205/0 | 49.3 ± 8.0 | Areca |
SCC, squamous cell carcinoma; ACC, adenocarcinoma; PCR-RFLP, polymerase chain reaction-restriction fragment length polymorphisms; (E,M), esophageal cancer of mixed histological types; (O,S), squamous cell carcinoma of the oropharynx; (L,S), squamous cell carcinoma of the larynx; (E,S), squamous cell carcinoma of the esophagus; (E,A), adenocarcinoma of the esophagus.
Studies | UATD cancer | Control | ||||||
---|---|---|---|---|---|---|---|---|
GG | GA | AA | GG | GA | AA | HWE | ||
Chiu CJ 2001 | 86 | 15 | 1 | 232 | 50 | 2 | 0.70 | |
Liu CJ 2005 | 175 | 16 | 1 | 120 | 24 | 2 | 0.53 | |
Chen WC 2005 | 125 | 12 | 0 | 88 | 14 | 0 | 0.46 | |
Gupta R 2008 | 61 | 23 | 10 | 114 | 19 | 0 | 0.38 | |
Oh SS 2010 (E,M) | 167 | 27 | 0 | 294 | 45 | 1 | 0.60 | |
Oh SS 2010 (O,S) | 160 | 49 | 5 | 641 | 195 | 13 | 0.68 | |
Oh SS 2010 (L,S) | 55 | 9 | 3 | 641 | 195 | 13 | 0.68 | |
Oh SS 2010 (E,S) | 19 | 8 | 0 | 641 | 195 | 13 | 0.68 | |
Oh SS 2010 (E,A) | 52 | 7 | 0 | 641 | 195 | 13 | 0.68 | |
Whiteman DC 2010 (E,A) | 157 | 84 | 12 | 842 | 403 | 48 | 0.98 | |
Whiteman DC 2010 (E,S) | 128 | 71 | 8 | 842 | 403 | 48 | 0.98 | |
Yang CM 2011 | 180 | 23 | 2 | 155 | 43 | 0 | 0.09 |
HWE, Hardy-Weinberg equilibrium; (E,M), esophageal cancer of mixed histological types; (O,S), squamous cell carcinoma of the oropharynx; (L,S), squamous cell carcinoma of the larynx; (E,S), squamous cell carcinoma of the esophagus; (E,A), adenocarcinoma of the esophagus.
As shown in Fig. 2, the heterogeneity of AA vs. GA +GG for all 12 studies was assessed and the P value was 0.47, suggesting the absence of heterogeneity. Thus, fixed-effects model was used to obtain the data. Overall OR was 1.54 (95% CI = 1.07-2.21), and the test for overall effect Z value was 2.34 (P = 0.02) for AA vs. GA +GG model. The results suggested that the AA genotype of the -308G/A polymorphism had a 54%-increased risk of UADT cancer, compared with the G carriers (GG and GA genotypes). Summary of the results of other genetic comparisons are listed in Table 3.
Forest plot of UADT cancer risk associated with the -308G/A polymorphism of the TNF-α gene in overall analysis (AA vs. GA+GG).
The squares and horizontal lines correspond to the study-specific OR and 95% CI. The area of the squares reflects the study-specific weight. The diamond represents the summary OR and 95% CI.
AA VS. GA +GG (recessive model) |
AA+GA VS. GG (dominant model) |
AA VS. GG (homozygote model) |
AA VS. GA (heterozygote model) |
A VS. G (allelic contrast) |
|||||||
---|---|---|---|---|---|---|---|---|---|---|---|
Na | OR(95%CI) | Pb | OR(95%CI) | Pb | OR(95%CI) | Pb | OR(95%CI) | Pb | OR(95%CI) | Pb | |
Total | 12 | 1.54 (1.07,2.21) | 0.48 | 0.92 (0.70,1.17) | < 0.01 | 0.71 (0.48,1.06) | 0.47 | 1.47 (0.99,2.15) | 0.49 | 0.94 (0.73,1.22) | < 0.01 |
Ethnicities | |||||||||||
East Asian | 5 | 1.10 (0.33,3.73) | 0.60 | 0.68 (0.52,0.88) | 0.21 | 1.03 (0.30,3.50) | 0.61 | 1.66 (0.49,5.63) | 0.57 | 0.70 (0.53,0.93) | 0.27 |
South Asian | 1 | 33.18 (1.92,573.62) | - | 3.25 (1.70,6.18) | - | 39.10 (2.25,678.56) | - | 17.43 (0.96,316.65) | - | 3.86 (2.16,6.87) | - |
Caucasian | 2 | 1.18 (0.72,1.92) | 0.68 | 1.15 (0.93,1.41) | 0.97 | 1.23 (0.75,2.02) | 0.70 | 1.08 (0.65,1.81) | 0.65 | 1.13 (0.95,1.34) | 0.91 |
Mixed | 4 | 1.60 (0.76,3.35) | 0.66 | 0.86 (0.66,1.12) | 0.12 | 1.53 (0.73,3.21) | 0.69 | 1.84 (0.85,4.01) | 0.48 | 0.88 (0.61,1.27) | 0.15 |
Cancer sites | |||||||||||
Oropharynx | 6 | 2.68 (1.34,5.35) | 0.13 | 0.86 (0.51,1.46) | < 0.01 | 2.70 (1.35,5.36) | 0.10 | 2.59 (1.23,5.46) | 0.37 | 0.93 (0.53,1.62) | < 0.01 |
Esophagus | 5 | 1.12 (0.70,1.79) | 0.96 | 1.05 (0.82,1.34) | 0.20 | 1.15 (0.72,1.86) | 0.94 | 1.05 (0.64,1.72) | 0.98 | 1.03 (0.83,1.29) | 0.20 |
Larynx | 1 | 3.01 (0.84,10.85) | - | 0.67 (0.35,1.28) | - | 2.69 (0.74,9.72) | - | 5.00 (1.21,20.73) | - | 0.84 (0.48,1.47) | - |
Histology type | |||||||||||
SCC | 9 | 1.81 (1.15,2.84) | 0.26 | 0.91 (0.64,1.28) | < 0.01 | 1.85 (1.18,2.92) | 0.26 | 1.70 (1.05,2.75) | 0.28 | 0.97 (0.69,1.36) | < 0.01 |
ACC | 2 | 1.21 (0.65,2.27) | 0.54 | 0.74 (0.28,1.99) | 0.02 | 1.23 (0.65,2.32) | 0.46 | 1.19 (0.62,2.29) | 0.89 | 0.74 (0.28,1.95) | 0.02 |
Mixed | 1 | 0.58 (0.02.14.35) | - | 1.03 (0.62,1.72) | - | 0.59 (0.02,14.47) | - | 0.55 (0.02,14.02) | - | 1.01 (0.62,1.65) | - |
Areca chewers | 2 | 2.53 (0.43,14.69) | < 0.01 | 0.63 (0.42,0.94) | 0.32 | 2.33 (0.40,13.71) | < 0.01 | 3.76 (0.66,21.49) | < 0.01 | 0.69 (0.47,1.01) | 0.07 |
a Number of comparisons
b P-value of Q-test for heterogeneity test
Chen WC 2005 was only included in the analysis of (AA+GA VS. GG) model and (A VS. G) model for lacking of TNF-α-308 AA genotype
Random-effects model was used when P-value for heterogeneity test < 0.05, otherwise, fixed-effects model was used.
Subgroup analyses were performed according to ethnicity, tumor site and tumor histology. For ethnicity (Fig. 3), all the studies were divided into 4 subgroups: East Asian, South Asian, Caucasian and mixed (mixed with different populations). In recessive model (AA vs. GA +GG), significantly increased risks were found among South Asian (OR = 33.18, 95%CI = 1.92-573.62, P = 0.02) but not in East Asian (OR = 1.10, 95%CI = 0.33-3.73), Caucasian (OR = 1.18, 95%CI = 0.72-1.92, P = 0.52) or Mixed group (OR = 1.60, 95%CI = 0.76-3.35). In the subgroup analysis of tumor site (Fig. 4), significantly increased risks were found in oropharynx cancer (OR = 2.68, 95%CI = 1.34-5.35, P = 0.01); no association was found with esophagus cancer (OR = 1.12, 95%CI = 0.70-1.79, P = 0.64), or larynx cancer (OR = 3.01, 95%CI = 0.84-10.85, P = 0.09). In the subgroup of cancer histology type (Fig. 5), significantly increased risks were found in squamous cell carcinoma (OR = 1.81, 95%CI = 1.15-2.84, P = 0.01); no association was found in other cancer histology type: adenocarcinoma (OR = 1.21, 95%CI = 0.65-2.27, P = 0.55) or mixed (OR = 0.58, 95%CI = 0.02-14.35, P = 0.74). Other genetic comparisons are summarized in Table 3.
Forest plot of UADT cancer risk associated with the -308G/A polymorphism (AA vs. GA+GG), stratified by ethnicity.
The squares and horizontal lines correspond to the study-specific OR and 95% CI. The area of the squares reflects the study-specific weight. The diamond represents the summary OR and 95% CI.
Forest plot of UADT cancer risk associated with the -308G/A polymorphism (AA vs. GA+GG), stratified by tumor site. The squares and horizontal lines correspond to the study-specific OR and 95% CI. The area of the squares reflects the study-specific weight. The diamond represents the summary OR and 95% CI.
Forest plot of UADT cancer risk associated with the -308G/A polymorphism (AA vs. GA+GG), stratified by histology type.
The squares and horizontal lines correspond to the study-specific OR and 95% CI. The area of the squares reflects the study-specific weight. The diamond represents the summary OR and 95% CI.
Cases and controls in two studies which investigated oropharynx cancer were all areca chewers. Compared with the GG homozygote group, the AA/AG genotypes were associated with a statistically significant decreased risk of cancer risk (OR = 0.63, 95%CI = 0.42-0.94) (Table 3).
Publication biasPublication bias of literature was assessed using Begg’s funnel plot and Egger’s test. The shapes of the funnel plots did not show any evidence of obvious asymmetry (Fig. 6). Egger’s test was also performed to provide statistical evidence of funnel plot asymmetry. The results did not suggest any evidence of publication bias (P = 0.53).
Begg’s funnel plot to detect publication bias test (AA vs. GA+GG).
Each point represents a separate study for the indicated association. Log (OR): natural logarithm of OR. Horizontal line represents size of effect.
Sensitivity analysis was performed by sequentially excluding each study to assess the stability of the current analysis. All the obtained results were similar and did not draw different conclusions.
UADT cancer contributes to approximately 5% of all new cancer cases worldwide (Ferlay et al. 2010). There is evidence that inflammation and cytokines play a crucial role in the development of UADT cancer (Marx 2004). TNF-α, a critical cytokine, has been demonstrated to play an important role in several solid tumors (Balkwill 2006). Previous studies showed that the -308G/A polymorphism was associated with UADT cancer risk; however, the results of these studies were controversial rather than conclusive (Chiu et al. 2001; Chen et al. 2005; Liu et al. 2005; Gupta et al. 2008; Oh et al. 2010; Whiteman et al. 2010; Yang et al. 2011a). In order to confirm the association between the -308G/A polymorphism and the risk of UADT cancer, we performed this meta-analysis.
In the overall meta-analysis, we found that in the recessive model, overall OR for AA versus GA + GG was 1.54 (95%CI = 1.07-2.21). TNF-α is a pleiotropic cytokine during the process of carcinogenesis. High doses of local TNF-α could suppress tumor via destructing tumor blood vessels and promoting lysis of tumor cells (Lejeune 2002). However, low doses of TNF-α in tumor microenvironment could promote tumor cell growth and tumor vascularization (Balkwill 2009). The reason that A allele carriers have higher UADT cancer risk might be that A allele could improve the transcription efficiency of TNF-α, and excessive expression of TNF-α is associated with cancer progression as mentioned above.
After ethnicity subgroup analysis, we found a significant association between the -308G/A polymorphism and UADT cancer risk in South Asian, but not in East Asian, Caucasian or mixed ethnic populations. There might be several reasons for this observation. First, different ethnic populations usually have different linkage disequilibrium (LD) patterns (Frisse et al. 2001). A polymorphism locus may link with another risk variant closely in certain ethnic population. Second, environmental context which strongly varies between populations may play a role. However, as only one case-control study in South Asian population was included, the results should be taken with caution. Additional studies are needed to validate the effect of -308G/A polymorphism on UADT cancer risk in different ethnic groups.
After subgroup analysis according to tumor site, results showed that the -308G/A polymorphism was associated with an increased risk of oropharyngeal cancer; but the association was not found in esophageal and laryngeal cancer subgroups. One possible explanation is that cancers at different sites might be exposed to different environmental factors, and that the -308G/A polymorphism may exert varying effects in different cancers (Pelucchi et al. 2008). It should be noticed that the sample size of larynx cancer was limited which might have insufficient statistical power to detect a slight effect (Zintzaras and Lau 2008).
When stratified by tumor histology, we found that the -308G/A polymorphism was associated with an increased risk of squamous cell carcinoma but the association was not found in the adenocarcinoma or the mixed group. One explanation for this discrepancy is that cancers of different histology type might have different expression profiles of TNF-α and thus the same polymorphism might play different roles in cancers of different origins.
Apparent gene-environment interaction was observed between the -308G/A polymorphism and areca chewing. Under dominant model, AA/AG carriers were associated statistically with significantly decreased risk of cancer compared with GG homozygotes. Previous studies have shown that the A allele could increase the production of TNF-α (Kroeger et al. 1997). TNF-α may induce the apoptosis process of tumor cells through binding to TNF receptors, which can lead to the activation of endonucleases and DNA fragmentation (van Horssen et al. 2006). TNF-α may also promote lysis of tumor cells by activating the anti-tumor T cell response (Xu et al. 2007). Furthermore, a previous study has shown that GG genotype could increase the risk of oral submucosa fibrosis (OSF), which is a pre-cancerous lesion with a high risk of malignant transformation (Chiu et al. 2001). Interestingly, the overall meta-analysis showed that AA genotype was associated with increased UADT cancer risk, which was contrary to the results of areca chewing subgroup. Several reasons may lead to the contradiction. First, we should notice that all of the cases and controls were East Asian in the areca chewing subgroup. Different ethnicities have different genetic background, different linkage disequilibrium patterns and different life style. These factors may interact with the -308G/A polymorphism and affect the results. Second, all participants in the areca chewing subgroup were male. Compared with the females, males have vital different hormone level which may interact with TNF-α and affect the pathogenesis of cancer (Calle and Kaaks 2004). Third, we should also notice that there were only two case-control studies in the areca chewing subgroup, which might have limited power to establish a reliable association. More case-control studies are required to validate the association.
Heterogeneity is very important in meta-analysis. There was no heterogeneity in overall comparisons in the recessive model. Even after stratified by ethnicity, tumor site and histological type of tumor, the heterogeneity was not significant in all the subgroups. Publication bias is another vital issue in meta-analysis. In this study, publication bias was assessed using Begg’s funnel plots and Egger’s test. No significant publication bias was found, suggesting the reliability of our results.
Some limitations of the present meta-analysis should be noticed. First, we only retrieved articles form PubMed, EMBASE and CNKI databases. It is possible that some related articles might be in some other databases and thus were not included. Secondly, although we did not find significant publication bias in this meta-analysis, it is possible that there were some unpublished results. Thirdly, data were not stratified by other factors, such as gender or life style as we could not extract sufficient information from the found documents. Moreover, in the stratified analysis, the participant number of some subgroups was small, which could lead to low statistic power to reveal a reliable association.
In conclusion, this meta-analysis provided evidence that the -308G/A polymorphism might contribute to an increased risk of UADT cancer susceptibility. Future case-control studies with larger sample size should be conducted to validate these findings.
This work was supported by grants from the National Natural Science Foundation of China (No. 30930100, 81072218, 30872872, 81001208, 81200791, 81061120531), International Science and Technology Cooperation Program of China (2012DFA31370), the Doctoral Program of the Ministry of Education of China (No. 20110181110055, 20120181120011) and the Specialized Research Fund for the Doctoral Program of Higher Education of China (20100181120057).
The authors have no conflict of interest to declare.