- Open Access
The Fasciola hepatica genome: gene duplication and polymorphism reveals adaptation to the host environment and the capacity for rapid evolution
© Cwiklinski et al.; licensee BioMed Central. 2015
- Received: 10 November 2014
- Accepted: 13 March 2015
- Published: 3 April 2015
The liver fluke Fasciola hepatica is a major pathogen of livestock worldwide, causing huge economic losses to agriculture, as well as 2.4 million human infections annually.
Here we provide a draft genome for F. hepatica, which we find to be among the largest known pathogen genomes at 1.3 Gb. This size cannot be explained by genome duplication or expansion of a single repeat element, and remains a paradox given the burden it may impose on egg production necessary to transmit infection. Despite the potential for inbreeding by facultative self-fertilisation, substantial levels of polymorphism were found, which highlights the evolutionary potential for rapid adaptation to changes in host availability, climate change or to drug or vaccine interventions. Non-synonymous polymorphisms were elevated in genes shared with parasitic taxa, which may be particularly relevant for the ability of the parasite to adapt to a broad range of definitive mammalian and intermediate molluscan hosts. Large-scale transcriptional changes, particularly within expanded protease and tubulin families, were found as the parasite migrated from the gut, across the peritoneum and through the liver to mature in the bile ducts. We identify novel members of anti-oxidant and detoxification pathways and defined their differential expression through infection, which may explain the stage-specific efficacy of different anthelmintic drugs.
The genome analysis described here provides new insights into the evolution of this important pathogen, its adaptation to the host environment and external selection pressures. This analysis also provides a platform for research into novel drugs and vaccines.
- Schistosoma Mansoni
- Adult Parasite
- Echinococcus Multilocularis
- Adult Fluke
- Hepatica Infection
The digenean trematode Fasciola hepatica is one of the most important pathogens of domestic livestock and has a global distribution [1-4]. The disease, fasciolosis, results in huge losses to the agricultural industry associated with poor food conversion, lower weight gains, impaired fertility and reduced milk (cattle) and wool (sheep) production. Heavy, acute infections can result in death, particularly in sheep and goats. Economic losses attributable to F. hepatica infection have been estimated at more than US$3 billion per annum worldwide [5,6], although even this estimate may be conservative as F. hepatica infection modulates its host’s immune response and its ability to resist or eliminate common microbial pathogens [7,8]. Fasciolosis is also an important zoonosis in regions where agricultural management practices are less advanced, particularly in South America and North Africa [3,9]. It is estimated that between 2.4 and 17 million people are infected with this liver fluke worldwide, with a further 91 million people living at risk, resulting in fasciolosis being included on the World Health Organization list of major neglected tropical diseases [1-3].
Here we provide a genome assembly for F. hepatica and assess genome-wide polymorphism and transcriptional profiles in order to identify key features of its genome that underlie its ability to migrate through different physiological environments, to parasitise different host species, and to respond rapidly to external selection pressures.
A large genome with high gene polymorphism
Fasciola hepatica assembly statistics
204 Kbp (REAPR 155 Kbpa)
Number of scaffolds ≥3 Kbp
Number of scaffolds ≥1 Kbp
Contig N50 (≥100 bp)
9.7 K bp
Number of contigs (≥100 bp)
Total assembly length
Total length of gaps
Number of RNAseq-supported gene models
Mean number of exons/gene
Mean exon size (95% range)
303 bp (36 bp – 1,369 bp)
Mean intron size (95% range)
3.7 Kbp (33 bp - 17.5 Kbp)
Proportion CEGMA core eukaryotic genes found
We investigated levels of polymorphism among F. hepatica genes by re-sequencing the genomes of individual fluke from each of five isolates, all from the UK. Substantial polymorphism among isolates was observed; 48% of genes exhibited at least one non-synonymous SNP and the level of non-synonymous nucleotide diversity, pi, averaged across 21.8 Mbp of coding sequence, was 5.2 × 10-4 (that is, two randomly sampled sequences differed approximately every 1,900 bp). By comparison, this figure is higher than in humans , similar to most vertebrates  and, on limited data, smaller than some parasitic nematode populations . Although F. hepatica is a self-fertilising hermaphrodite, and so has the potential to inbreed and lose genetic diversity, our data show that F. hepatica populations, as a whole, harbour substantial genetic variation. A likely explanation is that parasite populations are typically large, often larger than that of their hosts, which greatly slows any enhanced effects of genetic drift caused by self-fertilisation .
Expression patterns from multi-gene families reveal important developmental host-parasite interactions
To explore novel processes associated with developmental changes, genes were clustered based on expression pattern (Figure 3b). Several biological processes were greatly downregulated in maturing and adult parasites relative to metacercariae (Cluster 12, Figure 3), including molecules involved in cell adhesion (cadherins, integrins) and cytoskeletal proteins (talins) that could play an important role in sensing changes in the physiological environment and rapidly initiating excystment following ingestion by the definitive host. Clusters showing strongest patterns of differential expression (both up- and downregulation) were markedly over-represented by peptidases and terms associated with regulation of peptidase activity (Clusters 2, 8, 10 and 11, Figure 3). Genes associated with F. hepatica structure, particularly protein polymerisation and microtubule based movement, such as tubulin, dynein and surface tegumental genes, were highly upregulated in immature and adult fluke relative to earlier stages (Clusters 1 and 3, Figure 3b and Additional file 9: Table S7, Additional file 10: Table S8 and Additional file 11: Table S9). Our data show that, across the transcriptome as a whole, the strongest changes in expression were observed among different members of the multigenic protease and tubulin families, as detailed below.
Unlike the FhCL peptidases, the FhCB family consists of a single clade of seven members. Nevertheless, these peptidases were also temporally regulated (Figure 4c) such that three members (FhCB1, FhCB2 and FhCB3) exhibited parallel expression patterns to FhCL3 and were thus highly expressed in the NEJs and downregulated as the parasites matured. These data suggest a concerted role for FhCLs and FhCBs in the early infection stage. Also, the constitution and specific expression of a family of peptidases, asparaginyl endopeptidases or legumains that are responsible for the processing of the inactive FhCL and FhCB zymogens to functional enzymes  suggest specific and important developmental roles for these peptidases (Figure 4d). Thus, legumain 1, which was the most highly expressed member in NEJs could be required for activation of the FhCL3 and FhCB 1/2/3 peptides at the time the parasite emerges from its encysted stage in the intestine and initiates infection. Furthermore, legumain 3, which was switched on late in parasite development, is the prospective candidate for activating the FhCL1, FhCL2 and FhCL5 in mature blood-feeding adult parasites.
Anti-oxidant and detoxification systems
We investigated the evolution of anti-oxidant systems in F. hepatica, which are essential for adaptation to the host environment. Thus, as the parasite rapidly develops and enters different aerobic/anaerobic environments, anti-oxidant systems are critical not only for the detoxification of reactive oxygen and nitrogen (ROS, RNS) generated by endogenous cellular metabolism but also as a frontline defence against superoxide and nitric oxide radicals released by host immune effector cells such as macrophages, eosinophils and neutrophils [43,44]. Parasitic platyhelminths express genes encoding superoxide dismutase (SOD) which dismutates superoxide radicals to H2O2 but they do not possess catalase, the enzyme responsible for converting H2O2 to water and oxygen (although a catalase gene is present in the genome of free-living flatworms, such as Schmidtea mediterranea). A catalase gene was not found in the F. hepatica genome, which is consistent with other parasitic platyhelminths, where the function of peroxide detoxification has been supplanted by the newly discovered thiol-dependent peroxiredoxin and its reducing partner thioredoxin . The presence of the gene encoding the recently described thioredoxin glutathione reductase (TGR; Additional file 14: Table S11), together with the absence of distinct thioredoxin reductase and glutathione reductase genes verified that TGR is the sole reductive enzyme that links the thioredoxin-dependent and glutathione-dependent anti-oxidant defence systems in this parasite as observed in other platyhelminths [43,46,47]. The pivotal position of TGR between these two essential anti-oxidant systems makes it a promising target for the development of new anti-trematode drugs . Included in the repertoire of F. hepatica anti-oxidants are genes encoding SOD, glutathione transferases (GSTs) and three fatty acid binding proteins (FABP; Additional file 14: Table S11). With the exception of the GSTs, we found that in F. hepatica each of these components is encoded by a single gene, which is in stark contrast to the expanded anti-oxidant gene families in the closely related trematodes, O. viverrini and C. sinensis [13,15]. Peroxiredoxin, GST and FABP are secreted by F. hepatica via non-classical secretory pathways, perhaps as cargoes of exosomes , into the host circulation where they influence host immune responses by recruiting and activating M2 macrophages [7,48,49] and suppressing dendritic cell activity . These immunomodulatory effects contribute to the establishment of an immune suppressive environment, which aids parasite survival and the development of chronic disease.
The F. hepatica genome is one of the largest pathogen genomes sequenced to date but we found no evidence of genome duplication or repeat expansion to explain this. Why this large genome should have evolved is unclear, especially given that its replication may be energetically costly or slow cell division . For a parasite, such as Fasciola, that relies on the production of large numbers of eggs to facilitate transmission, one would expect strong selection against the accumulation of junk DNA if a large genome imposed a cost on egg production. It is possible that much of the non-coding portion of the F. hepatica genome is involved in gene regulation , and the size of F. hepatica’s genome may be related to its complex life cycle and variety of developmental morphs. If so, however, it would be difficult to explain why the F. hepatica genome is around three time the size of the Schistosoma genome , which has a similar life cycle. The large genome of F. hepatica therefore remains a paradox for which we may have to wait for comparative genome sequencing, currently underway , across other platyhelminth taxa to provide an answer.
The ability of F. hepatica to infect and survive in different tissue environments as it migrates from the intestine, through the liver and into the bile ducts is underpinned by gene duplication. Thus, our results show that, across the whole transcriptome, the strongest patterns of differential expression were observed among members of protease and tubulin gene families, and, in the case of proteases, these can be associated with changes in the active site and substrate specificity. While gene duplication appears to be a key process of adaptation to the parasitic life-style used by many helminths, it is notable that different helminth taxa have arrived at different evolutionary solutions. Comparison between F. hepatica and the bile-dwelling liver flukes C. sinensis and O. viverrini, shows the expansion of the anti-oxidant SOD families  and cathepsin F protease families [13,61] in C. sinensis and O. viverrini but not in F. hepatica, and conversely the expansion of the cathepsin L family in the latter species only, which suggests differences in how these related parasites tackle life within the same environment. Differences in the specificity or developmental expression of detoxification genes, such as members of the ABC and MATE families, or of the tubulin gene family, may be important in understanding why F. hepatica responds so markedly different to drugs at different stages of development compared to other digeneans. The exclusive activity of TCBZ against F. hepatica, and the potency of praziquantel to all these other digeneans except F. hepatica, points to a uniqueness in this parasite that needs to be resolved .
Fasciola hepatica is a highly adaptable parasite, evidenced by its ability to infect novel hosts, and it is notable that our results reveal high levels of polymorphism in genes specific to parasitic digeneans. Diversity within F. hepatica populations at genes important for the host-parasite interface may underpin a high evolutionary potential for F. hepatica to respond to changes in host availability or to other selection pressures. Similarly, the broad geographic range of F. hepatica would suggest that it is able to adapt to different climatic conditions and, in this respect, F. hepatica may also be able to respond to exploit changes in climate in temperate regions, such as the UK, where warmer and wetter winters favour transmission and increased prevalence . F. hepatica is seen to rapidly develop drug resistance to TCBZ  and the standing genetic diversity that we find across the genome suggests that it harbours the potential to evolve resistance to any novel drug treatment, which compounds the difficultly of controlling F. hepatica given the shortage of drugs effective against the juvenile stages. Nevertheless, the availability of a genome for F. hepatica that we provide, plus the characterisation of early molecular events in infection should help support the development of novel drugs and vaccines, particularly against the migrating juveniles that cause much of the pathology. For example, RNAi-mediated knockdown of cathepsin L and B expression has been shown to prevent NEJs from crossing the intestine  and vaccine trials with cathepsin L1 have provided partial protection against infection and pathology . The exploitation of a broader range of targets within the F. hepatica genome is now a priority given the widespread prevalence of resistance to TCBZ within F. hepatica populations. The commercial importance for agriculture to develop a replacement treatment to TCBZ may stimulate new treatments that could be translated to other important digenean parasites, including those of humans.
Source of parasite material
Adult parasites from each of five isolates were used for genome sequencing: (1) FhepLivSP, from the laboratory maintained Shrewsbury isolate (Ridgeway Research, UK); (2) FhepLivS1, a clonal line derived from the Shrewsbury population; (3) and (4) FhepLivR1 and R3, clonal lines derived from two isolates recovered from sheep in Northern England naturally infected with F. hepatica; and (5) FhepLivR2, a clonal line derived from a F. hepatica population from naturally infected sheep in South West Wales. For RNA sequencing, the following were used: (1) metacercariae and newly excysted juveniles (NEJ) at 1, 3 and 24 h post excystment from a North American isolate (Baldwin Aquatics Inc., Monmouth, OR, USA). Twenty-one–day-old juvenile flukes were recovered from mice infected with the same isolate; (2) an adult parasite recovered from the bile ducts of cattle naturally infected with F. hepatica in Uruguay. All animal work was conducted with ethical approval from the Universities of Liverpool (UK) and McGill (Canada).
Approximately 10 μg of DNA from a single, adult fluke taken from the FhepLivS1 strain was used to prepare Illumina TruSeq fragment libraries and 26 Gbp of 2 × 250 bp reads generated on an Illumina MiSeq (mean insert sizes 470 bp and 580 bp). Further individuals of FhepLivS1 were used to prepare Nextera Mate-Pair libraries (3 Kbp and 10 Kbp) and approximately 60 m 2× 100 bp Illumina reads from each library were generated. Approximately 200 μg of DNA prepared from several fluke of the FhepLivSP isolate was used to construct 2 Kbp, 5 Kbp and 8 Kbp mate-pair libraries, which were sequenced either on an Illumina GAII or HiSeq 2000. For each of the parasite isolates FhepLivSP, FhepLivS1, FhepLivR1, FhepLivR2 and FhepLivR3, DNA was isolated from a single adult; a fragment library was prepared and sequenced on a single lane of an Illumina HiSeq 2000 to yield approximately 24 Gbp of sequence (Centre for Genomic Research, Liverpool, UK). For RNA sequencing, Illumina TruSeq RNA libraries were prepared from biological replicates of metacercariae (3 replicates), NEJ 1 h (2 replicates), NEJ 3 h (2 replicates), NEJ 24 h (2 replicates), Juveniles 21 days (1 replicate) and Adult (1 replicate) (Genome Quebec, Montreal, Canada).
Assembly and annotation
Illumina MiSeq reads were trimmed to Q ≥30 and adaptors removed using Sickle and Perl and assembled using Newbler (Roche GS-Assembler v2.6) with flags set for large genome and a heterozygote sample. Mate-pair reads were first mapped to these contigs using Bowtie2  to remove duplicates and wrongly orientated reads, and scaffolded into contigs using SSPACE . Gap filling was achieved using GapFiller for 2× 250 bp and 2× 100 bp paired-end reads and run for three iterations (available as ENA accessions LN627018-LN647175). RNAseq data were mapped to scaffolds within the assembled genome greater than 3 Kbp using TopHat2 to identify transcribed regions and splice junctions. These, together with RNAseq data assembled using Trinity and S. mansoni protein sequence, were passed to the MAKER pipeline  to predict genes. Repeatmasker, Windowmasker and Dustmasker were used to identify repetitive regions. CEGMA v2.4 , which searches for 248 highly conserved genes, was used to assess the completeness of genome with the settings for vertebrates to allow long introns. REAPR  was used to assess the quality of scaffolding within the assembly and to produce an alternative, more conservative assembly by splitting scaffolds at locations with lower support (available as ENA accessions LN736597-LN774150). Homologs of F. hepatica predicted protein sequences were identified within UniProt using BLAST and functional domains identified using InterPro. InParanoid and MultiParanoid  were used to identify ortholog clusters from Schistosoma mansoni (v3.1.16), Clonorchis sinensis (v3.7), Schmidtea mediterranea and Echinococcus multilocularis (v29042013) predicted proteins [11-15,70].
Gene expression analysis
RNAseq libraries were mapped to MAKER gene models using TopHat2  and read counts extracted using htseq-count. Genes with a sum of at least five reads across all libraries were analysed for differential expression in edgeR  using a negative binomial model of successive developmental stages relative to metacercariae and with tagwise dispersion estimated from all samples. Hierarchical clustering, based on model coefficients, was used to group differentially expressed genes by similarity of expression into 12 clusters and hypergeometric tests used to test for over-representation of gene ontology terms within each cluster relative to the whole gene-set.
Genetic diversity analysis
Reads for each isolate were mapped to the genome using Bowtie2 and resulting bam files passed together to the GATK pipeline  for local realignment and SNP calling. SNP calls within genes were filtered by score >100, FS <60 and combined coverage between 10 and 250. Because gene duplicates collapsed within the assembly might give erroneously high diversity for some genes, genes were excluded with a median coverage greater than 213 and by heuristic scoring of SNPs appearing as heterozygotes in all samples. IGV was used to manually assess the success of filtering parameters. Levels of nucleotide diversity for different classes of SNPs were calculated for all genes. To identify the most polymorphic genes, a generalised linear model with a Poisson distribution was fitted to the number of non-synonymous SNPs as the response variable versus number of non-synonymous sites as a quadratic function. This accounts for the distribution of SNP counts within a gene, the fact that longer genes have the potential to have more SNPs and the possibility that genes of different lengths may evolve at different rates. Genes were classified according to their conservation across platyhelminths by the presence of orthologs retained across increasing taxonomic scales. This was preferred to assessing conservation on the basis of non-synonymous versus synonymous substitution ratios between species, since synonymous sites were saturated by multiple substitutions at such broad taxonomic scales. The robustness of the generalised linear model was tested by randomly sampling equal proportions of genes from each quartile of the length distribution for each class of conserved gene, to ensure that no bias was introduced if the length of a gene was correlated to its level of conservation. From the generalised linear model, the most polymorphic genes were identified having residuals in the top 1% quantile and hypergeometric tests were used to test for over-representation of GO terms within these highly polymorphic genes relative to the gene-set as a whole.
Discovery and characterisation of gene families
Discovery of F. hepatica gene families was carried out using BLAST analysis (NCBI v2.2.27 and v2.2.29), with available published gene sequences of interest (Additional file 16: Table S13), followed by manual annotation. Comparative analysis was carried out using the closely related trematode genome sequence datasets: Clonorchis sinensis and Schistosoma mansoni. Sequence alignment and phylogenetic analysis was carried out using Clustal Omega  and MEGA5 , respectively.
Data are freely available from WormBase ParaSite and the European Nucleotide Archive under accessions LN627018-LN647175 (assembly data), PRJEB6687 (genomic read data) and PRJEB6904 (transcriptomic read data).
We are grateful to Katherine Allen for producing clonal lines of liver fluke and Catherine Hartley, University of Liverpool and Paula Martin and Oliver Gladstone, Ridgeway Research, UK for technical assistance in producing parasite material. We are grateful for sequencing support provided by Margaret Hughes and Suzanne Kay at the Centre for Genomic Research, University of Liverpool, and Mathieu Bourgey, Genome Quebec, Canada. This study was funded by grants from the Biotechnology and Biological Sciences Research Council, UK (BB1002480/1) to JH, SP, DW and JLC, and a Ministère Économie, Innovation et Exportation (MEIE), Québec, award and a European Research Council Advanced Grant (HELIVAC, 322725) to JPD.
- Keiser J, Utzinger J. Food-borne trematodiases. Clin Microbiol Rev. 2009;22:466–83.View ArticlePubMed CentralPubMedGoogle Scholar
- Fürst T, Keiser J, Utzinger J. Global burden of human food-borne trematodiasis: a systematic review and meta-analysis. Lancet Infect Dis. 2012;12:210–21.View ArticlePubMedGoogle Scholar
- Robinson MW, Dalton JP. Zoonotic helminth infections with particular emphasis on fasciolosis and other trematodiases. Philos Trans R Soc Lond B Biol Sci. 2009;364:2763–76.View ArticlePubMed CentralPubMedGoogle Scholar
- Videnova K, Mackay DKJ. Availability of vaccines against major animal diseases in the European Union. Rev Sci Tech. 2012;31:971–8.PubMedGoogle Scholar
- Spithill TW, Carmona C, Piedrafita D, Smooker PM. Prospects for immunoprophylaxis against Fasciola hepatica (liver fluke). In: Targets, Screens, Drugs and Vaccines. Weinheim: Wiley-VCH Verlag GmbH & Co. KGaA; 2012. p. 465–84.View ArticleGoogle Scholar
- Piedraffita D, Spithill TW, Smith RE, Raadsma HW. Improving animal and human health through understanding liver fluke immunology. Parasite Immunol. 2010;32:572–81.Google Scholar
- Dalton JP, Robinson MW, Mulcahy G, O’Neill SM, Donnelly S. Immunomodulatory molecules of Fasciola hepatica: Candidates for both vaccine and immunotherapeutic development. Vet Parasitol. 2013;195:272–85.View ArticlePubMedGoogle Scholar
- Claridge J, Diggle P, McCann CM, Mulcahy G, Flynn R, McNair J, et al. Fasciola hepatica is associated with the failure to detect bovine tuberculosis in dairy cattle. Nat Commun. 2012;3:853.View ArticlePubMed CentralPubMedGoogle Scholar
- García HH, Moro PL, Schantz PM. Zoonotic helminth infections of humans: echinococcosis, cysticercosis and fascioliasis. Curr Opin Infect Dis. 2007;20:489–94.View ArticlePubMedGoogle Scholar
- Mas-Coma S, Bargues M-D, Valero MA. Fascioliasis and other plant-borne trematode zoonoses. Int J Parasitol. 2005;35:1255–78.View ArticlePubMedGoogle Scholar
- Berriman M, Haas BJ, LoVerde PT, Wilson RA, Dillon GP, Cerqueira GC, et al. The genome of the blood fluke Schistosoma mansoni. Nature. 2009;460:352–65.View ArticlePubMed CentralPubMedGoogle Scholar
- Protasio AV, Tsai IJ, Babbage A, Nichol S, Hunt M, Aslett MA, et al. A systematically improved high quality genome and transcriptome of the human blood fluke Schistosoma mansoni. PLoS Negl Trop Dis. 2012;6:e1455.View ArticlePubMed CentralPubMedGoogle Scholar
- Young ND, Nagarajan N, Lin SJ, Korhonen PK, Jex AR, Hall RS, et al. The Opisthorchis viverrini genome provides insights into life in the bile duct. Nat Commun. 2014;5:e4378.Google Scholar
- Zhou Y, Zheng H, Chen Y, Zhang L, Wang K, Guo J, et al. The Schistosoma japonicum genome reveals features of host-parasite interplay. Nature. 2009;460:345–51.View ArticlePubMed CentralGoogle Scholar
- Huang Y, Chen W, Wang X, Liu H, Chen Y, Guo L, et al. The carcinogenic liver fluke, Clonorchis sinensis: new assembly, reannotation and analysis of the genome and characterization of tissue transcriptomes. PLoS One. 2013;8:e54732.View ArticlePubMed CentralPubMedGoogle Scholar
- Young ND, Jex AR, Li B, Liu S, Yang L, Xiong Z, et al. Whole-genome sequence of Schistosoma haematobium. Nat Genet. 2012;44:221–5.View ArticlePubMedGoogle Scholar
- Sanderson AR. Maturation and probable gynogenesis in the liver fluke, Fasciola hepatica L. Nature. 1953;172:110–2.View ArticlePubMedGoogle Scholar
- Hirai H, Hirai Y, LoVerde PT. Evolution of sex chromosomes ZW of Schistosoma mansoni inferred from chromosome paint and BAC mapping analyses. Parasitol Int. 2012;61:684–9.View ArticlePubMedGoogle Scholar
- Komalamisra C. Chromosomes and C-banding of Opisthorchis viverrini. Southeast Asian J Trop Med Public Health. 1999;30:576–9.PubMedGoogle Scholar
- Park G-M, Im K-I, Huh S, Yong T-S. Chromosomes of the liver fluke, Clonorchis sinensis. Korean J Parasitol. 2000;38:201.View ArticlePubMed CentralPubMedGoogle Scholar
- Adelson DL, Raison JM, Edgar RC. Characterization and distribution of retrotransposons and simple sequence repeats in the bovine genome. Proc Natl Acad Sci. 2009;106:12855–60.View ArticlePubMed CentralPubMedGoogle Scholar
- Li W-H, Sadler LA. Low nucleotide diversity in man. Genetics. 1991;129:513–23.PubMed CentralPubMedGoogle Scholar
- Gayral P, Melo-Ferreira J, Glémin S, Bierne N, Carneiro M, Nabholz B, et al. Reference-free population genomics from next-generation transcriptome data and the vertebrate–invertebrate gap. PLoS Genet. 2013;9:e1003457.View ArticlePubMed CentralPubMedGoogle Scholar
- Skuce P, Stenhouse L, Jackson F, Hypša V, Gilleard J. Benzimidazole resistance allele haplotype diversity in United Kingdom isolates of Teladorsagia circumcincta supports a hypothesis of multiple origins of resistance by recurrent mutation. Int J Parasitol. 2010;40:1247–55.View ArticlePubMedGoogle Scholar
- Anderson TJC, Romero-Abel ME, Jaenike J. Genetic structure and epidemiology of Ascaris populations: patterns of host affiliation in Guatemala. Parasitol. 1993;107:319–34.View ArticleGoogle Scholar
- Scheiffele P. Cell-cell signaling during synapse formation in the CNS. Annu Rev Neurosci. 2003;26:485–508.View ArticlePubMedGoogle Scholar
- Zito K, Fetter RD, Goodman CS, Isacoff EY. Synaptic clustering of Fascilin II and Shaker: essential targeting sequences and role of Dlg. Neuron. 1997;19:1007–16.View ArticlePubMedGoogle Scholar
- Tuttle AM, Hoffman TL, Schilling TF. Rabconnectin-3a regulates vesicle endocytosis and canonical wnt signaling in zebrafish neural crest migration. Plos Biol. 2014;12:e1001852.View ArticlePubMed CentralPubMedGoogle Scholar
- Kalbe M, Haberl B, Haas W. Snail host finding by Fasciola hepatica and Trichobilharzia ocellata: compound analysis of “miracidia-attracting glycoproteins”. Exp Parasitol. 2000;96:231–42.View ArticlePubMedGoogle Scholar
- Andrews SJ. The life cycle of Fasciola hepatica. In: Dalton JP, editor. Fasciolosis. Wallingford: CABI; 1999.Google Scholar
- McVeigh P, Atkinson L, Marks NJ, Mousley A, Dalzell JJ, Sluder A, et al. Parasite neuropeptide biology: Seeding rational drug target selection? Int J Parasitol Drugs Drug Resist. 2012;2:76–91.View ArticlePubMed CentralPubMedGoogle Scholar
- Tielens AGM. Metabolism. In: Dalton JP, editor. Fasciolosis. Wallingford: CABI; 1999.Google Scholar
- McVeigh P, Maule AG, Dalton JP, Robinson MW. Fasciola hepatica virulence-associated cysteine peptidases: a systems biology perspective. Microbes Infect. 2012;14:301–10.View ArticlePubMedGoogle Scholar
- Robinson MW, Tort JF, Lowther J, Donnelly SM, Wong E, Xu W, et al. Proteomics and phylogenetic analysis of the cathepsin l protease family of the helminth pathogen Fasciola hepatica: expansion of a repertoire of virulence-associated factors. Mol Cell Proteomics. 2008;7:1111–23.View ArticlePubMedGoogle Scholar
- Turk V, Stoka V, Vasiljeva O, Renko M, Sun T, Turk B, et al. Cysteine cathepsins: from structure, function and regulation to new frontiers. Biochim Biophys Acta. 1824;2012:68–88.Google Scholar
- Corvo I, O’Donoghue AJ, Pastro L, Pi-Denis N, Eroy-Reveles A, Roche L, et al. Dissecting the active site of the collagenolytic cathepsin L3 protease of the invasive stage of Fasciola hepatica. PLoS Negl Trop Dis. 2013;7:e2269.View ArticlePubMed CentralPubMedGoogle Scholar
- Robinson MW, Menon R, Donnelly SM, Dalton JP, Ranganathan S. An integrated transcriptomics and proteomics analysis of the secretome of the helminth pathogen Fasciola hepatica: proteins associated with invasion and infection of the mammalian host. Mol Cell Proteomics. 2009;8:1891–907.View ArticlePubMed CentralPubMedGoogle Scholar
- Dalton JP, Brindley PJ, Donnelly S, Robinson MW. The enigmatic asparaginyl endopeptidase of helminth parasites. Trends Parasitol. 2009;25:59–61.View ArticlePubMedGoogle Scholar
- Ryan LA, Hoey E, Trudgett A, Fairweather I, Fuchs M, Robinson MW, et al. Fasciola hepatica expresses multiple α- and β-tubulin isotypes. Mol Biochem Parasitol. 2008;159:73–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Fuchs MA, Ryan LA, Chambers EL, Moore CM, Fairweather I, Trudgett A, et al. Differential expression of liver fluke β-tubulin isotypes at selected life cycle stages. Int J Parasitol. 2013;43:1133–9.View ArticlePubMedGoogle Scholar
- Devine C, Brennan GP, Lanusse CE, Alvarez LI, Trudgett A, Hoey E, et al. Inhibition of triclabendazole metabolism in vitro by ketoconazole increases disruption to the tegument of a triclabendazole-resistant isolate of Fasciola hepatica. Parasitol Res. 2011;109:981–95.View ArticlePubMedGoogle Scholar
- Sanabria R, Ceballos L, Moreno L, Romero J, Lanusse C, Alvarez L. Identification of a field isolate of Fasciola hepatica resistant to albendazole and susceptible to triclabendazole. Vet Parasitol. 2013;193:105–10.View ArticlePubMedGoogle Scholar
- Williams DL, Bonilla M, Gladyshev VN, Salinas G. Thioredoxin glutathione reductase-dependent redox networks in platyhelminth parasites. Antioxid Redox Signal. 2013;19:735–45.View ArticlePubMed CentralPubMedGoogle Scholar
- Lu J, Holmgren A. The thioredoxin antioxidant system. Free Radic Biol Med. 2014;66:75–87.View ArticlePubMedGoogle Scholar
- McGonigle L, Mousley A, Marks NJ, Brennan GP, Dalton JP, Spithill TW, et al. The silencing of cysteine proteases in Fasciola hepatica newly excysted juveniles using RNA interference reduces gut penetration. Int J Parasitol. 2008;38:149–55.View ArticlePubMedGoogle Scholar
- Prast-Nielsen S, Huang H-H, Williams DL. Thioredoxin glutathione reductase: Its role in redox biology and potential as a target for drugs against neglected diseases. Biochim Biophys Acta. 1810;2011:1262–71.Google Scholar
- Maggioli G, Silveira F, Martín-Alonso JM, Salinas G, Carmona C, Parra F. A recombinant thioredoxin-glutathione reductase from Fasciola hepatica induces a protective response in rabbits. Exp Parasitol. 2011;129:323–30.View ArticlePubMedGoogle Scholar
- Robinson MW, Hutchinson WF, Dalton JP, Donnelly S. Peroxiredoxin: a central player in immune modulation. Parasite Immunol. 2010;32:305–13.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- 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.View ArticlePubMed CentralPubMedGoogle Scholar
- Wilkinson R, Law CJ, Hoey EM, Fairweather I, Brennan GP, Trudgett A. An amino acid substitution in Fasciola hepatica P-glycoprotein from triclabendazole-resistant and triclabendazole-susceptible populations. Mol Biochem Parasitol. 2012;186:69–72.View ArticlePubMedGoogle Scholar
- Brennan GP, Fairweather I, Trudgett A, Hoey E, McCoy, McConville M, et al. Understanding triclabendazole resistance. Exp Mol Pathol. 2007;82:104–9.Google Scholar
- Virkel G, Lifschitz A, Sallovitz J, Ballent M, Scarcella S, Lanusse C. Inhibition of cytochrome P450 activity enhances the systemic availability of triclabendazole metabolites in sheep. J Vet Pharmacol Ther. 2009;32:79–86.View ArticlePubMedGoogle Scholar
- Devine C, Brennan GP, Lanusse CE, Alvarez LI, Trudgett A, Hoey E, et al. Potentiation of triclabendazole action in vivo against a triclabendazole-resistant isolate of Fasciola hepatica following its co-administration with the metabolic inhibitor, ketoconazole. Vet Parasitol. 2012;184:37–47.View ArticlePubMedGoogle Scholar
- Morphew RM, Eccleston N, Wilkinson TJ, McGarry J, Perally S, Prescott M, et al. Proteomics and in silico approaches to extend understanding of the glutathione transferase superfamily of the tropical liver fluke Fasciola gigantica. J Proteome Res. 2012;11:5876–89.PubMedGoogle Scholar
- Reed MB, Panaccio M, Strugnell RA, Spithill TW. Developmental expression of a Fasciola hepatica sequence homologous to ABC transporters. Int J Parasitol. 1998;28:1375–81.View ArticlePubMedGoogle Scholar
- Young ND, Hall RS, Jex AR, Cantacessi C, Gasser RB. Elucidating the transcriptome of Fasciola hepatica – A key to fundamental and biotechnological discoveries for a neglected parasite. Biotechnol Adv. 2010;28:222–31.View ArticlePubMedGoogle Scholar
- Pagel M, Johnstone RA. Variation across species in the size of the nuclear genome supports the junk-dna explanation for the c-value paradox. Proc Biol Sci. 1992;249:119–24.View ArticlePubMedGoogle Scholar
- Consortium TEP. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489:57–74.View ArticleGoogle Scholar
- Holroyd N, Sanchez-Flores A. Producing parasitic helminth reference and draft genomes at the Wellcome Trust Sanger Institute. Parasite Immunol. 2012;34:100–7.View ArticlePubMedGoogle Scholar
- Kang J-M, Bahk Y-Y, Cho P-Y, Hong S-J, Kim T-S, Sohn W-M, et al. A family of cathepsin F cysteine proteases of Clonorchis sinensis is the major secreted proteins that are expressed in the intestine of the parasite. Mol Biochem Parasitol. 2010;170:7–16.View ArticlePubMedGoogle Scholar
- Fox NJ, White PCL, McClean CJ, Marion G, Evans A, Hutchings MR. Predicting impacts of climate change on Fasciola hepatica risk. PLoS One. 2011;6:e16126.View ArticlePubMed CentralPubMedGoogle Scholar
- Zafra R, Pérez-Écija RA, Buffoni L, Moreno P, Bautista MJ, Martínez-Moreno A, et al. Early and late peritoneal and hepatic changes in goats immunized with recombinant cathepsin l1 and infected with Fasciola hepatica. J Comp Pathol. 2013;148:373–84.View ArticlePubMedGoogle Scholar
- Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Boetzer M, Henkel CV, Jansen HJ, Butler D, Pirovano W. Scaffolding pre-assembled contigs using SSPACE. Bioinformatics. 2011;27:578–9.View ArticlePubMedGoogle Scholar
- Cantarel BL, Korf I, Robb SMC, Parra G, Ross E, Moore B, et al. MAKER: An easy-to-use annotation pipeline designed for emerging model organism genomes. Genome Res. 2007;18:188–96.View ArticlePubMedGoogle Scholar
- Parra G, Bradnam K, Korf I. CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics. 2007;23:1061–7.View ArticlePubMedGoogle Scholar
- Hunt M, Kikuchi T, Sanders M, Newbold C, Berriman M, Otto TD. REAPR: a universal tool for genome assembly evaluation. Genome Biol. 2013;14:R47.View ArticlePubMed CentralPubMedGoogle Scholar
- Alexeyenko A, Tamas I, Liu G, Sonnhammer ELL. Automatic clustering of orthologs and inparalogs shared by multiple proteomes. Bioinformatics. 2006;22:e9–15.View ArticlePubMedGoogle Scholar
- Blythe MJ, Kao D, Malla S, Rowsell J, Wilson R, Evans D, et al. A dual platform approach to transcript discovery for the planarian Schmidtea mediterranea to establish RNAseq for stem cell and regeneration biology. PLoS One. 2010;5:e15617.View ArticlePubMed CentralPubMedGoogle Scholar
- Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14:R36.View ArticlePubMed CentralPubMedGoogle Scholar
- McCarthy DJ, Chen Y, Smyth GK. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 2012;40:4288–97.View ArticlePubMed CentralPubMedGoogle Scholar
- DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011;43:491–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol. 2011;7:539–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Centers for Disease Control and Prevention. Laboratory identification of parasitic diseases of public health concern. [http://www.cdc.gov/dpdx/fascioliasis/index.html]
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.