Skip to main content

Dynamic expression of cytokine and transcription factor genes during experimental Fasciola gigantica infection in buffaloes



Determining the mechanisms involved in the immune-pathogenesis of the tropical liver fluke, Fasciola gigantica, is crucial to the development of any effective therapeutic intervention. Here, we examined the differential gene expression of cytokines and transcription factors in the liver of F. gigantica-infected buffaloes, over the course of infection.


Water buffaloes (swamp type) were infected orally with 500 F. gigantica encysted metacercariae. Liver tissue samples were collected 3, 10, 28, 42, 70 and 98 days post-infection (dpi). Levels of gene expression of nine cytokines (IFN-γ, TGF-β, IL-1β, IL-4, IL-6, IL-10, IL-12B, IL-13 and IL-17A) and four transcription factors (T-bet, GATA-3, Foxp3 and ROR-γτ) were determined using quantitative real-time PCR (qRT-PCR). We evaluated any correlation between gene expression of these immune-regulatory factors and the severity of liver pathology.


Histopathological examination revealed that cellular infiltration, hemorrhage and fibrosis without calcification in the liver parenchyma of infected buffaloes, increased over the course of infection. This progressive pathology was attributed to dysregulated and excessive inflammatory responses induced by infection. The early infection phase (3–10 dpi) was marked by a generalized immunosuppression and elevated TGF-β expression in order to facilitate parasite colonization. A mixed Th1/Th2 immune response was dominant from 28 to 70 dpi, to promote parasite survival while minimizing host tissue damage. During late infection (98 dpi), the response was biased towards Th1/Treg in order to inhibit the host’s Th2 protective response and promote chronic infection. Both IL-10 and IL-17A and the Th17/Treg balance, played key roles in mediating the inflammatory and immunoregulatory mechanisms in the liver during chronic fasciolosis.


Our data showed distinct CD4+ T helper (Th) polarization and cytokine dysregulation in response to F. gigantica infection in water buffaloes over the course of infection. Characterizing the temporal expression profiles for host immune genes during infection should provide important information for defining how F. gigantica adapts and survives in the liver of buffaloes and how host immune responses influence F. gigantica pathogenicity.


Fasciolosis is a zoonotic parasitic disease caused by infection with the digenetic trematode flukes of the genus Fasciola. While Fasciola hepatica is prevalent in temperate regions, F. gigantica is more widespread in Africa and Asia [1, 2]. Migration of these flukes inside the body of the host causes severe damage to the liver parenchyma and gall-bladder [3,4,5]. Buffaloes are economically important animals for the farming communities in developing countries. Infection of buffaloes with F. gigantica is common in southern China and other geographic regions of the world [6]. Infection can cause poor animal health and significant loss of meat and milk production, with considerable financial implications [3, 7]. Fasciola gigantica flukes specifically target the liver of their definitive host. Effective and balanced local immunity is therefore essential for detecting and controlling these hepatotropic parasites, and for limiting hepatic damage.

Liver flukes are, however, efficient immune-modulators and produce many effectors in order to exploit the host immune response to ensure their survival. A recent study in experimentally infected buffaloes reported a modest increase in the level of Th2-type immune cytokines during early F. gigantica colonization and immunosuppression during chronic F. gigantica infection [8]. Other studies reported a pro-inflammatory or a mixed Th1/Th2 immune response during early infection, and heightened Th2 and Treg responses during chronic infection. This heightened response was assumed to play roles in restoring the host tissue integrity by damping excessive inflammatory response [9, 10]. How inflammation contributes to the pathogenesis of F. gigantica is a complex and multi-faceted story that is still unfolding.

Hepatic immune-inflammatory mechanisms are essential to maintain liver homeostasis and, if dysregulated (e.g. due to parasite infection), can lead to liver pathology and dysfunction. The abnormal production of cytokines and/or transcription factors can lead to inadequate control of Fasciola infection [11, 12]. CD4+ T-cells are subdivided into Th1, Th2, Th17 and regulatory T-cells (Treg) subsets, based on their pattern of cytokine production [13]. Transcription factors T-bet, GATA-3, Foxp3 and ROR-γτ play important roles in the differentiation of Th1, Th2, Treg and Th17 cells respectively, and mediate the production of cytokines in these cells [14,15,16]. Although immunological impairment and polarization of the Th1/Th2 balance are major consequences of F. gigantica-induced liver pathology, the expression profile and dynamic changes of Th1/Th2 cytokines during F. gigantica infection has not been completely elucidated. Also, the role of Th17 and Treg cells in the pathogenesis of F. gigantica infection is still not well-defined.

In the present study, we hypothesized that F. gigantica infection impairs the balance of Th subsets (Th1/Th2/Th17) and Treg, thus contributing to the immune-pathogenesis of fasciolosis. A temporal study of gene expression of nine cytokines (IFN-γ, IL-1β, IL-12B, IL-4, IL-6, IL-10, IL-13, IL-17A and TGF-β) and four transcription factors (T-bet, GATA-3, Foxp3 and ROR-γτ) in the livers from buffaloes infected with F. gigantica was conducted using quantitative real-time PCR (qRT-PCR). Data analyses revealed a large number of differentially regulated genes, which exhibited temporal profiles of expression across the time course study. Our results showed evidence of a step-change in gene expression from an ‘early’ TGF-β-associated immune-suppersive response (3–10 dpi), to a mixed Th1/Th2 immune response (28–70 dpi) and a ‘late’ predominantly Th1/Treg-driven response (98 dpi). These results provide new insights into the dynamic immune response of buffaloes to F. gigantica over the course of experimental infection.


Parasite strain

Eggs of F. gigantica were obtained from the gall bladder and faeces of naturally infected buffaloes slaughtered for human consumption at local abattoirs (Nanning, Guangxi, P.R. China). Protocols used for the preparation of F. gigantica eggs, snail infection with miracidia and harvesting of encysted metacercariae (EM), were performed as previously described [16]. EM were stored in sterile phosphate buffered solution (PBS) at 4 °C. We employed PCR amplification and sequencing of the second internal transcribed spacer (ITS-2) of ribosomal DNA (rDNA) to genotype EM, as described previously [17]. Species identity was confirmed as F. gigantica based on absolute homology to the known ITS-2 sequence of F. gigantica from Guangxi province (GenBank: AJ557569). The viability of EM was examined microscopically and only those with viability greater than 90% were used.


