Association of Common Polymorphisms in TNFA, NFkB1 and NFKBIA with Risk and Prognosis of Esophageal Squamous Cell Carcinoma

Background Tumour necrosis factor-alpha (TNF-α) and nuclear factor of kappa light chain gene enhancer in activated B cells (NF-κB) play critical role in carcinogenesis processes like tumour initiation, proliferation, migration and invasion. Single nucleotide polymorphisms in TNF-α, NF-κB and its inhibitor IκB genes were shown to be associated with susceptibility and prognosis of several cancers; however, their role in esophageal squamous cell carcinoma (ESCC) is not well recognised. Therefore, in present study, we aimed to investigate association of common polymorphisms in TNFA, NFkB1 and NFKBIA with risk and prognosis of ESCC in northern Indian population. Methods We genotyped 290 ESCC patients (including 162 followed up cases) and 311 mean age, gender and ethnicity matched controls for TNFA -308G>A, NFkB1 -94ATTG ins/del and NFKBIA (-826C>T and 3’UTRA>G) polymorphisms using PCR alone or followed by RFLP and TaqMan assay. Results TNFA -308GA genotype was associated with increased risk of ESCC specifically in females and in patients with regional lymph node involvement, while, NFKBIA -826CT+TT genotype conferred decreased risk of ESCC in females. Haplotypes of NFKBIA -826C>T and 3’UTRA>G polymorphisms, C-826G3’UTR and T-826A3’UTR, were associated with reduced risk of ESCC. No independent role of NFkB1 -94ATTG ins/del polymorphism in susceptibility of ESCC was found. Multi-dimensionality reduction analysis showed three factor model TNFA-308, NFKBIA-826, NFKBIA 3’UTR as better predictor for risk of ESCC. Furthermore, combined risk genotype analysis of all studied polymorphisms showed increased risk of ESCC in patients with 1-3 risk genotype compared to ‘0’ risk genotype. Survival analysis did not show any significant prognostic effect of studied polymorphisms. However, in stepwise multivariate analysis, metastasis was found to be independent prognostic predictor of ESCC patients. Conclusion TNFA-308 and NFKBIA (-826C>T and 3’UTRA>G) polymorphisms may play role in susceptibility but not in prognosis of ESCC patients in northern Indian population.


