Associations Between Milk Performance Traits in Holstein Cows and 16 Candidate SNPs Identified by Arrayed Primer Extension (APEX) Microarray

An oligonucleotide microarray—which allows for parallel genotyping of many SNPs in genes involved in cow milk protein biosynthesis—was used to identify which of the 16 candidate SNPs are associated with milk performance traits in Holstein cows. Four hundred cows were genotyped by the developed and validated microarray. Significant associations were found between four single SNPs, namely DGAT1 (acyloCoA:diacylglycerol acyltransferase), LTF (lactoferrin), CSN3 (kappa-casein), and GHR (growth hormone receptor) and with fat and protein yield and percentage. Many significant associations between combined genotypes (two SNPs) and milk performance traits were found. The associations between the combined genotypes DGAT1/LTF and DGAT1/LEPTIN analyzed traits are presented as examples.

The microarray based on APEX (Arrayed Primer Extension) is a fast and reliable method for multiple SNP analysis of potential application in marker-assisted selection. After further development, the chip may prospectively be used for dairy cattle paternity analysis and evolutionary studies.     

Promoter variants at AP2 box region of Hsp70.1 affect thermal stress response and milk production traits in Frieswal cross bred cattle

Heat shock proteins (Hsp) are known to play major role in protection of cells from thermal stress. Nucleotide polymorphisms within the promoter of Hsp affect degree of expression and inducibility of Hsp mRNA. The present study aimed to investigate the effect of polymorphism within promoter region on the cellular expression of Hsp70.1 mRNA and association of identified polymorphisms with the physiological parameters during summer stress and milk production traits in dairy cattle. Two hundred Frieswal cows were genotyped using double PCR-RFLP to identify deletion of cytosine within the Hsp70.1 promoter AP2 box at base position 895. Homozygous wild type genotypes (CC) were found in lower frequency (39.29, n = 78) than heterozygous cytosine deletion mutant genotypes (C −) (60.71, n = 122). In the observed physiological parameters (rectal temperature, respiration rate and heat tolerance coefficient), cows that were homozygous wild types had better significant (P < 0.05) summer tolerance than the heterozygous deletion genotypes. Cytosine deletion mutation in the promoter region negatively affected (P < 0.01) the expression of Hsp70.1 mRNA in peripheral bovine mononuclear cells (PBMC) subjected to in vitro heat stress. Further association of observed polymorphism with the milk production traits was significant as the heterozygous cytosine deletion cows had lower total milk yield, peak yield, yield at 300 days, protein% (P < 0.01) and fat% (P < 0.05) than the native wild type promoter cows. The results from the present study suggest that the promoter region of bovine hsp70.1 gene is polymorphic and may be useful in selection of dairy cows for relatively better thermotolerance and higher milk production.     

牛α1-AT基因5’侧翼区新SNPS的鉴定及其与产奶性状的相关性

以中国荷斯坦奶牛(Chinese Holstein dairy cattle)为对象,以α1 抗胰蛋白酶基因(α1-AT)为候选基因,扩增5′侧翼区668 bp和999 bp的片段并进行测序.首次发现,在+3 142 bp处P1和+4 408 bp处P2分别发生C-T、T-C突变.随后采用PCR-RFLP方法对随机采自6个牛场,共计294头牛进行了检测,遗传变异和产奶性状分析结果显示：2个位点的等位基因在群体中都有分布,且处于中度多态.P1位点A和B等位基因的频率分别为50.34%和49.66%; AA、AB和BB基因型频率分别为23.81%、53.06%和23.13%;P2位点E和F等位基因的频率分别为30.61%和69.39%, EE、EF和FF基因型频率分别为11.90%、37.41%和50.68%. χ2适合性检验表明,该群体在P1位点的突变达到Hardy-Weinberg平衡状态（P>0.05）,在P2位点未达到平衡.基因与产奶性状关联分析表明,AB基因型个体产奶量与脂蛋比显著高于AA基因型个体(P＜0.05);FF基因型个体乳蛋白率显著高于EF基因型个体(P＜0.05);9种单倍型组合与乳脂率、乳蛋白率、体细胞数、产奶量及脂蛋比均存在不同程度相关性.     

Single nucleotide polymorphisms, haplotypes and combined genotypes of LAP3 gene in bovine and their association with milk production traits