Eight to ten month-old (80–100 kg body weight) water buffaloes (n = 35) were obtained from a local breeder and were identified as swamp type by karyotypic analysis. Buffaloes were kept in separate concrete floor pens. Commercial feed and clean water were provided ad libitum for all animals during the entire study period. None of the buffaloes had been used previously for any experimental procedure. Animals were confirmed as negative in terms of prior infection with liver flukes, by negative fecal examination and negative serum F. gigantica-specific IgG-antibody-based ELISA prior to the start of the study. All animals were treated with a single dose of triclabendazole (5% w/v) in order to eliminate any potential existing fluke infection that may have been missed on laboratory examination. Following triclabendazole treatment, buffaloes were allowed to acclimatize for 30 days to avoid any residual efficacy of the treatment on the establishment of experimental F. gigantica infection.

Animal inoculation and tissue collection

Thirty-five buffaloes were assigned randomly to seven different groups (5 buffaloes/group). Group I was composed of 5 buffaloes that were mock-incoulated with PBS only. Buffaloes in Groups II-VII were each infected with 500 viable metacercariae by oral gavage. Control buffaloes were euthanized at the start of the experiment to obtain baseline values for hepatic tissue pathology and gene expression. Animals from each of the six infected groups were sacrificed and their livers were harvested at 3, 10, 28, 42, 70 and 98 days post infection (dpi), for histopathological and molecular studies. Group III-VII buffaloes were examined clinically on a weekly basis for the development of clinical signs of fasciolosis. The control group served as a baseline point of reference for monitoring the progressive changes in gene expression over the course of infection. However, the inclusion of matched control groups receiving placebo (treated with vehicle only) and euthanized at the same time points as infected groups would have strengthened the power of the study. Alternatively, control buffaloes could have been kept alive (rather than killing them at the start) and liver punch biopsy samples obtained from them at each of the above time points. Despite the early sacrifice of the control animals for reasons of economy and resources, we believe the effects we have observed to be due to the experimental manipulation.

Gross examination and histopathological evaluation

At six time points after infection (indicated above), animals in each infected group were sacrificed, their livers were harvested and examined for pathological lesions and the presence of the flukes. Parasite eggs were recovered by filtering bile fluid through a 0.15 mm pore size mesh. Fasciola gigantica infection was confirmed by observing gross pathological lesions, associated with flukes in the livers and/or by the presence of flukes and eggs in the bile ducts. Samples of liver tissue (~8 g) showing pathological lesions were collected from each animal. Tissue samples were resuspended in 10% PBS-buffered formalin solution overnight, then dehydrated in alcohol, rinsed in xylene, and embedded in paraffin. 3 μm sections of paraffin-embedded tissue were mounted onto glass slides, and stained with hematoxylin and eosin (H&E). Stained tissue sections were examined microscopically at 400× magnification and imaged using a Zeiss Axio Imager manual upright research microscope. Additional liver tissue samples for RNA extraction were collected and kept in RNA store buffer (Tiangen Biotech, Beijing, China), snap frozen in liquid nitrogen and stored at -80 °C.

RNA isolation

Total RNA was extracted from frozen liver tissue samples by RNAprep Pure Tissue Kit (Tiangen Biotech, Beijing, China) following the manufacturer’s instructions. RNA integrity was examined by 2% agarose gel electrophoresis and quantified by NanoDrop 2000/2000c Spectrophotometer analysis (Thermo Scientific, Waltham, US).

Quantification of cytokine and transcription factor gene expression

Quantitative gene expression analysis was performed on liver samples obtained from uninfected control animals, and from infected animals on 3, 10, 28, 42, 70 and 98 dpi. Levels of mRNA expression of nine cytokines (IFN-γ, IL-1β, IL-12B, IL-4, IL-6, IL-10, IL-13, IL-17A and TGF-β) and four transcription factors (T-bet, GATA-3, Foxp3 and ROR-γτ) were determined using quantitative real-time PCR (qRT-PCR). All qRT-PCR primers used in the study are described in Table 1. Complementary DNA (cDNA) was synthesized from 500 ng RNA samples using a PrimerScript™ RT reagent kit (TaKaRa Bio, Dalian, China). qRT-PCR was performed using SYBR®Premix Ex Taq™ II (Tli RNaseH Plus, TaKaRa Bio) and a CFX96 real-time PCR instrument (Bio-Rad, Hercules, US). To determine the specificity of amplification, melting curve analysis was applied to all final PCR products. The efficiency of qRT-PCR, and relative quantification (RQ) of gene expression, were analyzed using the comparative 2–ΔΔCq method [18]. The level of expression of each gene was normalized using glyceraldehyde-3-phosphate dehydrogenase (GAPDH) as the reference housekeeping gene.

Table 1 List of primers used in the SYBR green-based qRT-PCR analysis

Data analysis

Statistical analysis and graph production were performed using GraphPad Prism (GraphPad Software Inc., La Jolla, CA, USA, version 6.02.). Levels of cytokine and transcription factor mRNA expression between uninfected and infected groups were compared at different time points after infection using one-way analysis of variance (ANOVA) with post-hoc LSD multiple comparison tests. Results were presented as F(DFn, DFd) and P-value. Pearson’s correlation coefficients (r-value) were used to detect any correlation between the measured level of Th1, Th2, Treg and Th17 immune cytokines and transcription factors (T-bet, GATA-3, Foxp3 and ROR-γτ) gene expression between infected groups, followed by a two-tailed post-hoc test and presented as a P-value. Data shown represent the mean ± SEM of results from five buffaloes. The level of significance for all analyses was evaluated with a confidence interval > 95% (P < 0.05).


Gross and histopathological attributes

