Transcriptomic profiling of three-dimensional cholangiocyte spheroids long term exposed to repetitive Clonorchis sinensis excretory-secretory products.

BACKGROUND: Biliary tract infection with the carcinogenic human liver fluke, Clonorchis sinensis, provokes chronic inflammation, epithelial hyperplasia, periductal fibrosis, and even cholangiocarcinoma. Complications are proportional to the intensity and duration of the infection. In addition to mechanical irritation of the biliary epithelia from worms, their excretory-secretory products (ESPs) cause chemical irritation, which leads to inflammation, proliferation, and free radical generation. METHODS: A three-dimensional in vitro cholangiocyte spheroid culture model was established, followed by ESP treatment. This allowed us to examine the intrinsic pathological mechanisms of clonorchiasis via the imitation of prolonged and repetitive in vivo infection. RESULTS: Microarray and RNA-Seq analysis revealed that ESP-treated cholangiocyte H69 spheroids displayed global changes in gene expression compared to untreated spheroids. In ESP-treated H69 spheroids, 185 and 63 probes were found to be significantly upregulated and downregulated, respectively, corresponding to 209 genes (p < 0.01, fold change > 2). RNA-Seq was performed for the validation of the microarray results, and the gene expression patterns in both transcriptome platforms were well matched for 209 significant genes. Gene ontology analysis demonstrated that differentially expressed genes were mainly classified into immune system processes, the extracellular region, and the extracellular matrix. Among the upregulated genes, four genes (XAF1, TRIM22, CXCL10, and BST2) were selected for confirmation using quantitative RT-PCR, resulting in 100% similar expression patterns in microarray and RNA-Seq. CONCLUSIONS: These findings broaden our understanding of the pathological pathways of liver fluke-associated hepatobiliary disorders and suggest a novel therapeutic strategy for this infectious cancer.

High-quality reference genome for Clonorchis sinensis.

The Chinese liver fluke, Clonorchis sinensis, causes the disease clonorchiasis, affecting ~35 million people in regions of China, Vietnam, Korea and the Russian Far East. Chronic clonorchiasis causes cholangitis and can induce a malignant cancer, called cholangiocarcinoma, in the biliary system. Control in endemic regions is challenging, and often relies largely on chemotherapy with one anthelmintic, called praziquantel. Routine treatment carries a significant risk of inducing resistance to this anthelmintic in the fluke, such that the discovery of new interventions is considered important. It is hoped that the use of molecular technologies will assist this endeavour by enabling the identification of drug or vaccine targets involved in crucial biological processes and/or pathways in the parasite. Although draft genomes of C. sinensis have been published, their assemblies are fragmented. In the present study, we tackle this genome fragmentation issue by utilising, in an integrated way, advanced (second- and third-generation) DNA sequencing and informatic approaches to build a high-quality reference genome for C. sinensis, with chromosome-level contiguity and curated gene models. This substantially-enhanced genome provides a resource that could accelerate fundamental and applied molecular investigations of C. sinensis, clonorchiasis and/or cholangiocarcinoma, and assist in the discovery of new interventions against what is a highly significant, but neglected disease-complex.

Long Non-Coding RNA and mRNA Expression Analysis in Liver of Mice With Clonorchis sinensis Infection.

