Association between five types of Tumor Necrosis Factor-α Gene Polymorphism and Hepatocellular Carcinoma Risk: A Meta-Analysis

DOI: https://doi.org/10.21203/rs.3.rs-26916/v1

Abstract

Background: Research focusing on the relationship between five types of tumor necrosis factor-alpha (TNF-α) SNPs and the risk of hepatocellular carcinoma (HCC) were still controversial. Hereby, we performed a meta-analysis to determine the association between TNF-α promoter SNPs: -1031 T/C, -863 C/A, -857 C/T, -308 G/A, and -238 G/A with HCC risk.

Methods: We interrogated articles from journal database: PubMed, Pro-Quest, EBSCO, Science Direct, and Springer to determine the relationship between five types of SNPs in TNF-α gene with HCC risk. RevMan 5.3 software was used for analysis in fixed/random effect models.

Results: This meta-analysis included 23 potential articles from 2004-2018 with 3,237 HCC cases and 4,843 controls. We found that SNP -863 C/A were associated with a significantly increased HCC risk (A vs C, OR=1.31, 95% CI=1.03-1.67; CA/AA vs CC, OR=1.19, 95% CI=1.03-1.36). Similar results were obtained in -857 C/T (TT/CT vs CC, OR=1.31, 95% CI=1.06-1.62), -308 G/A (G vs A, OR=1.98, 95% CI=1.62-2.42; AA/GA vs GG, OR=1.95, 95% CI=1.53-2.49; GG/GA vs AA, OR=2.52, 95% CI=1.69-3.76; AA vs GG, OR=3.14, 95% CI=2.06-4.79; and GA vs GG, OR=2.07, 95% CI=1.60-2.68), and -238 G/A (A vs G, OR=1.50, 95% CI=1.16-1.94; AA vs GG, OR=3.87, 95% CI=1.32-11.34; GA/GG vs AA, OR=2.67, 95% CI=1.17-6.10). While no associations were observed between SNP TNF-α -1031 T/C and HCC risk.

Conclusions: The present meta-analysis showed that TNFα SNPs -863C/A, -857 C/T, -308 G/A, and -238 G/A were associated with the risk of HCC.

Background

Hepatocellular carcinoma (HCC) or liver cancer accounts for the cancer with the fourth highest mortality rate in the world. In 2018, there are estimated to be 841,080 of HCC new cases worldwide [1]. Patients with HCC usually have poor prognosis and high mortality rates, even in developed countries [2]. Until now, the pathogenesis of HCC has not been fully understood, but it is known so far that it is influenced by hepatitis B and C virus infections and also influenced by environmental factors (smoking, alcohol, aflatoxin B1) [3]. It is known that there are differences in the risk of HCC within each person, in which case the host factor has an important role [4].

Tumor Necrosis Factor-α (TNF-α) is an important inflammatory cytokine in the development of liver disease. This cytokine can cause hepatic injury, cirrhosis and eventually promote hepatocellular carcinoma [5, 6]. Several previous studies have identified some Single Nucleotide Polymorphism (SNP)s in TNF-α gene, especially in the promoter region. SNPs of TNF-α -1031T/C, -863C/A, -857C/T, -308G/A, and − 238G/A are SNPs in the TNF-α promoter site that have often been investigated regarding their association with HCC in several previous studies [7, 8]. It was also said that those SNPs could affect TNF-α production at the transcription level [9, 10].

High production of TNF-is related to the increase of pro-inflammatory cytokine secretion, the activation of proto oncogenes and several genes associated with cell growth, invasion, and cancer cells metastasis [11, 12]. Excessive production of TNFα can also induce the generation of free radicals in the form of Reactive Oxygen Species which can cause further liver damage and genomic instability [13]. It is also said that high TNF-α expression is an independent predictor of poor survival in HCC patients [14].

The results of various previous studies regarding the relationship between TNF-α polymorphism and HCC across various ethnicities and populations show mixed results, related [1517] or there is no relationship [18, 19]. Research on TNF-α gene SNPs in patients with Hepatitis B Virus (HBV) infection is specific for each ethnicity and shows different results in each population [20]. As the results regarding the role of these five SNPs of TNF-α genes against HCC are still controversial, we conducted this meta-analysis to determine the relationship between those TNF-α gene polymorphisms and HCC.

Methods

Database searching

We conducted an electronic database searching to identify all previously published cohort or case-control studies investigating the association between five types of TNF-α gene SNPs (-1031T/C, -863C/A, -857C/T, -308G/A, and − 238G/A) and risk of HCC. We searched data from PubMed, ProQuest, EBSCO, Science Direct, and Springer by using MeSH terms: “Tumor Necrosis Factor-alpha” or “TNF-alpha” and “polymorphism” or “SNP” or “single nucleotide polymorphism” or “variant” and “hepatocellular carcinoma” or “HCC” or “liver cancer”. The search was conducted in September 2019 - January 2020. We also performed a manual search was also performed to obtain potential sources cited in other meta-analysis.

Criteria For Inclusion And Exclusion

All included studies should meet the following criteria: (1) investigate the association between any of the five SNPs (-1031T/C, -863C/A, -857C/T, -308G/A, and − 238G/A) in TNF-α gene and risk of HCC; (2) have cohort or case-control study design; (3) risk of HCC was reported as RR (relative risk) and/or OR (odds ratio) with 95% Confidence Interval (CI) or provide sufficient data to extract RR and/or OR with 95% CI data; (4) include human subjects; (6) in English language. We excluded studies with the following criteria: (1) studies with design other than case-control or cohort; (2) duplicated studies; (3) studies with unmeasurable population and not qualified data; (4) using non-English language; (5) Studies in which the full text or main data could not be obtained.

Study Selection

