Subsets of Inflammatory Cytokine Gene Polymorphisms are Associated with Risk of Carcinogenic Liver Fluke Opisthorchis viverrini-Associated Advanced Periductal Fibrosis and Cholangiocarcinoma
Article information
Abstract
Opisthorchis viverrini infection induces chronic inflammation, and a minor proportion of infected individuals develop advanced periductal fibrosis (APF) and cholangiocarcinoma (CCA). Inflammatory cytokines and/or their gene polymorphisms may link to these biliary pathologies. We therefore investigated associations among cytokine gene polymorphisms and cytokine production in 510 Thai cases infected with O. viverrini who presented with APF+ or APF−, as established by abdominal ultrasonography as well as in patients diagnosed with CCA. Levels of pro-inflammatory and anti-inflammatory cytokines were determined in culture supernatants after stimulation of peripheral blood mononuclear cells (PBMCs) with O. viverrini excretory-secretory (ES) products. Pro-inflammatory cytokines, IL-1β, IL-6, IFN-γ, LT-α, and TNF-α were significantly increased in CCA patients compared with non-CCA (APF− and APF+) cases. Polymorphisms in genes encoding IL-1β-511C/T, IL-6-174G/C, IFN-γ +874T/A, LT-α +252A/G, and TNF-α −308G/A were then investigated by using PCR-RFLP or allele specific-PCR (AS-PCR) analyses. In the CCA cases, LT-α +252A/G and TNF-α −308G/A heterozygous and homozygous variants showed significantly higher levels of these cytokines than the wild type. By contrast, levels of cytokines in wild type of IFN-γ +874T/A were significantly higher than the variants in CCA cases. IFN-γ +874T/A polymorphisms were associated with advanced periductal fibrosis, whereas IL-6 −174G/C polymorphisms were associated with CCA. To our knowledge, these findings provide the first demonstration that O. viverrini infected individuals carrying several specific cytokine gene polymorphisms are susceptible to develop fibrosis and CCA.
INTRODUCTION
Chronic infection with the carcinogenic liver fluke Opisthorchis viverrini results in continuous irritation, chronic inflammation, restlessness tissue damage, and repair of the bile tract, including the second order bile ducts [1]. Consequently, infection with O. viverrini is strongly associated with chronic hepatobiliary disease, including cholangitis, cholecystitis, obstructive jaundice, advanced periductal fibrosis, and cholangiocarcinoma (CCA) [1,2]. Ultrasonographic assessment of hepatobiliary abnormalities of otherwise asymptomatic persons stool positive for infection with O. viverrini revealed that about 25% of these cases developed advanced periductal fibrosis, and CCA was suspected in up to 0.5% of them based on abdominal ultrasound examinations [3,4]. The mechanism underlying the different outcomes of O. viverrini infection is unknown. Our extensive research experience with opisthorchiasis led us the suspect that host inflammatory responses to the liver fluke might play pivotal roles in the outcome [1,2]. Inflammation mediated by cytokines has been described in infectious diseases and malignancies [5,6]. Specifically, the liver fluke infection induces marked inflammation characterized by infiltration by neutrophils and mononuclear cells into the vicinity of the infected bile ducts. Proteins and other metabolites released by the O. viverrini worms activate Toll-like receptor 4 and induce expression of the potent chemo-attractant IL-8 and pro-inflammatory IL-6 cytokines [7].
Gene polymorphisms affect the cytokine expression and influence inflammatory processes in response to infectious diseases and cancer [8]. In the context of liver fluke infection, hepatobiliary disease and associated bile duct cancer, we have reported marked elevation in expression of IL-6 in response to O. viverrini excretory-secretory (ES) antigens during chronic opisthorchiasis, particularly where biliary periductal fibrosis is marked [9]. In addition, we observed significant increases in both pro-inflammatory and anti-inflammatory cytokines during cholangiocarcinoma associated with chronic liver fluke infection [10]. We now hypothesize that cytokine gene polymorphisms influence cytokine expression profiles, which in turn may result in diverse pathogenic sequelae during infection with O. viverrini. This study described below investigated associations among cytokine gene polymorphisms and cytokine production during chronic opisthorchiasis and associated biliary tract diseases, including CCA.
MATERIALS AND METHODS
Study participants
This study was approved by the Khon Kaen University Human Ethics Committee, reference nos. HE480528 and HE 581272. Written informed consents were obtained from all participants or patients. Four hundred participants positive for infection with O. viverrini, including 200 without APF− and 200 with APF+, were recruited from sites endemic for O. viverrini infection in Khon Kaen or neighboring provinces in northeastern Thailand, as described [9]. This study also included 110 histologically proven cases of cholangiocarcinoma (CCA); these latter participants had undergone hepatectomy at Srinagarind Hospital, Faculty of Medicine, Khon Kaen University, Khon Kaen, Thailand.
Opisthorchis viverrini-ES product preparation
ES products of O. viverrini (Ov-ES) were prepared as described [11,12]. Briefly, viable liver flukes recovered at necropsy from experimentally infected hamsters were cultured in RPMI-1640 containing antibiotics and the protease inhibitors (0.1 mM phenylmethanesulfonyl fluoride (PMSF), 1 mM leupeptin and 0.1 mM N-[N-(L-3-trans-carboxyoxiran-2-carbonyl)-L-leucine]-agmatine, E-64). The cultures were maintained in vitro at 37°C and supernatants containing the Ov-ES were collected twice daily for 7 days. The culture supernatants were subjected to centrifugation at 3,000 g for 10 min to pellet the eggs. Clarified supernatants were pooled, dialyzed against PBS, and concentrated by membrane ultrafiltration (PM 10, Amicon, Danvers, Masssacusetts, USA). Lipopolysaccharide (LPS) was removed using Triton-X114 extraction [13], followed by passage through Bio-Beads SM2 (Bio-Rad, Hercules, CA, USA) to remove Triton-X114. Finally, OvES was filtered through a 0.2 μm membrane and then aliquoted and stored at −80°C. The LPS concentration was determined using the Limulus amebocyte assay (less than 100 ng/ml failed to affect cell proliferation or cytokine production, unpublished) [14].
Peripheral blood mononuclear cells (PBMCs) culture
PBMC were isolated from heparinized blood by density gradient centrifugation using Ficoll-Paque (Sigma Diagnostic, St. Louis, Missouri, USA), counted on a hematocytometer and viability assessed using 0.5% Trypan blue dye. PBMC were cultured at 2×104 cells/ml in RPMI-1640 supplemented with antibiotics and 5% fetal bovine serum with or without Ov-ES, 20 μg/ml. Phytohemaglutinin (PHA)-L (Difco Laboratories, Detroit, Michigan, USA) at 10 μg/ml was used as positive control of cytokine stimulation. The cells were cultured at 37°C in a humidified 5% CO2 incubator for 48 hr. Supernatants were harvested and stored at −80°C for cytokine level determination.
Cytokine detection
A 25–50 μl of supernatant was used for cytokine assays. The secreted cytokines were quantified using the FlowCytomix bead-based multiplexing assays kit, human Th1/Th2 (IFN-γ, IL-1β, IL-2, IL-4, IL-5, IL-6, IL-8, IL-10, IL-12p70, TNF-α, and LT-α) 11-plex kit (BMS810FF, BenderMedSystems, Vienna, Austria) and quantified using flow cytometry (FC500, Beckman-Coulter, Brea, California, USA). Cytokine levels were established by interpolation from standard curves established using human recombinant cytokines (BMS8400FF, BenderMedSystems), after which the findings were expressed as pg/ml.
DNA extraction and genotyping
Genomic DNA was extracted from PBMC from each participant using the Puregene DNA Isolation Kit (Gentra systems, Minneapolis, Minnesota, USA), quantified using the NanoDrop ND-1000 and stored at −20°C. Single nucleotide polymorphisms (SNPs) for IL-1β-511C/T, IL-6 −174G/C, LT-α +252A/G, and TNF-α −308G/A polymorphisms were analyzed by restriction fragment length polymorphism (RFLP), as described [15–18]. Allele specific-PCR (AS-PCR) was used to analyze for IFN-γ +874T/A polymorphism [19]. Table 1 lists the PCR thermal cycling conditions and oligonucleotide primers specific for 5 cytokine gene polymorphisms, IL-1β −511C/T, IL-6 −174G/C, IFN-γ +874T/A, LT-α +252A/G, and TNF-α −308G/A. Amplified products were evaluated by 1.5% (IFN-γ), 3.0% (IL-6, LT-α and TNF-α), and 8.0% (IL-1β) agarose gel electrophoresis containing ethidium bromide (0.5 μg/ml). After electrophoresis, the ethidium-stained fragments were visualized and documented (Multi Doc-ItTM Digital Imaging System, UVP, Upland, California, USA) (Fig. 1).

