Tumor Necrosis Factor-α Induced Protein 8 Polymorphism and Risk of Non-Hodgkin’s Lymphoma in a Chinese Population: A Case-Control Study

Abstract

Background

Non-Hodgkin’s lymphoma (NHL) has been reported to be associated with autoimmune and pro-inflammatory response, and genetic polymorphisms of candidate genes involved in autoimmune and pro-inflammatory response may influence the susceptibility to NHL. To evaluate the role of such genetic variations in risk of NHL, we conducted a case-control study of 514 NHL patients and 557 cancer-free controls in a Chinese population.

Method

We used the Taqman assay to genotype six potentially functional single nucleotide polymorphisms (SNPs) in six previously reported inflammation and immune-related genes (TNF rs1799964T>C, LTA rs1800683G>A, IL-10 rs1800872T>G, LEP rs2167270G>A, LEPR rs1327118C>G, TNFAIP8 rs1045241C>T). Logistic regression models were used to estimate odds ratios (ORs) and 95% confidence intervals (95% CI).

Results

We observed a significantly increased risk of NHL associated with the TNFAIP8 rs1045241C>T polymorphism (adjusted OR = 3.03; 95% CI = 1.68–5.45 for TT vs. CC and adjusted OR = 2.03; 95% CI = 1.53–2.69 for CT/TT vs. CC). The risk associated with the T allele was more evident in subgroups of 40–60 year-old, non-smokers or light-smokers (less than 25 pack-years), and subjects with normal weight or overweight. Risk for both B and T cell non-Hodgkin’s lymphoma was elevated for CT/TT genotypes (adjusted OR = 1.95, 95% CI = 1.41–2.70 for B cell NHL and adjusted OR = 2.22, 95% CI = 1.49–3.30 for T cell NHL), particularly for DLBCL (adjusted OR = 2.01, 95%CI = 1.41–2.85) and FL (adjusted OR = 2.53, 95% CI = 1.17–5.45). These risks were not observed for variant genotypes of other five SNPs compared with their common homozygous genotypes.

Conclusions

The polymorphism of TNFAIP8 rs1045241C>T may contribute to NHL susceptibility in a Chinese population. Further large-scale and well-designed studies are needed to confirm these results.

About the Authors:

Yan Zhang

Affiliations Department of Medical Oncology, Fudan University Shanghai Cancer Center, Shanghai, China, Department of Oncology, Shanghai Medical College, Fudan University, Shanghai, China

Meng-Yun Wang

Affiliations Department of Oncology, Shanghai Medical College, Fudan University, Shanghai, China, Cancer Research Laboratory, Fudan University Shanghai Cancer Center, Shanghai, China

Jing He

Affiliations Department of Oncology, Shanghai Medical College, Fudan University, Shanghai, China, Cancer Research Laboratory, Fudan University Shanghai Cancer Center, Shanghai, China

Jiu-Cun Wang

Affiliations Ministry of Education Key Laboratory of Contemporary Anthropology and State Key Laboratory of Genetic Engineering, School of Life Sciences, Fudan University, Shanghai, China, Fudan-Taizhou Institute of Health Sciences, Taizhou, Jiangsu, China

Ya-Jun Yang

Affiliations Ministry of Education Key Laboratory of Contemporary Anthropology and State Key Laboratory of Genetic Engineering, School of Life Sciences, Fudan University, Shanghai, China, Fudan-Taizhou Institute of Health Sciences, Taizhou, Jiangsu, China

Li Jin

Affiliations Ministry of Education Key Laboratory of Contemporary Anthropology and State Key Laboratory of Genetic Engineering, School of Life Sciences, Fudan University, Shanghai, China, Fudan-Taizhou Institute of Health Sciences, Taizhou, Jiangsu, China

Zhi-Yu Chen

Affiliations Department of Medical Oncology, Fudan University Shanghai Cancer Center, Shanghai, China, Department of Oncology, Shanghai Medical College, Fudan University, Shanghai, China

Xue-Jun Ma

Affiliations Department of Oncology, Shanghai Medical College, Fudan University, Shanghai, China, Department of Radiation Oncology, Fudan University Shanghai Cancer Center, Shanghai, China

Meng-Hong Sun

Affiliations Department of Oncology, Shanghai Medical College, Fudan University, Shanghai, China, Department of Pathology, Fudan University Shanghai Cancer Center, Shanghai, China

Kai-Qin Xia

Affiliations Department of Oncology, Shanghai Medical College, Fudan University, Shanghai, China, Cancer Research Laboratory, Fudan University Shanghai Cancer Center, Shanghai, China, Department of Pathology, Fudan University Shanghai Cancer Center, Shanghai, China

Xiao-Nan Hong

* E-mail: [email protected] (XYZ); [email protected] (XNH)

Affiliations Department of Medical Oncology, Fudan University Shanghai Cancer Center, Shanghai, China, Department of Oncology, Shanghai Medical College, Fudan University, Shanghai, China

Qing-Yi Wei

Affiliations Cancer Research Laboratory, Fudan University Shanghai Cancer Center, Shanghai, China, Department of Epidemiology, the University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America

Xiao-Yan Zhou

* E-mail: [email protected] (XYZ); [email protected] (XNH)

Affiliations Department of Oncology, Shanghai Medical College, Fudan University, Shanghai, China, Cancer Research Laboratory, Fudan University Shanghai Cancer Center, Shanghai, China, Department of Pathology, Fudan University Shanghai Cancer Center, Shanghai, China

Introduction