Even though buffaloes did not exhibit clear clinical signs, F. gigantica induced a wide range of pathological lesions over the course of infection (Table 2). During the early stage of infection (3–70 dpi), immature flukes migrated through the intestinal wall, abdominal cavity, liver capsule and liver parenchyma, to reach the bile ducts where they attain sexual maturity. Migration of the juvenile flukes through the host’s tissues was characterized by classical inflammatory signs and accumulation of fibroblasts towards the end of this stage, but without obvious fibrosis. The total number (and average length) of flukes recovered from the livers of infected groups at 28, 42, 70 and 98 dpi were: 3 flukes (1.5 mm), 22 flukes (2.5 mm), 65 flukes (8 mm), and 36 flukes (14 mm), respectively. No fluke was observed in the liver or bile duct before 28 dpi, probably because they were too small to be visible to the naked-eye. Hepatic hemorrhage, swelling, necrosis and viscous bile were first detected at 10 dpi. Tissue fibrosis was observed at 70 dpi, along with disappearance of hemorrhage. Seven parameters were used to describe the histopathological features in buffalo livers during infection (Table 3). Histopathological characteristics of the infected liver tissues (Fig. 1) included an infiltration of inflammatory cells, such as neutrophils and lymphocytes at 3 dpi. As infection progressed, an accumulation of eosinophils, monocytes and red blood cells (RBCs) was detected (10–28 dpi), while fibroblast formation and bile duct hyperplasia were observed later (42 dpi). During the chronic stage (70–98 dpi), significant fibrosis and necrosis, without calcification, were observed. Fluke eggs were detected in the bile and feces of infected animals at 98 dpi. Granulomas with necrotic centers, heavy infiltrates of lymphocytes, RBCs and eosinophils were present at 98 dpi.

Table 2 The presence (+) and absence (−) of gross pathological lesions, flukes and fluke eggs in the liver of buffaloes experimentally infected with F. gigantica
Table 3 The presence (+) and absence (−) of the histopathological changes observed in the liver of buffaloes experimentally infected with F. gigantica
Fig. 1

Histopathological characteristics of the livers of buffaloes infected with Fasciola gigantica. a At 3 dpi, there was local hyperemia associated with mild filtration of lymphocytes and neutrophils. b At 10 dpi, there was scattered vacuolation of hepatocytes, consistent with fat, along with mild to moderate focal necrosis. c At 28 dpi, diffuse intravascular coagulation, severe infiltration of inflammatory cells mainly neutrophils and lymphocytes, and granular degeneration of the cytoplasm were observed. d At 42 dpi, moderate to severe multifocal hemorrhages and necrosis, infiltration of eosinophils, RBCs and monocytes, and accumulation of fibroblasts were detected. e At 70 dpi, there were mild to moderate multifocal bile duct hyperplasia and focal coagulative necrosis associated with collagen deposition. f At 98 dpi, severe periportal fibrosis associated with multifocal inflammatory infiltrate, cellular debris, moderate multifocal hemorrhages, and granulomas with necrotic centres were detected. g Liver tissue from uninfected animal showed normal histological architecture of the hepatic tissue. h Adult flukes in the intrahepatic bile duct, along with epithelial hyperplasia of the duct. In all figures, tissue sections were stained with H&E and arrows point at the corresponding morphological features described above. Scale-bars: a-g, 50 μm; h, 100 μm

Gene expression of cytokines

To determine how immune response changed over the course of infection, gene expression of nine cytokines was assessed in the liver tissue of 29 infected buffaloes and five uninfected buffaloes using qRT-PCR. RNA from one animal from group VII did not pass the quality control check and was therefore excluded from the analysis. Results showed that infection had a significant impact on gene expression of Th1, Th2, Th17, and Treg cytokines (Fig. 2a and Table 4). The transcriptional profile of cytokine genes showed an immunosuppressive state during early infection, a mixed Th1/Th2 response as the infection progressed, and a shifting to Th1/Treg response, associated with greater histopathological changes and fibrosis during late stage infection.

Fig. 2

Temporal changes of the mRNA expression of cytokines and transcription factors in the liver of buffaloes experimentally infected with Fasciola gigantica. The X axis represents days post infection (dpi), control animals are represented by empty bars and Y axis represents the mRNA relative expression of target genes relative to GAPDH based on 2–ΔΔCq calculation. Columns show the means and error bars show SEMs. The log number of mRNA relative expression of a IFN-γ, IL-1β, IL-4, IL-6, IL-10, IL-12B, IL-13, IL-17A, TGF-β and b T-bet, GATA-3, Foxp3, and ROR-γτ, are shown. c Trends in the ratios between each pair of transcription factors during the course of infection, which was used as indicators of the balance between Th1/Th2 and Treg/Th17. Significant differences of each time-point compared with control non-infected (NC) group: *P < 0.05; **P < 0.01; ***P < 0.001 or ****P < 0.0001 (analyzed by one-way ANOVA, post-hoc LSD test)

Table 4 Relative gene expression in the liver of buffaloes experimentally infected with F. gigantica compared to the control group. The relative quantities were obtained according to the comparative Cq method (2–ΔΔCq). P values indicate statistical significance (P < 0.05)

The gene expression of Th1 cytokines (IFN-γ and IL-12B) and IL-10 showed no significant change during early infection. IFN-γ showed significantly higher expression at 70 dpi (P = 0.0079) and 98 dpi (P = 0.0003), when compared with the control. Relative gene expression of both of IL-12B (P < 0.0001) and IL-10 (P = 0.0036) peaked at 98 dpi. Fluke infection also upregulated IL-4 expression at 28 dpi (P = 0.0028), whereas IL-13 significantly increased at 28 dpi (P = 0.0259) and 70 dpi (P = 0.0426), peaking at 98 dpi (P < 0.0001). TGF-β mRNA was rapidly elevated at 3 dpi (P = 0.0003), followed by a decrease to the basal level and remained low thereafter. Upregulation of IL-17A was noted at 28 dpi (P = 0.0114), however expression was relatively low based on a high Cq value, suggesting it to be of little diagnostic value as a marker for Th17 cytokine in this study. The kinetics of the cytokines IL-1β and IL-6 mRNA were similar, both peaking at 28 dpi (P = 0.0007 for IL-1β; P = 0.0002 for IL-6) and 98 dpi (P < 0.0001 for IL-1β; P = 0.0008 for IL-6).