Primer sequences, restriction enzymes, and size of digestion products in the PCR-RFLP and AS-PCR analyses

Representative amplicons and digestion fragments corresponding to gene polymorphisms. (A) 190 and 114 bp fragments specific for C allele, 304 bp fragment specific for T allele of IL-1b-511C/T polymorphism. (B) 233 bp and 54 bp bands specific for G allele, 122 bp, 111 bp, and 54 bp bands specific for C allele of IL-6-174G/C polymorphism. C) 262 bp fragment specific for T and A alleles of IFN-γ +874T/A polymorphism. (D) 427 bp fragment specific for A allele, 231 and 196 bp fragments specific for G allele of LT-α +252A/G polymorphism. (E) 107 bp fragment specific for A allele, 87 bp fragments specific for G allele of TNF-α −308G/A polymorphism. M, 50 bp DNA marker.
Statistical analysis
The Kolmogorov-Smirnov test was used to test whether the data were distributed normally. Statistical analysis was performed with non-parametric tests (Kruskal-Wallis and Mann-Whitney), with assistance of the IBM SPSS Statistics for Windows package, version 19 (SPSS, Chicago, Illinois, USA). Hardy-Weinberg equilibrium (HWE) analyses (chi-square tests, χ2) were used to compare the expected and observed genotype frequencies in the non-CCA groups of study participants. Allele frequency differences between cases (CCA) and controls (non-CCA) were analyzed, and the odd ratios (OR) and 95% confidence intervals (95% CI) were calculated.
RESULTS
Demographic and clinicopathological characteristics
Table 2 summarizes the clinicopathological characteristics of participants in this study who had been diagnosed with CCA as well as the liver fluke infected cases without CCA. The APF− group included 109 females and 91 males, mean age of 46.8 years (range, 20–60 years) and the APF+ group included 108 females and 92 males, mean age of 47.4 years (range, 25–60 years). The CCA group was composed of 72 men and 38 women, and the mean age was 55.9 years (range, 40–76 years).
Cytokine production in Ov-ES stimulated PBMCs
Fig. 2 summarizes the net production of all cytokines quantified following Ov-ES stimulation of PBMC from the APF−, APF+, and CCA cases. Levels of IL-1β, IFN-γ, and LT-α were significantly higher in CCA patients compared with the APF− and APF+ groups. TNF-α production was significantly higher in CCA and APF+ groups compared with APF−. Notably, IL-6 was significantly higher in APF+ cases compared with the CCA and APF− groups. Differences were not apparent among the groups for 6 other cytokines, IL-2, IL-4, IL-5, IL-8, IL-10, and IL-12p70.
Gene polymorphism and cytokine production
Overall comparison regardless of CCA or non-CCA groups did not reveal significant differences in levels of cytokines among gene polymorphisms except for LT-α, where net cytokine production in LT-α +252A/G heterozygous (AG) and homozygous (GG) variants were significantly elevated compared to the wild type (Fig. 3). However, when subdividing into CCA and non-CCA groups (APF− and APF+), LT-α +254A/G heterozygous (AG), and homozygous (GG) variants showed significantly higher levels of cytokines than the wild type in CCA group; 2.35-fold and 2.24-fold, respectively. Similarly, TNF-α-308G/A heterozygous (GA) and homozygous (AA) variants showed significantly higher levels of cytokines than the wild type in CCA group with fold changes of 12.49 and 4.76, respectively. Interestingly, IFN-γ+874T>A heterozygous (TA) and homozygous (AA) was significantly lower (0.56- and 0.2-fold, respectively) than the wild type in the CCA group.
Distribution of cytokine gene polymorphisms among patient groups
Table 4 presents the genotypes and allele distribution of IL-1β −511C/T, IL-6 −174G/C, IFN-γ+874T/A, LT-α +254A/G, and TNF-α−308G/A polymorphisms in the participants, those with APF−, APF+, and CCA. Comparison among groups did not reveal significant differences in the genotype or allele frequencies for IL-1β−511C/T, LT-α +252A/G, and TNF-α−308G/A. However, comparing the IL6 −174G/C genotypes between the CCA and non-CCA (APF− or APF+) groups revealed that the presence of heterozygous (GC) genotype associated with CCA (OR=2.35 and 2.95 for APF− vs CCA and APF+ vs CCA, respectively). Moreover, a stronger association with CCA was found with the homozygous variant (CC) (OR=5.50 and 7.52 for APF− vs CCA and APF+ vs CCA, respectively). In addition, C alleles at this locus also associated with CCA (OR=3.36) when the CCA and APF+ groups were compared. Comparing between APF+ and APF− groups, IFN-γ+874T>A heterozygous (TA) and homozygous (AA) variants associated with biliary periductal fibrosis (OR=3.10 and 2.74, respectively). The IFN-γ+874T>A homozygous variant (AA) genotype was also associated with CCA (OR=2.20) as compared between APF− and CCA groups, whereas the A allele associated with both fibrosis (OR=1.60, APF− vs APF+) and CCA (OR=1.71, APF− vs CCA). In the case of TNF-α −308G/A polymorphism, only the TNF-α −308 A allele appeared to be associated with CCA (OR=1.61 and OR=1.5 in APF− vs CCA and in APF+ vs CCA, respectively).