Clonorchiasis is recognized as an important zoonotic parasitic disease worldwide. However, the roles of host long non-coding RNAs (lncRNAs) and messenger RNAs (mRNAs) in the response to Clonorchis sinensis (C. sinensis) infection remain unknown. Here we compared the expression of lncRNAs and mRNAs in the liver tissue of mice infected with C. sinensis, in order to further understand the molecular mechanisms of clonorchiasis. A total of 388 lncRNAs and 1,172 mRNAs were found to be differentially expressed with absolute value of fold change (FC) ≥ 2.0 and p < 0.05 by microarray. Compared with controls, Gm6135 and 4930581F22Rik were the most over- and under-expressed lncRNAs; flavin-containing monooxygenase 3 (Fmo3) and deleted in malignant brain tumors 1 (Dmbt1) were the most over- and under-expressed mRNAs. Moreover, functional annotation showed that the significantly different mRNAs were related with "FOXO signaling pathway", "Wnt signaling pathway", and "AMPK signaling pathway". Remarkably, lncRNA Gm8801 were significantly correlated with mRNA glycerol-3-phosphate acyltransferase mitochondrial (Gpam), insulin receptor substrate 2 (Irs2), and tumor necrosis factor receptor superfamily member 19 (Tnfrsf19) in ceRNA networks. These results showed that the expression profiles of lncRNAs and mRNAs in the liver changed after C. sinensis infection. Our results provided valuable insights into the lncRNAs and mRNAs involved in clonorchiasis pathogenesis, which may be useful for future control strategies.

New insights from Opisthorchis felineus genome: update on genomics of the epidemiologically important liver flukes.

BACKGROUND: The three epidemiologically important Opisthorchiidae liver flukes Opisthorchis felineus, O. viverrini, and Clonorchis sinensis, are believed to harbour similar potencies to provoke hepatobiliary diseases in their definitive hosts, although their populations have substantially different ecogeographical aspects including habitat, preferred hosts, population structure. Lack of O. felineus genomic data is an obstacle to the development of comparative molecular biological approaches necessary to obtain new knowledge about the biology of Opisthorchiidae trematodes, to identify essential pathways linked to parasite-host interaction, to predict genes that contribute to liver fluke pathogenesis and for the effective prevention and control of the disease. RESULTS: Here we present the first draft genome assembly of O. felineus and its gene repertoire accompanied by a comparative analysis with that of O. viverrini and Clonorchis sinensis. We observed both noticeably high heterozygosity of the sequenced individual and substantial genetic diversity in a pooled sample. This indicates that potency of O. felineus population for rapid adaptive response to control and preventive measures of opisthorchiasis is higher than in O. viverrini and C. sinensis. We also have found that all three species are characterized by more intensive involvement of trans-splicing in RNA processing compared to other trematodes. CONCLUSION: All revealed peculiarities of structural organization of genomes are of extreme importance for a proper description of genes and their products in these parasitic species. This should be taken into account both in academic and applied research of epidemiologically important liver flukes. Further comparative genomics studies of liver flukes and non-carcinogenic flatworms allow for generation of well-grounded hypotheses on the mechanisms underlying development of cholangiocarcinoma associated with opisthorchiasis and clonorchiasis as well as species-specific mechanisms of these diseases.

[Genomics and transcriptomics of the Chinese liver fluke Clonorchis sinensis (Opisthorchiidae, Trematoda)].

The review summarizes the results of first genomic and transcriptomic investigations of the liver fluke Clonorchis sinensis (Opisthorchiidae, Trematoda). The studies mark the dawn of the genomic era for opisthorchiids, which cause severe hepatobiliary diseases in humans and animals. Their results aided in understanding the molecular mechanisms of adaptation to parasitism, parasite survival in mammalian biliary tracts, and genome dynamics in the individual development and the development of parasite-host relationships. Special attention is paid to the achievements in studying the codon usage bias and the roles of mobile genetic elements (MGEs) and small interfering RNAs (siRNAs). Interspecific comparisons at the genomic and transcriptomic levels revealed molecular differences, which may contribute to understanding the specialized niches and physiological needs of the respective species. The studies in C. sinensis provide a basis for further basic and applied research in liver flukes and, in particular, the development of efficient means to prevent, diagnose, and treat clonorchiasis.

Dysregulation of hepatic microRNA expression profiles with Clonorchis sinensis infection.