Spearman’s correlation analysis between each pair of cytokines was performed based on the RQ values from infected groups (Fig. 3a). This analysis revealed a significant positive correlation between the expression of IL-12B mRNA and Th2 cytokines (IL-4: r(29) = 0.44, P = 0.0336; IL-13: r(29) = 0.41, P = 0.0542), and between IFN-γ and IL-12B (r(29) = 0.62, P = 0.0017) or IL-13 (r(29) = 0.51, P = 0.0124). The gene expression of IL-10 was significantly correlated with both Th1 and Th2 cytokines (IFN-γ: r(29) = 0.49, P = 0.0171; IL-12B: r(29) = 0.50, P = 0.014; IL-4: r(29) = 0.80, P < 0.0001; IL-13: r(29) = 0.51, P = 0.014) and with IL-17A (r(29) = 0.73, P < 0.0001). Positive correlation was also found between IL-4 and IL-13 (r(29) = 0.44, P = 0.0336) gene expression. There was only a statistically significant relationship between IL-17A expression and IL-4 r(29) = 0.87, P < 0.0001) and IL-10 (r(29) = 0.73, P < 0.0001), but not with any of the other cytokines. Although not significant, the relative mRNA expression of TGF-β showed an inverse relationship to any other cytokine.

Fig. 3

Pearson’s correlation analysis between mRNA expressions of Th1/Th2/Th17/Treg type a cytokines, b transcription factors, and c cytokines vs transcription factors in the liver of buffaloes infected with Fasciola gigantica (tested samples, n = 29) by pairwise comparison. Connecting lines illustrate positive (r > 0) and negative (r < 0) correlation indices

Expression profiles of transcription factors

Fasciola gigantica infection altered the expression of T-bet, GATA-3 and ROR-γτ, with the expression of Foxp3 being least affected (Fig. 2b and Table 4). Compared with the control, mRNA of T-bet was highly expressed at 98 dpi (P = 0.0113), while the expression of the GATA-3, Foxp3 and ROR-γτ genes was upregulated at 28 dpi (P = 0.0035, P = 0.0424 and P = 0.0002, respectively). Spearman’s correlation analysis indicated that the expression of T-bet and Foxp3 was correlated (r(29) = 0.63, P = 0.0011). Positive correlation was also found between the expression of ROR-γτ and GATA-3 (r(29) = 0.43, P = 0.039) or Foxp3 (r(29) = 0.42, P = 0.0439) (Fig. 3b). As shown in Fig. 3c, the mRNA expression of GATA-3 correlated significantly with that of IL-4 (r(29) = 0.54, P = 0.0072), but was less correlated with another Th2 cytokine IL-13 (r(29) = 0.017, P = 0.9387), whereas, the expression of ROR-γτ was correlated with that of IL-17A (r(29) = 0.54, P = 0.0072). There was a tendency of correlation between the expression of T-bet and Th1 cytokine IFN-γ (r(29) = 0.22, P = 0.305) or IL-12B (r(29) = 0.33, P = 0.1265) and between the expression of Foxp3 and IL-10 (r(29) = 0.38, P = 0.0718).

Th1/Th2/Th17/Treg balance

To better understand the temporal trends in gene expression of Th1, Th2, Th17 and Treg in the liver post-infection, the ratios between the expression of transcription factors were evaluated using the pairwise comparison method (Fig. 2c). The RQ ratio of T-bet/GATA-3 was employed as an indicator of Th1/Th2 cytokine gene expression pattern. This ratio was higher in the liver at 98 dpi when compared to the uninfected group (P = 0.0049). This analysis also revealed that Th1/Th17 (P = 0.0246) and Treg/Th17 (P = 0.0126) were both upregulated at 98 dpi, when compared with the samples tested at other time points (3–70 dpi). However, no significant change was found in either Th1/Treg, Th2/Treg, or Th2/Th17 expression ratios at any of the examined time points. These data suggest a predominance of the Th1 and Treg cellular immune response at 98 dpi.


Our findings provide an increased understanding of the immuno-inflammatory response of buffaloes to F. gigantica infection through examining the correlation between gene expression of nine cytokines and four transcription factors in Th1, Th2, Th17 and Treg cells, together with the liver pathology during early, mid and late stages of infection. We initially characterized gross and histopathological changes in the liver at six time points following infection compared to uninfected control. Although infected buffaloes did not exhibit clinical signs, they developed hepatic gross pathologies similar to what have been observed in other studies [15, 19]. Histopathological features also agreed with previous studies on F. gigantica [9, 20, 21] and F. hepatica [22]. An unexpected finding was that none of the examined liver from infected animals showed calcium deposits. The difference in liver calcification between our study and previous reports may relate to variations in the experimental conditions, such as using different infectious doses, different parasite strains, the permissibility of the host species, the duration of the experimental infection, and using single infection vs repeated infections [21,22,23].

Next, we have shown that buffaloes showed a substantial liver immune response subsequent to oral challenge with 500 F. gigantica encysted metacercariae. Certain genes exhibited shifting temporal expression patterns as the infection proceeded. Expression of cytokine and transcription factor genes were decreased at 3 and 10 dpi, suggesting a general immunosuppression state to allow parasite colonization. Cytokine reduction during the early stage of F. gigantica infection has been reported previously [24, 25]. A previous study has documented the downregulation of MHC-II related genes, as well as the suppression of the host pro-inflammatory (Th1) immune response during early F. gigantica infection [16]. This immunosuppression might be attributed to an increased expression of the immunosuppressive cytokine TGF-β. While TGF-β was upregulated at 3 dpi, Foxp3 (Treg transcription factor) was upregulated at a later stage of infection (28 dpi), indicating that Tregs (also known as Foxp3-expressing cells) play a more diverse role than causing early immunosuppression in regulating immune responses. Although increased TGF-β at 3 dpi may have promoted the generation of Foxp3-expressing Tregs, increased IL-6 at 28 dpi may have abrogated the suppressive function of Tregs [26] by antagonizing TGF-β-induced Treg generation and stimulating Treg differentiation to Th17, as reported previously [27].

The inflammatory cascades triggered by F. gigantica must be tightly coordinated in order to avoid severe liver pathology. Th2 cells mediate humoral responses via induction of IL-4, IL-5, IL-9, IL-10 and IL-13 cytokines, which are important for controlling extracellular parasites [28]. In agreement with this, our results at 28 dpi showed a significant increase in the expression of the anti-inflammatory Th2-associated cytokines (IL-4, IL-13) and transcription factor (GATA-3), indicating the induction of Th2 response in liver tissue, probably to limit the fluke development and to balance the increased levels of inflammatory cytokines IL-1β, IL-6 and IL-17A. This response was expected because IL-4 is known to have the greatest effect in inducing Th2 differentiation [29] and the GATA-3 transcription factor is a key regulator of Th2 differentiation [30]. IL-4 stimulation may have activated STAT6, which is the major signal transducer in IL-4-mediated Th2 differentiation [31]. STAT6 is known to promote Th2 differentiation by stimulating GATA-3 [32]. Therefore, this coordinated response may be mediated by the parasite and/or the host in order to promote parasite survival while minimizing host hepatic damage. Th2-predominant response together with the suppression of Th1/Th17 response has been reported previously [33,34,35,36]. Interestingly, some molecules secreted by Fasciola species can suppress the differentiation of Th17 cells independently of Th2 cells differentiation, by altering the function of dendritic cells [37].