Leucine aminopeptidase 3 (LAP3) is an aminopeptidase which catalyses the removal of N-terminal amino acids and is involved in protein maturation and degradation. In this study, we detected the polymorphisms of LAP3 gene by polymerase chain reaction-single-strand conformation polymorphism (PCR-SSCP) and DNA sequencing methods in 916 individuals from three Chinese cattle breeds including Chinese Holstein, Luxi Yellow and Bohai Black. One novel single nucleotide polymorphism (SNP) (g.24564G>A ss196003366) and four previously deposited SNPs in the GenBank database (g.24794T>G, g.24803T>C, g.24846T>C, g.25415T>C) were detected. Three of the SNPs (g.24794T>G, g.24803T>C, g.24846T>C) were firstly found to be linked completely and regarded as a SNP g.24794M>N by PCR-SSCP and DNA sequencing in the tested breeds. The allelic frequencies and genetic indices of the SNPs were different in three Chinese cattle populations. The SNPs and their genetic effects on milk production traits in Chinese Holsteins were evaluated. Least squares analysis showed that cows with genotype MM had higher fat percentage and protein percentage than genotype NN (P<0.05); and the cows with g.25415T>C-CC genotype had higher protein rate than ones with TT genotype (P<0.05). In addition, eight haplotypes and 23 combined genotypes were identified based on the nine genotypes and the association between combined genotypes and milk production traits were analyzed. Statistic results showed that the cows with genotype combination MAT/MGC have higher protein and fat rate and lower SCS. Our finding demonstrated that the LAP3 gene possibly contributed to conducting association analysis and can be used as molecular marker in milk production traits and other performance for animal breeding.     

Functional polymorphisms of caveolin-1 variants as potential biomarkers of esophageal squamous cell carcinoma

Abstract

Objective: To investigate the association of caveolin-1 (CAV1) genetic variants (C239A (rs1997623), G14713A (rs3807987), G21985A (rs12672038), T29107A (rs7804372)) with esophageal squamous cell carcinoma (ESCC) susceptibility.

Methods: A total of 427 patients with ESCC and 427 healthy controls were genotyped using the polymerase chain reaction and restriction fragment length polymorphism (PCR-RFLP) method.

Results: There were significant differences between patients and controls in distributions of their genotypes and allelic frequencies in G14713A and T29107A polymorphisms. Furthermore, haplotype analysis revealed that haplotypes CAAT and CAGT were associated with high risk for ESCC, while haplotype CGGA was protective against ESCC. Stratified analysis showed the associations between the SNPs (G14713A and T29107A) and ESCC risk were noteworthy among female patients and patients who never smoke or drank alcohol.

Conclusions: Genetic polymorphisms of CAV1 G14713A and T29107A might affect an individual’s susceptibility in developing ESCC, making them efficient potential genetic biomarkers for early detection of ESCC.     

Associations between newly discovered polymorphisms in the Bos taurusgrowth hormone receptor gene and performance traits in Holstein–Friesian dairy cattle

Variations in the growth hormone receptor (GHR) gene sequence are associated with performance traits in cattle. For example, the single nucleotide polymorphism (SNP) F279Y in transmembrane exon 8 has a strong association with milk yield. In this study, 32 previously unreported, putative novel SNPs (31 in the 5′ non‐coding region) were identified by resequencing ～19 kb of the GHR gene in genomic DNA from 22 cattle of multiple breeds. A population of 848 Holstein–Friesian AI sires was subsequently genotyped for the 32 putative novel SNPs and seven published SNPs (including F279Y, one in exon 1A promoter and five in exon 10). Associations between each segregating SNP and genetic merit for performance were quantified in the 848 Holstein–Friesians using weighted animal linear mixed models. Six of the published SNPs and seven of the novel SNPs were associated with at least one of the traits – milk yield, fat yield, protein yield, fat percentage, protein percentage, somatic cell score, calving interval, survival and growth and size traits. Even when the allelic substitution effect (P < 0.001) of F279Y was accounted for, the allelic substitution effect of one of the novel SNPs (GHR4.2) in the 5′ non‐coding region of GHR was associated with a lactation milk yield of 37.46 kg (P < 0.001). GHR4.2 and F279Y were not in linkage disequilibrium (r² = 0.00, D’ = 0.04) in the 848 Holstein–Friesians, indicating that their association with milk yield was independent.     

Identification of single nucleotide polymorphisms in the bovine follicle‐stimulating hormone receptor and effects of genotypes on superovulatory response traits