Non-Hodgkin’s lymphoma (NHL) incidence rates have been increasing in both developed and developing countries with about 355,900 new cases in the world annually [1]. In China, the most common subtype of NHL is diffuse large B cell lymphoma (DLBCL), whereas follicular lymphoma (FL) is less common than in Western countries. Extranodal lesions and T/NK cell NHLs (eg. Extranodal NK/T-cell lymphoma) appear to be more common in China [2]. However, the exact causes of NHL remain largely unknown. Some evidence has showed that immune dysfunction may be one of the risk factors [3], and single nucleotide polymorphisms (SNPs) in immune and inflammatory response genes may play an important role in lymphomagenesis [4], [5], [6], [7].

TNF/LTA and IL-10 genes code for immunoregulatory cytokines that can mediate inflammation, apoptosis and Th1/Th2 balance [6], [7], and they may be good candidate genes for studying lymphomagenesis. The cytokines TNF-α and LT-α are thought to influence lymphomagenesis through up-regulation of pro-inflammatory and anti-apoptotic signals, possibly via the NK-κB pathway [5]. Some evidence also showed that polymorphisms in IL-10 may regulate TNF-α levels and thus contribute to activation of the NK-κB pathway [8]. A pooled analysis including 7,999 NHL cases and 8,452 controls from 14 case-control studies was carried out by InterLymph Consortium, which showed that LTA 252A?G (rs909253), IL-1-3575T?A (rs1800890), and particularly TNF-308G?A (rs1800629) were associated with an increase risk of DLBCL in non-Hispanic white populations [9]. Purdue et al. reported that IL-10-3575T?A and TNF-863C?A (rs1800630) were associated with an elevated risk of DLBCL in an Australian case-control study [5]. A recent genome-wide association study (GWAS) of FL has identified additional two variants in the 6p21 chromosomal region [10], which is the TNF gene location, suggesting that genetic variants in these regions may influence NHL susceptibility.

The tumor necrosis factor-α induced protein 8 (TNFAIP8) family are newly identified proteins that are important for inflammation and immune homeostasis [11]. They play roles in anti-inflammation by negatively regulating T cell receptor (TCR) and Toll-like receptor (TLR) signaling [12]. But the association between TNFAIP8 polymorphisms and NHL risk has not been reported so far, particularly in Chinese populations.

The circulating levels of adipocytokines, including adiponectin, resistin and leptin may also alter immune system function and chronic inflammatory response. Leptin has pro-inflammatory properties and stimulates the growth of certain cancer cells as well as circulating pro-inflammatory cytokines, such as TNF-α and interleukin [13]. Associations between NHL and polymorphisms in the leptin (LEP) and leptin receptor (LEPR) gene have also been reported. Skibola et al. found that the LEP19G allele was associated with an increased risk of NHL, particularly FL [14]. A similar result was reported by a UK study, in which the LEPR 223ArgArg genotype was shown to be associated with an increased risk of FL among women [15].

To test the hypothesis that polymorphisms in inflammation and immune-related genes (such as TNF, LTA, IL-10, LEP, LEPR and TNFAIP8) may be associated with susceptibility to NHL, we conducted a case-control study in a Chinese population and genotyped six potentially functional SNPs in the above-mentioned candidate genes.

Materials and Methods

Ethics Statement

The study was approved by the Institutional Review Board of Fudan University Shanghai Cancer Center. Participation was voluntary. All participants singed a written informed consent, and all clinical investigation was conducted according to the principles expressed in the Declaration of Helsinki consent.

Study Population

The study population was identified from histologically confirmed NHL cases diagnosed and treated between June 2005 and September 2011 at Fudan University Shanghai Cancer Center. All patients came from the Eastern China, including Shanghai, Jiangsu Province and the surrounding regions. Enrollment criteria included the following: Chinese Han ethnicity, HIV antibody negative, and without a diagnosis of second primary malignancy. All cases were classified and reviewed according to the 2008 WHO classification of tumors of haematopoietic and lymphoid tissues [16]. Of the 732 eligible cases, 514 (70.2%) consented to participate in the study and provided blood samples. Additionally, 557 cancer-free control subjects were all selected from Taizhou Longitudinal Study (TZL) at the same time period. All the control subjects were frequency matched to the cases on age (according to cases’ age groups by every 5 years) and sex. The cancer-free controls were genetically unrelated ethnic Han Chinese who were not selected from family members of patients and had no blood relationship. All participants were all from Eastern China, and there was a major geographical overlap between the cases and controls. TZL was initiated in July, 2007, in Taizhou, Jiangsu province of China, as described previously [17].

Data Collection

After signing a written informed consent, all cases and controls were asked to provide a blood sample of about 2 ml and completed a questionnaire. Cases were asked about information including age, sex, ethnicity, status of smoking and drinking, the height and weight before diagnosis. The controls were also asked to recall the same questions prior to the recruitment date. A total of 514 cases and 557 controls provided blood samples, but 121 (23.5%) cases and 22 (3.9%) controls did not provide information about smoking and drinking, and 91 (17.7%) cases failed to give the exact information of height and weight before diagnosis.

SNP Selection