BACKGROUND: Clonorchiasis remains an important zoonotic parasitic disease worldwide. The molecular mechanisms of host-parasite interaction are not fully understood. Non-coding microRNAs (miRNAs) are considered to be key regulators in parasitic diseases. The regulation of miRNAs and host micro-environment may be involved in clonorchiasis, and require further investigation. METHODS: MiRNA microarray technology and bioinformatic analysis were used to investigate the regulatory mechanisms of host miRNA and to compare miRNA expression profiles in the liver tissues of control and Clonorchis sinensis (C. sinensis)-infected rats. RESULTS: A total of eight miRNAs were downregulated and two were upregulated, which showed differentially altered expression profiles in the liver tissue of C. sinensis-infected rats. Further analysis of the differentially expressed miRNAs revealed that many important signal pathways were triggered after infection with C. sinensis, which were related to clonorchiasis pathogenesis, such as cell apoptosis and inflammation, as well as genes involved in signal transduction mechanisms, such as pathways in cancer and the Wnt and Mitogen-activated protein kinases (MAPK) signaling pathways. CONCLUSIONS: The present study revealed that the miRNA expression profiles of the host were changed by C. sinensis infection. This dysregulation in miRNA expression may contribute to the etiology and pathophysiology of clonorchiasis. These results also provide new insights into the regulatory mechanisms of miRNAs in clonorchiasis, which may present potential targets for future C. sinensis control strategies.

Peroxiredoxin 6 expression is inversely correlated with nuclear factor-κB activation during Clonorchis sinensis infestation.

Clonorchis sinensis is a carcinogenic human liver fluke. Its infection promotes persistent oxidative stress and chronic inflammation environments in the bile duct and surrounding liver tissues owing to direct contact with worms and their excretory-secretory products (ESPs), provoking epithelial hyperplasia, periductal fibrosis, and cholangiocarcinogenesis. We examined the reciprocal regulation of two ESP-induced redox-active proteins, NF-κB and peroxiredoxin 6 (Prdx6), during C. sinensis infection. Prdx6 overexpression suppressed intracellular free-radical generation by inhibiting NADPH oxidase2 and inducible nitric oxide synthase activation in the ESP-treated cholangiocarcinoma cells, substantially attenuating NF-κB-mediated inflammation. NF-κB overexpression decreased Prdx6 transcription levels by binding to two κB sites within the promoter. This transcriptional repression was compensated for by other ESP-induced redox-active transcription factors, including erythroid 2-related factor 2 (Nrf2), hypoxia inducible factor 1α (HIF1α), and CCAAT/enhancer-binding protein ß (C/EBPß). Distribution of immunoreactive Prdx6 and NF-κB was distinct in the early stages of infection in mouse livers but shared concomitant localization in the later stages. The intensity and extent of their immunoreactive staining in infected mouse livers are proportional to lesion severity and infection duration. The constitutive elevations of Prdx6 and NF-κB during C. sinensis infection may be associated with more severe persistent hepatobiliary abnormalities mediated by clonorchiasis.

Characterization and identification of differentially expressed microRNAs during the process of the peribiliary fibrosis induced by Clonorchis sinensis.

Clonorchis sinensis (C. sinensis) infection can lead to biliary fibrosis. MicroRNAs (miRNAs) play important roles in regulation of genes expression in the liver diseases. However, the differential expression of miRNAs that probably regulates the portal fibrogenesis caused by C. sinensis has not yet been investigated. Hepatic miRNAs expression profiles from C. sinensis-infected mice at different time-points were analyzed by miRNA microarray and validated by quantitative real-time PCR (qRT-PCR). 349 miRNAs were differentially expressed in the liver of the C. sinensis-infected mice at 2, 8 or 16weeks post infection (p.i.), compared with those at 0week p.i., and there were 143 down-regulated and 206 up-regulated miRNAs among them. These all dysregulated miRNAs were potentially involved in the pathological processes of clonorchiasis by regulation of cancer-related signaling pathway, TGF-ß signaling pathway, MAPK signaling pathway, Toll-like receptor signaling pathway, PI3K /AKT signaling pathway, etc. 169 of these dysregulated miRNAs were predicted to be involved in the TGF/Smads signaling pathway which plays an important role in the biliary fibrosis caused by C. sinensis. Additionally, miRNA-32, miRNA-34a, miRNA-125b and miRNA-497 were negatively correlated with Smad7 expression, indicating these miRNAs may specifically down-regulate Smad7 expression and participate in regulation of biliary fibrosis caused by C. sinensis. The results of the present study for the first time demonstrated that miRNAs were differentially expressed in the liver of mice infected by C. sinensis, and these miRNAs may play important roles in regulation of peribiliary fibrosis caused by C. sinensis, which may provide possible therapeutic targets for clonorchiasis.