In dairy cows, there is evidence that failure to respond to superovulation protocols is a heritable trait. In women, genotyping for the p.N680S single nucleotide polymorphism (SNP) in the follicle‐stimulating hormone receptor (FSHR) gene may help identify poor responders before ovarian stimulation is initiated. Our objectives were to identify SNPs in the coding region of the bovine FSHR gene and to investigate the effect of FSHR genotypes on superovulatory response in Holstein cattle. Sequencing of FSHR exons 1–10 revealed seven SNPs. Three were non‐synonymous mutations (c.337C>G, c.871A>G and c.1973C>G). SNP c.337C>G encodes for a proline‐to‐alanine (p.Pro113Ala) amino acid replacement in the extracellular ligand‐binding domain of the receptor. PCR‐RFLP analyses showed that homozygous GG Holstein cows present a higher percentage of viable embryos, whereas GG and CG animals have less unfertilised oocytes. SNP c.871A>G results in an isoleucine‐to‐valine (p.Ile291Val) modification, and homozygous AA animals present lower embryo yield after superovulatory treatments. SNP c.1973C>G corresponds to a threonine‐to‐serine (p.The658Ser) modification in the intracellular carboxyl‐terminal domain of the FSHR protein, and homozygous GG Holstein cows were associated with a lower embryo yield and a higher percentage of unfertilised oocytes. Our results suggest that specific alleles of the bovine FSHR gene are associated with variations in embryo yield and in the number of unfertilised oocytes.     

SNP exploring in the middle and terminal regions of the IGF-1 gene and association with production and reproduction traits in Holstein cattle

Five primer sets were designed in order to identify single nucleotide polymorphisms (SNPs) in middle and terminal exons (2 to 6) and in some flanking intronic regions of the bovine insulin-like growth factor 1 (IGF-1) gene. Sequencing results of PCR products for 10% of animals showed no variant in exons but a SNP at intron 4 was occurred. Both polymerase chain reaction–restriction fragment length polymorphism (PCR–RFLP) and high resolution melting (HRM) methods were developed to genotype samples. The PCR–RFLP results showed the presence of three fragments on agarose gel for the C allele due to two cleavage sites while two fragments for the T allele were observed. Melting curves of 123 bp fragments in HRM analysis showed a difference between temperature melting (T_m) of two homozygous genotypes as the CC genotypes had higher T_m than the TT genotypes. Melting curve of the CT genotype was different and crossed two parallel patterns of homozygous genotypes. The frequencies of the CC, CT and TT genotypes were 0.6, 0.37 and 0.03, respectively. Also, the estimated allele frequencies were 0.785 and 0.215 for the C and T alleles, respectively. Results showed higher accuracy of the HRM analysis compared to the PCR–RFLP method. Least square means (LSMs) comparison of the different genotypes in the SNP showed significant association with milk fat yield trait in the first lactation and open days after the second calving. The polymorphism did not have a significant effect on other milk production or reproduction traits. It seems that other variants or QTLs known in this region underlie genetic variation in the production and reproduction of dairy cattle.     

Haplotype distribution in the GLI3 gene and their associations with growth traits in cattle

The glioma-associated oncogene family zinc finger 3 gene (GLI3) mediates in all vertebrates hedgehog (Hh) signaling that plays an essential role in the induction and patterning of numerous cell types during invertebrate and vertebrate development. In this study, a total of 6 single nucleotide polymorphisms (SNPs: 1–6) were identified by polymerase chain reaction–single stranded conformational polymorphism (PCR–SSCP) and DNA pool sequencing, including all 13 exons and 12 exon–intron boundaries within the bovine GLI3 gene. 16 haplotypes and 13 combined genotypes were revealed and the linkage disequilibrium was assessed in 708 individuals representing three main cattle breeds from China. The statistical analyses indicated that the SNP2, 3 and 4 are associated with the body weight at birth and 6 months in Nanyang cattle population (P < 0.05). No significant association was detected between 11 combined genotypes and body weight at five different ages. Our results provide evidence that polymorphisms in the GLI3 gene are associated with growth traits, and may be used for marker-assisted selection in beef cattle breeding program.     

SYNE1-QK1 SNPs,G × G and G × E interactions on the risk of hyperlipidaemia