Two investigators independently performed the electronic search and retrieved the articles that matched with our searched terms. Any disagreement was settled by discussion and consensus with all the authors. Final decision was merely based on the agreements of all authors.

Data Extraction

A standardized reporting form was used to extract the data from each article which included the first author’s name, year of publication, population country, TNF-α SNP type, study design, etiology of HCC, SNP genotyping method, controls, and frequencies of SNP. Hardy–Weinberg equilibrium test was performed and its significance of the control groups was calculated when the original information was not provided.

Quality Assessment

Newcastle - Ottawa quality assessment scale (NOQS) was used for measuring the quality of the included studies. This scale was designed through collaboration between the Universities of Newcastle, Australia and Ottawa, Canada. The purpose of this scale was to assess the quality of observational studies for producing good meta-analysis. The studies were qualified as high quality (9 stars), medium quality (7–8 stars), and low quality (less than 7 stars) [21].

Data analysis

Analysis were conducted using Review Manager 5.3 software (The Cochrane Collaboration, UK). Hardy-Weinberg Equilibrium was examined by Chi Square test when the original information was not provided. We evaluate 5 genetic models (allele, dominant, recessive, codominant major vs minor homozygote, and codominant heterozygote vs major homozygote) for each SNP type separately. Heterogeneity assumption was assessed with Cochrane Q statistic and I2 statistic. The pool estimated ORs was calculated with either fixed or random effects model assumptions. If Q test showed significant result (p-value < 0.05), we used a random effects model. Otherwise, if Q test showed insignificant result (p-value > 0.05), we used a fixed effect model. We also calculated the 95% confidence interval (CI) of pool estimated OR. Inverted funnel plots were conducted to find any presence of publication bias.

Results

Study Selection

According to the PRISMA flow diagram (Fig. 1), we initially obtained 29,182 studies through primary database searching and 12 through manual searching. After screening the titles/abstracts we selected 374 potentially relevant articles. Among them, 260 were excluded due to non-English language (19 studies) and unsuitable study designs (review papers/case reports/ cross sectional/meta-analysis/experimental; 241 studies). Then, 114 full text studies were checked for their eligibility. Some studies were excluded due to duplicates or irrelevant study design/insufficient information/unqualified articles until we finally obtained 23 included studies.

Eventually, as much as 23 potential articles were included which recruited 9,792 participants and consisted of 3,237 HCC cases and 4,843 controls (Table 1). The papers were published 2004 to 2018. Each study has sample size ranged from 45 to 1,624 participants. Of all included studies in this meta-analysis, most studies came from Asia, in which 7 studies were from China [7, 2227], 4 from Taiwan [15, 2830], 2 from India [18, 31], 2 from Korea [8, 32], 1 was Japanese [33], 1 was Thai [34], and 1 was Turkey [35]. There were also some studies with non-Asian countries, including 2 from Egypt [36, 37], 1 from Brazil [38], 1 from Italy [39], and 1 from Tunisia [40]. From the 23 included studies, the aetiology of HCC was mostly caused by HBV (12 studies), followed by mixed cause (8 studies), HCV (2 studies), and alcohol/smoking (1 study). Some studies observed 1 locus, and others observed more than 1 loci. Studies observing polymorphism at -308 was the most frequent (19 studies), while SNP − 1031 was the less frequently observed (5 studies).

Table 1

The basic characteristics of the included studies

No

First author

Year

TNF-α SNPs

Study design

Sample size

Population

SNP genotyping method

HCC etiology

Control

NOQS

Case (n = 3,237)

Control (n = 4,843)

1

Heneghan

2004

-308 G/A, -238 G/A

Case control

98

75

China

PCR-RFLP

Mixed

Healthy subjects

8

2

Ho

2004

-308 G/A

Case control

74

289

Taiwan

PCR-RFLP

Mixed

Healthy subjects

8

3

Migita

2005

-308 G/A

Case control

48

188

Japan

PCR-SSP

HBV

HBV without HCC

7

4

Niro

2005

-1031 T/C, -863 C/A, -308 G/A, -238 G/A

Case control

30

96

Italy

Sequencing

HBV

SR

7

5

Jeng

2007

-308 G/A

Case control

108

108

Taiwan

Sequencing

Mixed

Healthy subjects

9

6

Kummee

2007

-863 C/A, -308 G/A, -238 G/A

Case control

50

150

Thailand

PCR-RFLP

HBV

Healthy subjects

7

7

Akkiz

2009

-308 G/A

Case control

110

110

Turkey

PCR-RFLP

Mixed

Healthy subjects

9

8

Wang

2010

-308 G/A, -238 G/A

Case control

230

158

China

Sequencing

HBV

SR

8

9

Chen

2011

-1031 T/C, -863 C/A, -857 C/T, -238 G/A

Case control

126

126

China

Sequencing

HBV

Healthy subjects

9

10

Shi

2011

-308 G/A

Case control

88

88

China

PCR-RFLP

HBV

Healthy subjects

9

11

Qiu

2012

-863 C/A, -857 C/T

Case control

195

189

China

PCR-RFLP

HBV

SR

7

12

Radwan

2012

-308 G/A

Case control

128

160

Egypt

PCR-RFLP

HCV

Healthy subjects

8

13

Shi

2012

-308 G/A

Case control

73

116

China

PCR-RFLP

Mixed

Healthy subjects

8

14

Yang

2012

-863 C/A

Case control

772

852

China

RT-PCR

Mixed

Healthy subjects

7

15

Teixeira

2013

-308 G/A, -238 G/A

Case control

111

202

Brazil

PCR-SSP

Mixed

Healthy subjects

5

16

Panigrahi

2014

-863 C/A, -857 C/T, -238 G/A

Case control

14

85

India

PCR-RFLP

HBV

Healthy subjects

6

17

Saxena

2014

-308 G/A

Case control

59

139

India

PCR-RFLP

HBV