Oxidative stress-mediated mouse liver lesions caused by Clonorchis sinensis infection.

Clonorchis sinensis is a high-risk pathogenic helminth that strongly provokes inflammation, epithelial hyperplasia, periductal fibrosis, and even cholangiocarcinoma in chronically infected individuals. Chronic inflammation is associated with an increased risk of various cancers due to the disruption of redox homeostasis. Accordingly, the present study was conducted to examine the time course relationship between histopathological changes and the appearance of oxidative stress markers, including lipid peroxidation, enzymes involved in lipid peroxidation, and mutagenic DNA adducts in the livers of mice infected with C. sinensis, as well as proinflammatory cytokines in infected mouse sera. Histopathological phenotypes such as bile duct epithelial hyperplasia, periductal fibrosis, edema and inflammatory infiltration increased in infected livers in a time-dependent manner. Intense immunoreactivity of lipid peroxidation products (4-hydroxy-2-nonenal; malondialdehyde), cyclooxygenase-2, 5-lipoxygenase and 8-oxo-7,8-dihydro-2'-deoxyguanosine were concomitantly observed in these injured regions. We also found elevated expressions of cyclooxygenase-2 and 5-lipoxygenase in C. sinensis excretory-secretory product-treated cholangiocarcinoma cells. Moreover, the levels of proinflammatory cytokines such as TNF-α, ILß-1 and IL-6 were differentially upregulated in infected sera. With regard to oxidative stress-mediated carcinogenesis, our findings suggest that C. sinensis infestation may disrupt host redox homeostasis, creating a damaging environment that favors the development of advanced hepatobiliary diseases such as clonorchiasis-associated cholangiocarcinoma.

Involvement of PSMD10, CDK4, and Tumor Suppressors in Development of Intrahepatic Cholangiocarcinoma of Syrian Golden Hamsters Induced by Clonorchis sinensis and N-Nitrosodimethylamine.

BACKGROUND: Clonorchis sinensis is a group-I bio-carcinogen for cholangiocarcinoma (CCA). Although the epidemiological evidence links clonorchiasis and CCA, the underlying molecular mechanism involved in this process is poorly understood. In the present study, we investigated expression of oncogenes and tumor suppressors, including PSMD10, CDK4, p53 and RB in C. sinensis induced hamster CCA model. METHODS: Different histochemical/immunohistochemical techniques were performed to detect CCA in 4 groups of hamsters: uninfected control (Ctrl.), infected with C. sinensis (Cs), ingested N-nitrosodimethylamine (NDMA), and both Cs infected and NDMA introduced (Cs+NDMA). The liver tissues from all groups were analyzed for gene/protein expressions by quantitative PCR (qPCR) and western blotting. PRINCIPAL FINDINGS: CCA was observed in all hamsters of Cs+NDMA group with well, moderate, and poorly differentiated types measured in 21.8% ± 1.5%, 13.3% ± 1.3%, and 10.8% ± 1.3% of total tissue section areas respectively. All CCA differentiations progressed in a time dependent manner, starting from the 8th week of infection. CCA stroma was characterized with increased collagen type I, mucin, and proliferative cell nuclear antigen (PCNA). The qPCR analysis showed PSMD10, CDK4 and p16INK4 were over-expressed, whereas p53 was under-expressed in the Cs+NDMA group. We observed no change in RB1 at mRNA level but found significant down-regulation of RB protein. The apoptosis related genes, BAX and caspase 9 were found downregulated in the CCA tissue. Gene/protein expressions were matched well with the pathological changes of different groups except the NDMA group. Though the hamsters in the NDMA group showed no marked pathological lesions, we observed over-expression of Akt/PKB and p53 genes proposing molecular interplay in this group which might be related to the CCA initiation in this animal model. CONCLUSIONS/SIGNIFICANCE: The present findings suggest that oncogenes, PSMD10 and CDK4, and tumor suppressors, p53 and RB, are involved in the carcinogenesis process of C. sinensis induced CCA in hamsters.