This study aimed to assess the relationship of 3 spectrin repeat containing nuclear envelope protein 1 (SYNE1) and 4 KH domain containing RNA binding (QK1) single nucleotide polymorphisms (SNPs), their haplotypes, gene-gene (G × G), gene-environment (G × E) interactions and hypercholesterolaemia (HCH) and hypertriglyceridaemia (HTG) in the Chinese Maonan minority. The genetic make-up of the SYNE1-QK1 SNPs in 1932 unrelated subjects (normal, 641; HCH, 649; and HTG, 642) was obtained by next-generation sequencing technologies. The genotypic frequencies of following SNPs were suggestively distinctive between the control and HCH groups (rs2623963, rs7745725, rs9459317, rs16897566), or between the control and HTG groups (rs2623963, rs1358317, rs7745725, rs1923608, rs16897566 SNPs; P < .05, respectively). Multiple-locus linkage disequilibrium analysis indicated that the identified SNPs were not inherited independently. Several haplotypes and gene-gene interaction haplotypes among the detected SNPs may be related with an increased morbidity of HCH (C-G-A, C-G-G and C-G-G-T-C-A-T) and HTG (C-G-G, G-T-G-C, C-G-G-G-T-G-C and C-G-G-T-C-A-T), whereas others may be related with an decreased risk of HCH (G-A-A, G-C-A-T, C-A-A-T-C-A-T and G-A-A-G-C-A-T) and HTG (G-A-A, G-C-A-T, C-A-A-T-C-A-T and G-A-A-G-C-A-T). The association evaluation based on haplotypes and gene-gene interactions could improve the power of detecting the risk of dyslipidaemia than anyone of SNP alone. There was significant three-locus model involving SNP-SNP, haplotype-haplotype/environment and G × G interactions (P < .05-0.001) that were detected by GMDR in HCH and HTG groups. Different interactions between genetic and environmental factors would produce different redundancy or synergy effects on the morbidity of HCH and/or HTG.     

Identification of splice variants,expression analysis and single nucleotide polymorphisms of the PRMT2 gene in dairy cattle

Protein arginine N-methyltransferase 2 (PRMT2), also named HRMT1L1, belongs to the Bovine Protein arginine N-methyltransferase (PRMT) genes which are involved in the immune response. To explore the variability of the PRMT2 gene and resistance to mastitis in cows, splice variant (SV), and single nucleotide polymorphisms (SNPs) were identified in this study. A SV (PRMT2-SV) lacking exon 7 (98-bp) of the PRMT2 gene was found in healthy and mastitis-infected mammary gland tissues. Two of four SNPs were significantly associated with bovine milk yield and protein content. Further, we estimated the relative expression of PRMT2-SV in the mammary gland tissue of dairy cattle by using quantitative real-time polymerase chain reaction. The result showed that expression of the PRMT2-SV mRNA was significantly upregulated 4.02-fold (p < 0.05) in infected mammary tissues (n = 5) compared to healthy tissues (n = 5). Our findings reveal that PRMT2-SV may play an important role in mastitis resistance in dairy cattle. The SNPs may be used as a possible candidate SNPs for marker-assisted selection and management in Chinese Holstein cattle.     

Combined effects of p53 and MDM2 polymorphisms on susceptibility and surgical prognosis in hepatitis B virus-related hepatocellular carcinoma

The p53 signaling pathway works as a potent barrier to tumor progression. Two single nucleotide polymorphisms (SNPs) in the gene loci of p53 pathway, p53 codon 72 Arg72Pro and MDM2 SNP309 (T>G), have been shown to cause perturbation of p53 function, but the effect of the two SNPs on the risk of hepatocellular carcinoma (HCC) remains inconsistent. This study investigated the influence of combined p53 Arg72Pro and MDM2 SNP309 on the risk of developing HCC in patients with chronic hepatitis B virus infection, and evaluated the significance of the two combined SNPs on patient prognosis. In total, 350 HCC patients, 230 non-HCC patients, and 96 healthy controls were genotyped for the p53 Arg72Pro and MDM2 SNP309. The combined p53 Pro/Pro and MDM2 G/G genotype was significantly associated with HCC risk (P= 0.047). Multivariate analysis indicated that combined p53 Pro/Pro and MDM2 G/G genotype was an independent factor affecting recurrence and survival (P<0.05). Patients with combined p53 Pro/Pro and MDM2 G/G genotypes had a poorer prognosis than other genotypes, P<0.01 for both disease-free survival (DFS) and overall survival (OS). DFS and OS rates also differed significantly between Barcelona Clinic Liver Cancer (BCLC) stage A patients with combined p53 Pro/Pro and MDM2 G/G and other genotypes (P<0.05). Thus, the combined p53 Pro/Pro and MDM2 G/G genotype is associated with increased risk of developing HCC and is an independent adverse prognostic indicator in early stage HCC.     