Healthy subjects

8

18

Jin

2015

-1031 T/C, -857 C/T, -238 G/A

Case control

224

206

Korea

Single base primer extension assay

HBV

LC

7

19

Sghaier

2015

-308 G/A, -238 G/A

Case control

15

200

Tunisia

PCR-RFLP

HBV

Healthy subjects

8

20

Shin

2015

-1031 T/C, -863 C/A, -857 C/T, -308 G/A, -238 G/A

Case control

157

201

Korea

PCR-RFLP

Mixed

Healthy subjects

8

21

Yang

2015

-1031 T/C, -863 C/A, -857 C/T, -308 G/A

Case control

298

889

Taiwan

PCR-RFLP

Smoking and alcohol

Healthy subjects

7

22

Tsai

2017

-308 G/A

Case control

200

200

Taiwan

Sequencing

HBV

LC

6

23

Tharwat

2018

-308 G/A

Retrospective cohort

29

16

Egypt

PCR-RFLP

HCV

Healthy subjecys

8

ASC = Asymptomatic Carrier; LC = Liver Cirrhosis; HCC = Hepatocellular Carcinoma; PCR-RFLP = polymerase chain reaction-restriction fragment length polymorphism; PCR-SSP = polymerase chain reaction-sequence-specific primers; RT-PCR = Real-time polymerase chain reaction; NOQS = Newcastle-Ottawa quality assessment scale.

Association Between Snp Tnf-α -1031 And Hcc Risk

Only five studies regarding the association between SNP TNF-α -1031 and HCC risk with 825 cases and 1518 controls were available. The number of cases for CC and TC genotypes was reported together in the studies by Niro et al [39] and Jin et al [32] which could only be used for dominant-model analysis (CC/TC vs TT; Table 2). As the heterogeneity among studies for all models (I2) was less than 60% and p-value for the heterogeneity was more than 0.05, fixed-effects models were applied. However, we did not obtain any significant association between SNP TNF-α -1031 and HCC risk in all model analyses. For the estimation of publication bias, we found no visual asymmetry in Funnel Plot analysis.

Table 2

Pooled risk estimates for SNPs TNF-α -1031 T/C, -863 C/A, -857 C/T, -308 G/A, and − 238 G/A with HCC risk

SNP TNF-α

N

OR (95% CI)

P value for Z test

I2 for heterogenity

P value for heterogenity

-1031 T/C

         

C vs T

3

1.06 [0.89–1.26]

0.52

0.77

0

CC + TC vs TT

5

1.04 [0.86–1.27]

0.69

0.97

0

CC vs TT

3

1.39 [0.88–2.18]

0.15

0.99

0

CC vs TC + TT

3

1.41 [0.91–2.21]

0.13

1

0

TC vs TT

3

0.96 [0.77–1.19]

0.70

0.64

0

-863 C/A

         

A vs C

7

1.31 [1.03–1.67]

0.03*

0.007

66

AA + CA vs CC

8

1.19 [1.03–1.36]

0.02*

0.35

10

AA vs CC

7

1.43 [0.98–2.10]

0.07

0.50

0

AA vs CA + CC

7

1.28 [0.89–1.86]

0.19

0.84

0

CA vs CC

7

1.41 [1.00-1.99]

0.05

0.0007

74

-857 C/T

         

T vs C

4

1.16[0.94–1.43]

0.17

0.45

0

TT + CT vs CC

5

1.31 [1.06–1.62]

0.01*

0.24

27

TT vs CC

4

0.75 [0.42–1.36]

0.35

0.87

0

TT vs CT + CC

4

0.77 [0.44–1.35]

0.36

0.77

0

CT vs CC

4

1.02 [0.77–1.34]

0.89

0.54

0

-308 G/A

         

A vs G

17

1.98 [1.62–2.42]

< 0.001*

0.03

44

AA + GA vs GG

19

1.95 [1.53–2.49]

< 0.001*

0.003

54

AA vs GG

13

3.14 [2.06–4.79]

< 0.001*

0.69

0

AA vs GA + GG

13

2.52 [1.69–3.76]

< 0.001*

0.81

0

GA vs GG

17

2.07 [1.60–2.68]

< 0.001*

0.006

52

-238 G/A

         

A vs G

8

1.50 [1.16–1.94]

0.002*

0.07

46

AA + GA vs GG

9

1.39 [0.87–2.24]

0.17

0.03

52

AA vs GG

5

3.87 [1.32–11.34]

0.01*

0.70

0

AA vs GA + GG

8

2.67 [1.17–6.10]

0.02*

0.87

0

GA vs GG

8

1.28 [0.72–2.28]

0.39

0.01

61

Association Between Snp Tnf-α -863 And Hcc Risk

We included eight studies with 1642 cases and 2746 controls to determine the association between SNP TNF-α -863 and HCC risk. The numbers of cases for CA and AA genotypes were reported together in the study by Niro et al [39], thus it could only be used for dominant-model analysis (AA/CA vs CC; Table 2). In studies with heterogeneity among studies (I2) more than 60% and p-value for the heterogeneity was less than 0.05, we applied random-effects models. We found significant association between the allele model of A versus C of TNF-α C/A SNP with the risk of HCC (OR = 1.31, 95% CI = 1.03–1.67, p = 0.03). The dominant-model analysis (CA/AA vs CC) also showed significant association between SNP TNF-α 863 C/A and HCC risk (OR = 1.19, 95% CI = 1.03–1.36, p = 0.02; Fig. 2). As heterogeneity was found in the statistical analyses, we did sensitivity analyses to evaluate the sources of heterogeneity. We found that heterogeneity between studies was mainly caused by the study by Kummee et al [34], as after this study was excluded, no significant heterogeneity was found.

Association Between Snp Tnf-α -857 And Hcc Risk