The expression dynamics of transforming growth factor-ß/Smad signaling in the liver fibrosis experimentally caused by Clonorchis sinensis.

BACKGROUND: Liver fibrosis is a hallmark of clonorchiasis suffered by millions people in Eastern Asian countries. Recent studies showed that the activation of TGF-ß/Smad signaling pathway can potently regulate the hepatic fibrogenesis including Schistosoma spp. and Echinococcus multilocularis-caused liver fibrosis. However, little is known to date about the expression of transforming growth factor-ß (TGF-ß) and other molecules in TGF-ß/Smad signaling pathway which may play an important role in hepatic fibrosis caused by C. sinensis. METHODS: A total of 24 mice were individually infected orally with 45 metacercariae, both experimental mice and mocked-infected control mice were anesthetized at 4 week post-infection (wk p.i.), 8 wk p.i. and 16 wk p.i., respectively. For each time-point, the liver and serum from each animal were collected to analyze histological findings and various fibrotic parameters including TGF-ß1, TGF-ß receptors and down-stream Smads activation, as well as fibrosis markers expression. RESULTS: The results showed that collagen deposition indicated by hydroxyproline content and Masson's trichrome staining was increased gradually with the development of infection. The expression of collagen type α1 (Col1a) mRNA transcripts was steadily increased during the whole infection. The mRNA levels of Smad2, Smad3 as well as the protein of Smad3 in the liver of C. sinensis-infected mice were increased after 4 wk p.i. (P < 0.05, compared with normal control) whereas the TGF-ß1, TGF-ß type I receptor (TGFßRI) and TGF-ß type II receptor (TGFßRII) mRNA expression in C. sinensis-infected mice were higher than those of normal control mice after 8 wk p.i. (P < 0.05). However, the gene expression of Smad4 and Smad7 were peaked at 4 wk p.i. (P < 0.05), and thereafter dropped to the basal level at 8 wk p.i., and 16 wk p.i., respectively. The concentrations of TGF-ß1 in serum in the C. sinensis-infected mice at 8 wk p.i. and 16 wk p.i (P < 0.05) were significantly higher than those in the control mice. CONCLUSIONS: The results of the present study indicated for the first time that the activation of TGF-ß/Smad signaling pathway might contribute to the synthesis of collagen type I which leads to liver fibrosis caused by C. sinensis.

Pathological Lesions and Inducible Nitric Oxide Synthase Expressions in the Liver of Mice Experimentally Infected with Clonorchis sinensis.

The nitric oxide (NO) formation and intrinsic nitrosation may be involved in the possible mechanisms of liver fluke-associated carcinogenesis. We still do not know much about the responses of inducible NO synthase (iNOS) induced by Clonorchis sinensis infection. This study was conducted to explore the pathological lesions and iNOS expressions in the liver of mice with different infection intensity levels of C. sinensis. Extensive periductal inflammatory cell infiltration, bile duct hyperplasia, and fibrosis were commonly observed during the infection. The different pathological responses in liver tissues strongly correlated with the infection intensity of C. sinensis. Massive acute spotty necrosis occurred in the liver parenchyma after a severe infection. The iNOS activity in liver tissues increased, and iNOS-expressing cells with morphological differences were observed after a moderate or severe infection. The iNOS-expressing cells in liver tissues had multiple origins.