Differential response to stocking rates and feeding by two genotypes of Holstein-Friesian cows in a pasture-based automatic milking system

The throughput of automatic milking systems (AMS) is likely affected by differential traffic behavior and subsequent effects on the milking frequency and milk production of cows. This study investigated the effect of increasing stocking rate and partial mixed ration (PMR) on the milk production, dry matter intake (DMI), feed conversion efficiency (FCE) and use of AMS by two genotypes of Holstein-Friesian cows in mid-lactation. The study lasted 8 weeks and consisted in a factorial arrangement of two genotypes of dairy cattle, United States Holstein (USH) or New Zealand Friesian (NZF), and two pasture-based feeding treatments, a low stocking rate system (2 cows/ha) fed temperate pasture and concentrate, or a high stocking rate system (HSR; 3 cows/ha) fed same pasture and concentrate plus PMR. A total of 28 cows, 14 USH and 14 NZF, were used for comparisons, with 12 cows, six USH and six NZF, also used for tracking of animal movements. Data were analyzed by repeated measure mixed models for a completely randomized design. No differences (P>0.05) in pre- or post-grazing herbage mass, DMI and FCE were detected in response to increases in stocking rate and PMR feeding in HSR. However, there was a significant (P<0.05) grazing treatment×genotype×week interaction on milk production, explained by differential responses of genotypes to changes in herbage mass over time (P<0.001). A reduction (P<0.01) in hours spent on pasture was detected in response to PMR supplementation in HSR; this reduction was greater (P=0.01) for USH than NZF cows (6 v. 2 h, respectively). Regardless of the grazing treatment, USH cows had greater (P=0.02) milking frequency (2.51 v. 2.26±0.08 milkings/day) and greater (P<0.01) milk yield (27.3 v. 16.0±1.2 kg/day), energy-corrected milk (24.8 v. 16.5±1.0 kg/day), DMI (22.1 v. 16.6±0.8 kg/day) and FCE (1.25 v. 1.01±0.06 kg/kg) than NZF cows. There was also a different distribution of milkings/h between genotypes (P<0.001), with patterns of milkings/h shifting (P<0.001) as a consequence of PMR feeding in HSR. Results confirmed the improved FCE of grazing dairy cows with greater milk production and suggested the potential use of PMR feeding as a tactical decision to managing HSR and milkings/day in AMS farms.     

Polymorphisms of <Emphasis Type="Italic">STAT5A</Emphasis> gene and their association with milk production traits in Holstein cows

The STAT5A gene was studied as a candidate gene for five milk production traits (milk yield at 305 days, protein percentage, fat percentage, lactose percentage and dry matter percentage) in Holstein cows. According to the sequence of bovine STAT5A gene, two pairs of primers (P1 and P2) were designed to detect polymorphisms of STAT5A gene in 401 Holstein cows by PCR-RFLP and PCR-SSCP. The results showed that the products amplified by primers P1 and P2 displayed polymorphisms. For P1, three genotypes (AA, AG, and GG) were detected, and the frequency of AA/AG/GG was 0.252/0.486/0.262, respectively. Sequence analysis revealed a single nucleotide substitution A–G at 14217 bp (GenBank NC_007317) of bovine STAT5A gene while compared GG genotype with AA genotype. The differences of the least squares means for the four milk production traits (milk yield at 305 days, fat percentage, lactose percentage and dry matter percentage) between AA, AG and GG were not significant (P > 0.05). Least squares mean of protein percentage for AG or GG was significantly higher than that for AA (P < 0.05); the difference of the least squares mean for protein percentage was not significant between AG and GG (P > 0.05). For P2, three genotypes (CC, CT, and TT) were detected in Holstein cows, and the frequency of CC/CT/TT was 0.751/0.234/0.015, respectively. Sequencing revealed an insertion CCT at 17266 (NC_007317) of bovine STAT5A gene while compared CC genotype with TT genotype. The differences of the least squares means for the three milk production traits (protein percentage, lactose percentage and dry matter percentage) between CC, CT and TT were not significant (P > 0.05). Least squares mean of milk yield at 305 days for TT or CT was significantly higher than that for CC (P < 0.05); the difference of the least squares mean for milk yield at 305 days was not significant between TT and CT (P > 0.05). Least squares mean of fat percentage for CC or CT was significantly higher than that for TT (P < 0.05); the difference of the least squares mean for fat percentage was not significant between CC and CT (P > 0.05). The results preliminarily indicated that allele G of A14217G polymorphic site of STAT5A gene is a potential DNA marker for improving protein percentage in dairy cattle, 17266indelCCT polymorphic site of STAT5A gene is a potential DNA marker for improving milk yield at 305 days and fat percentage in dairy cattle.     