Comparison of genotypes and allele frequencies for IL-1β, IL-6, IFN-γ, LT-α, and TNF-α in CCA cases and the participants in the APF+ and the APF− groups
DISCUSSION
Infection with the liver fluke O. viverrini frequently leads to chronic inflammation of the bile ducts and results in unremitting tissue damage that, in turn, can cause hepatobiliary diseases including CCA. A subset of people with opisthorchiasis develop overt hepatobiliary disease, with about 25% of them displaying marked biliary tract fibrosis and as many as 0.5% of infected cases exhibiting incipient bile duct cancer upon ultrasound examination [4,20,21]. Why infection leads to different disease outcomes among the population in high liver fluke infection endemic sites remains to be established. Divergent immunological responses among infected persons may contribute to the outcome, in addition to other co-factors [22]. The pro-inflammatory cytokine IL-6 is significantly elevated in infected individuals with marked and chronic hepatobiliary fibrosis [9,23] and cholangiocytes upregulate expression of IL- 6 and IL-8 in response to the Ov-ES components [24,25]. Moreover, the recent study revealed both pro-inflammatory and anti-inflammatory cytokines were increased in the O. viverrini-associated CCA compared to uninfected normal controls [10].
Here, we investigated pro-inflammatory and anti-inflammatory cytokine profiles among different groups of patients infected with O. viverrini: liver fluke infected individuals exhibiting biliary tract disease (advanced periductal fibrosis, APF+) or without overt biliary periductal fibrosis (APF−), and liver fluke infection-associated CCA cases. We observed that pro-inflammatory cytokines, IL-1β, IL-6, IFN-γ, LT-α, and TNF-α, though no other cytokines were increased significantly during CCA in comparison with participants infected with O. viverrini but without CCA, i.e., the APF− and APF+ groups. These findings indicated that marked inflammatory responses during opisthorchiasis increased the risk for cholangiocarcinogenesis.
Gene polymorphisms are known to affect cytokine levels [8,26,27]. Single nucleotide polymorphism (SNP) can affect gene function in many ways [28]. If it occurs at the non-coding region, it can affect the binding of transcriptional factors, interfere RNA splicing, and increase or decrease messenger RNA degradation. SNP in coding region can result in 2 consequences; synonymous and nonsynonymous SNPs. Synonymous SNPs do not affect the protein sequence, but can still affect the protein conformation such as the silent polymorphism in the MDR1 gene changes substrate specificity [29]. Nonsynonymous SNPs change the amino acid sequence of protein and subsequent protein function as shown in cystic fibrosis patients [30]. For SNPs in our liver fluke associated CCA, we report here that certain cytokine gene polymorphisms affect proinflammatory cytokine production. However, detail molecular mechanisms need to be further investigated.
Accordingly, we hypothesized that cytokine gene polymorphisms may underline variations in inflammatory cytokine response in response to liver fluke infection and downstream biliary tract inflammation and cholangiocarcinogenesis. The gene polymorphisms IL-1β −511C/T, IL6 −174G/C, LT-α +252A/G, and TNF-α−308G/A generated the high-producer phenotypes, whereas IFN-γ+874T/A exhibited a low-production phenotype. In agreement with these reports, both LT-α +252A/G and TNF-α−308G/A heterozygous and homozygous variants showed significant higher levels of cytokine than the wild type. Similarly, both IL-1β −511C/T and IL-6 −174G/C variants produced more cytokines than the wild type, even though the difference was not statistically significant. By contrast, the variants of IFN-γ+874T/A produced less cytokines than the wild type. We note, however, that these cytokine production profiles were observed only in the CCA cases, and not among the non-CCA groups. The difference of cytokine production from the same gene polymorphism between non-CCA and CCA patients might indicate divergence in cytokine regulation/networking among the groups in this study.
Subsequently we investigated whether gene polymorphisms of these pro-inflammatory cytokine were related to the development of fibrosis or CCA. An IL-6 −174G>C polymorphism represents a potential modulator of risk for malignancy with ethnicity and tumor site dependent trends [31]. The present findings revealed that both heterozygous (GC) and homozygous (CC) variant of IL6 −174G>C polymorphisms were associated with cholangiocarcinogenesis in this Thai population. Given that the heterozygous (GC) and homozygous (CC) variants of IL6 −174G>C are high producer phenotypes, these gene polymorphisms might explain elevated levels of IL-6 during in CCA, as reported earlier. The A allele of the TNF-α −308G/A polymorphism appears to associate with the disease progression of hepatitis B viral infection and hepatocellular carcinoma [32]. In likely fashion, the present findings revealed that the A allele of TNF-α −308G/A polymorphism associated with CCA.
Interferon gamma plays key roles in the immune responses and in antitumor mechanisms [33,34]. IFN-γ+874T>A polymorphisms are associated with viral and bacterial infections, which result in chronic inflammation and carcinogenic [35–37]. Here, IFN-γ+874T>A heterozygous (TA) and homozygous (AA) variants associated with fibrosis, and the IFN-γ+874T>A homozygous (AA) variant was also associated with CCA. As noted, these variants produce less IFN-γ than the TT (wild type) genotype. Together these findings indicate that in a subset of liver fluke infected individuals carrying these alleles, the phenotype of which is reduced levels of IFN-γ, are at elevated risk for biliary periductal fibrosis and CCA. Mechanisms underlying the association are not known. However, for anti-tumor activity, IFN-γ plays a central role in the recognition and elimination of transformed cells through its cytostatic, pro-apoptotic, and immune-provoking effects [33].
To summarize, cytokine and gene polymorphism profiles of Thais infected with the liver fluke O. viverrini were characterized. To our knowledge, this is the first demonstration that persons infected with O. viverrini inherited with certain gene polymorphisms that result in aberrant cytokine production are susceptible to develop fibrosis and CCA. Pro-inflammatory cytokines were increased in persons with cholangiocarcinoma (CCA). Different cytokine gene polymorphisms were found among patients with or without fibrosis and CCA. Interestingly, a subset of patients with the IL6 −174G>C polymorphisms, the IFN-γ+874T>A heterozygous (TA) and homozygous (AA) variants, and the A allele of TNF-α −308G/A associated with CCA. Further studies on these gene polymorphisms and in relation to other cytokines are warranted since they may identify genotypes and phenotypes predictive of risk for CCA.
ACKNOWLEDGMENTS
This work was supported by the Higher Education Research Promotion and National Research University Project of Thailand, Office of the Higher Education Commission, through the Health Cluster (SHeP-GMS), Khon Kaen University, the Thailand Research Fund (TRF) grant no. RTA 5680006, the National Institute of Allergy and Infectious Diseases (NIAID), NIH, award no. P50AI098639, and National Cancer Institute (NCI), National Institutes of Health (NIH) award no. CA164719. B. Sripa is a TRF Senior Research Scholar. The content is solely the responsibility of the authors and does not necessarily represent the official views of the TRF, NIAID, NCI, or the NIH or the funders.
Notes
CONFLICT OF INTEREST
We have no conflict of interest related to this study.