Comparison of two serpins of Clonorchis sinensis by bioinformatics, expression, and localization in metacercaria.

Clonorchiasis, which has been an important public health problem in China, is caused by ingestion of raw or undercooked fish contaminated by live metacercaria. Therefore, preventing fish from infecting is of great significance for controlling the disease. SERPINs (serine protease inhibitors) are well known as negative regulators of hemostasis, thrombolysis, and innate immune responses. In the present study, two full-length sequences encoding SERPIN were identified from metacercaria cDNA library of Clonorchis sinensis (C. sinensis) and were denominated as CsSERPIN and CsSERPIN3, respectively. Bioinformatics analysis showed that the two sequences shares 35.9% identity to each other. Both of the sequences have SERPIN domain and the greatest difference between the two domains is the reactive centre loop. Transmembrane region was found in CsSERPIN3 while not in CsSERPIN. The expression of the two CsSERPINs was significantly higher at the life stage of metacercaria than that of adult. The transcription levels of CsSERPIN and CsSERPIN3 at metacercaria stage were 3.249- and 11.314-fold of that at adult stage, respectively. Furthermore, the expression of CsSERPIN was 4.32-fold of that of CsSERPIN3 at metacercaria stage. Immunobiochemistry revealed that CsERPIN was dispersed at subtegument and oral sucker of metacercaria, while CsSERPIN3 localized intensely in the tegument of metacercaria of C. sinensis inside of the cyst wall. All these indicated that the CsSERPINs play important roles at metacercaria stage of the parasite. CsSERPIN may take part in regulation of endogenous serine proteinase and CsSERPIN3 may be involved in immune evasion and be a potential candidate for vaccine and drug target for clonorchiasis.

Biochemical characterization and functional analysis of fructose-1,6-bisphosphatase from Clonorchis sinensis.

Fructose-1,6-bisphosphatase (FBPase), a key regulatory enzyme of gluconeogenesis, plays an essential role in metabolism and development of most organisms. To the wealth of available knowledge about FBPase from Clonorchis sinensis (CsFBPase), in this study, the characteristics of CsFBPase and its potential role in pathogenesis of clonorchiasis were investigated. The Km value of CsFBPase was calculated to be 41.9 uM. The optimal temperature and pH of CsFBPase were 37 °C and pH 7.5-8.0, respectively. In addition, Mg(2+) or K(+) played a regulatory role in enzyme activity of CsFBPase. Both transcriptional and translational level of CsFBPase were higher in metacercariae (one of larva stages) than those in adult worm (P < 0.05). CsFBPase were observed to extensively express in the intestine, vitellaria and tegument of adult worms and ubiquitously in metacercariae. Moreover, CsFBPase was confirmed as a component of excretory/secretory products. Consequently, the translocation of CsFBPase could be detected on epithelial cells of bile duct in liver of C. sinensis infected rat. Recombinant CsFBPase can specifically bind to the membrane of human hepatic stellate cell line LX-2 by immunofluorescence analysis and stimulated proliferation and activation of LX-2 which demonstrated by Cell Counting Kit-8 and upregulation of key fibrosis-related factors, such as α-smooth muscle actin, collagen I and collagen III using qRT-PCR. Thus, we predicated that CsFBPase might be a multifunctional enzyme which played as both regulatory enzyme and virulence factor in pathogenesis of C. sinensis infection.

Influencing factors for cure of clonorchiasis by praziquantel therapy: infection burden and CYP3A5 gene polymorphism.