Introduction
Chronic inflammation, a critical component of tumour microenvironment, is involved in pathogenesis of approximately 25% of all human cancers including esophageal cancer (EC) [1,2]. Tumour necrosis factor-alpha (TNF-α) and nuclear factor of kappa light chain gene enhancer in activated B cells (NF-κB) are two major mediators of inflammation in cancer and they are intricately linked to malignant processes like tumour initiation, proliferation, invasion and angiogenesis [3,4].
TNF-α gene (TNFA) spans about 3 kilo-base pair on chromosome 6p21.3 and contains 4 exons. It encodes a proinflammatory pleiotropic cytokine that plays role in cell differentiation, proliferation, apoptosis and immunity [5]. Deregulation of TNF-α is implicated in wide spectrum of diseases like diabetes, osteoporosis, autoimmune diseases and cancers [6]. Also, expression of TNF-α was found to progressively increase from precancerous to cancerous lesions of esophageal carcinoma pointing its critical role in EC [7]. Several functional polymorphisms are present in the promoter region of TNFA; however the most documented single nucleotide polymorphism (SNP) is located at 308 nucleotide position resulting in substitution of G to A (rs1800629). Epidemiological studies exploring association of TNFA-308 G>A polymorphism with cancer risk are inconsistent. While, some studies have shown increased risk of carcinoma of breast [8], colon [9], oral cavity [10] and cervix [11], other studies did not find any role of the polymorphism in cervical [12] and prostate cancers [13]. In esophageal cancer, there are two reports which failed to find any association of the polymorphism [14,15]. TNFA-308 polymorphism may also have prognostic implication as it was found to confer adverse outcome to head & neck cancer and gastro-esophageal patients [15,16].
NF-κB is a family of transcription factors that are activated by TNF-α in classical canonical pathway [17]. There are five members of NF-κB: RelA, RelB, c-Rel, p50/105 (encoded by NFKB1 gene; chromosomal location: 4q23-q24), and p52/p100. The dimeric form of NF-κB, p50/RelA, is the most common form [18]. In un-stimulated cell, NF-κB remains sequestered in cytoplasm by its inhibitor IκB. Following activating stimuli, IκBs are phosphorylated and degraded, so NF-κB is activated and is translocated to the nucleus to initiate the target gene expression [19]. The IκB family also constitutes several members among which IκBα (encoded by NFKBIA, located on 14q13) is classical form that can be found in the cytoplasm and nuclei [20]. NF-κB was reported to be constitutively activated in esophageal squamous cell carcinoma (ESCC) tissues and ESCC cell lines (Eca109 and EC9706) [21,22]. While, a study has also shown amplification and overexpression of NFKBIA gene in KYSE series EC cell lines [23]. Several polymorphisms are present in NFKB1 (1900 SNPs) and NFKBIA (158 SNPs) according to dbSNP database (www.ncbi.nlm.nih.gov/snp), however, previous studies have extensively explored role of common polymorphic variants in promoter region of NFKB1 (-94 ATTG ins/del; rs28720239) and NFKBIA (-826 C>T; rs2233406) and 3'UTR region of NFKBIA (3'UTR A>G; rs696) in susceptibility and prognosis of various cancers [24][25][26][27][28][29]. Literature exploring association of these NFKB1 and NFKBIA variants in EC are missing till now. So, in the present study, we investigated the association of TNFA-308 G>A, NFKB1 -94ATTG ins/del and NFKBIA (-826 C>T and 3'UTR A>G) polymorphisms with susceptibility to esophageal squamous cell carcinoma (ESCC) or its clinical phenotypes, their interaction with environmental risk factors and their role in survival outcome of ESCC patients.

Ethics statement
The study was approved by the ethical committee of Sanjay Gandhi Postgraduate Institute of Medical Sciences, Lucknow (28 th Sanjay Gandhi Post Graduate Institute of Medical Sciences Ethics Committee Meeting on May 14, 2003) and written informed consent was obtained from all recruited individuals. All the clinical and demographical details were recorded anonymously by one individual (data collector), while experiments and analysis was performed by another individual (investigator). The study subjects were coded in random numbers and identity of patients (name and their hospital identity number) was never revealed to investigator at any stage of the study. Blood samples and other details from all study subjects were collected after taking their informed written consent.

Subject recruitment
During the period of 7 years from 2005 to 2012, 311 EC cases were recruited from out-patients clinic of Gastroenterology and Radiotherapy departments, Sanjay Gandhi Post Graduate Institute of Medical Sciences, Lucknow. The diagnosis of patient was confirmed by pathologist. EC patients with two major histopathologies: squamous cell carcinoma (SCC) and adenocarcinoma (ADC) were recruited in the study. However, the frequency of ADC cases (6.8%) was low in our study and etiologies of SCC and ADC are entirely different; therefore, we carried our analysis in 290 incident ESCC cases only. Patients with pre-malignant conditions like corrosive esophageal injury, Achalasia injury, Barrett's esophagus or Gastroesophageal reflux diseass and prior treatment before enrolment were excluded from the study. Data on demography, clinical parameters (like tumour location and regional lymph node involvements) and environmental exposures (alcohol consumption, smoking of cigarette or 'bidi' and smokeless tobacco usage) were recorded only for the patients in a standard proforma from medical records and personal interview as described previously [30]. During the same time frame, 311 mean age, gender and ethnicity matched cancer free controls were recruited in the study. The exclusion criteria of controls were presence of any malignancy, chronic disease or pre-malignant conditions. For survival analysis, patients who had received radiotherapy (RT) with or without concurrent cisplatin based chemotherapy (CT) were followed up during 2005-2012. Many surgically resected cases and patients who discontinued their treatment could not be followed up, so prognostic evaluation was limited to 162 cases. The inclusion criteria of patients for survival analysis was patients who had received RT based treatment regime i.e. patients who had received RT or RT along with CT. However, patients who were referred for surgery or for palliative treatment or with previous history of surgical resection were excluded in the survival analysis. The maximum follow up time was 72.6 months and the median follow up time was 9.3 months. Dosage of RT and CT, dysphagia grade and duration, tumour stage, Karnofsky performance score and other parameters were recorded for each followed patients. The patients were prospectively followed up every six months from the time of enrolment through telephonic conversations or their visit to out-patient clinic until death or end of the study. The regional lymph node involvement was recorded as 'present' or 'absent'. Around, 47% (120/255) of cases had regional lymph node involvement in overall-group of patients, while, in followed up patients group, regional lymph node positivity was recorded in only 39.5% (64/162) of cases.

DNA extraction and SNP genotyping
A 3 ml venous blood was collected from each subjects in sterile EDTA vacutainer tubes and stored at -40°C until DNA extraction. Genomic DNA was extracted from peripheral blood leukocytes using standard salting out method [31]. The genotyping of TNFA-308 G>A and NFKB1 -94ins/del ATTG polymorphisms were carried out through PCR based methods using primer sequences as described previously [32,33]. NFKBIA -826 C>T and 3'UTR A>G polymorphisms were genotyped by TaqMan assay and PCR RFLP [34] respectively. The details of genotyping methods are described in Table S1.
To improve genotyping quality and validation, 20% of samples were re-genotyped by other laboratory personnel and genotyping results were reproducible with no discrepancy. Lab personals were also blinded to case control status during genotyping to eliminate bias.

Statistical analysis
The effective sample size for case-control study was calculated using Quanto 1.1 ver. software [35], while minimum sample size for case-only survival analysis for obtaining 80% power was estimated using PS version 2.1.3.1 software [36]. Genotypic frequencies of each studied polymorphisms among control subjects were checked for Hardy-Weinberg equilibrium (HWE) using goodness-of-fit χ 2 test. Binary logistic regression was applied to calculate Odds ratio (OR) and 95% Class interval (CI) for various predictors after adjustment of covariates like age and gender. Haplotype analysis of NFKBIA polymorphisms was conducted using SNPAnalyzer version 1.0 [37]. Since response rate was low, case only analysis was performed for gene-environment interaction. In case of multiple comparisons, False discovery rate (FDR) test was applied to avoid type 1 error and the threshold value was taken as 0.10. Multi-factor Dimensionality Reduction (MDR) analysis was performed to evaluate the high order interaction between the polymorphisms using MDR 3.0.2 software (www.multifactordimensionalityreduction.org). MDR software gives number of output parameters like cross validation consistency (CVC), testing accuracy (TA), balanced training accuracy for different interactions and single best model is identified as interaction that had maximum CVC and TA. Statistical significance of the model was evaluated using a 1000-fold permutation test. Kaplan Meier and Log rank tests were carried out to estimate the difference in survival times according to genotypes and clinical/demographical characteristics. Survival time was calculated from date of ESCC diagnosis to death of patients or date of last follow up. Univariate Cox regression analysis was done to determine predictive factor of ESCC survival by estimating Hazard ratio (HRs) and 95% CI. Multivariate analysis was also performed, in which all variables were first entered together in single step and after that in stepwise manner also, to identify independent prognostic predictor of ESCC. Two models (forward selection and backward elimination) were employed in stepwise Cox regression analysis. All statistical analyses were performed with SPSS software version 15.0 (SPSS, Chicago, Illinos, USA) and differences were taken as significant when two sided Pvalue was less than 0.05.

Results
The power calculation analysis showed that at minimum minor allelic frequency (MAF) of 5.4% (as reported for Gujarati Indian population in Hapmap database for TNFA-308 G>A polymorphism) and genetic effect of 2.0, the case control pair of 240 was sufficient to achieve 80% power.

Characteristics of subjects included in susceptibility analysis
The demographic and clinical characteristics of cases and control are shown in Table 1. There was no significant difference in median age and gender distribution between cases (Median age = 57 years and Males: 72.8%) and controls (Median age: 55 years and Males: 71.1%). Majority of patients (59.3%) had tumour at middle third location and 47% of cases had regional lymph node involvement. Data on environmental risk factors showed that 79% of cases used tobacco in some form (chewing, smoking or snuff) and 28.5% had alcohol drinking habit.

Association of the TNFA-308 G>A, NFKB1 -94ATTG ins/del and NFKBIA (-826 C>T and 3' UTR A>G) polymorphisms with overall risk of ESCC
Chi square test showed genotypic frequencies of polymorphisms in TNFA, NFKB1 and NFKBIA were in accordance with HWE in controls (P>0.05 in each case). When genotypic distribution of TNFA -308 G>A polymorphism was compared between cases and controls, 1.7 fold increased risk of ESCC was observed with TNFA -308 GA genotype compared to GG genotype (OR = 1.73, 95% CI = 1.13-2.67, P = 0.013). Also, TNFA -308 A allele was associated with enhanced risk of ESCC compared to G allele (OR = 1.62, 95% CI = 1.08-2.42, P = 0.019) ( Table 2). No associations of NFKB1 -94 ins/del, and NFKBIA (-826 C>T and 3'UTR A>G) polymorphisms were observed with the overall risk of ESCC.

Survival analysis
Univariate and multivariate Cox regression analysis showed no significant effect of TNFA-308 G>A, NFKB1 -94ATTG ins/ del, NFKBIA (-826 C>T and 3'UTR A>G) polymorphisms on the survival outcome of ESCC patients (Table 3 and 4). However, when all of the variables in the study (genetic polymorphisms and clinico-pathological variables) were entered in single step in multivariate analysis, presence of metastasis conferred poor survival outcome to ESCC patients (HR = 3.39, 95% CI = 1.44-7.98, P = 0.005). Similarly, stepwise regression analysis (in both forward selection and backward elimination model) indicated metastasis as an independent predictor of prognosis of ESCC patients (HR = 2.72, 95% CI =1.44-5.14, P = 0.002) ( Table 4).

Interaction of the polymorphisms with environmental risk factors
A case only analysis was performed to evaluate interaction of the polymorphisms with environmental risk factors, however, no significant modulation in risk of ESCC was found in tobacco users, smokers and alcohol users (Data not shown).

Combined risk genotype analysis of TNFA-308 G>A, NFKB1 -94ATTG ins/del and NFKBIA (-826 C>T and 3'UTR A>G) polymorphisms
We also performed combined risk genotype analysis to study the effect of risk genotypes of all four polymorphisms on susceptibility to ESCC. For this, we pooled risk genotype (OR >1) of all four SNPs into new variable according to the number of risk genotypes, ranging from 0-4 (in case of protective association, we reversed the reference group). A significant dose dependent risk of ESCC was observed with combined risk genotype of all four SNPs (P trend = 0.001). Moreover, in dichotomized analysis, patients with 1-3 risk genotype had 1.67 fold higher risk of ESCC compared to patients with '0' risk genotype (OR =1.67, 95% CI = 1.21-2.31, P =0.002, q value = 0.005) ( Table 8). Role of combined risk genotype in the prognosis of ESCC patients was also analyzed. However, no significant difference in median survival was found among patients carrying '0', '1' '2' or '3'risk genotypes (Log rank P value = 0.882). Cox regression analysis also did not show any significant hazard of death with any of the risk genotypes i.e., '1' '2' or '3'risk genotypes (Data not shown). Furthermore, there was no effect of dichotomized risk genotype (1-3 vs. 0 risk genotype) on survival outcome of ESCC patients (Table 8).

Multi-dimensionality reduction (MDR) analysis
MDR analysis showed that TNFA -308 polymorphism was the best one factor model and TNFA -308 and NFKBIA 3'UTR were the best model for two factors and TNFA-308, NFKBIA -826, NFKBIA 3'UTR polymorphisms as best model for three factors. However, P value of permutation testing was significant only for three factor model (P value = 0.006, 0.264 and 0.214 for three, two and one factor model respectively) (Table S2). Thus, three factor model TNFA-308, NFKBIA-826, NFKBIA 3'UTR with testing accuracy of 0.57 and 10/10 CVC is better predictor risk of ESCC compared to one/two factors.

Discussion
Inter-individual variation in ESCC may be partly attributed to genetic variants in inflammatory and immune-responsive genes [38,39]. In the present study, we investigated the role of common polymorphic variants in TNFA, NFKB1 and NFKBIA (the major mediators of inflammatory and immune response in malignancy) with risk and prognosis of ESCC in northern Indian population. We observed that TNFA-308 G>A polymorphism was associated with enhanced risk of ESCC especially in females and in patients with regional lymph node involvements. No independent role of either NFKB1 -94ATTG ins/del or NFKBIA (-826 C>T and 3'UTR A>G) polymorphisms in ESCC susceptibility was found, however, two haplotypes of NFKBIA polymorphisms (C -826 G 3 ' UTR and T -826 A 3 ' UTR) seem to have protective role in ESCC. MDR analysis showed TNFA -308, NFKBIA -826 and NFKBIA 3'UTR polymorphisms as better predictor for risk of ESCC. Moreover, when we pooled risk genotype of all the polymorphisms, a significant increased risk of ESCC was observed in subjects with ≥1 risk genotypes compared to subjects with '0' risk genotype.
The MAF of TNFA -308 G>A polymorphism in present study was 7.1% which is similar to Hapmap Guajarati Indian (GIH) and Han Chinese (CHB) population (http// hapmap.ncbi.nlm.nih.gov/). We found a significant 1.7 fold increased risk of ESCC with TNFA -308 GA genotype and A allele in comparison to GG genotype to G allele respectively. Literature suggests strong association of TNFA -308 A allele with MHC haplotype HLA-A1-B8-DR3, which is in turn linked with increased production of TNF-α [40,41]. Also, TNFA -308A allele is stronger transcriptional activator than G allele and production of TNF-α was reported to be higher in monocyte culture of individuals with TNFA-308 GA compared to that with GG genotype [42,43]. Thus, increased risk of ESCC with TNFA -308 GA genotype and A allele may be due to higher TNF-α production, which promotes malignant processes through induction of premalignant chemokines, angiogenic mediators, reactive oxygen intermediates and inflammatory mediators [4]. In contrast to our finding, two previous studies did not find independent association of TNFA -308 G>A polymorphism with susceptibility to gastro-esophageal cancer or ESCC [14,15]. However, their results should be interpreted cautiously as TNFA-308 polymorphism was not in HWE in controls in these studies (P HWE were less than 0.05). For this reason, we could not perform meta-analysis of TNFA -308 polymorphism in esophageal cancer as one of exclusion criteria for studies in Table 3. Univariate survival analysis of selected gene polymorphisms in ESCC. meta-analysis is disagreement of polymorphism with HWE in controls. NF-κB is a major transcription regulator of immune response, apoptosis and cell-cycle genes [18]. A four base pair deletion polymorphism in NFKB1 i.e., -94ATTG del is extensively explored in various cancers. The deletion allele of the polymorphism abolishes binding site for nuclear protein resulting in reduced promoter activity [44]. Although, prior studies have shown significant association of the NFKB1 polymorphism with carcinoma of urinary bladder, prostate, cervix and nasopharynx [33,[45][46][47], no role of the polymorphism in ESCC susceptibility was found in the present study. Similar to our study, Riemann et al. did not find association of the polymorphism with bladder cancer, colorectal cancer, renal cell carcinoma and B cell chronic lymphocytic leukemia in German subjects [48,49]. The different findings of NFKB1 -94ATTG del polymorphism may be due to its cancer specific or population specific nature of association. We also examined association of two functional polymorphisms in NFKBIA (-826C>T and 3'UTRA>G) with ESCC susceptibility, however no significant role in ESCC was observed.

Genotypes N (%) Median survival (in months) Log Rank P value HR (95% CI) P value
Univariate survival analysis through Kaplan Meier test and Cox regression showed no difference in median survival or death hazard of ESCC patients with any of the selected gene polymorphism. Also, in multivariate analysis, no effect of studied gene polymorphisms in prognosis of ESCC was found. However, presence of metastasis was associated with worse survival outcome of ESCC patients in multivariate analysis, which is well known fact. A previous study showed that TNFA-308 AA genotype was associated with poor prognosis of tumour at gastro-esophageal (GE) junction study [15]. This could be due to entirely different etiology of tumour located at GE junction compared to esophageal squamous cell carcinoma. The lack of prognostic role of the selected gene polymorphisms in ESCC patients in the present study may be due to small number of alive cases (15.4%) compared to dead or lost to follow up cases (84.6%) or short median follow-up. Stratification based on gender showed that TNFA-308 GA was associated with increased risk of ESCC in females, in contrast, NFKBIA -826 CT genotype conferred female-specific decreased risk of ESCC. These is sexual dimorphism in immune/inflammatory response and females exhibits more vigorous cellular/ humoral immune reactions compared to males [50]. Also, female hormone progesterone was found to increases TNF-α secretion in activated monocytes [51]. So, female specific association of TNFA-308 A allele may explained by higher TNF-α level in females leading to aggressive inflammatory response. However, these findings (female specific associations) need to be confirmed in larger cohort and should be interpreted cautiously due to low number of female patients in the present study.
ESCC shows considerable heterogeneity in clinical phenotypes and different genetic factors/pathways are implicated in development of specific esophageal tumour phenotypes [52,53]. In line with this, NFKBIA -826 CT+TT and 3'UTR AG+GG genotypes seem to be associated with lower risk of upper and lower third esophageal tumour respectively. Furthermore, TNFA-308 GA genotype conferred significant increased risk of regional lymph node involvement in ESCC patients. Similar to this, Kim et al. had shown association of the TNFA -308 polymorphism with lymph node metastasis in gastric cancer patients [54].
Regional lymph node involvement usually associated with poor clinical outcome of the cancer patients. However, in the present study no modification in prognosis of esophageal cancer patients was observed with lymph node status. This may be due to low number (64/162, 39.5%) of followed up patients had regional lymph node involvement. So, despite association TNFA-308 GA with regional lymph node involvement in ESCC patients, no role of TNFA-308 polymorphism in clinical outcome of ESCC patients was found.
We next carried out the haplotype analysis of two NFKBIA polymorphisms which showed significant decreased risk of ESCC with C -826 G 3 ' UTR and T-826 A 3 ' UTR haplotypes. Bioinformatics analysis showed that NFKBIA -826 C>T polymorphism may affect binding of transcription factor like BRAC1/2, TBP and MYB, while 3'UTRA>G polymorphism may affect binding of mir-196a (www.SNPinfo.com). Few studies have also examined the effect of NFKBIA 3'UTRA>G polymorphism on expression of the protein, however, finding are inconsistent. While, the expression of NFKBIA was reported to be higher in peri-tumour tissues from colorectal cancer patients with 3'UTR AA+AG genotypes than those with 3'UTR GG genotype [55], another study did not find difference in expression of NFKBIA in melanoma patients with different 3'UTR genotypes [56]. Data demonstrating exact functional role of NFKBIA -826 and 3'UTR polymorphisms in ESCC are lacking, however, it is reasonable to assume that differential expression of IκB due to NFKBIA specific haplotypes may result in different NF-κB activation, further leading to differential risk for ESCC.
Cancer is a multi-genic disease in which single SNP may only have a modest independent effect on disease phenotype and multiple SNPs may provide a more accurate representation of the risk. So we carried a high order Table 5. Gender specific association of selected polymorphism with risk of ESCC.

Genotypes
Males Females  Table 6. Association of selected polymorphisms with clinical characteristics (tumour location and regional lymph node involvement) and risk of ESCC.  interaction analysis to find out the best predictive model for the risk of ESCC. A three factors model TNFA-308, NFKBIA -826, NFKBIA 3'UTR polymorphisms was predicted to be best model with maximum testing accuracy and CV consistency. Furthermore, when we pooled the risk genotype (with OR>1) of all the polymorphisms, significant increased risk of ESCC of was observed with combined risk genotype in dose responsive manner. Also, individuals with 1-3 risk genotypes were found to have higher risk of ESCC compared to those with '0' risk genotype. However, no role of combined risk genotype on survival outcome of ESCC patients was found. These findings suggest the joint effect of studied gene polymorphisms in susceptibility but not in prognosis of ESCC. Interaction of selected polymorphisms with environmental risk factors were also analyzed, which did not show any significant outcome. This implies that there may be other classes of genes which might show interaction with environmental risk factor for developing risk of ESCC in northern Indian population. The present study has several strengths like well defined set of cases and controls, agreement of genotypic data with HWE in controls, consistencies of MAF of polymorphisms with Hapmap GIH data and adoption of stringent quality control measures. Limitation of study is low sample size in subgroups, absence of qualitative environmental exposure data in cases and short median follow-up in survival analysis.

Genotypes
In summary, our results suggest independent role of TNFA -308 G>A polymorphism and combined effect of TNFA and NFKBIA gene polymorphisms in susceptibility of ESCC in northern Indian population. However, none of the genetic variants seem to have implications in prognosis of ESCC. Table S1. Details of genotyping methods of TNFA-308 G>A, NFKB1 -94ATTG ins/del and NFKBIA (-826 C>T and 3'UTR A>G) polymorphisms. (DOCX)