We only obtained five studies about the association between SNP TNF-α -857 and HCC risk with 716 cases and 1005 controls. The number of cases for TT and CT genotypes was reported together in the studies by Jin et al [32] which could only be used for dominant-model analysis (TT/CT vs CC; Table 2). As the heterogeneity among studies for all models (I2) was less than 60% and p-value for the heterogeneity was more than 0.05, fixed-effects models were applied. A significant association between SNP TNF-α -857 C/T and HCC risk was found in dominant-model analyses (OR = 1.31, 95% CI = 1.06–1.62, p = 0.01; Fig. 3). We did not observe any visual asymmetry in Funnel Plot analysis regarding the publication bias.

Association Between Snp Tnf-α -308 And Hcc Risk

Interestingly, we found most included studies analyzing the association between SNP TNF-α -308 G/A and HCC risk. The number of cases for AA and GA genotypes was reported together in the studies by Niro et al [39] and Jin et al [32] which could only be used for dominant-model analysis (AA/GA vs GG; Table 2). In studies with heterogeneity among studies (I2) more than 60% and p-value for the heterogeneity was less than 0.05, random-effects models were applied. All five genetic models showed significant association with risk of HCC with OR for allele model = 1.98, 95% CI = 1.62–2.42, p < 0.001; OR for dominant model = 1.95, 95% CI = 1.53–2.49, p < 0.001; OR for recessive model = 2.52, 95% CI = 1.69–3.76, p < 0.001; OR for codominant major vs minor homozygote model = 3.14, 95% CI = 2.06–4.79, p < 0.001; and OR for codominant heterozygote vs major homozygote model = 2.07, 95% CI = 1.60–2.68, p < 0.001; Fig. 4). As heterogeneity was found in the statistical analyses, we did sensitivity analyses to evaluate the sources of heterogeneity. We found that heterogeneity between studies was mainly caused by the studies by Ho et al [28] and Akkiz et al [35], as after these studies were excluded, no significant heterogeneity was found.

Association Between Snp Tnf-α -238 And Hcc Risk

Nine studies with 831 cases and 1293 controls were included to determine the association between SNP TNF-α -238 and HCC risk. Similar to the TNF-a -863 SNP, the numbers of cases for GA and AA genotypes were reported together in the study by Niro et al [39], thus it could only be used for dominant-model analysis (GA/AA vs GG; Table 2). In studies with heterogeneity among studies (I2) more than 60% and p-value for the heterogeneity was less than 0.05, we applied random-effects models. The allele model of A versus G of SNP TNF-α -238 G/A was significantly associated with risk of HCC (OR = 1.50, 95% CI = 1.16–1.94, p = 0.002). The codominant-model analysis (AA vs GG) also showed significant association between SNP TNF-α 238 G/A and HCC risk (OR = 3.87, 95% CI = 1.32–11.34, p = 0.01). The recessive model analysis (AA vs GA + GG) proved significant association between SNP TNF-α 238 G/A and HCC risk as well (OR = 2.67, 95% CI = 1.17–6.10, p = 0.02; Fig. 5). As heterogeneity was found in the statistical analyses, we did sensitivity analyses to evaluate the sources of heterogeneity. We found that heterogeneity between studies was mainly caused by the study by Teixeira et al [38], as after this study was excluded, no significant heterogeneity was found.

Sensitivity Analysis

In this meta-analysis, there were no significant change in ORs by deleting a particular study, which indicated that no single study influenced the statistical significance of the overall results.

Discussion

The development of HCC depends on some factors such as viral infection, environmental, behavioral, metabolism, and genetics [41, 42]. The contribution of SNP as a genetic factor is widely studied related to its role in the development of HCC, in which the prevalence of SNP is different in each population [5, 43, 44]. From the host factor, TNF-α is thought to play an important role in hepatocarcinogenesis through the induction of fibrogenic factors. Tumor Necrosis Factor-α is associated with increased necroinflammatory activity and liver fibrosis. Necroinflammation of hepatocytes triggers mutagenesis and activation of oncogens from protooncogenes in host cells, causing HCC [45].

Five biallelic SNPs in the promoter region Tumor Necrosis Factor-α gene were known: -1031T/C, -863C/A, -857C/T, -308G/A, and − 238G/A [7]. There have been many studies conducted to determine the association between these five TNF-α SNPs and HCC risk. Several meta-analyses have also been performed to analyze the association between TNF-α SNP and HCC. However, to date, only few meta-analyses have analysed all these five SNPs with HCC risk. Hereby, we performed this updated meta-analysis to clarify the independent role of each of these five SNPs on HCC risk.

As research on SNP TNF-α -1031 T/C and HCC risk is still limited, we only found five studies investigating this SNP, even though many studies have shown significant association between this SNP with other diseases such as polycystic ovary syndrome and endometriosis [46, 47]. In a study conducted by Shin et al survival of HCC cases with TNF-α -1031 wild type (TT) genotype or SNP TC genotype was significantly better than those with the SNP CC genotype [8]. However, in this meta-analysis, no significant relationship was found between the SNP and HCC risk in all genetic models. Meta-analysis conducted by Wei et al also shows that there is no significant relationship between this SNP and HCC risk [48].

This study shows a significant relationship between SNP TNF-α -863 C/A with HCC risk in allele models and dominant model analysis. A meta-analysis conducted by Wei et al also showed a significant relationship betweehis SNP and HCC in dominant and codominant (CA vs CC) model analysis [48]. Polymorphism of TNF-α -863 C/A in the promoter region can influence TNF-α expression, however, the result is still conflicting. Some research suggest that it may increase TNF-α expression [49, 50], while other research show the opposite result in which it may decrease TNF-α expression [27, 51].