Although CD4+ Th17 cells play a role in the control of a variety of parasitic infections, the role of IL-17 produced by CD4+ Th17 cells in immunity to F. gignatica has not been clearly defined [38, 39]. Our current and previous studies [8], demostrated that IL-17 may play a role in the inflammatory process during early F. gigantica infection. Th17 cells are a newly-identified class of effector T cell, which produces IL-17A and IL-17F. Th17 cells can contribute to resistance to many intracellular [40] and extracellular parasites [8]. In addition to controlling infection, IL-17 expression has been associated with inflammatory and allergic responses [41]. TGF-β and IL-6 act cooperatively and non-redundantly to promote Th17 activity [27]. In agreement with this, TGF-β was upregulated at 3 dpi, and both IL-6 and IL-17A were upregulated at 28 dpi. IL-6 was found to be important in suppressing Treg generation in order to promote Th17 differentiation as discussed above. The expression of retinoic acid-related orphan receptors ROR-γt, a key transcription factor in Th17 differentiation [42, 43], was positively correlated with the RQ level of IL-17A at 28 dpi.

At 70 and 98 dpi a mixed Th1/Th2 type profile, supported by co-dominance of IFN-γ, IL-1β, IL-12B, IL-6, IL-4 and IL-13 cytokines, was observed as consistent with the result obtained by Kumar [44]. The high expression of Th1 cytokines (IL-1β and IL-6), suggesting acute inflammatory activity, correlated with hemorrhage, necrosis and severe infiltration of inflammatory cells in the liver. Th2 cytokines, on the other hand, often lead to chronic inflammatory response by promoting fibrosis and tissue remodeling [45, 46]. The accumulation of fibroblasts during chronic infection suggests that Th2 cytokine-mediated tissue repair was taking place. This state of immune homeostasis is probably required for parasite persistence and was correlated with an increased number of flukes. The T-bet/GATA-3 ratio, a marker of Th1/Th2 [47], was not altered during acute infection, but was significantly increased when the infection became chronic, suggesting more bias towards a Th1-mediated inflammatory response during late infection. This observation contradicts previous studies in mice, cattle and buffaloes infected with F. gigantica that reported a mixed Th1/Th2 response with a predominance of a Th2-biased pattern [48, 49]. Other studies have reported a Th2 response in early infection and an increased Th0-type response during chronic F. gigantica infection in cattle [35] and immunosuppression during chronic F. gigantica infection in buffaloes [8]. These differing findings between studies may result from different experimental conditions.

A synergism between IL-1β and IL-17 has previously been reported [50]. In our study, a strong correlation was observed between IL-17A and IL-4 and IL-10, as well as between ROR-γτ and GATA-3 and Foxp3, implying that induction of Th17 was paralleled by induction of both Th2 and Treg. This refutes the antagonistic relationship previously reported between a polarized Th2/Treg immune response and suppression of Th1/Th17 cytokines [11, 51,52,53]. A gradual upregulation of IL-10 mRNA expression was observed during early infection, in agreement with our previous finding in buffalo’s serum [8]. IL-10 levels significantly increased at 98 dpi, which was also reported in a study on F. hepatica-infected sheep [54]. IL-10 is considered to be a regulatory cytokine rather than a Th2 type cytokine [55, 56]. High levels of IL-10 can suppress excessive inflammation induced by the parasite [57, 58]. In our study, elevated levels of IL-10 and the concurrent upregulation of Th1 and Th2 cytokine genes in the liver suggest that IL-10 was required for the Th2 response to F. gigantica infection, instead of serving an anti-inflammatory function. Our data also showed that the RQ level of Foxp3 increased at 28 dpi. This level remained relatively high over the course of infection, which may have contributed to the Th1/Th2 balance in order to enhance parasite survival in the liver and to protect the host from over-inflammation [59]. Foxp3-expressing Treg cells are critical to maintaining immune homeostasis by minimizing tissue pathology, while modulating host immune response against helminth infections [60, 61]. The monitoring of Th17 and Treg cell offers a promising new perspective on the pathogenesis of F. gigantica infection and deserves further exploration.


This is the first study to characterize the correlation between the expression of Th1, Th2, Th17, and Treg cytokines and transcription factors with liver pathology in buffaloes, during the course of experimental F. gigantica infection. The expression patterns for the examined genes indicated that there were periods of differential regulation during F. gigantica infection, which may suggest either a mechanism of immune evasion based on modulation of transcription or a mechanism used by host tissue to limit the infection and tissue damage. Gene expression profiles showed a significant T-cell adaptive immune suppression between 3 and 10 dpi to facilitate parasite colonization and a more substantial transcript differential expression change (mixed Th1/Th2 immune response) from 28 to 70 dpi, which may contribute to the parasite survival while minimizing host tissue damage. During late infection (98 dpi), the response was biased towards Th1/Treg in order to limit the host’s Th2 protective response and promote chronic infection, which conincided with an increase in the number of recovered flukes and greater histopathological changes and fibrosis. Clearly, regulation of immune cytokine and transcription factor gene expression during F. gigantica infection is a complex mechanism involving a variety of parasite-specific and host-specific factors, perhaps depending on the parasite’s needs for survival during various phases of infection. Research into the mechanisms governing differential expression of these immune-related genes, may shed light on the actual role of these immune regulatory factors in the pathogenesis of F. gigantica infection in buffaloes. Much remains to be elucidated about the role of various types of inflammatory responses, and the cross-talk between types of T-helper responses, in fasciolosis.



Analysis of variance


Complementary deoxyribonucleic acid

Cq :

Quantification cycle


Degrees of freedom denominator


Degrees of freedom numerator


Deoxyribonucleic acid


Days post-infection


Enzyme linked immunosorbent assay


Encysted metacercariae


Forkhead box P3


Glyceraldethde-3-phosphate dehydrogenase