Chemotherapy of clonorchiasis with praziquantel (PZQ) is effective but about 15% of treated cases have been reported uncured. The present study investigated correlation of single nucleotide polymorphisms (SNPs) of the cytochrome P450 gene, CYP3A5 and cure of clonorchiasis. A total of 346 egg passing residents were subjected and treated by 3 doses of 25 mg/kg PZQ. Reexamination recognized 33 (9.5%) uncured and 313 cured. Numbers of eggs per gram of feces (EPGs) before treatment were significantly lower in the cured group than in the uncured group (2,011.2±3,600.0 vs 4,998.5±7,012.0, P<0.001). DNAs of the subjects were screened for SNPs at 7 locations of CYP3A5 using PCR. In the uncured group, the SNP frequencies at g.-20555G>A and g.27526C>T of CYP3A5 were 15.2% and 9.1% while those were 3.8% and 1.0%, respectively, in the cured group. The cure rate was significantly lower in the cases with SNP at g.27526C>T and EPGs≥1,000. In conclusion, EPGs and SNPs of CYP3A5 are factors which influence cure of clonorchiasis by PZQ therapy. It is strongly suggested to recommend 2-day medication for individuals with high EPGs≥1,000.

Transcriptional induction of minichromosome maintenance protein 7 (Mcm7) in human cholangiocarcinoma cells treated with Clonorchis sinensis excretory-secretory products.

Clonorchiasis is an infection associated with bile duct malignancy and subsequent development of cholangiocarcinoma. This disease is mainly caused by Clonorchis sinensis worms and their excretory-secretory products (ESP). However, the precise molecular mechanisms of carcinogenesis remain to be determined. Previously, we established differential gene expression profiles from microarrays containing 23,920 human genes of known function in a human cholangiocarcinoma cell line, HuCCT1, treated with ESP. Among the upregulated genes, we focused on minichromosome maintenance protein 7 (Mcm7), which is implicated in various cancer types, and analyzed transcriptional regulation mediated by ESP to further elucidate its role in cholangiocarcinoma development. Global histone acetylation levels were increased in ESP-treated cells, along with histone acetyltransferase (HAT) protein expression. Detailed promoter analysis using reporter and chromatin immunoprecipitation assays revealed that transcriptional activation of Mcm7 is mediated by HAT recruitment to the promoter region upon C. sinensis ESP treatment. These findings contribute to clarification of the intrinsic mechanism underlying the cellular carcinogenesis process stimulated by Mcm7 in C. sinensis-treated host cells.

Single nucleotide polymorphisms of cytokine genes are associated with fibrosis of the intrahepatic bile duct wall in human clonorchiasis.

This study examined the association of cytokine gene polymorphisms with intrahepatic bile duct wall fibrosis in human clonorchiasis. A total of 240 residents in Heilongjiang, China underwent ultrasonography, blood sampling, and stool examination. Single nucleotide polymorphism (SNP) sites for IFN-gamma (+874 T/A), IL-10 (-1,082 G/A, -819 C/T, -592 C/A), TNF-alpha (-308 G/A), and TGF-beta1 (codon 10 T/C, codon 25 G/C) genes were observed with the TaqMan allelic discrimination assay. No significant correlation was observed between individual cytokine gene polymorphisms and intrahepatic duct dilatation (IHDD). Among individuals with clonorchiasis of moderate intensity, the incidence of IHDD was high in those with IFN-gamma intermediate-producing genotype, +874AT (80.0%, P = 0.177), and in those with TNF-alpha low-producing genotype, -308GG (63.0%, P = 0.148). According to the combination of IFN-gamma and TNF-alpha genotypes, the risks for IHDD could be stratified into high (intermediate-producing IFN-gamma and low producing TNF-alpha), moderate, and low (low-producing IFN-gamma and high producing TNF-alpha) risk groups. The incidence of IHDD was significantly different among these groups (P = 0.022): 88.9% (odds ratio, OR = 24.0) in high, 56.5% (OR = 3.9) in moderate, and 25.0% (OR = 1) in low risk groups. SNP of IFN-gamma and TNF-alpha genes may contribute to the modulation of fibrosis in the intrahepatic bile duct wall in clonorchiasis patients.