In the present study, we also found a significant relationship between SNP TNF-α-857 C/T with dominant model analysis. Limited studies were present regarding the relationship between this polymorphism with HCC risk, thus there were only 7 studies that came within our inclusion criteria. This is different from the meta-analysis conducted by Wei et al which showed no significant relationship between this SNP and HCC risk, however in that meta-analysis there were only 3 included studies [48].

There have been many previous studies investigating SNP TNF-α -308 G/A with HCC risk, thus we got 19 included studies. All five genetic analysis models showed a significant relationship between the SNP and HCC risk. This is in line with several previous meta-analysis studies, including those conducted by Hu et al and Tavakoulpour and Sali on allele models and dominant model analyses [16, 52], Wei et al on codominant dan dominant model analyses [48], and Xiao et al in all except recessive model analysis [17].

Studies on SNP TNF-α -238 G/A have also been extensively performed. Nevertheless, meta-analyses of these SNPs are still limited and yield conflicting results. The present study showed a significant relationship between this SNP and HCC risk. This is in line with a meta-analysis conducted by Xiao et al showing that this SNP is associated with HCC risk, although in that study only HBV-related HCC was studied [17]. Another meta-analysis conducted by Hu et al, however, shows no significant relationship between this SNP with HCC risk in Asian population [52].

This meta-analysis still had several limitations. To date, there has rarely been any meta-analysis discussing the relationship between SNPs of TNF-α -1031 T/C, -863 C/A, -857 C/T, and − 238 G/A with HCC. Most present meta-analysis only focus on SNP − 308 G/A. We also got a small number of studies on SNPs − 1031 T/C, -863 C/A, -857 C/T, due to the limited number of studies. Research on these three SNPs, both observational and meta-analysis, is further needed in the future to find out the role of each of the three SNPs on HCC so that a stronger power of study will be obtained.

The second limitation of this study was that we only included studies in English language so that it does not rule out the exclusion of any good research in non-English languages. Third, HCC is a multifactorial process involving host, agent, and environment factors. The role of other host genomes, viruses, and the environment can influence the course and risk of HCC. Variations in genes adjacent to the TNF-α gene are said to also be able to regulate TNF-α expression [53]. Epigenetic factors which include histone modification, non-coding RNA, and gene methylation can also affect TNF-α expression and contribute to HCC risk [54]. Fourth, the study populations in our included studies came from various ethnicities (Asian, African, European and South American). Further research is required to determine the effects of these various TNF-α SNPs in each different population.

Conclusion

In conclusion, this meta-analysis determined the association between SNP TNF-α -1031 T/C, -863 C/A, -857 C/T, -308 G/A, and − 238 G/A and HCC risk. Our novel data demonstrated that all these SNPs but − 1031 T/C might increase HCC risk. Therefore, genetic predisposition, especially polymorphism of TNF-α gene may play important role in the pathogenesis of HCC. Further studies in larger population and analysing gene-environment interaction are required to provide better understanding between TNF-α polymorphisms and the risk of HCC.

Abbreviations

HCC = Hepatocellular carcinoma; TNF-α = Tumor Necrosis Factor-α; SNP = Single Nucleotide Polymorphism; HBV = Hepatitis B Virus; RR = relative risk; OR = odds ratio; CI = confidence interval; NOQS = Newcastle - Ottawa quality assessment scale.

Declarations

Ethics approval and consent to participate

Not applicable 

Consent for publish

Not applicable 

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request. 

Competing interests

The authors declare that they have no conflict of interest. 

Funding

All publication funds from the beginning of data collection to manuscript preparation and our publication costs are supported by Dato ’Sri Prof. Dr. Tahir through the Tahir Professorship Program, Indonesia. 

Authors’ Contributions

This study was designed by CDKW and RH, CDKW and FCA conducted the literature search, S performed statistical analysis, GIP and RH wrote the manuscript. All authors read and approved the final manuscript. 

Acknowledgments

Not applicable.

Authors’ Information

Citrawati Dyah Kencono Wungu: https://orcid.org/0000-0001-5180-957X

Fis Citra Ariyanto: https://orcid.org/0000-0001-6883-195X

Gwenny Ichsan Prabowo: -

Soetjipto: -

Retno Handajani: -