GATA binding protein


Hematoxylin and eosin






Internal transcribed spacer


Least significant difference


Messenger ribonucleic acid


Quantitative reverse transcription-polymerase chain reaction


Red blood cells


Ribosomal DNA


Ribonucleic acid


Retinoid-related orphan nuclear receptor


Relative quantitation


Standard error of the mean


T-box expressed in T cells


Transforming growth factor


T helper


Regulatory T cell


Weeks post-infection


  1. 1.

    Mas-Coma S, Funatsu I, Bargues M. Fasciola hepatica and lymnaeid snails occurring at very high altitude in South America. Parasitology. 2001;123:115–27.

    Article  Google Scholar 

  2. 2.

    Cywinska A. Epidemiology of fascioliasis in human endemic areas. J Helminthol. 2005;79:207–16.

    Article  Google Scholar 

  3. 3.

    Yadav S, Sharma R, Kalicharan A, Mehra U, Dass R, Verma A. Primary experimental infection of riverine buffaloes with Fasciola gigantica. Vet Parasitol. 1999;82:285–96.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Ahmedullah F, Akbor M, Haider M, Hossain M, Khan M, Hossain M, et al. Pathological investigation of liver of the slaughtered buffaloes in Barisal district. Ban J Vet Med. 2007;5:81–5.

    Google Scholar 

  5. 5.

    Raadsma H, Kingsford N, Spithill T, Piedrafita D. Host responses during experimental infection with Fasciola gigantica or Fasciola hepatica in merino sheep: I. Comparative immunological and plasma biochemical changes during early infection. Vet Parasitol. 2007;143:275–86.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Yuan W, Liu JM, Lu K, Li H, Duan MM, Feng JT, et al. Molecular identification and seasonal infections of species of Fasciola in ruminants from two provinces in China. J Helminthol. 2016;90:359–63.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Chai JY. Intestinal flukes. In: Murrell KD, Fried B, editors. Food-borne parasitic zoonoses. World class parasites, vol. 11. Boston: Springer; 2007. p. 53–115.

    Google Scholar 

  8. 8.

    Zhang FK, Guo AJ, Hou JL, Sun MM, Sheng ZA, Zhang XX, et al. Serum levels of cytokines in water buffaloes experimentally infected with Fasciola gigantica. Vet Parasitol. 2017;244:97–101.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Dalton JP (Editor). Fasciolosis. Wallingford: CABI Publishing; 1999.

  10. 10.

    Mas-Coma S, Bargues MD, Valero MA. Fascioliasis and other plant-borne trematode zoonoses. Int J Parasitol. 2005;35:1255–78.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    O'neill SM, Brady MT, Callanan JJ, Mulcahy G, Joyce P, Mills KH, et al. Fasciola hepatica infection downregulates Th1 responses in mice. Parasite Immunol. 2000;22:147–55.

    Article  PubMed  Google Scholar 

  12. 12.

    Zhang W, Moreau E, Hope J, Howard C, Huang W, Chauvin A. Fasciola hepatica and Fasciola gigantica: comparison of cellular response to experimental infection in sheep. Exp Parasitol. 2005;111:154–9.

  13. 13.

    Sallusto F, Lanzavecchia A. Heterogeneity of CD4+ memory T cells: functional modules for tailored immunity. Eur J Immunol. 2009;39:2076–82.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Infante-Duarte C, Kamradt T. Th1/Th2 balance in infection. Semin Immunopathol. 1999;21:317–38.

    CAS  Article  Google Scholar 

  15. 15.

    Zhang WY, Moreau E, Yang BZ, Li ZQ, Hope JC, Howard CJ, et al. Humoral and cellular immune responses to Fasciola gigantica experimental infection in buffaloes. Res Vet Sci. 2006;80:299–307.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Zhang FK, Zhang XX, Elsheikha HM, He JJ, Sheng ZA, Zheng WB, et al. Transcriptomic responses of water buffalo liver to infection with the digenetic fluke Fasciola gigantica. Parasit Vectors. 2017;10:56.

  17. 17.

    Huang WY, He B, Wang CR, Zhu XQ. Characterisation of Fasciola species from mainland China by ITS-2 ribosomal DNA sequence. Vet Parasitol. 2004;120:75–83.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008;3:1101–8.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Molina EC, Skerratt LF. Cellular and humoral responses in liver of cattle and buffaloes infected with a single dose of Fasciola gigantica. Vet Parasitol. 2005;131:157–63.

    Article  PubMed  Google Scholar 

  20. 20.

    Shaikh AA, Bilqees F, Khan MM. Bile duct hyperplasia and associated abnormalities in the buffaloes infected with Fasciola gigantica. Pak J Zool. 2004;36:231–8.

    Google Scholar 

  21. 21.

    Okaiyeto SO, Salami OS, Dnbirni SA, Allam L, Onoja II. Clinical, gross and histopathological changes associated with chronic fasciolosis infection in a dairy farm. J Vet Adv. 2012;2:444–8.

    Google Scholar 

  22. 22.

    Marcos LA, Yi P, Machicado A, Andrade R, Samalvides F, Sánchez J, et al. Hepatic fibrosis and Fasciola hepatica infection in cattle. J Helminthol. 2007;81:381–6.

    Article  PubMed  Google Scholar 

  23. 23.

    Grange D, Dhumeaux D, Couzineau P, Bismuth H, Bader J-P. Hepatic calcification due to Fasciola gigantica. Arch Surg. 1974;108:113–5.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Gironènes N, Valero MA, García-Bodelón MA, Chico-Calero I, Punzón C, Presno M, Mas-Coma S. Immune suppression in advanced chronic fascioliasis: an experimental study in a rat model. J Infect Dis. 2007;195:1504–12.

    Article  Google Scholar 

  25. 25.

    Rodríguez E, Noya V, Cervi L, Chiribao ML, Brossard N, Chiale C, et al. Glycans from Fasciola hepatica modulate the host immune response and TLR-induced maturation of dendritic cells. PLoS Negl Trop Dis. 2015;9:e0004234.

    Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Pasare C, Medzhitov R. Toll pathway-dependent blockade of CD4+ CD25+ T cell-mediated suppression by dendritic cells. Science. 2003;299:1033–6.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Mangan PR, Harrington LE, O'quinn DB, Helms WS, Bullard DC, Elson CO, et al. Transforming growth factor-β induces development of the TH17 lineage. Nature. 2006;441:231–4.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Sher A, Coffman RL. Regulation of immunity to parasites by T cells and T cell-derived cytokines. Annu Rev Immunol. 1992;10:385–409.

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Swain SL, Weinberg AD, English M, Huston G. IL-4 directs the development of Th2-like helper effectors. J Immunol. 1990;145:3796–806.

    CAS  PubMed  Google Scholar 

  30. 30.

    Zhu J, Min B, Hu-Li J, Watson CJ, Grinberg A, Wang Q, et al. Conditional deletion of Gata3 shows its essential function in TH1-TH2 responses. Nat Immunol. 2004;5:1157.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Takeda K, Kamanaka M, Tanaka T, Kishimoto T, Akira S. Impaired IL-13-mediated functions of macrophages in STAT6-deficient mice. J Immunol. 1996;157:3220–2.

    CAS  PubMed  Google Scholar 

  32. 32.

    Zhu J, Guo L, Watson CJ, Hu-Li J, Paul WE. Stat6 is necessary and sufficient for IL-4’s role in Th2 differentiation and cell expansion. J Immunol. 2001;166:7276–81.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Mulcahy G, Dalton JP, Cathepsin L. Proteinases as vaccines against infection with Fasciola hepatica (liver fluke) in ruminants. Res Vet Sci. 2001;70:83–6.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Waldvogel AS, Lepage MF, Zakher A, Reichel MP, Eicher R, Heussler VT. Expression of interleukin 4, interleukin 4 splice variants and interferon gamma mRNA in calves experimentally infected with Fasciola hepatica. Vet Immunol Immunopathol. 2004;97:53–63.

  35. 35.

    Ingale SL, Singh P, Raina OK, Mehra UR, Verma AK, Gupta SC, et al. Interferon-gamma and interleukin-4 expression during Fasciola gigantica primary infection in crossbred bovine calves as determined by real-time PCR. Vet Parasitol. 2008;15:158–61.

    Article  Google Scholar 

  36. 36.

    Ingale SL, Singh P, Raina OK, Verma AK, Channappanavar R, Mehra UR. Interleukin-2 and interleukin-10 gene expression in calves experimentally infected with Fasciola gigantica. Livest Sci. 2010;131:141–3.

    Article  Google Scholar 

  37. 37.

    Dowling DJ, Hamilton CM, Donnelly S, La Course J, Brophy PM, Dalton J, et al. Major secretory antigens of the helminth Fasciola hepatica activate a suppressive dendritic cell phenotype that attenuates Th17 cells but fails to activate Th2 immune responses. Infect Immun. 2010;78:793–801.

  38. 38.

    Noya V, Brossard N, Rodríguez E, Dergan-Dylon LS, Carmona C, Rabinovich GA, et al. A mucin-like peptide from Fasciola hepatica instructs dendritic cells with parasite specific Th1-polarizing activity. Sci Rep. 2017;7:40615.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Valero MA, Perez-Crespo I, Chillón-Marinas C, Khoubbane M, Quesada C, Reguera-Gomez M, et al. Fasciola hepatica reinfection potentiates a mixed Th1/Th2/Th17/Treg response and correlates with the clinical phenotypes of anemia. PLoS One. 2017;12:e0173456.

    Article  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Lieberman LA, Cardillo F, Owyang AM, Rennick DM, Cua DJ, Kastelein RA, et al. IL-23 provides a limited mechanism of resistance to acute toxoplasmosis in the absence of IL-12. J Immunol. 2004;173:1887–93.

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    Yen D, Cheung J, Scheerens H, Poulet F, McClanahan T, Mckenzie B, et al. IL-23 is essential for T cell–mediated colitis and promotes inflammation via IL-17 and IL-6. J Clin Invest. 2006;116:1310–6.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  42. 42.

    Yang XO, Pappu BP, Nurieva R, Akimzhanov A, Kang HS, Chung Y, et al. T helper 17 lineage differentiation is programmed by orphan nuclear receptors RORα and RORγ. Immunity. 2008;28:29–39.

    CAS  Article  PubMed  Google Scholar 

  43. 43.

    Wan YY, Regulatory T. Cells: immune suppression and beyond. Cell Mol Immunol. 2010;7:204–10.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Kumar N, Raina OK, Nagar G, Prakash V, Jacob SS. Th1 and Th2 cytokine gene expression in primary infection and vaccination against Fasciola gigantica in buffaloes by real-time PCR. Parasitol Res. 2013;112:3561–8.

    Article  PubMed  Google Scholar 

  45. 45.

    Allen JE, Sutherland TE. Host protective roles of type 2 immunity: parasite killing and tissue repair, flip sides of the same coin. Semin Immunol. 2014;26:329–40.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  46. 46.

    Cortés A, Muñoz-Antoli C, Esteban JG, Toledo R. Th2 and Th1 responses: clear and hidden sides of immunity against intestinal helminths. Trends Parasitol. 2017;33:678–93.

    Article  PubMed  Google Scholar 

  47. 47.

    Chakir H, Wang H, Lefebvre DE, Webb J, Scott FW. T-Bet/GATA-3 ratio as a measure of the Th1/Th2 cytokine profile in mixed cell populations: predominant role of GATA-3. J Immunol Methods. 2003;278:157–69.

    CAS  Article  PubMed  Google Scholar 

  48. 48.

    Chantree P, Phatsara M, Meemon K, Chaichanasak P, Changklungmoa N, Kueakhai P, et al. Vaccine potential of recombinant cathepsin B against Fasciola gigantica. Exp Parasitol. 2013;135:102–9.

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    Molina EC. Serum interferon-gamma and interleukins-6 and-8 during infection with Fasciola gigantica in cattle and buffaloes. J Vet Sci. 2005;6:135–9.

    PubMed  Google Scholar 

  50. 50.

    Ilarregui JM, van Beelen AJ, Fehres CM, Bruijns SC, García-Vallejo JJ, van Kooyk Y. New roles for CD14 and IL-β linking inflammatory dendritic cells to IL-17 production in memory CD4+ T cells. Immunol Cell Biol. 2016;94:907–16.

    CAS  Article  PubMed  Google Scholar 

  51. 51.

    Donnelly S, O'Neill SM, Sekiya M, Mulcahy G, Dalton JP. Thioredoxin peroxidase secreted by Fasciola hepatica induces the alternative activation of macrophages. Infect Immun. 2005;73:166–73.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  52. 52.

    Flynn RJ, Mannion C, Golden O, Hacariz O, Mulcahy G. Experimental Fasciola hepatica infection alters responses to tests used for diagnosis of bovine tuberculosis. Infect Immun. 2007;75:1373–81.

    CAS  Article  PubMed  Google Scholar 

  53. 53.

    Donnelly S, Stack CM, O'Neill SM, Sayed AA, Williams DL, Dalton JP. Helminth 2-Cys peroxiredoxin drives Th2 responses through a mechanism involving alternatively activated macrophages. FASEB J. 2008;22:4022–32.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  54. 54.

    Hacariz O, Sayers G, Flynn RJ, Lejeune A, Mulcahy G. IL-10 and TGF-β1 are associated with variations in fluke burdens following experimental fasciolosis in sheep. Parasite Immunol. 2009;31:613–22.

    CAS  Article  PubMed  Google Scholar 

  55. 55.

    Anderson CF, Oukka M, Kuchroo VJ, Sacks D. CD4+ CD25 Foxp3 Th1 cells are the source of IL-10–mediated immune suppression in chronic cutaneous leishmaniasis. J Exp Med. 2007;204:285–97.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  56. 56.

    Jankovic D, Kullberg MC, Feng CG, Goldszmid RS, Collazo CM, Wilson M, et al. Conventional T-bet+Foxp3 Th1 cells are the major source of host-protective regulatory IL-10 during intracellular protozoan infection. J Exp Med. 2007;204:273–83.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  57. 57.

    Mahanty S, Mollis SN, Ravichandran M, Abrams JS, Kumaraswami V, Jayaraman K, et al. High levels of spontaneous and parasite antigen-driven interleukin-10 production are associated with antigen-specific hyporesponsiveness in human lymphatic filariasis. J Infect Dis. 1996;173:769–72.

    CAS  Article  PubMed  Google Scholar 

  58. 58.

    Anuradha R, Munisankar S, Dolla C, Kumaran P, Nutman TB, Babu S. Parasite antigen-specific regulation of Th1, Th2, and Th17 responses in Strongyloides stercoralis infection. J Immunol. 2015;195:2241.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  59. 59.

    Scalzo K, Plebanski M, Apostolopoulos V. Regulatory T-Cells: immunomodulators in health and disease. Curr Top Med Chem. 2006;6:1759–68.

    CAS  Article  PubMed  Google Scholar 

  60. 60.

    Zhou S, Jin X, Chen X, Zhu J, Xu Z, Wang X, et al. Heat shock protein 60 in eggs specifically induces Tregs and reduces liver immunopathology in mice with schistosomiasis japonica. PLoS One. 2015;10:e0139133.

    Article  PubMed  PubMed Central  Google Scholar 

  61. 61.

    Escamilla A, Zafra R, Pérez J, McNeilly TN, Pacheco IL, Buffoni L, et al. Distribution of Foxp3+ T cells in the liver and hepatic lymph nodes of goats and sheep experimentally infected with Fasciola hepatica. Vet Parasitol. 2016;230:14.

    CAS  Article  PubMed  Google Scholar 

  62. 62.

    Castillo-Velázquez U, Aranday-Cortés E, Gutiérrez-Pabello JA. Alternative activation modifies macrophage resistance to Mycobacterium bovis. Vet Microbiol. 2011;151:51–9.

    Article  PubMed  Google Scholar 

  63. 63.

    Mingala CN, Konnai S, Venturina FA, Onuma M, Ohashi K. Quantification of water buffalo (Bubalus bubalis) cytokine expression in response to inactivated foot-and-mouth disease (FMD) vaccine. Res Vet Sci. 2009;87:213–7.

    CAS  Article  PubMed  Google Scholar 

  64. 64.

    Puech C, Dedieu L, Chantal I, Rodrigues V. Design and evaluation of a unique SYBR green real-time RT-PCR assay for quantification of five major cytokines in cattle, sheep and goats. BMC Vet Res. 2015;11:65.

    Article  PubMed  PubMed Central  Google Scholar 

  65. 65.

    Pleasance J, Wiedosari E, Raadsma HW, Meeusen E, Piedrafita D. Resistance to liver fluke infection in the natural sheep host is correlated with a type-1 cytokine response. Parasite Immunol. 2011;33:495–505.

    CAS  Article  PubMed  Google Scholar 

  66. 66.

    Maeda Y, Ohtsuka H, Tomioka M. Effect of progesterone on Th1/Th2/Th17 and regulatory T cell-related genes in peripheral blood mononuclear cells during pregnancy in cows. Vet Res Commun. 2013;37:43–9.

    Article  PubMed  Google Scholar 