We first screened five frequently reported (TNF, LTA, IL-10, LEP, LEPR) and one newly identified (TNFAIP8) inflammation and immune-related genes from published papers. We then searched SNPs in these six genes by NCBI dbSNP database (http://www.ncbi.nlm.nih.gov/) for all populations, and used HapMap database (http://hapmap.ncbi.nlm.nih.gov/) to focus on CHB (Han Chinese in Beijing, China) population. SNPinfo (http://snpinfo.niehs.nih.gov/) website was used as a tool to predict SNP functions. Potentially functional SNPs were defined to fit at least three of following four criteria: (1) SNPs located at the two ends of these genes, such as 5′ near gene, 5′ untranslated region (UTR), 3′ near gene or 3′UTR; (2) the minor allele frequency (MAF) was≥5% in the Hapmap CHB population. (3) SNP variation may affect transcription factor binding sites (TFBS) activity in the putative promoter region (here defined as 2-kb upstream from the first exon). (4) SNP may affect the microRNA-binding sites activity. Then eight SNPs were identified (TNF rs1799964T>C, LTA rs1800683G>A, IL-10 rs1800872T>G, LEP rs2167270G>A, LEP rs4728096 T>C, LEPR rs1327118C>G, TNFAIP8 rs1045241C>T, TNFAIP8 rs11064A>G). However, the linkage disequilibrium (LD) analysis revealed that LEP rs2167270G>A and LEP rs4728096 T>C were in quite high LD (r2 = 0.93), and TNFAIP8 rs1045241C>T and rs11064A>G were also in LD with r2 = 0.85. As a result, we selected following six potentially functional and representative SNPs (TNF rs1799964T>C, LTA rs1800683G>A, IL-10 rs1800872T>G, LEP rs2167270G>A, LEPR rs1327118C>G, TNFAIP8 rs1045241C>T) for further genotyping.

Genotyping

Genomic DNA was extracted from each blood sample by using the Qiagen Blood DNA Mini Kit (Qiagen Inc., Valencia, CA) according to the manufacturer’s instructions. DNA purity and concentration were determined by spectrophotometric measurement of absorbance at 260 and 280 nm by a UV spectrophotometer (Nano Drop Technologies, Inc., Wilmington, DE) and all are qualified.

All TaqMan assays for this study including the pre-designed SNP-genotyping assay mix containing PCR primers and probes were purchased from ABI (Applied Biosystems, Foster City, CA). Genotyping were conducted on the ABI 7900HT detection system (Applied Biosystems). To ensure the accuracy of genotyping results, four negative controls (no DNA) and four duplicated samples were included in each of the 384-well plates for the quality control. The assays were repeated for 5% of the samples, and the results were 100% concordant. The analyzed fluorescence results were then auto-called into the genotypes using the built-in SDS2.2 software of the system.

Statistical Analysis

BMI was calculated as weight (kg) divided by the square of the height (m). In this study, we used the BMI cutoff points as suggested by the Cooperative Meta-Analysis Group of Working Group on Obesity in China [18]. If BMI<18.5 kg/m2, the individuals were defined as lower than normal weight, 18.5≤BMI≤24.0 kg/m2 was defined as normal weight and BMI>24.0 kg/m2 was defined as overweight. Smoking status was divided into smokers and non-smokers by whether or not they had smoked for more than one year. Those who drank alcoholic beverages at least once a week for one year or more were defined as alcohol users, while the others were non-users. Differences in the distributions of the alleles and genotypes as well as demographic characters, smoking status, alcohol use and BMI between the cases and controls were evaluated by the Chi-square test. The Hardy–Weinberg equilibrium (HWE) of genotype distribution in the controls was tested by a goodness-of-fit Chi-square test. Unconditional univariate and multivariate logistic regression models were applied to calculate crude and adjusted odds ratios (ORs) and 95% confidence intervals (95% CI). We used the multiple imputation (MI) method by SAS 9.1 software to handle missing covariate information. All covariates were imputed, when calculated adjusted ORs, according to the distributions of observed values of either cases or controls. All statistical tests were two-sided, and P<0.05 was considered statistically significant. All analyses were performed using SAS Software, version 9.1 (SAS Institute, Cary, NC).

Results

Characteristics of the Study Population

There were 514 NHL cases and 557 cancer-free controls included in this study, whose DNA samples were available. The frequency distributions of demographic and some selected characteristics of the participants are shown in Table 1. There was no statistical difference in the distributions of age and sex between cases and controls because of frequency matching by design. The mean age was 49.3 years for the cases (±14.1; range, 15–85) and 49.6 years for the controls (±13.5; range, 20–85; P = 0.895), and 64.0% of the cases and 63.4% of the controls were male (P = 0.830). The controls were more likely to be smokers and alcohol users (P<0.0001 and P = 0.0001, respectively) and have higher BMI (P<0.0001) compared with the cases. Consequently, smoking status, alcohol use and BMI were adjusted for in the subsequent multivariate logistic regression analyses. Of the 514 cases, 336 (65.4%) were B cell lymphoma, and 178 (34.6%) were T cell and natural killer (NK) cell lymphoma. After stratified, 233 (45.3%), 52 (10.1%), 51 (9.9%), 146 (28.4%), 32 (6.2%) were DLBCL, FL, other B cell lymphoma, NK/T cell lymphoma, and other T cell lymphoma, respectively. Of all cases, 322 (62.6%) had an Ann Arbor Stage of I–II and 192 (37.4%) had a later Ann Arbor Stage of III–IV.

[Figure omitted. See PDF.]

Table 1. Characteristics of Non-Hodgkin’s lymphoma cases and cancer-free controls.

https://doi.org/10.1371/journal.pone.0037846.t001

Association between Selected SNPs and Risk of NHL

Genotype distributions of the selected six SNPs in cases and controls and their associations with NHL risk are presented in Table 2. All observed genotype distributions among controls were in agreement with the Hardy–Weinberg equilibrium (P = 0.285 for rs1799964, P = 0.777 for rs180683, P = 0.872 for rs1800872, P = 0.733 for rs2167270, P = 0.559 for rs1327118, P = 0.420 for rs1045241). Significant difference in the genotype frequencies was observed between the cases and controls for TNFAIP8 rs1045241C>T (P<0.0001). When the rs1045241CC genotype was used as the reference, the CT heterozygous, TT homozygous and combined CT/TT genotypes were associated with significantly increased risk of NHL (adjusted OR = 1.88; 95% CI = 1.39–2.53 for CT heterozygous genotype, adjusted OR = 3.03; 95% CI = 1.68–5.45 for TT homozygous genotype, adjusted OR = 2.03; 95% CI = 1.53–2.69 for CT/TT genotype) after adjustment for age, sex, BMI, smoking and drinking status. However, no significantly altered NHL risk was found for variant genotypes of the other five SNPs compared with their common genotypes. We also evaluated the combined effect of all the six SNPs. We had divided the subjects into seven groups according to the number of combined variant genotypes. The “0” risk genotype group was used as the reference, and unconditional logistic regression models were applied to calculate OR and 95%CI for each group. But we did not find any significant association between the combined effect of these six SNPs and risk of NHL (data not shown).

[Figure omitted. See PDF.]

Table 2. Genotypes distributions of the selected functional polymorphisms among NHL cases and cancer-free controls and their associations with NHL risk.

https://doi.org/10.1371/journal.pone.0037846.t002

Stratified Analysis

We further evaluated the association between the six candidate SNPs and risk of NHL by subgroups of age, sex, smoking status, alcohol use, BMI, common subtypes and Ann Arbor stage. No statistical significances were found for the SNPs except TNFAIP8 rs1045241 C>T. The stratified analysis results of TNFAIP8 rs1045241 C>T are presented in Table 3. In general, an increased risk associated with rs1045241 CT/TT genotypes was more evident in subgroups of 40–60 year-old individuals (adjusted OR = 2.99, 95% CI = 2.03–4.43), non-smokers (adjusted OR = 1.90, 95% CI = 1.35–2.68) or smoked less than 25 pack-years (adjusted OR = 2.45, 95% CI = 1.32–4.54), and normal weight (adjusted OR = 1.76, 95% CI = 1.23–2.52) or overweight groups (adjusted OR = 2.07, 95% CI = 1.32–3.25). Moreover, the patients with rs1045241 CT/TT genotypes were associated with risk of both B and T cell lymphoma (adjusted OR = 1.95, 95% CI = 1.41–2.70 and adjusted OR = 2.22, 95% CI = 1.49–3.30, respectively). After stratifying by histological subtypes, we observed increased risk for DLBCL (adjusted OR = 2.01, 95%CI = 1.41–2.85) and FL (adjusted OR = 2.53, 95%CI = 1.17–5.45), but no significant association was found for NK/T cell lymphoma (adjusted OR = 0.79, 95%CI = 0.27–2.33). However, after we verified homogeneity assumption by using a Chi square-based Q-test. The results indicated that an increased NHL risk associated with CT/TT genotypes was particularly more pronounced only in subgroups of 40–60 year-old (P = 0.03), and patients with a later Ann Arbor stage of III–IV (P = 0.04).

[Figure omitted. See PDF.]

Table 3. Stratification analysis of the association between TNFAIP8 rs1045241C>T and NHL risk.

https://doi.org/10.1371/journal.pone.0037846.t003

Discussion

Genetic polymorphisms in immune-related genes that regulate the immune and inflammation response may play an important role in the incidence of NHL [19], [20]. In this case-control study, we reported that TNFAIP8 rs1045241C>T was significantly associated with an increased risk of NHL in a Chinese population. Stratified analyses revealed that subgroups of 40–60 years, non-smokers or light-smokers (≤25 pack-years), subjects with normal weight or overweight were more likely to have been diagnosed with NHL, especially the subtypes of DLBCL and FL. According to epidemiologic data about NHL in China, the main subtypes of lymphoma are DLBCL, FL and NK/T lymphoma, and other subtypes have a small sample size, and therefore we only analyzed the dominant subtypes of the cases. These results support the hypothesis that some polymorphism in inflammation and immune-related genes may be associated with risk of NHL and its common subtypes.

TNFAIP8, also known as GG2-1, MDC-3.13, SCC-S2, is located on chromosome 5 (5q23.1). It was first identified in a human head and neck squamous cell carcinoma (HNSCC) cell line [21]. Recently, it has been reported to be associated with oncogenesis, immunity, and inflammation in several studies [11], [12], [22]. TNFAIP8 mRNA over-expression has been found in various malignant cell lines, such as breast cancer [23], non-small cell lung cancer [24], and esophageal squamous cell carcinoma [25]. Evidence showed that TNFAIP8 expression was upregulated by TNF-α induced NF-κB pathway activation in cancer cell lines, which can inhibit caspase-8 and reduce apoptosis [26], [27]. Tumor necrosis factor-α induced protein 8-like 2 (TIPE2), a member of the TNFAIP8 family, was originally identified as a gene abnormally expressed in the inflamed spinal cord of mice with experimental autoimmune encephalomyelitis [12]. TIPE2-deficent mice were more likely to suffer from chronic inflammation diseases and multiple organ inflammation [12]. In humans, the abnormal expression of TIPE2 was associated with systemic autoimmunity [28], diabetic nephropathy [29], and hepatitis B [30]. These studies supported that TIPE2 plays an important role in maintaining immune homeostasis.

However, fewer studies have focused on the association between polymorphisms of TNFAIP8 and NHL risk. The SNP TNFAIP8 rs1045241C>T included in the present study is located on the 3′UTR of TNFAIP8, and the SNP function prediction shows that it may affect the microRNA-binding sites activity (http://snpinfo.niehs.nih.gov/snpfunc.htm). Some evidence has indicted that a genetic polymorphism in a microRNA target site influence transcriptional and post-transcriptional gene expression in cancers [31]. Though we found that rs1045241C>T was associated with the risk of NHL, the exact mechanisms of this relationship remained unknown. Further functional studies have been conducting to confirm our results.

In the present study, we did not find any significant association between the investigated polymorphisms of TNF rs1799964T>C, LTA rs1800683G>A, IL-10 rs1800872T>G, LEP rs2167270G>A, LEPR rs1327118C>G and NHL risk in a Chinese population. However, many pieces of evidence have shown that other SNPs in these genes like TNF–308 G>A (rs1800629), LTA 252A>G (rs909253), IL10–3575T>A (rs1800890), LEPR 233Q>R (rs1137101) increased the risk of NHL, particularly DLBCL in non-Hispanic white populations [9], [14], [15], [32], [33]. Our studies did not cover these SNPs because their MAF was <5% in CHB population according to Hapmap database (http://hapmap.ncbi.nlm.nih.gov/). But Xiao et al. reported that no association between TNF–308, LTA 252 polymorphisms and histological subtypes, disease stages in Chinese NHL patients [34]. We believe that this inconsistency may be due to ethnic and demographic differences between Chinese populations and white populations as well as the smaller sample size that may also have led to insufficient statistical power. So our results need to be replicated in additional studies in different populations with large sample sizes.

Fewer studies reported an association between LEP 19 A>G, LEPR 233Q>R and risk of NHL in Chinese populations, though some previous published studies revealed that obesity was a risk factor of NHL in developed countries [35]. Obesity may play the role of inducing lymphomagenesis by the co-effects of LEP 19 A>G, LEPR 233Q>R polymorphisms and immune dysfunction [14], [15], [36]. Inconsistently, we found the NHL cases were more likely to have lower BMI compared with the cancer-free controls (P<0.0001). One possible reason may due to ethnic difference in BMI and its contribution to lymphoma. It is likely that populations of developing and developed countries have different body size, and Asians have lower BMI than Westerners. Thus the Chinese standard of obesity was applied to the subjects in this study, but obesity did not appear closely related to lymphomagenesis. However, another possibility was that obesity may not a direct risk factor for NHL in Chinese populations. There may be some intricate interaction among obesity, diet, hormone, genetics, and susceptibility to NHL but the present study did not have enough power to detect such interactions. Further larger, transcontinental studies are needed to fully explore and elaborate the extent of such interactions.

In the stratified analyses, the risk effect of the TNFAIP8 rs1045241 CT/TT genotype was more evident in subgroups of 40–60 year-old, non-smoker or light-smoker, and normal BMI or overweight. Stratifying by subtypes and stage, the risk effect of rs1045241 CT/TT genotype was evident for both B and T cell lymphoma, especially for DLBCL and FL as well as later Ann Arbor Stage of ?–?. But we noticed that after verification of homogeneity assumption, only age and stage showed marginally statistical significance. So, there was no strong evidence of an effect between the selected covariates and NHL risk. The exact association between selected covariates and NHL risk needs to be further investigated by large studies. Up to now, there is still inconsistence about the effect of smoking, alcohol use, and BMI on NHL risk, and the mechanisms are multiple. Wang et al. once reported that autoimmune conditions, obesity, and later birth order could contribute to lymphomagensis through an alteration of the proinflammatory pathway, specifically involving common genetic variants in TNF and IL-10 [8]. Given the results of our study, NHL development may be of complex process that includes a great number of events. A single factor would have a limited effect on the susceptibility. Once the exposure of environmental factors was accumulated to a certain level, alteration took place in immune and proinflammatory pathways, and genetic variants or other unknown events may also participate in the process. On the other hand, our study was not large enough, and possible selection bias may exist due to a higher-than expected non-response rate in case group. Well-designed large studies are needed to test the hypothesis that TNFAIP8 polymorphisms may influence NHL susceptibility, particularly by the possible interactions with environmental risk factors.

Several limitations exist in the present study. First, the sample size was not large enough, so selection bias may be inevitable, and the study did not provide enough statistical power to detect gene-gene and gene-environment interaction. Second, we selected only one functional SNP for each candidate gene, which restricted further analysis to identify other potentially important associations. Moreover, different histological subtypes of NHL, such as DLBCL, FL, and NK/T cell lymphoma may have different risk factors and etiologies. We did not have enough information about environmental exposure from NHL patients, such as use of hair dye, family history of NHL, occupation and EB-virus infection to be included in stratified analyses.

In summary, we identified TNFAIP8 rs1045241C>T to be associated with risk of NHL in a Chinese population, particularly DLBCL and FL. To the best of our knowledge, this is the first study to report a TNFAIP8 polymorphism associated with NHL susceptibility in a Chinese population, which provides additional evidence that inflammation and immune-related genetic polymorphisms play an important role in lymphomagenesis. Further large-scale well-designed population-based studies are needed to validate our findings, define the high-risk population and deepen our understanding of NHL pathogenesis.

Acknowledgments

We thank Yan-Zi Gu, Huan Chen, Ting-Yan Shi, and Mei-Ling Zhu for their help in collecting blood samples and DNA extraction, and Yu-Hu Xin for his assistance in technical support.

Author Contributions

Conceived and designed the experiments: XYZ. Performed the experiments: YZ KQX. Analyzed the data: YZ MYW JH. Contributed reagents/materials/analysis tools: JCW YJY LJ ZYC XJM MHS QYW. Wrote the paper: YZ XYZ QYM. Supplied clinical information: XNH.

Citation: Zhang Y, Wang M-Y, He J, Wang J-C, Yang Y-J, Jin L, et al. (2012) Tumor Necrosis Factor-α Induced Protein 8 Polymorphism and Risk of Non-Hodgkin’s Lymphoma in a Chinese Population: A Case-Control Study. PLoS ONE7(5): e37846. https://doi.org/10.1371/journal.pone.0037846

References

1. Jemal A, Bray F, Center MM, Ferlay J, Ward E, et al. (2011) Global cancer statistics. CA Cancer J Clin 61: 69–90.A. JemalF. BrayCenter MM, Ferlay J, Ward E, et al.2011Global cancer statistics.CA Cancer J Clin616990

2. Li X, Li G, Gao Z, Zhou X, Zhu X (2011) The relative frequencies of lymphoma subtypes in China: A nationwide Study of 10002 Cases by the Chinese lymphoma study group. Annals of Oncology 22: iv141.X. LiG. LiZ. GaoX. ZhouX. Zhu2011The relative frequencies of lymphoma subtypes in China: A nationwide Study of 10002 Cases by the Chinese lymphoma study group.Annals of Oncology22iv141

3. Fernberg P, Chang ET, Duvefelt K, Hjalgrim H, Eloranta S, et al. (2010) Genetic variation in chromosomal translocation breakpoint and immune function genes and risk of non-Hodgkin lymphoma. Cancer Causes Control 21: 759–769.P. FernbergET ChangK. DuvefeltH. HjalgrimS. Eloranta2010Genetic variation in chromosomal translocation breakpoint and immune function genes and risk of non-Hodgkin lymphoma.Cancer Causes Control21759769

4. Dong LM, Potter JD, White E, Ulrich CM, Cardon LR, et al. (2008) Genetic Susceptibility to Cancer: The Role of Polymorphisms in Candidate Genes. JAMA: The Journal of the American Medical Association 299: 2423–2436.LM DongJD PotterE. WhiteCM UlrichLR Cardon2008Genetic Susceptibility to Cancer: The Role of Polymorphisms in Candidate Genes.JAMA: The Journal of the American Medical Association29924232436

5. Purdue MP, Lan Q, Kricker A, Grulich AE, Vajdic CM, et al. (2006) Polymorphisms in immune function genes and risk of non-Hodgkin lymphoma: findings from the New South Wales non-Hodgkin Lymphoma Study. Carcinogenesis 28: 704–712.MP PurdueQ. LanA. KrickerAE GrulichCM Vajdic2006Polymorphisms in immune function genes and risk of non-Hodgkin lymphoma: findings from the New South Wales non-Hodgkin Lymphoma Study.Carcinogenesis28704712

6. Lan Q (2006) Cytokine polymorphisms in the Th1/Th2 pathway and susceptibility to non-Hodgkin lymphoma. Blood 107: 4101–4108.Q. Lan2006Cytokine polymorphisms in the Th1/Th2 pathway and susceptibility to non-Hodgkin lymphoma.Blood10741014108

7. Lan Q, Wang SS, Menashe I, Armstrong B, Zhang Y, et al. (2011) Genetic variation in Th1/Th2 pathway genes and risk of non-Hodgkin lymphoma: a pooled analysis of three population-based case-control studies. British Journal of Haematology 153: 341–350.Q. LanSS WangI. MenasheB. ArmstrongY. Zhang2011Genetic variation in Th1/Th2 pathway genes and risk of non-Hodgkin lymphoma: a pooled analysis of three population-based case-control studies.British Journal of Haematology153341350

8. Wang SS, Cozen W, Cerhan JR, Colt JS, Morton LM, et al. (2007) Immune Mechanisms in Non-Hodgkin Lymphoma: Joint Effects of the TNF G308A and IL10 T3575A Polymorphisms with Non-Hodgkin Lymphoma Risk Factors. Cancer Research 67: 5042–5054.SS WangW. CozenJR CerhanJS ColtLM Morton2007Immune Mechanisms in Non-Hodgkin Lymphoma: Joint Effects of the TNF G308A and IL10 T3575A Polymorphisms with Non-Hodgkin Lymphoma Risk Factors.Cancer Research6750425054

9. Rothman N, Skibola CF, Wang SS, Morgan G, Lan Q, et al. (2006) Genetic variation in TNF and IL10 and risk of non-Hodgkin lymphoma: a report from the InterLymph Consortium. The Lancet Oncology 7: 27–38.N. RothmanCF SkibolaSS WangG. MorganQ. Lan2006Genetic variation in TNF and IL10 and risk of non-Hodgkin lymphoma: a report from the InterLymph Consortium.The Lancet Oncology72738

10. Conde L, Halperin E, Akers NK, Brown KM, Smedby KE, et al. (2010) Genome-wide association study of follicular lymphoma identifies a risk locus at 6p21.32. Nature Genetics 42: 661–664.L. CondeE. HalperinNK AkersKM BrownKE Smedby2010Genome-wide association study of follicular lymphoma identifies a risk locus at 6p21.32.Nature Genetics42661664

11. Lou Y, Liu S (2011) The TIPE (TNFAIP8) family in inflammation, immunity, and cancer. Molecular Immunology. Y. LouS. Liu2011The TIPE (TNFAIP8) family in inflammation, immunity, and cancer.Molecular Immunology

12. Sun H, Gong S, Carmody RJ, Hilliard A, Li L, et al. (2008) TIPE2, a Negative Regulator of Innate and Adaptive Immunity that Maintains Immune Homeostasis. Cell 133: 415–426.H. SunS. GongRJ CarmodyA. HilliardL. Li2008TIPE2, a Negative Regulator of Innate and Adaptive Immunity that Maintains Immune Homeostasis.Cell133415426

13. Tilg H, Moschen AR (2006) Adipocytokines: mediators linking adipose tissue, inflammation and immunity. Nature Reviews Immunology 6: 772–783.H. TilgAR Moschen2006Adipocytokines: mediators linking adipose tissue, inflammation and immunity.Nature Reviews Immunology6772783

14. Skibola CF, Holly EA, Forrest MS, Hubbard A, Bracci PM, et al. (2004) Body Mass Index, Leptin and Leptin Receptor Polymorphisms, and Non-Hodgkin Lymphoma.pdf>. Cancer Epidemiol Biomarkers Prev 13: 779–786.CF SkibolaEA HollyMS ForrestA. HubbardPM Bracci2004Body Mass Index, Leptin and Leptin Receptor Polymorphisms, and Non-Hodgkin Lymphoma.pdf>.Cancer Epidemiol Biomarkers Prev13779786

15. Willett EV, Skibola CF, Adamson P, Skibola DR, Morgan GJ, et al. (2005) Non-Hodgkin's lymphoma, obesity and energy homeostasis polymorphisms. British Journal of Cancer 93: 811–816.EV WillettCF SkibolaP. AdamsonDR SkibolaGJ Morgan2005Non-Hodgkin's lymphoma, obesity and energy homeostasis polymorphisms.British Journal of Cancer93811816

16. Swerdlow SH, Campo E, NL H (2008) World Health Organisation classification of tumours of haematopoietic and lymphoid tissues. Lyon: IARC. SH SwerdlowE. CampoH. NL2008World Health Organisation classification of tumours of haematopoietic and lymphoid tissues.LyonIARC

17. Wang X, Lu M, Qian J, Yang Y, Li S, et al. (2009) Rationales, design and recruitment of the Taizhou Longitudinal Study. BMC Public Health 9: 223.X. WangM. LuJ. QianY. YangS. Li2009Rationales, design and recruitment of the Taizhou Longitudinal Study.BMC Public Health9223

18. Zhou BF (2002) Predictive values of body mass index and waist circumference for risk factors of certain related diseases in Chinese adults–study on optimal cut-off points of body mass index and waist circumference in Chinese adults. Biomed Environ Sci 15: 83–96.BF Zhou2002Predictive values of body mass index and waist circumference for risk factors of certain related diseases in Chinese adults–study on optimal cut-off points of body mass index and waist circumference in Chinese adults.Biomed Environ Sci158396

19. Wang SS, Cerhan JR, Hartge P, Davis S, Cozen W, et al. (2006) Common Genetic Variants in Proinflammatory and Other Immunoregulatory Genes and Risk for Non-Hodgkin Lymphoma. Cancer Research 66: 9771–9780.SS WangJR CerhanP. HartgeS. DavisW. Cozen2006Common Genetic Variants in Proinflammatory and Other Immunoregulatory Genes and Risk for Non-Hodgkin Lymphoma.Cancer Research6697719780

20. Liang X, Caporaso N, McMaster ML, Ng D, Landgren O, et al. (2009) Common genetic variants in candidate genes and risk of familial lymphoid malignancies. British Journal of Haematology 146: 418–423.X. LiangN. CaporasoML McMasterD. NgO. Landgren2009Common genetic variants in candidate genes and risk of familial lymphoid malignancies.British Journal of Haematology146418423

21. Patel S, Wang FH, Whiteside TL, Kasid U (1997) Identification of seven differentially displayed transcripts in human primary and matched metastatic head and neck squamous cell carcinoma cell lines: implications in metastasis and/or radiation response. Oral Oncol 33: 197–203.S. PatelFH WangTL WhitesideU. Kasid1997Identification of seven differentially displayed transcripts in human primary and matched metastatic head and neck squamous cell carcinoma cell lines: implications in metastasis and/or radiation response.Oral Oncol33197203

22. Freundt EC, Bidere N, Lenardo MJ (2008) A Different TIPE of Immune Homeostasis. Cell 133: 401–402.EC FreundtN. BidereMJ Lenardo2008A Different TIPE of Immune Homeostasis.Cell133401402

23. Kumar D, Gokhale P, Broustas C, Chakravarty D, Ahmad I, et al. (2004) Expression of SCC-S2, an antiapoptotic molecule, correlates with enhanced proliferation and tumorigenicity of MDA-MB 435 cells. Oncogene 23: 612–616.D. KumarP. GokhaleC. BroustasD. ChakravartyI. Ahmad2004Expression of SCC-S2, an antiapoptotic molecule, correlates with enhanced proliferation and tumorigenicity of MDA-MB 435 cells.Oncogene23612616

24. Dong Q-Z, Zhao Y, Liu Y, Wang Y, Zhang P-X, et al. (2010) Overexpression of SCC-S2 correlates with lymph node metastasis and poor prognosis in patients with non-small-cell lung cancer. Cancer Science 101: 1562–1569.Q-Z DongY. ZhaoY. LiuY. WangP-X Zhang2010Overexpression of SCC-S2 correlates with lymph node metastasis and poor prognosis in patients with non-small-cell lung cancer.Cancer Science10115621569

25. Hadisaputri YE, Miyazaki T, Suzuki S, Yokobori T, Kobayashi T, et al. (2011) TNFAIP8 Overexpression: Clinical Relevance to Esophageal Squamous Cell Carcinoma. Annals of Surgical Oncology. YE HadisaputriT. MiyazakiS. SuzukiT. YokoboriT. Kobayashi2011TNFAIP8 Overexpression: Clinical Relevance to Esophageal Squamous Cell Carcinoma.Annals of Surgical Oncology

26. You Z, Ouyang H, Lopatin D, Polver PJ, Wang CY (2001) Nuclear factor-kappa B-inducible death effector domain-containing protein suppresses tumor necrosis factor-mediated apoptosis by inhibiting caspase-8 activity. J Biol Chem 276: 26398–26404.Z. YouH. OuyangD. LopatinPJ PolverCY Wang2001Nuclear factor-kappa B-inducible death effector domain-containing protein suppresses tumor necrosis factor-mediated apoptosis by inhibiting caspase-8 activity.J Biol Chem2762639826404

27. Kumar D, Whiteside TL, Kasid U (2000) Identification of a Novel Tumor Necrosis Factor-a-inducible Gene,SCC-S2, Containing the Consensus Sequence of a Death Effector Domain of Fas-associated Death Domain-like Interleukin-b-converting Enzyme-inhibitory Protein. J Biol Chem 275: 2973–2978.D. KumarTL WhitesideU. Kasid2000Identification of a Novel Tumor Necrosis Factor-a-inducible Gene,SCC-S2, Containing the Consensus Sequence of a Death Effector Domain of Fas-associated Death Domain-like Interleukin-b-converting Enzyme-inhibitory Protein.J Biol Chem27529732978

28. Li A, Chen Y, Zhao X, Niu Y, Cong P, et al. (2010) Characterization and transcriptional regulation analysis of the porcine TNFAIP8L2 gene. Molecular Genetics and Genomics 284: 185–195.A. LiY. ChenX. ZhaoY. NiuP. Cong2010Characterization and transcriptional regulation analysis of the porcine TNFAIP8L2 gene.Molecular Genetics and Genomics284185195

29. Zhang S, Zhang Y, Wei X, Zhen J, Wang Z, et al. (2010) Expression and regulation of a novel identified TNFAIP8 family is associated with diabetic nephropathy. Biochimica et Biophysica Acta (BBA) - Molecular Basis of Disease 1802: 1078–1086. pp. 1078–1086.S. ZhangY. ZhangX. WeiJ. ZhenZ. Wang2010Expression and regulation of a novel identified TNFAIP8 family is associated with diabetic nephropathy.Biochimica et Biophysica Acta (BBA) - Molecular Basis of Disease 18021078–108610781086

30. Xi W, Hu Y, Liu Y, Zhang J, Wang L, et al. (2011) Roles of TIPE2 in hepatitis B virus-induced hepatic inflammation in humans and mice. Molecular Immunology 48: 1203–1208.W. XiY. HuY. LiuJ. ZhangL. Wang2011Roles of TIPE2 in hepatitis B virus-induced hepatic inflammation in humans and mice.Molecular Immunology4812031208

31. Yu Z, Li Z, Jolicoeur N, Zhang L, Fortin Y, et al. (2007) Aberrant allele frequencies of the SNPs located in microRNA target sites are potentially associated with human cancers. Nucleic Acids Research 35: 4535–4541.Z. YuZ. LiN. JolicoeurL. ZhangY. Fortin2007Aberrant allele frequencies of the SNPs located in microRNA target sites are potentially associated with human cancers.Nucleic Acids Research3545354541

32. Skibola CF, Bracci PM, Nieters A, Brooks-Wilson A, de Sanjose S, et al. (2010) Tumor Necrosis Factor (TNF) and Lymphotoxin- (LTA) Polymorphisms and Risk of Non-Hodgkin Lymphoma in the InterLymph Consortium. American Journal of Epidemiology 171: 267–276.CF SkibolaPM BracciA. NietersA. Brooks-WilsonS. de Sanjose2010Tumor Necrosis Factor (TNF) and Lymphotoxin- (LTA) Polymorphisms and Risk of Non-Hodgkin Lymphoma in the InterLymph Consortium.American Journal of Epidemiology171267276

33. Wang SS, Purdue MP, Cerhan JR, Zheng T, Menashe I, et al. (2009) Common Gene Variants in the Tumor Necrosis Factor (TNF) and TNF Receptor Superfamilies and NF-kB Transcription Factors and Non-Hodgkin Lymphoma Risk. PLoS One 4: e5360.SS WangMP PurdueJR CerhanT. ZhengI. Menashe2009Common Gene Variants in the Tumor Necrosis Factor (TNF) and TNF Receptor Superfamilies and NF-kB Transcription Factors and Non-Hodgkin Lymphoma Risk.PLoS One4e5360

34. Xiao H, Zhang K (2010) Genetic polymorphisms of tumor necrosis factor-α and lymphotoxin-α in Chinese patients with non-Hodgkin lymphoma. Annals of Hematology 90: 725–727.H. XiaoK. Zhang2010Genetic polymorphisms of tumor necrosis factor-α and lymphotoxin-α in Chinese patients with non-Hodgkin lymphoma.Annals of Hematology90725727

35. Larsson SC, Wolk A (2007) Obesity and risk of non-Hodgkin's lymphoma: A meta-analysis. International Journal of Cancer 121: 1564–1570.SC LarssonA. Wolk2007Obesity and risk of non-Hodgkin's lymphoma: A meta-analysis.International Journal of Cancer12115641570

36. Chen Y, Zheng T, Lan Q, Foss F, Kim C, et al. (2010) Cytokine polymorphisms in Th1/Th2 pathway genes, body mass index, and risk of non-Hodgkin lymphoma. Blood 117: 585–590.Y. ChenT. ZhengQ. LanF. FossC. Kim2010Cytokine polymorphisms in Th1/Th2 pathway genes, body mass index, and risk of non-Hodgkin lymphoma.Blood117585590