References

  1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A (2018) Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA: a cancer journal for clinicians 68 (6):394-424. doi:10.3322/caac.21492
  2. Fact Sheets: Liver Cancer (2016) http://gco.iarc.fr/today/data/pdf/fact-sheets/cancers/cancer-fact-sheets-7.pdf.
  3. Brar TS, Hilgenfeldt E, Soldevila-Pico C (2018) Etiology and Pathogenesis of Hepatocellular Carcinoma. In: Liu C (ed) Precision Molecular Pathology of Liver Cancer. Springer International Publishing, Cham, pp 1-15. doi:10.1007/978-3-319-68082-8_1
  4. Nahon P, Zucman-Rossi J (2012) Single nucleotide polymorphisms and risk of hepatocellular carcinoma in cirrhosis. Journal of Hepatology 57 (3):663-674. doi:https://doi.org/10.1016/j.jhep.2012.02.035
  5. Mathew S, Abdel-Hafiz H, Raza A, Fatima K, Qadri I (2016) Host nucleotide polymorphism in hepatitis B virus-associated hepatocellular carcinoma. World journal of hepatology 8 (10):485-498. doi:10.4254/wjh.v8.i10.485
  6. Ma L, Chen S, Mao X, Lu Y, Zhang X, Lao X, Qin X, Li S (2018) The association between TNFR gene polymorphisms and the risk of Hepatitis B Virus-Related Liver Diseases in Chinese population. Scientific Reports 8 (1):9240. doi:10.1038/s41598-018-27623-7
  7. Chen X, Zhang L, Chang Y, Shen T, Wang L, Zhuang H, Lu F (2011) Association of TNF-alpha genetic polymorphisms with hepatocellular carcinoma susceptibility: a case-control study in a Han Chinese population. The International journal of biological markers 26 (3):181-187. doi:10.5301/jbm.2011.8580
  8. Shin SP, Kim NK, Kim JH, Lee JH, Kim JO, Cho SH, Park H, Kim MN, Rim KS, Hwang SG (2015) Association between hepatocellular carcinoma and tumor necrosis factor alpha polymorphisms in South Korea. World J Gastroenterol 21 (46):13064-13072. doi:10.3748/wjg.v21.i46.13064
  9. El-Tahan RR, Ghoneim AM, El-Mashad N (2016) TNF-α gene polymorphisms and expression. SpringerPlus 5 (1):1508. doi:10.1186/s40064-016-3197-y
  10. Cheong JY, Cho SW, Hwang IL, Yoon SK, Lee JH, Park CS, Lee JE, Hahm KB, Kim JH (2006) Association between chronic hepatitis B virus infection and interleukin-10, tumor necrosis factor-α gene promoter polymorphisms. Journal of Gastroenterology and Hepatology 21 (7):1163-1169. doi:10.1111/j.1440-1746.2006.04304.x
  11. Martins GR, Gelaleti GB, Moschetta MG, Maschio-Signorini LB, Zuccari DA (2016) Proinflammatory and Anti-Inflammatory Cytokines Mediated by NF-κB Factor as Prognostic Markers in Mammary Tumors. Mediators of inflammation 2016:9512743. doi:10.1155/2016/9512743
  12. Méndez-García LA, Nava-Castro KE, Ochoa-Mercado TL, Palacios-Arreola MI, Ruiz-Manzano RA, Segovia-Mendoza M, Solleiro-Villavicencio H, Cázarez-Martínez C, Morales-Montor J (2019) Breast Cancer Metastasis: Are Cytokines Important Players During Its Development and Progression? Journal of interferon & cytokine research : the official journal of the International Society for Interferon and Cytokine Research 39 (1):39-55. doi:10.1089/jir.2018.0024
  13. Wang Z, Li Z, Ye Y, Xie L, Li W (2016) Oxidative Stress and Liver Cancer: Etiology and Therapeutic Targets. Oxidative medicine and cellular longevity 2016:7891574. doi:10.1155/2016/7891574
  14. Wang H, Liu J, Hu X, Liu S, He B (2016) Prognostic and Therapeutic Values of Tumor Necrosis Factor-Alpha in Hepatocellular Carcinoma. Medical science monitor : international medical journal of experimental and clinical research 22:3694-3704. doi:10.12659/msm.899773
  15. Jeng J-E, Tsai J-F, Chuang L-Y, Ho M-S, Lin Z-Y, Hsieh M-Y, Chen S-C, Chuang W-L, Wang L-Y, Yu M-L, Dai C-Y, Chang J-G (2007) Tumor necrosis factor-alpha 308.2 polymorphism is associated with advanced hepatic fibrosis and higher risk for hepatocellular carcinoma. Neoplasia 9 (11):987-992. doi:10.1593/neo.07781
  16. Cheng K, Zhao YJ, Liu L, Wan JJ (2013) Tumor necrosis factor-α 238 G/A polymorphism and risk of hepatocellular carcinoma: evidence from a meta-analysis. Asian Pacific journal of cancer prevention : APJCP 14 (5):3275-3279. doi:10.7314/apjcp.2013.14.5.3275
  17. Xiao Q, Fu B, Chen P, Liu ZZ, Wang W, Ye Q (2016) Three polymorphisms of tumor necrosis factor-alpha and hepatitis B virus related hepatocellular carcinoma: A meta-analysis. Medicine (Baltimore) 95 (50):e5609-e5609. doi:10.1097/MD.0000000000005609
  18. Saxena R, Chawla YK, Verma I, Kaur J (2014) IFN-γ (+874) and not TNF-α (−308) is associated with HBV-HCC risk in India. Molecular and Cellular Biochemistry 385 (1):297-307. doi:10.1007/s11010-013-1838-9
  19. Somi MH, Najafi L, Noori BN, Alizadeh AH, Aghah MR, Shavakhi A, Ehsani MJ, Aghazadeh R, Masoodi M, Amini S, Baladast M, Zali MR (2006) Tumor necrosis factor-alpha gene promoter polymorphism in Iranian patients with chronic hepatitis B. Indian journal of gastroenterology : official journal of the Indian Society of Gastroenterology 25 (1):14-15
  20. Heidari Z, Moudi B, Mahmoudzadeh Sagheb H, Moudi M (2016) Association of TNF-α Gene Polymorphisms with Production of Protein and Susceptibility to Chronic Hepatitis B Infection in the South East Iranian Population. Hepat Mon 16 (11):e41984-e41984. doi:10.5812/hepatmon.41984
  21. Wells GA, Shea B, O’connell D, Peterson J, Welch V, Losos M, Tugwell P (2016) The Newcastle-Ottawa Scale (NOS) for assessing the quality of nonrandomised studies in meta-analyses. 2013.
  22. Heneghan MA, Johnson PJ, Clare M, Ho S, Harrison PM, Donaldson PT (2003) Frequency and nature of cytokine gene polymorphisms in hepatocellular carcinoma in Hong Kong Chinese. International Journal of Gastrointestinal Cancer 34 (1):19-26. doi:10.1385/IJGC:34:1:19
  23. Wang B, Wang J, Zheng Y, Zhou S, Zheng J, Wang F, Ma X, Zeng Z (2010) A study of TNF-alpha-238 and -308 polymorphisms with different outcomes of persistent hepatitis B virus infection in China. Pathology 42 (7):674-680. doi:10.3109/00313025.2010.523696
  24. Shi HZ, Ren P, Lu QJ, Niedrgethmnn M, Wu GY (2012) Association between EGF, TGF-β1 and TNF-α gene polymorphisms and hepatocellular carcinoma. Asian Pacific journal of cancer prevention : APJCP 13 (12):6217-6220. doi:10.7314/apjcp.2012.13.12.6217
  25. Qiu B, Wang X, Zhang P, Shi C, Zhang J, Qiu W, Wang W, Li D (2012) Association of TNF-α promoter polymorphisms with the outcome of persistent HBV infection in a northeast Chinese Han population. Acta biochimica et biophysica Sinica 44 (8):712-718. doi:10.1093/abbs/gms046
  26. Shi Z, Du C (2011) Tumor necrosis factor alpha 308 G/A polymorphism and hepatocellular carcinoma risk in a Chinese population. Genetic testing and molecular biomarkers 15 (7-8):569-572. doi:10.1089/gtmb.2011.0008
  27. Yang Y, Qiu X-Q, Yu H-P, Zeng X-Y, Bei C-H (2012) TNF-α -863 polymorphisms and the risk of hepatocellular carcinoma. Exp Ther Med 3 (3):513-518. doi:10.3892/etm.2011.418
  28. Ho S-Y, Wang Y-J, Chen H-L, Chen C-H, Chang C-J, Wang P-J, Chen HHW, Guo H-R (2004) Increased Risk of Developing Hepatocellular Carcinoma Associated with Carriage of the TNF2 Allele of the −308 Tumor Necrosis Factor-α Promoter Gene. Cancer Causes & Control 15 (7):657-663. doi:10.1023/B:CACO.0000036173.99930.75
  29. Yang MD, Hsu CM, Chang WS, Yueh TC, Lai YL, Chuang CL, Wang SC, Jeng LB, Ji HX, Hsiao CL, Wu CN, Tsai CW, Chung JG, Bau DT (2015) Tumor Necrosis Factor-α Genotypes Are Associated with Hepatocellular Carcinoma Risk in Taiwanese Males, Smokers and Alcohol Drinkers. Anticancer research 35 (10):5417-5423
  30. Tsai J-F, Chen S-C, Lin Z-Y, Dai C-Y, Huang J-F, Yu M-L, Chuang W-L (2017) Independent and additive interaction between polymorphisms of tumor necrosis factor α−308 and lymphotoxin α+252 on risk of hepatocellular carcinoma related to hepatitis B. The Kaohsiung Journal of Medical Sciences 33 (9):453-457. doi:https://doi.org/10.1016/j.kjms.2017.04.009
  31. Panigrahi R, Sarkar N, Biswas A, Pal A, Saha D, Singh SP, Panigrahi MK, Bandopadhyay M, Chakrabarti S, Chakravarty R (2014) Association of TNF-α Promoter Polymorphism with HBV Associated Disease Outcome Among HBV Infected Patients from Orissa, Southern Part of East India. J Clin Exp Hepatol 4 (3):202-208. doi:10.1016/j.jceh.2014.08.002
  32. Jin YJ, Lee D, Chung YH, Kim JA, Kim SE, Lee YS, Shin ES, Ryu SH, Jang MK, Lee JE, Park NH (2015) Tumor Necrosis Factor-Alpha Gene Polymorphism Associated With Development of Hepatitis B Virus-associated Hepatocellular Carcinoma. Journal of clinical gastroenterology 49 (8):e76-81. doi:10.1097/mcg.0000000000000261
  33. Migita K, Miyazoe S, Maeda Y, Daikoku M, Abiru S, Ueki T, Yano K, Nagaoka S, Matsumoto T, Nakao K, Hamasaki K, Yatsuhashi H, Ishibashi H, Eguchi K (2005) Cytokine gene polymorphisms in Japanese patients with hepatitis B virus infection--association between TGF-beta1 polymorphisms and hepatocellular carcinoma. J Hepatol 42 (4):505-510. doi:10.1016/j.jhep.2004.11.026
  34. Kummee P, Tangkijvanich P, Poovorawan Y, Hirankarn N (2007) Association of HLA-DRB1*13 and TNF-alpha gene polymorphisms with clearance of chronic hepatitis B infection and risk of hepatocellular carcinoma in Thai population. Journal of viral hepatitis 14 (12):841-848. doi:10.1111/j.1365-2893.2007.00880.x
  35. Akkiz H, Bayram S, Bekar A, Ozdil B, Akgöllü E, Sümbül AT, Demiryürek H, Doran F (2009) G-308A TNF-alpha polymorphism is associated with an increased risk of hepatocellular carcinoma in the Turkish population: case-control study. Cancer epidemiology 33 (3-4):261-264. doi:10.1016/j.canep.2009.06.001
  36. Radwan MI, Pasha HF, Mohamed RH, Hussien HI, El-Khshab MN (2012) Influence of transforming growth factor-β1 and tumor necrosis factor-α genes polymorphisms on the development of cirrhosis and hepatocellular carcinoma in chronic hepatitis C patients. Cytokine 60 (1):271-276. doi:10.1016/j.cyto.2012.05.010
  37. Tharwat E, Gad GFM, Nazmy MH, Mohamed HI, Hamza N, Wahid A, Ibrahim ARN (2019) Impact of IL-27p28 (rs153109) and TNF-α (rs1800629) Genetic Polymorphisms on the Progression of HCV Infection in Egyptian Patients. Immunological investigations 48 (3):255-267. doi:10.1080/08820139.2018.1510958
  38. Teixeira AC, Mendes CT, Jr., Marano LA, Deghaide NH, Secaf M, Elias J, Jr., Muglia V, Donadi EA, Martinelli AL (2013) Alleles and genotypes of polymorphisms of IL-18, TNF-α and IFN-γ are associated with a higher risk and severity of hepatocellular carcinoma (HCC) in Brazil. Human immunology 74 (8):1024-1029. doi:10.1016/j.humimm.2013.04.029
  39. Niro GA, Fontana R, Gioffreda D, Valvano MR, Lacobellis A, Facciorusso D, Andriulli A (2005) Tumor necrosis factor gene polymorphisms and clearance or progression of hepatitis B virus infection. Liver international : official journal of the International Association for the Study of the Liver 25 (6):1175-1181. doi:10.1111/j.1478-3231.2005.01166.x
  40. Sghaier I, Zidi S, Mouelhi L, Dabbech R, Ghazouani E, Brochot E, Stayoussef M, Yacoubi-Loueslati B (2015) The relationship between TNF alpha gene polymorphisms (−238/−308), TNF RII VNTR (p75) and outcomes of hepatitis B virus infection in Tunisian population. Gene 568 (2):140-145. doi:https://doi.org/10.1016/j.gene.2015.05.029
  41. Zamor PJ, deLemos AS, Russo MW (2017) Viral hepatitis and hepatocellular carcinoma: etiology and management. J Gastrointest Oncol 8 (2):229-242. doi:10.21037/jgo.2017.03.14
  42. An P, Xu J, Yu Y, Winkler CA (2018) Host and Viral Genetic Variation in HBV-Related Hepatocellular Carcinoma. Front Genet 9:261-261. doi:10.3389/fgene.2018.00261
  43. Pontoriero AC, Trinks J, Hulaniuk ML, Caputo M, Fortuny L, Pratx LB, Frías A, Torres O, Nuñez F, Gadano A, Argibay P, Corach D, Flichman D (2015) Influence of ethnicity on the distribution of genetic polymorphisms associated with risk of chronic liver disease in South American populations. BMC Genetics 16 (1):93. doi:10.1186/s12863-015-0255-3
  44. Thursz M, Yee L, Khakoo S (2011) Understanding the host genetics of chronic hepatitis B and C. Seminars in liver disease 31 (2):115-127. doi:10.1055/s-0031-1276642
  45. Aroucha DCBL, do Carmo RF, Moura P, Silva JLA, Vasconcelos LRS, Cavalcanti MSM, Muniz MTC, Aroucha ML, Siqueira ERF, Cahú GGOM, Pereira LMMB, Coêlho MRCD (2013) High tumor necrosis factor-α/interleukin-10 ratio is associated with hepatocellular carcinoma in patients with chronic hepatitis C. Cytokine 62 (3):421-425. doi:https://doi.org/10.1016/j.cyto.2013.03.024
  46. Abutorabi R, Baradaran A, Sadat Mostafavi F, Zarrin Y, Mardanian F (2015) Evaluation of Tumor Necrosis Factor Alpha Polymorphism Frequencies in Endometriosis. Int J Fertil Steril 9 (3):329-337. doi:10.22074/ijfs.2015.4548
  47. Yun J-H, Choi J-W, Lee K-J, Shin J-S, Baek K-H (2011) The promoter -1031(T/C) polymorphism in tumor necrosis factor-alpha associated with polycystic ovary syndrome. Reproductive Biology and Endocrinology 9 (1):131. doi:10.1186/1477-7827-9-131
  48. Wei Y, Liu F, Li B, Chen X, Ma Y, Yan L, Wen T, Xu M, Wang W, Yang J (2011) Polymorphisms of Tumor Necrosis Factor-Alpha and Hepatocellular Carcinoma Risk: A HuGE Systematic Review and Meta-Analysis. Digestive Diseases and Sciences 56 (8):2227-2236. doi:10.1007/s10620-011-1617-y
  49. Higuchi T, Seki N, Kamizono S, Yamada A, Kimura A, Kato H, Itoh K (1998) Polymorphism of the 5′-flanking region of the human tumor necrosis factor (TNF)-α gene in Japanese. Tissue Antigens 51 (6):605-612. doi:10.1111/j.1399-0039.1998.tb03002.x
  50. Soga Y, Nishimura F, Ohyama H, Maeda H, Takashiba S, Murayama Y (2003) Tumor necrosis factor-alpha gene (TNF-α) −1031/−863, −857 single-nucleotide polymorphisms (SNPs) are associated with severe adult periodontitis in Japanese. Journal of Clinical Periodontology 30 (6):524-531. doi:10.1034/j.1600-051X.2003.00287.x
  51. Laddha NC, Dwivedi M, Begum R (2012) Increased Tumor Necrosis Factor (TNF)-α and its promoter polymorphisms correlate with disease progression and higher susceptibility towards vitiligo. PloS one 7 (12):e52298. doi:10.1371/journal.pone.0052298
  52. Hu Q, Lou G-G, Liu Y-C, Qian L, Lv B-D (2014) The Tumor Necrosis Factor-α-308 and -238 Polymorphisms and Risk of Hepatocellular Carcinoma for Asian Populations: A Meta-Analysis. Current Therapeutic Research 76:70-75. doi:https://doi.org/10.1016/j.curtheres.2014.04.001
  53. Vázquez-Huerta DI, Alvarez-Rodríguez BA, Topete-Reyes JF, Muñoz-Valle JF, Parra-Michel R, Fuentes-Ramírez F, Salazar-López MA, Valle Y, Reyes-Castillo Z, Cruz-González A, Brennan-Bourdon LM, Torres-Carrillo N (2014) Tumor necrosis factor alpha -238 G/A and -308 G/A polymorphisms and soluble TNF-α levels in chronic kidney disease: correlation with clinical variables. Int J Clin Exp Med 7 (8):2111-2119
  54. Wang Q, Roy B, Turecki G, Shelton RC, Dwivedi Y (2018) Role of Complex Epigenetic Switching in Tumor Necrosis Factor-α Upregulation in the Prefrontal Cortex of Suicide Subjects. Am J Psychiatry 175 (3):262-274. doi:10.1176/appi.ajp.2017.16070759