Download references


We are grateful to Xiao-Xuan Zhang, Wen-Bin Zheng and Jian-Gang Ma from Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences, and Yang-Qun Kang, Xue-Fang Mei, Chang-Hong Hu, Yu Zhang, Yao-Yao Zhang, and Yi-Ying Liang from Guangxi University for their technical assistance.


Project financial support was provided by National Key Basic Research Program (973 Program) of China (Grant No. 2015CB150300) and by National Natural Science Foundation of China (Grant No. 31260605).

Availability of data and materials

The data supporting the findings of this article are included within the article.

Author information




WS, WYH and XQZ conceived and designed the study, and critically revised the paper. WS, ZYW, FKZ and KJL prepared the experimental samples. WS, ZYW and ZAS performed the experiments. WS, ZYW, FKZ and DYW analyzed the data. WS drafted the manuscript. HME helped in the data analysis and critical revision of the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Wei-Yi Huang or Xing-Quan Zhu.

Ethics declarations

Ethics approval

The study design was reviewed and approved by the Animal Ethics Committee of Guangxi University. Animals used in the study were handled in accordance with good animal practices as required by the Animal Ethics Procedures and Guidelines of the People’s Republic of China.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Shi, W., Wei, Z., Elsheikha, H.M. et al. Dynamic expression of cytokine and transcription factor genes during experimental Fasciola gigantica infection in buffaloes. Parasites Vectors 10, 602 (2017).

Download citation


  • Fasciola gigantica
  • Buffaloes
  • Liver
  • Cytokines
  • Transcription factors
  • Gene expression