Figures
Abstract
Fasciola hepatica is a parasitic trematode that infects a wide range of mammalian hosts, including livestock and humans, in temperate and tropical regions globally. This trematode causes the disease fascioliasis, which consists of an acute phase (≤ 12 weeks) during which juvenile parasites migrate through the host liver tissues, and a chronic phase (> 12 weeks) following the establishment of adult parasites in the liver bile ducts. Few studies have explored the progression of the host response over the course of Fasciola infection in the same animals. In this study, we characterized transcriptomic changes in peripheral blood mononuclear cells (PBMCs) collected from sheep at three time points over the first eight weeks of infection relative to uninfected controls. In total, 183 and 76 genes were found to be differentially transcribed at two and eight weeks post-infection respectively. Functional and pathway analysis of differentially transcribed genes revealed changes related to T-cell activation that may underpin a Th2-biased immune response against this parasite. This first insight into the dynamics of host responses during the early stages of infection improves the understanding of the pathogenesis of acute fascioliasis, informs vaccine development and presents a set of PBMC markers with diagnostic potential.
Citation: Alvarez Rojas CA, Scheerlinck J-P, Ansell BRE, Hall RS, Gasser RB, Jex AR (2016) Time-Course Study of the Transcriptome of Peripheral Blood Mononuclear Cells (PBMCs) from Sheep Infected with Fasciola hepatica. PLoS ONE 11(7): e0159194. https://doi.org/10.1371/journal.pone.0159194
Editor: Ulrike Gertrud Munderloh, University of Minnesota, UNITED STATES
Received: April 20, 2016; Accepted: June 28, 2016; Published: July 20, 2016
Copyright: © 2016 Alvarez Rojas et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: Data is accessible from the NCBI under BioProject ID: PRJNA327701.
Funding: This work was supported by the Australian Research Council and (ARC) grant number DP110103175 and the National Health and Medical Research Council (NHMRC) of Australia as well as by a Victorian Life Sciences Computation Initiative (VLSCI) grant number VR0007 on its Peak Computing Facility at the University of Melbourne, an initiative of the Victorian Government (to RB Gasser).
Competing interests: The authors have declared that no competing interests exist.
Introduction
Fasciola hepatica, commonly known as liver fluke, is responsible for acute and chronic liver disease (= fascioliasis) in various mammals worldwide, including humans. The disease has important consequences in infected livestock due to losses arising from reduced productivity and fertility, mortalities, failure to thrive and the condemnation of livers [1]. The infection establishes via oral ingestion of metacercariae, usually encysted on aquatic vegetation. Upon passage through the digestive tract, the juvenile stage of the parasite emerges in the small intestine, penetrates the gut wall and proceeds to migrate to and through the liver of the host animal [2].
During host invasion, F. hepatica secretes a potent mix of digestive proteases, allowing it to undertake migration through and feed on host tissues, particularly liver parenchyma, before reaching its final destination, the bile ducts. This developmental process takes approximately six to 12 weeks [2]. At sufficient infection intensity, symptoms of acute fascioliasis appear, especially in sheep, as severe abdominal pain, blood-loss and associated anaemia/hypoproteinaemia. Even at low infection intensities, livers of infected sheep may show signs of inflammation, haemorrhage and tissue necrosis. Once established in the bile duct, F. hepatica can survive as adults for several years, leading to a chronic immunological response, severe cholangitis and fibrosis of the bile ducts. At sufficient intensity, such chronic infections are often associated with blockage of the bile ducts and concomitant distension of the gall bladder, as well as atrophy of liver parenchyma, resulting from tissues unable to recover from the initial juvenile assault [2, 3].
Fasciola infection is associated with a non-inflammatory Th2-biased response [4–6], partially regulated by “parasite-driven” immunomodulation [7–10]. Previous investigations have followed humoral and cellular responses to Fasciola infection in sheep and cattle using histopathological and immunological tools [7, 11–16]. Other investigations have focused on the detection of specific interleukins or other immune proteins (i.e. IFN-γ, IL-4, IL-10 and TGF-β1) from PBMCs in response to excretory/secretory (ES) products from the parasite in sheep [7, 8], cattle [10, 17] and mice [10]. Recently, Rojas-Caraballo et al. [18] followed the progression of the host response to Fasciola infection using RNA microarray analysis; however, though informative, this work was conducted in mice, an experimental host. On other hand, investigations on the natural host have provide evidence of the involvement of Galectins interacting with the parasite [19, 20]. Despite this research, our understanding of how the host response changes over time, particularly with respect to acute and chronic stages of infection, remains limited.
In a previous study, we examined the transcriptomic responses of liver tissue in sheep against those experimentally infected with F. hepatica for 8 weeks, and defined the involvement of specific genes associated with the host’s metabolism, immune response and tissue repair/regeneration, highlighting an apparent overlapping function of many genes involved in these processes [21]. In the present study, we explore transcriptomic changes in PBMCs taken from the same sheep before infection, and 2 and 8 weeks following experimental infection, and compare these changes to uninfected controls. We aimed to provide a detailed analysis of the transcription of genes involved in the complex immune networks in PBMCs and to gain a better understanding of the progression of the immune response of sheep to this parasite which might assist efforts to develop vaccines against F. hepatica.
Material and Methods
Animals and experimental design
All animal experimental procedures were approved by the Ethics Committee of the Faculty of Veterinary and Agricultural Sciences, The University of Melbourne (Ethics application ID 1312938.1). Four female Merino sheep (< 6months of age) received a single oral individual dose of 180 metacercariae of F. hepatica purchased from Baldwin Aquatics Inc., Oregon, USA. Four additional female animals of the same age and breeding stock were maintained as controls and kept in a separate raised pen in the same facility. All animals were verified helminth free by faecal inspection prior to the start of the study. Individual sheep were clinically monitored on a weekly basis checking for general condition and alert status. No signs of illness nor mortality was observed prior to the experimental endpoint. All food, water and animal management practices were consistent for infected and control animals throughout the study. Eight weeks after infection, each animal was euthanized by intravenous injection of pentobarbitone sodium (172.5 mg/kg) and immediately necropsied.
PBMC isolation
Forty ml of blood were obtained weekly from the jugular vein of each animal (control and infected groups), starting from one day prior to the experimental infection. PBMCs were isolated before infection and 2 and 8 weeks post infection (WPI), as described previously [22]. Briefly, blood was collected in the presence of 100 μl of heparin (5000 units/ml) and centrifuged for 20 min at 600g, the buffy coat was collected and then diluted with 10 ml of phosphate-buffered saline (PBS) before being under-laid with 2 ml of Ficoll (1.077 specific gravity; Pharmacia Biotech, Sweden) and centrifuged for 30 min at 100g. The lymphocytes were collected and washed four times in PBS and stored in RNAlater® (Life Technologies, United States) following the manufacturer’s instructions.
RNA extraction and library preparation
PBMCs from each animal were collected at the above mentioned time-points and centrifuged at 600g for 5 min. An aliquot of the pelleted cells (200 μl/sample) was washed once in PBS, and TriPure reagent (Roche Diagnostics, Switzerland) was added. RNA was isolated following the manufacturer’s instructions. Purified RNA was resuspended in RNAse-free water (Life technologies). RNA quality was assessed using a Bioanalyzer 2100 (Agilent, USA); samples with an RNA quality indicator (RQI) of ≥8 were used for sequencing. Library preparation and RNA sequencing were performed as described previously [21].
Bioinformatic analysis
Paired-end RNA-seq reads for each sample/replicate were filtered for quality and adapter-trimmed using Trimmomatic software [23] (sliding window: 4 bp, leading and trailing: 3 bp, minimum read length: 100 bp; Phred quality: 25). Filtered reads were subjected to K-mer correction using the program corrector [24]. Data have been deposited in the Sequence Read Archive (SRA) from the National Center for Biotechnology Information (NCBI; SubmissionID: SUB1665337, BioProject ID: PRJNA327701). We used two strategies to define differentially transcribed genes (DTGs) in NOISeq. (1) First, because the immune response can vary significantly from individual to individual, we analyzed variability within the same animal two (T1) and eight weeks (T2) post infection, compared with the pre-infection timepoint (for control and infected animals) using NOISeq-sim (k = 0.5, q = 0.8). We filtered these DTGs and retained the genes that were differentially transcribed in at least 3 of 4 animals at each time point. Finally, we discarded any DTGs from infected animals found to be differentially transcribed over the same time-points in the uninfected controls. (2) In a second analysis, we treated infected and uninfected animals as biological replicates of the experimental and control groups, respectively. For this, we NOISeq-bio (k = 0.5, q = 0.9) comparing transcripts per million (TPM) values at 2 and 8 WPI in infected animals to the uninfected control animals. We kept for further analysis only those genes identified as being differentially transcribed at statistical significance using both analytical strategies.
Gene Ontology and KEGG analysis
All gene functional enrichment analyses relating to DTGs defined for infected animals were assessed based on the gene models provided in the Ovis aries genome assembly version Oar_v3.1 (GCA_000298735.1) hosted by Ensembl [25], and gene names and descriptions as they were annotated in Biomart application [26] as of 01/04/2016. Over-representation of GO terms in up- and down-regulated DTGs at each time-point and with respect to the control animals were identified using PANTHER [27] (adjusted P<0.05). DTGs were assigned to conserved biological pathways in the Kyoto Encyclopedia of Genes and Genomes (KEGG) using Kobas [28], using an expected value of 10−5. Pathway-enrichment analysis was conducted using the Fisher hyper-geometric mean tests to identify significantly enriched KEGG pathways and enzyme hierarchies for DTGs relative to all KEGG-annotated genes [29, 30].
Results and Discussion
Experimental infection was successful in all the four animals showing a similar response with an elevated number of eosinophils since the 4th week after infection (as previously shown in a previous study [21]. A detailed account of the RNA-seq data generated in this study can be found in S1 Table. In total more than 100 million reads were mapped to the Ovis aries genome; a mean of 58% of these reads mapped to recognized coding gene models. Because immune responses within individual animals are dependent in part on the genetics of the animal [31], we considered the possibility that it may not be appropriate to assess each animal as a biological replicate. We, therefore, adopted two strategies for our RNA-seq analyses; (Strategy1) in which we assessed each animal as an independent experiment using Noiseq-sim (which was developed for single replicate RNA-seq; see methods for details) and (Strategy 2) in which we treated infected and uninfected animals as biological replicates of an experimental versus control group using Noiseqbio.
Based on this approach, with strategy 1, we identified 216 differentially transcribed genes (DTGs; 149 up- and 67 down-regulated) at 2 WPI and 95 DTGs (78 up- and 17 down-regulated) in at least 3 of the 4 infected animals and 0 of the 4 control animals at contemporary time-points (S2 Table). With strategy 2, we identified 4,819 (2,187 up- and 2,632 down-regulated) and 2,537 (1,627 up and 910 down-regulated) DTGs at 2 and 8 WPI, respectively (S2 Table). In filtering these results for DTGs identified using both strategies 1 and 2, we identified a consensus set of 183 DTGs (135 up- and 48 down-regulated; S1 Fig) and 76 DTGs (63 up- and 13 down-regulated; S1 Fig) at 2 and 8 WPI respectively (Table 1). Among the consensus set of DTGs (henceforth referred to only as ‘DTGs’) at 2 and 8 WPI, we found 31 up-regulated at both time-points (Fig 1; i.e., continuously up-regulated over the course of the infection), with a single gene down-regulated at both time points (transcript ENSOART00000001087); this gene encodes the exocyst complex component 6 (EXOC6) involved in the docking of exocytic vesicles with fusion sites on the plasma membrane. Detailed information on the genes shown to be differentially transcribed, with their respective TPM values, can be found in S3 Table and S4 Table for up- and down- regulated genes at 2 WPI respectively and S5 Table and S6 Table for up- and down- regulated genes at 8 WPI, respectively.
Analysis was performed using the package NOISeq.
We explored the functional roles of the DTGs at 2 and 8 WPI. Table 2 shows overrepresented GO terms (biological processes) for up- and down-regulated genes at 2 and 8 WPI. Similarly, Table 3 shows the overrepresented GO terms (molecular function) for up- and down-regulated genes at 2 WPI. There was no overrepresentation of GO terms related to molecular function at 8 WPI. Fig 2 shows a schematic representation of significantly (p<0.05) enriched KEGG pathways in up-and down-regulated genes from ovine PBMCs at 2 and 8 WPI.
Gene transcription in PBMCs at 2 WPI
At 2 WPI, the juvenile flukes are between 1 to 2.5 mm in length, and they can be found penetrating the liver capsule and burrowing into the parenchyma, causing sub-capsular haemorrhage with eosinophil and macrophage infiltration [2, 32, 33].
Overall, from the total number of DTGs at 2 WPI, 115 genes were enriched in transcription at relative to other time-points and/or uninfected controls. We found no evidence of specific biological pathways up-regulated at 2 WPI that were not also up-regulated at 8 WPI; however, we did see a number of DTGs that were specifically up-regulated at this earlier stage of infection. “T-cell activation” was the most notable GO term (biological process) over-represented for DTGs at 2 WPI. A number of proteins have been identified to play a role in the T-cell activation in the literature [34]. Among these are the galectins, which are part of a family of animal lectins that bind N-acetyllactosamine-containing glycans. Galectins also have roles in pathogen recognition, shaping the course of adaptive immune responses and fine-tuning the inflammatory response [35]. In our study, two isoforms of galectin-9 (transcripts ENSOART00000017447 and ENSOART00000001365) were found to be up-regulated at 2 and 8 WPI. There is evidence galectins are part of the host-parasite interaction, for example promoting adhesion and changing mucus properties during parasite infection [20] or directly interacting with the parasite, as is the case of galectin-14 and galectin-11 with F. hepatica in the bile ducts [20]. We have previously demonstrated the up-regulation of the same genes in ovine liver parenchyma infected with F. hepatica [21]. An important role for galectins is their involvement in T-cell differentiation: galectin-9 together with TGF-beta1 have synergistic effects on the rate of conversion of inflammatory to regulatory T-cells in vitro inducing phosphorylation of SMAD2/3, ERK1/2 and formation of the SMAD2/3-SMAD4 complex [36]. The activation of SMAD3 has also been involved in the expression of Galectin-9 together with TGF-β [37]. Interestingly, we found Smad2 and tgf-β1 were up-regulated at 2 WPI. Smad2 and Smad3 are essential for the TGF-beta-mediated regulation of T-cell plasticity [38]. Elevated levels of peripheral blood mononuclear cell (PBMC)-derived transforming growth factor (TGF)-β1 have been observed in early phases of the infection with F. hepatica in cattle [17]. F. hepatica-induced TGF-β plays a critical role in bystander suppression of autoantigen-specific Th1 and Th17 responses that mediate autoimmune diseases [39].
Aside from galectins, other lectins also play an important role in the immune response [40, 41]. GO analysis revealed an over-representation of “positive regulation of cytokine production” and “innate immune response” involving (amongst others) three genes encoding C-type lectin receptors (CLRs; cleC6A, cleC4E and cleC3B). These receptors are crucial for tailoring immune responses to different pathogens, triggering distinct signalling pathways that induce the expression of specific cytokines which determine T-cell polarization [42, 43]. CLEC4E (also known as MINCLE) induces inflammatory cytokine production and drives the infiltration of neutrophils into damaged tissue [44]. This is very interesting since the migration of F. hepatica in liver produces extensive damage and subsequent cellular infiltration; our date suggests that CLEC4E may be an important driver of in the extensive liver damage and subsequent cellular infiltration associated with the juvenile migration of F. hepatica. Considering the comparatively moderate damage caused by Fasciola juveniles in cattle during the acute phase [2, 3], it would also be interesting to explore the transcriptional profile of the cleC4E (mincle) ortholog in cattle during the migration of juvenile flukes through liver tissue. A role for different C-type lectins receptors from macrophages (mannose receptor and Dectin-1,) has been described in response to F. hepatica antigens [45]. CLRs clearly merit further study with respect to the progression of the response to F. hepatica infection in sheep and other natural hosts.
Another important group of proteins with immune function that was shown to be up-regulated at 2 WPI are the S100 proteins, which are involved in “leukocyte migration” and the “innate immune response” (Table 2). The S100 protein family consists of 24 members comprising a group of proinflammatory molecules released by phagocytes [46, 47]. Six members of this family were up-regulated in PBMCs at 2 WPI (S100A4, A8, A9, A10, A11 and A12). S100A8, A9 and A12 S100A12 belong to the calgranulin subfamily and are known to be involved in host-responses to, for example, Echinococcus granulosus [48], Plasmodium falciparum[49], Leishmania major [50], Eimeria brunetti and E. acervulina [51] and F. hepatica [21]. Our analysis of GO terms for molecular function shows significant overrepresentation of “RAGE receptor binding” (Table 3). RAGE (i.e. receptor for advanced glycation end products) is a multi-ligand receptor of the immunoglobulin family that plays a role in the activation and differentiation of T-cells [52, 53] and has been shown to transduce the extracellular effects of S100A4, S100A11, S100A12 and other S100A proteins [54]. Further studies are necessary to elucidate the role of S100 proteins and RAGE in T-cell differentiation and other aspects of the host immune response to parasitic infection.
Also interesting among genes up-regulated specifically at 2 WPI is gata-1. GATA-1 expression suppresses the development of the Th1 response and transcription of associated genes (e.g., ifnɣ and cxcr3) in human T-cells [55]. Other up-regulated genes at 2 WPI, include four isoforms of cd300. CD300s are type I transmembrane glycoproteins with a single IgV-like extracellular domain [56]. They are crucial to remodelling host tissue and contribute to tuning, directing or terminating active immune responses. This is highly relevant in a disease, like fascioliasis, that involves the remodelling of the liver tissue, especially in the initial stages of the infection.
By comparison, Fasciola infection related to relatively few down-regulated genes generally. Among down-regulated genes at 2 WPI, we identified enrichment of genes associated with “T cell receptor signaling pathway”, “NTK cell cytotoxicity”, “B-cell receptor” and “Leukocyte transendothelial migration” pathways. For the T-cell receptor and NTK cytotoxicity signaling pathway the gene encoding the ZP70 was down-regulated, this is an essential kinase in initiating T-cell responses by the antigen receptor [57]. For the “B cell receptor pathway” the gene encoding CD22 was found down-regulated. CD22 is a member of the Siglec (sialic-acid-binding immunoglobulin-like lectin) family, this is an accessory co-receptor expressed on the B-cell surface that negatively modulates the B cell antigen receptors [58]. This suggests that this down-regulation of CD22 might result in a positive regulation of B cell antigens receptors in F. hepatica infection at 2 WPI. The infection with F. hepatica brings as a consequence an increase in transendothelial migration towards the damaged tissue in liver after migration of the liver fluke. Therefore a down-regulation in this mechanism should be elicited by the immune system [59]. It has been described that fatty acid from F. hepatica (Fh12) binds protein inhibiting TLR4 activation and suppresses the inflammatory cytokines [60]. Other components of the F. hepatica tegument also have anti-inflammatory functions in vitro suppressing dendritic cell maturation and function [61, 62].
Gene transcription in PBMCs at 8 WPI
At 8 WPI, migrating juvenile flukes should be about 10 mm in length and most will be entering the bile ducts, or migrating through the left liver lobe [2]. In the liver it is possible to find haemorrhagic areas, hepatic hyperplasia and fibrosis as part of the tissue repair process. Infiltrating lymphocytes and eosinophils will still be present at 8 WPI, suggesting that the parasites will still have close contact with immune cells, however this is more restricted compared with the situation at 2 WPI. By 8 WPI, relative to early-stages of infection, we found enrichment among up-regulated DTGs for biological processes including: “lipoxin metabolic process”, “leukotriene metabolic process” and “negative regulation of immune response” (Table 2). In the first two processes, the proteins ALOX5, ALOX5AP and ALOX15 play an important role in the production of lipoxins, which are members of a family of bioactive products generated from arachidonic acid [63]. Lipoxins have been described as playing a role in "immunomodulation" in Toxoplasma gondii infection [64] and suggested to influence T-lymphocyte effector functions in the setting of polarized T-helper cell responses (Th1 and Th2) [63]. In a previous study [21] we identified 572 gens up-regulated in liver tissue 8 WPI post infection with F. hepatica, forteen of these genes were identified to be up-regulated in PBMCs at this time point. None of the forty two down-regulated genes identified in liver at 8 WPI [21] were identified to be down-regulated in PBMCs at this time point in this study.
Gene transcription in PBMC is constitutively altered with Fasciola infection
In addition to DTGs found at 2 or at 8 WPI, we described a large number of genes differentially transcribed throughout (i.e. constitutively) F. hepatica infection in sheep. An overall appraisal of these DTGs showed KEGG and GO enrichment related to the immune system (Table 2), including “positive regulation of T-cell activation”, “regulation of interferon-gamma production”, “Leukocyte migration”, “positive regulation of cytokine production” and the “innate immune response”. Of particular interest are the up-regulation, at 2 and 8 WPI, of the same genes involved in processes including: “complement and coagulation cascades”, “chemokine” and the “TGF-beta signaling pathways” (Fig 2).
Three genes encoding proteins related with the complement system were up-regulated at 2 WPI: complement component 1, q subcomponent, B chain (c1qb) and q subcomponent, C chain (c1qc) and complement component 4 binding protein, alpha (c4bpa). These representatives of the complement system exert a variety of functions over immune cells. For example, complement activation fragments such as C4b interact with complement receptors and regulators modulating T-cell lineage and cytokine production [65]. C4BP also induces B-cell proliferation [66]. Complement inhibition is crucial for parasite survival within the host tissue or, were relevant, during blood-feeding [67]. Previous studies suggest that extracts from F. hepatica adults inhibit both the classical and alternative pathways of complement activation in normal bovine and human sera [68]. Similarly, a reduction in the classical pathway of complement-mediated lysis by normal rat sera has been described using extracts from metacercariae, eggs and E/S products [69]. Our results suggest that the parasite has an influence over the complement activation pathways and this deserves future investigations.
It is well known that the F. hepatica infection will attract numerous inflammatory cells to the liver [2]. Among the genes involved in these processes, ccl2 and ccl23 were differentially transcribed at both 2 and 8 WPI. CCL2, also known as monocyte chemoattractant protein-1 (MCP-1), is a member of the C-C chemokine family and a highly potent chemo-attractor of monocytes and CD41 T-cells. Monocytes/macrophages are found to be the major source of CCL2 in the peripheral blood [70]. Unlike the other C-C chemokines, CCL2 regulates Th-cell differentiation by polarizing Th0 cells toward a Th2 phenotype [71]. CCL2 is released locally in response to gastrointestinal nematode infection in mice and has been implicated in resistance to Trichuris muris by steering the host towards a Th2-type response [72]. CCL23 has been related to suppression of the production and release of neutrophils from the bone marrow [73]. Considering the important role that neutrophils play in the host response to many metazoan parasites [74–76], including F. hepatica [77, 78], it maybe that up-regulation of CCL23 is a consequence of an immunomodulatory influence by F. hepatica to block neutrophil production; this is the first suggestion of this kind of interaction with this parasite in sheep.
Constitutive changes in the metabolism of PBMCs in relation to “arachidonic acid” are also notable. At 2 WPI, we found up-regulation of arachidonate 5-lipoxygenase (alox5), arachidonate 5-lipoxygenase-activating protein (alox5ap) and alox15 in PBMCs. Interestingly, alox5ap and alox5 were also found to be up-regulated in liver tissue at 8 weeks after infection [21]. ALOX5 oxidizes arachidonic acid to produce 5-hydroperoxy-6,8,11,14-eicosatetraenoic acid, which is an intermediate for various leukotrienes of lipid mediators [79]. It has a pivotal role in the regulation of specific antibody responses by managing primary B-cells and Tfh (follicular) cells in vivo [80]. It has also suggested that 5-Lipoxygenase negatively regulates Th1 response during Brucella abortus infection in mice [81]; an up-regulation of 5-lipoxygenase-activating protein here appears to point to another pathway supporting Th1 suppression in F. hepatica infected sheep.
“Fc gamma R-mediated phagocytosis” was the only biological pathway reduced at both 2 and 8 WPI. This pathway is associated with phagocytosis of antibody-coated pathogens. Phagocytosis is commonly reduced with alternative activation of macrophages [82], which is a common hallmark of helminth-induced immunomodulation [83] and mediated by Fasciola hepatica, specifically, at least in part through the parasite secretion of thioredoxin peroxidase [84].
Surprisingly, a high number of genes encoding interleukins did not show any transcriptional activity at any of the time points examined in infected or control animals. This suggests that transcription of interleukins is reduced or non-existent in peripheral blood in animals and/or interleukins are only transcribed in cells that are directly interacting with the parasite. Interestingly, it has been recently described that tegumental antigens from F. hepatica induce unresponsive or anergic-like T cells in a mouse model [85], which might be related with our findings. In the case of interleukins considered to be typical of a Th2 response (IL4, IL5 and IL13) transcription levels were reported only at some time points and insufficiently to detect statistically-significant differences in transcription here. IL10 transcription levels were reported at all the time points studied, however no significant differences in transcription were found.
Conclusions
To our knowledge this is the first RNA-seq study to examine gene transcription in PBMCs from a natural host (sheep) infected with F. hepatica. We previously reported the exploration of the transcriptome of ovine liver tissue infected with the liver fluke [21]. The current study focused on the temporal progression of the host response to Fasciola, examining changes in the transcriptional profile of peripheral blood mononuclear cells (PBMCs) through the first 8 weeks of infection in the same animals used in our previous study. The results presented identify key genes involved in the immune processes activated by the presence of the parasite in PBMCs. Important biological events seem to be modified for the presence of the parasite, including the complement system, the chemokine signaling, T-cell activation and metabolic processes that influence the immune response. Important groups of genes appear to interact with, or at least be stimulated in response to the parasite including lectins, S100 and CD300. We also identify a number of genes that are differentially transcribed at the 2 time points examined suggesting that they can be used as markers of infection in PBMCs. Most changes in PBMC are largely uniform across the 8 weeks of infection. Some changes, e.g., positive regulation of T-cell activation and leukocyte migration, are specific to early infection. Others, relating more to lipoxin metabolic process and Fc gamma R-mediated phagocytosis appear directed towards later stage infection. The next steps in the investigation of host/parasite interaction between sheep and the liver fluke should involve the development of tools to measure the expression of proteins from differentially transcribed genes. These tools should contribute to a better knowledge of the interaction between host and parasite and also to a better understanding of the kind of immune response elicited by the parasite.
Supporting Information
S1 Fig. Venn diagrams showing consensus of found using both strategy 1 and 2 from material and methods.
https://doi.org/10.1371/journal.pone.0159194.s001
(TIF)
S1 Table. Number of RNA-seq reads from control and infected animals.
https://doi.org/10.1371/journal.pone.0159194.s002
(XLSX)
S2 Table.
A shows the differentially transcribed genes in each infected animal using NOISeq-sim comparing transcription levels at pre infection (T1) with 2 weeks post infection (T2). Also it shows the number of DTGs that were found in at least three animals and number of DTGs that were found uniquely in infected and not in control animals. B shows the results for the comparison between pre infection (T0) and 8 weeks post infection (T8). C shows the number of DTGs after comparing control and infected animals at transcription level at 2 weeks post infection (T1) with 8 weeks post infection (T2)
https://doi.org/10.1371/journal.pone.0159194.s003
(XLSX)
S3 Table. List of up-regulated genes and TPM (transcripts per million) values in PBMCs from sheep at 2 weeks post infection with Fasciola hepatica.
https://doi.org/10.1371/journal.pone.0159194.s004
(XLSX)
S4 Table. List of DOWN-regulated genes and TPM (transcripts per million) values in PBMCs from sheep at 2 weeks post infection with Fasciola hepatica.
https://doi.org/10.1371/journal.pone.0159194.s005
(XLSX)
S5 Table. List of up-regulated genes and TPM (transcripts per million) values in PBMCs from sheep at 8 weeks post infection with Fasciola hepatica.
https://doi.org/10.1371/journal.pone.0159194.s006
(XLSX)
S6 Table. List of DOWN-regulated genes and TPM (transcripts per millions) values in PBMCs from sheep at 8 weeks post infection with Fasciola hepatica.
https://doi.org/10.1371/journal.pone.0159194.s007
(XLSX)
Acknowledgments
This work was supported by the Australian Research Council and (ARC) grant number DP110103175 and the National Health and Medical Research Council (NHMRC) of Australia as well as by a Victorian Life Sciences Computation Initiative (VLSCI) grant number VR0007 on its Peak Computing Facility at the University of Melbourne, an initiative of the Victorian Government (to RB Gasser).
Author Contributions
Conceived and designed the experiments: ARJ JPS RBG. Performed the experiments: CAAR. Analyzed the data: CAAR ARJ BREA RSH. Contributed reagents/materials/analysis tools: JPS RBG RSH BREA. Wrote the paper: CAAR ARJ.
References
- 1. Schweizer G, Braun U, Deplazes P, Torgerson PR. Estimating the financial losses due to bovine fasciolosis in Switzerland. Vet Rec. 2005;157(7):188–93. pmid:16100368
- 2.
Behm CA. Pathology, pathophysiology and clinical aspects. In: Dalton JP, editor. Fasciolosis. Wallingford, UK: CAB International; 1999. p. 185–217.
- 3. Boray JC. Experimental fascioliasis in Australia. Adv Parasitol. 1969;7:95–210. Epub 1969/01/01. pmid:4935272.
- 4. O'Neill SM, Brady MT, Callanan JJ, Mulcahy G, Joyce P, Mills KH, et al. Fasciola hepatica infection downregulates Th1 responses in mice. Parasite immunology. 2000;22(3):147–55. Epub 2000/02/15. pmid:10672196.
- 5. Vukman KV, Adams PN, Metz M, Maurer M, O’Neill SM. Fasciola hepatica Tegumental Coat Impairs Mast Cells’ Ability To Drive Th1 Immune Responses. The Journal of Immunology. 2013;190(6):2873–9. pmid:23418624
- 6. Mendes EA, Mendes TA, dos Santos SL, Menezes-Souza D, Bartholomeu DC, Martins IV, et al. Expression of IL-4, IL-10 and IFN-gamma in the liver tissue of cattle that are naturally infected with Fasciola hepatica. Veterinary parasitology. 2013;195(1–2):177–82. Epub 2013/05/08. pmid:23648284.
- 7. Zhang WY, Moreau E, Hope JC, Howard CJ, Huang WY, Chauvin A. Fasciola hepatica and Fasciola gigantica: Comparison of cellular response to experimental infection in sheep. Exp Parasitol. 2005;111(3):154–9. http://dx.doi.org/10.1016/j.exppara.2005.06.005. pmid:16126202
- 8. Moreau E, Chauvin A, Boulard C. IFN gamma and IL-10 production by hepatic lymph node and peripheral blood lymphocytes in Fasciola hepatica infected sheep. Parasite (Paris, France). 1998;5(4):307–15. Epub 1999/01/08. pmid:9879553.
- 9. Prowse RK, Chaplin P, Robinson HC, Spithill TW. Fasciola hepatica cathepsin L suppresses sheep lymphocyte proliferation in vitro and modulates surface CD4 expression on human and ovine T cells. Parasite immunology. 2002;24(2):57–66. Epub 2002/03/05. pmid:11874560.
- 10. 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(1–2):53–63. Epub 2004/01/01. pmid:14700537.
- 11. Ferre I, Ortega-Mora LM, Rojo-Vázquez FA. Serum and bile antibody responses (IgG and IgA) during subclinical Fasciola hepatica infection in sheep. Veterinary parasitology. 1997;68(3):261–7. http://dx.doi.org/10.1016/S0304-4017(96)01070-9. pmid:9066071
- 12. Santiago N, Hillyer GV. Antibody Profiles by EITB and ELISA of Cattle and Sheep Infected with Fasciola hepatica. The Journal of Parasitology. 1988;74(5):810–8. pmid:3138397
- 13. Chauvin A, Moreau E, Boulard C. Responses of Fasciola hepatica infected sheep to various infection levels. Vet Res. 2001;32(1):87–92. pmid:11254181
- 14. Phiri IK, Phiri AM, Harrison LJS. Serum antibody isotype responses of Fasciola-infected sheep and cattle to excretory and secretory products of Fasciola species. Veterinary parasitology. 2006;141(3–4):234–42. http://dx.doi.org/10.1016/j.vetpar.2006.05.019. pmid:16797844
- 15. Raadsma HW, Kingsford NM, Suharyanta, Spithill TW, 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. Veterinary parasitology. 2007;143(3–4):275–86. http://dx.doi.org/10.1016/j.vetpar.2006.09.008. pmid:17069976
- 16. Chauvin A, Bouvet G, Boulard C. Humoral and cellular immune responses to Fasciola hepatica experimental primary and secondary infection in sheep. International Journal for Parasitology. 1995;25(10):1227–41. http://dx.doi.org/10.1016/0020-7519(95)00039-5. pmid:8557470
- 17. Flynn RJ, Mulcahy G. The roles of IL-10 and TGF-β in controlling IL-4 and IFN-γ production during experimental Fasciola hepatica infection. International Journal for Parasitology. 2008;38(14):1673–80. pmid:18597757
- 18. Rojas-Caraballo J, López-Abán J, Fernández-Soto P, Vicente B, Collía F, Muro A. Time-course investigation of the gene expression profile during Fasciola hepatica infection: A microarray-based study. Genomics Data. 2015;6:89–91. pmid:PMC4664735.
- 19.
Preston S, Dunphy J, Beddoe T, Meeusen E, Young A. Evaluation of the Role of Galectins in Parasite Immunity. In: Stowell RS, Cummings DR, editors. Galectins: Methods and Protocols. New York, NY: Springer New York; 2015. p. 371–95.
- 20. Young AR, Barcham GJ, Kemp JM, Dunphy JL, Nash A, Meeusen EN. Functional characterization of an eosinophil-specific galectin, ovine galectin-14. Glycoconjugate journal. 2009;26(4):423–32. Epub 2008/09/24. pmid:18810635.
- 21. Alvarez Rojas CA, Ansell BR, Hall RS, Gasser RB, Young ND, Jex AR, et al. Transcriptional analysis identifies key genes involved in metabolism, fibrosis/tissue repair and the immune response against Fasciola hepatica in sheep liver. Parasites & vectors. 2015;8(1):715. Epub 2015/03/18. pmid:25779253.
- 22. Scheerlinck JP, Casey G, McWaters P, Kelly J, Woollard D, Lightowlers MW, et al. The immune response to a DNA vaccine can be modulated by co-delivery of cytokine genes using a DNA prime-protein boost strategy. Vaccine. 2001;19(28–29):4053–60. Epub 2001/06/28. pmid:11427282.
- 23. Lohse M, Bolger AM, Nagel A, Fernie AR, Lunn JE, Stitt M, et al. RobiNA: a user-friendly, integrated software solution for RNA-Seq-based transcriptomics. Nucleic acids research. 2012;40(Web Server issue):W622–7. Epub 2012/06/12. pmid:22684630; PubMed Central PMCID: PMCPMC3394330.
- 24.
Brown C, Howe A, Zhang Q, Pyrkosz A, Brom T. A Reference-Free Algorithm for Computational Normalization of Shotgun Sequencing Data. Available at: http://arxiv.org/abs/1203.4802v2201210/11/2014. Available from: http://arxiv.org/abs/1203.4802v2.
- 25. Yates A, Akanni W, Amode MR, Barrell D, Billis K, Carvalho-Silva D, et al. Ensembl 2016. Nucleic acids research. 2016;44(D1):D710–D6. pmid:26687719
- 26. Yates A, Akanni W, Amode MR, Barrell D, Billis K, Carvalho-Silva D, et al. Ensembl 2016. Nucleic acids research. 2016;44(D1):D710–6. Epub 2015/12/22. pmid:26687719; PubMed Central PMCID: PMCPMC4702834.
- 27. Mi H, Poudel S, Muruganujan A, Casagrande JT, Thomas PD. PANTHER version 10: expanded protein families and functions, and analysis tools. Nucleic acids research. 2016;44(D1):D336–42. Epub 2015/11/19. pmid:26578592; PubMed Central PMCID: PMCPMC4702852.
- 28. Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, et al. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic acids research. 2011;39(suppl 2):W316–W22.
- 29. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic acids research. 2000;28(1):27–30. Epub 1999/12/11. pmid:10592173; PubMed Central PMCID: PMCPMC102409.
- 30. Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic acids research. 2012;40(Database issue):D109–14. Epub 2011/11/15. pmid:22080510; PubMed Central PMCID: PMCPMC3245020.
- 31. Ardia DR, Parmentier HK, Vogel LA. The role of constraints and limitation in driving individual variation in immune response. Functional Ecology. 2011;25(1):61–73.
- 32. Golbar HM, Izawa T, Juniantito V, Ichikawa C, Tanaka M, Kuwamura M, et al. Immunohistochemical characterization of macrophages and myofibroblasts in fibrotic liver lesions due to Fasciola infection in cattle. The Journal of veterinary medical science / the Japanese Society of Veterinary Science. 2013;75(7):857–65. Epub 2013/02/16. pmid:23411482.
- 33. Chauvin A, Boulard C. Local immune response to experimental Fasciola hepatica infection in sheep. Parasite (Paris, France). 1996;3(3):209–15.
- 34. Malissen B, Gregoire C, Malissen M, Roncagalli R. Integrative biology of T cell activation. Nature immunology. 2014;15(9):790–7. 10.1038/ni.2959. pmid:25137453
- 35. Rabinovich GA, Toscano MA. Turning 'sweet' on immunity: galectin-glycan interactions in immune tolerance and inflammation. Nat Rev Immunol. 2009;9(5):338–52. pmid:19365409
- 36. Lv K, Zhang Y, Zhang M, Zhong M, Suo Q. Galectin-9 promotes TGF-beta1-dependent induction of regulatory T cells via the TGF-beta/Smad signaling pathway. Molecular medicine reports. 2013;7(1):205–10. Epub 2012/10/16. pmid:23064791.
- 37. Wu C, Thalhamer T, Franca Rafael F, Xiao S, Wang C, Hotta C, et al. Galectin-9-CD44 Interaction Enhances Stability and Function of Adaptive Regulatory T Cells. Immunity. 2014;41(2):270–82. pmid:25065622
- 38. Takimoto T, Wakabayashi Y, Sekiya T, Inoue N, Morita R, Ichiyama K, et al. Smad2 and Smad3 are redundantly essential for the TGF-beta-mediated regulation of regulatory T plasticity and Th1 development. J Immunol. 2010;185(2):842–55. Epub 2010/06/16. pmid:20548029.
- 39. Walsh KP, Brady MT, Finlay CM, Boon L, Mills KHG. Infection with a Helminth Parasite Attenuates Autoimmunity through TGF-β-Mediated Suppression of Th17 and Th1 Responses. The Journal of Immunology. 2009;183(3):1577–86. pmid:19587018
- 40. Drickamer K, Taylor ME. Recent insights into structures and functions of C-type lectins in the immune system. Current opinion in structural biology. 2015;34:26–34. Epub 2015/07/15. pmid:26163333; PubMed Central PMCID: PMCPMC4681411.
- 41. Rabinovich GA, van Kooyk Y, Cobb BA. Glycobiology of immune responses. Annals of the New York Academy of Sciences. 2012;1253:1–15. Epub 2012/04/25. pmid:22524422; PubMed Central PMCID: PMCPMC3884643.
- 42. Geijtenbeek TBH, Gringhuis SI. Signalling through C-type lectin receptors: shaping immune responses. Nat Rev Immunol. 2009;9(7):465–79. pmid:19521399
- 43. Richardson MB, Williams SJ. MCL and mincle: C-type lectin receptors that sense damaged self and pathogen associated molecular patterns. Frontiers in Immunology. 2014;5.
- 44. Yamasaki S, Ishikawa E, Sakuma M, Hara H, Ogata K, Saito T. Mincle is an ITAM-coupled activating receptor that senses damaged cells. Nature immunology. 2008;9(10):1179–88. Epub 2008/09/09. pmid:18776906.
- 45. Guasconi L, Serradell MC, Garro AP, Iacobelli L, Masih DT. C-type lectins on macrophages participate in the immunomodulatory response to Fasciola hepatica products. Immunology. 2011;133(3):386–96. pmid:21595685
- 46. Donato R, Cannon BR, Sorci G, Riuzzi F, Hsu K, Weber DJ, et al. Functions of S100 Proteins. Current molecular medicine. 2013;13(1):24–57. pmid:PMC3707951.
- 47. Roth J, Vogl T, Sorg C, Sunderkötter C. Phagocyte-specific S100 proteins: a novel group of proinflammatory molecules. Trends in Immunology. 24(4):155–8. pmid:12697438
- 48. Basika T, Muñoz N, Casaravilla C, Irigoín F, Batthyány C, Bonilla M, et al. Phagocyte-specific S100 proteins in the local response to the Echinococcus granulosus larva. Parasitology. 2012;139(02):271–83.
- 49. Lee H-W, Kim T-S, Kang Y-J, Kim J-Y, Lee S, Lee W-J, et al. Up-regulated S100 calcium binding protein A8 in Plasmodium-infected patients correlates with CD4+CD25+Foxp3 regulatory T cell generation. Malaria Journal. 2015;14(1):1–11.
- 50. Guerfali FZ, Laouini D, Guizani-Tabbane L, Ottones F, Ben-Aissa K, Benkahla A, et al. Simultaneous gene expression profiling in human macrophages infected with Leishmania major parasites using SAGE. BMC Genomics. 2008;9:238. Epub 2008/05/23. pmid:18495030; PubMed Central PMCID: PMCPMC2430024.
- 51. Lopez-Bernad F, del Cacho E, Gallego M, Quilez J, Sanchez-Acedo C. Immunohistochemical study of S100-like protein in Eimeria brunetti and Eimeria acervulina. Veterinary parasitology. 2000;88(1–2):1–6. Epub 2000/02/19. pmid:10681017.
- 52. Chen Y, Akirav EM, Chen W, Henegariu O, Moser B, Desai D, et al. RAGE ligation affects T cell activation and controls T cell differentiation. J Immunol. 2008;181(6):4272–8. Epub 2008/09/05. pmid:18768885; PubMed Central PMCID: PMCPMC2643976.
- 53. Akirav EM, Preston-Hurlburt P, Garyu J, Henegariu O, Clynes R, Schmidt AM, et al. RAGE expression in human T cells: a link between environmental factors and adaptive immune responses. Plos One. 2012;7(4):e34698. Epub 2012/04/18. pmid:22509345; PubMed Central PMCID: PMCPMC3324532.
- 54. Donato R. RAGE: a single receptor for several ligands and different cellular responses: the case of certain S100 proteins. Current molecular medicine. 2007;7(8):711–24. Epub 2008/03/12. pmid:18331229.
- 55. Sundrud MS, Vancompernolle SE, Eger KA, Bruno TC, Subramaniam A, Mummidi S, et al. Transcription factor GATA-1 potently represses the expression of the HIV-1 coreceptor CCR5 in human T cells and dendritic cells. Blood. 2005;106(10):3440–8. Epub 2005/08/11. pmid:16091457; PubMed Central PMCID: PMC1895046.
- 56. Clark GJ, Ju X, Tate C, Hart DNJ. The CD300 family of molecules are evolutionarily significant regulators of leukocyte functions. Trends in Immunology. 2009;30(5):209–17. pmid:19359216
- 57. Wang H, Kadlecek TA, Au-Yeung BB, Goodfellow HES, Hsu L-Y, Freedman TS, et al. ZAP-70: An Essential Kinase in T-cell Signaling. Cold Spring Harbor Perspectives in Biology. 2010;2(5):a002279. pmid:PMC2857167.
- 58. Nitschke L. The role of CD22 and other inhibitory co-receptors in B-cell activation. Current Opinion in Immunology. 2005;17(3):290–7. http://dx.doi.org/10.1016/j.coi.2005.03.005. pmid:15886119
- 59. Schnoor M. Endothelial actin-binding proteins and actin dynamics in leukocyte transendothelial migration. J Immunol. 2015;194(8):3535–41. Epub 2015/04/08. pmid:25848070.
- 60. Martin I, Caban-Hernandez K, Figueroa-Santiago O, Espino AM. Fasciola hepatica fatty acid binding protein inhibits TLR4 activation and suppresses the inflammatory cytokines induced by lipopolysaccharide in vitro and in vivo. J Immunol. 2015;194(8):3924–36. Epub 2015/03/18. pmid:25780044; PubMed Central PMCID: PMCPMC4390499.
- 61. Hamilton CM, Dowling DJ, Loscher CE, Morphew RM, Brophy PM, O'Neill SM. The Fasciola hepatica tegumental antigen suppresses dendritic cell maturation and function. Infect Immun. 2009;77(6):2488–98. Epub 2009/04/01. pmid:19332532; PubMed Central PMCID: PMCPMC2687350.
- 62. 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(2):793–801. Epub 2009/11/18. pmid:19917714; PubMed Central PMCID: PMCPMC2812196.
- 63. Parkinson JF. Lipoxin and synthetic lipoxin analogs: an overview of anti-inflammatory functions and new concepts in immunomodulation. Inflammation & allergy drug targets. 2006;5(2):91–106. Epub 2006/04/15. pmid:16613568.
- 64. Machado FS, Aliberti J. Role of lipoxin in the modulation of immune response during infection. International immunopharmacology. 2008;8(10):1316–9. Epub 2008/08/09. pmid:18687293.
- 65. Kemper C, Atkinson JP. T-cell regulation: with complements from innate immunity. Nat Rev Immunol. 2007;7(1):9–18. pmid:17170757
- 66. Brodeur SR, Angelini F, Bacharier LB, Blom AM, Mizoguchi E, Fujiwara H, et al. C4b-binding protein (C4BP) activates B cells through the CD40 receptor. Immunity. 2003;18(6):837–48. Epub 2003/06/24. pmid:12818164.
- 67. Schroeder H, Skelly PJ, Zipfel PF, Losson B, Vanderplasschen A. Subversion of complement by hematophagous parasites. Developmental & Comparative Immunology. 2009;33(1):5–13. http://dx.doi.org/10.1016/j.dci.2008.07.010.
- 68. Montgomery TD, Leid RW, Wescott RB. Interaction of bovine complement with Fasciola hepatica. Veterinary parasitology. 1986;19(1–2):55–65. Epub 1986/01/01. pmid:3962163.
- 69. Baeza E, Poitou I, Villejoubert C, Boulard C. Complement depletion in rats infected with Fasciola hepatica: in vivo and in vitro studies. Veterinary parasitology. 1994;51(3–4):219–30. Epub 1994/02/01. pmid:8171824.
- 70. Ansari AW, Bhatnagar N, Dittrich-Breiholz O, Kracht M, Schmidt RE, Heiken H. Host chemokine (C-C motif) ligand-2 (CCL2) is differentially regulated in HIV type 1 (HIV-1)-infected individuals. International Immunology. 2006;18(10):1443–51. pmid:16916890
- 71. Gu L, Tseng S, Horner RM, Tam C, Loda M, Rollins BJ. Control of TH2 polarization by the chemokine monocyte chemoattractant protein-1. Nature. 2000;404(6776):407–11. Epub 2000/04/04. pmid:10746730.
- 72. deSchoolmeester ML, Little MC, Rollins BJ, Else KJ. Absence of CC chemokine ligand 2 results in an altered Th1/Th2 cytokine balance and failure to expel Trichuris muris infection. J Immunol. 2003;170(9):4693–700. Epub 2003/04/23. pmid:12707348.
- 73. Shih CH, van Eeden SF, Goto Y, Hogg JC. CCL23/myeloid progenitor inhibitory factor-1 inhibits production and release of polymorphonuclear leukocytes and monocytes from the bone marrow. Experimental hematology. 2005;33(10):1101–8. Epub 2005/10/13. pmid:16219532.
- 74. Hermosilla C, Caro TM, Silva LM, Ruiz A, Taubert A. The intriguing host innate immune response: novel anti-parasitic defence by neutrophil extracellular traps. Parasitology. 2014;141(11):1489–98. Epub 2014/04/12. pmid:24721985.
- 75. Chuah C, Jones MK, Burke ML, McManus DP, Owen HC, Gobert GN. Defining a pro-inflammatory neutrophil phenotype in response to schistosome eggs. Cellular microbiology. 2014;16(11):1666–77. Epub 2014/06/06. pmid:24898449.
- 76. Chen F, Wu W, Millman A, Craft JF, Chen E, Patel N, et al. Neutrophils prime a long-lived effector macrophage phenotype that mediates accelerated helminth expulsion. Nature immunology. 2014;15(10):938–46. pmid:PMC4479254.
- 77. Jefferies JR, Turner RJ, Barrett J. Effect of Fasciola hepatica excretory-secretory products on the metabolic burst of sheep and human neutrophils. International Journal for Parasitology. 1997;27(9):1025–9. http://dx.doi.org/10.1016/S0020-7519(97)00067-2. pmid:9363484
- 78. Jefferies JR, Corbett E, Barrett J, Turner RJ. Polarization and chemokinesis of ovine and human neutrophils in response to Fasciola hepatica excretory-secretory products. Int J Parasitol. 1996;26(4):409–14. Epub 1996/04/01. pmid:8773528.
- 79. Radmark O, Werz O, Steinhilber D, Samuelsson B. 5-Lipoxygenase: regulation of expression and enzyme activity. Trends Biochem Sci. 2007;32(7):332–41. Epub 2007/06/20. pmid:17576065.
- 80. Nagashima T, Ichimiya S, Kikuchi T, Saito Y, Matsumiya H, Ara S, et al. Arachidonate 5-Lipoxygenase Establishes Adaptive Humoral Immunity by Controlling Primary B Cells and Their Cognate T-Cell Help. The American Journal of Pathology. 2011;178(1):222–32. pmid:PMC3070564.
- 81. Fahel JS, de Souza MB, Gomes MT, Corsetti PP, Carvalho NB, Marinho FA, et al. 5-Lipoxygenase negatively regulates Th1 response during Brucella abortus infection in mice. Infect Immun. 2015;83(3):1210–6. Epub 2015/01/15. pmid:25583526; PubMed Central PMCID: PMCPMC4333460.
- 82. Varin A, Mukhopadhyay S, Herbein G, Gordon S. Alternative activation of macrophages by IL-4 impairs phagocytosis of pathogens but potentiates microbial-induced signalling and cytokine secretion. Blood. 2010;115(2):353–62. Epub 2009/11/03. pmid:19880493; PubMed Central PMCID: PMCPMC2808158.
- 83. Hewitson JP, Grainger JR, Maizels RM. Helminth immunoregulation: The role of parasite secreted proteins in modulating host immunity. Molecular and Biochemical Parasitology. 2009;167(1):1–11. pmid:19406170
- 84. 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(1):166–73. Epub 2004/12/25. pmid:15618151; PubMed Central PMCID: PMCPMC538930.
- 85. Aldridge A, O'Neill SM. Fasciola hepatica tegumental antigens induce anergic-like T cells via dendritic cells in a mannose receptor-dependent manner. Eur J Immunol. 2016;46(5):1180–92. Epub 2016/03/05. 10.1002/eji.201545905. 26931640. pmid:26931640