Assessment of association between variants and haplotypes of the IGF2 gene in beef cattle

Insulin-like growth factor 2 (IGF2) is a fetal growth and differentiation factor that plays an important role in muscle growth and in myoblast proliferation and differentiation. The aim of this study was to examine the association of the IGF2 polymorphism with growth traits in beef cattle breed. Four single nucleotide polymorphisms (SNPs: 1–4) were identified in the bovine IGF2 by sequencing pooled DNA samples (Pool-Seq) and forced polymerase chain reaction–restriction fragment length polymorphism (Forced PCR–RFLP) methods. The result of haplotype analysis of four SNPs showed that eight haplotypes and eighteen combined genotypes were revealed, and the linkage disequilibrium and evolutionary relationship were assessed in 1522 individuals representing four purebred cattle breeds from China. The statistical analyses indicated that the 4 SNPs and 18 combined genotypes or haplotypes are associated with the body weight at 18 and 24 months in Jiaxian cattle population (P < 0.05 or P < 0.01). Our results provide evidence that polymorphisms in the IGF2 gene are associated with growth traits, and may be used for marker-assisted selection in beef cattle breeding program.     

Leptin −2548 G/A polymorphisms are associated to clinical progression of oral cancer and sensitive to oral tumorization in nonsmoking population

Oral cancer is causally associated with environmental carcinogens, and the susceptibility to carcinogen-mediated tumorigenesis is proposed to be genotype-dependent. Leptin (LEP) and leptin receptor (LEPR) both play a crucial role in the mediation of physiological reactions and carcinogenesis and may serve as a candidate biomarker of oral cancer. The current case-control study aimed to examine the effects of LEP −2548 G/A (rs7799039), LEPR K109R (rs1137100), and LEPR Q223R (rs1137101) single-nucleotide polymorphisms (SNPs) with or without interacting to environmental carcinogens on the risk for oral squamous cell carcinoma. The SNPs of three genetic allele, from 567 patients with oral cancer and 560 healthy controls in Taiwan were analyzed. The results shown that the patients with polymorphic allele of LEP −2548 have a significant low risk for the development of clinical stage (A/G: adjusted odds ratio [AOR] = 0.670, 95% confidence interval [CI] = 0.454-0.988, P < 0.05; A/G + G/G: AOR = 0.676, 95% CI = 0.467-0.978, P < 0.05) compared to patients with ancestral homozygous A/A genotype. In addition, an interesting result was found that the impact of LEP −2548 G/A SNP on oral carcinogenesis in subjects without tobacco consumption is higher than subjects with tobacco consumption. These results suggest that the genetic polymorphism of LEP −2548 G/A (rs7799039), LEPR K109R (rs1137100), and LEPR Q223R (rs1137101) were not associated to the susceptibility of oral cancer; SNP in LEP −2548 G/A showed a poor clinicopathological development of oral cancer; population without tobacco consumption and with polymorphic LEP −2548 G/A gene may significantly increase the risk to have oral cancer.     

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-γ (+874 T/A), IL-10 (-1,082 G/A, -819 C/T, -592 C/A), TNF-α (-308 G/A), and TGF-β1 (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-γ intermediate-producing genotype, +874AT (80.0%, P = 0.177), and in those with TNF-α low-producing genotype, -308GG (63.0%, P = 0.148). According to the combination of IFN-γ and TNF-α genotypes, the risks for IHDD could be stratified into high (intermediate-producing IFN-γ and low producing TNF-α), moderate, and low (low-producing IFN-γ and high producing TNF-α) 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-γ and TNF-α genes may contribute to the modulation of fibrosis in the intrahepatic bile duct wall in clonorchiasis patients.