- Research article
- Open Access
Differences in transcription between free-living and CO2-activated third-stage larvae of Haemonchus contortus
BMC Genomicsvolume 11, Article number: 266 (2010)
The disease caused by Haemonchus contortus, a blood-feeding nematode of small ruminants, is of major economic importance worldwide. The infective third-stage larva (L3) of this gastric nematode is enclosed in a cuticle (sheath) and, once ingested with herbage by the host, undergoes an exsheathment process that marks the transition from the free-living (L3) to the parasitic (xL3) stage. This study explored changes in gene transcription associated with this transition and predicted, based on comparative analysis, functional roles for key transcripts in the metabolic pathways linked to larval development.
Totals of 101,305 (L3) and 105,553 (xL3) expressed sequence tags (ESTs) were determined using 454 sequencing technology, and then assembled and annotated; the most abundant transcripts encoded transthyretin-like, calcium-binding EF-hand, NAD(P)-binding and nucleotide-binding proteins as well as homologues of Ancylostoma-secreted proteins (ASPs). Using an in silico-subtractive analysis, 560 and 685 sequences were shown to be uniquely represented in the L3 and xL3 stages, respectively; the transcripts encoded ribosomal proteins, collagens and elongation factors (in L3), and mainly peptidases and other enzymes of amino acid catabolism (in xL3). Caenorhabditis elegans orthologues of transcripts that were uniquely transcribed in each L3 and xL3 were predicted to interact with a total of 535 other genes, all of which were involved in embryonic development.
The present study indicated that some key transcriptional alterations taking place during the transition from the L3 to the xL3 stage of H. contortus involve genes predicted to be linked to the development of neuronal tissue (L3 and xL3), formation of the cuticle (L3) and digestion of host haemoglobin (xL3). Future efforts using next-generation sequencing and bioinformatic technologies should provide the efficiency and depth of coverage required for the determination of the complete transcriptomes of different developmental stages and/or tissues of H. contortus as well as the genome of this important parasitic nematode. Such advances should lead to a significantly improved understanding of the molecular biology of H. contortus and, from an applied perspective, to novel methods of intervention.
Parasitic nematodes of livestock animals are of major socio-economic importance worldwide due to the diseases and production losses that they cause. Haemonchus contortus (order Strongylida) is one of the most important nematodes of the abomasum (stomach) of small ruminants. The disease (= haemonchosis) caused by this parasite represents ~15% of all gastrointestinal diseases of small ruminants worldwide and causes major financial losses http://www.fao.org. Haemonchus contortus is a blood-feeding worm that causes anaemia and associated complications, leading to death in severely affected animals . In the abomasum of the ruminant host, the oviparous adult females release eggs via the faeces into the environment, in which first-stage larvae (L1s) develop and then hatch (within ~1 day, depending on conditions). The L1s develop into second-stage larvae (L2s), which then moult to become infective third-stage larvae (L3s). The cuticle from the L2 is retained as a sheath around the L3 and protects this stage from environmental pressures. Infective L3s are ingested with herbage by the host, pass through the forestomachs and undergo an exsheathment process to then establish, via the parasitic fourth-stage larvae (L4s), to dioecious haematophagous adults within ~3 weeks .
The exsheathment process marks the transition from the free-living (L3) to the parasitic (xL3) larval stage, at which H. contortus grows, develops and commences feeding on the blood of the host . Early studies examined aspects of this process in nematodes in response to a host stimulus. Evidence indicated that at least two separate pathways control these responses [3–5]. These pathways are stimulated by gaseous CO2 and mediated by carbonic anhydrase, leading to the secretion of the neurotransmitter noradrenaline . Noradrenaline then mediates the activation of genes required for further development and the release of exsheathment fluid to induce moulting of the L2 cuticle through a separate pathway. In the absence of CO2 stimulation, both pathways appear to be inhibited in L3s by an analogue of the "insect juvenile hormone", which is implicated in the regulation of development in insects [3–6]. It has been suggested that the CO2 cue is sensed by chemosensory neurons of the amphids [7, 8], which are sensory organs located in the head of L3s and are exposed to the external environment via a pore . The exsheathment process can be induced in vitro  and is recognised to suitably represent the process in vivo in the host animal [10, 11]. In spite of knowledge of the exsheathment process, the role of carbonic anhydrase and the chemoreception are poorly understood in H. contortus, and the exact nature of the regulation and components of the "downstream" molecular pathways are not yet known .
Despite these substantial knowledge gaps for H. contortus, Rogers and Petronijevic  did hypothesize early that there is likely to be a specific pattern of gene expression during the development of nematodes, whereby genes which control continuous "house-keeping" processes and maintain survival are expressed on a constitutive basis, whereas genes which have specific functions in one or more stages of development are regulated for expression exclusively at specific developmental time points. Although studies of H. contortus have shown that the pattern of transcription differs between free-living and parasitic stages [12–14], there has been no detailed study of the molecular alterations occurring during the early phase of transition to parasitism in this parasite, although there is some information for the canine hookworm, Ancylostoma caninum [15–18]. Insights into the transcriptome of H. contortus during this critical phase of establishment in its host would enhance knowledge of developmental processes at the molecular level and might identify new intervention targets.
Advances in genomic, proteomic and bioinformatic technologies [19–22] now provide opportunities for exploring the molecular basis of developmental switches in H. contortus and other nematodes. In particular, the combined use of next-generation sequencing, such as 454 technology, SOLiD and Illumina/Solexa [19, 23–25] and improved bioinformatic algorithms for the analysis and annotation of expressed sequence tag (EST) datasets  are suited to elucidate molecular changes at the transcriptomic level. In the present study, we (i) provide the first detailed insights into the transcriptome of H. contortus during its transition from the L3 to the exsheathed, parasitic xL3 stage using a next generation sequencing-based approach, (ii) predict, employing a number of bioinformatic approaches, the functional roles of these molecules in larval development and the metabolic pathways linked to this transition, and (iii) discuss the implications of the findings in relation to the fundamental, developmental biology of nematodes as well as applied aspects of developing new methods of nematode control.
Totals of 101,305 (L3) and 105,553 (xL3) expressed sequence tags (ESTs; average of 210-216 bases in length) were generated using 454 technology. The vast majority (95%) of these ESTs matched publicly available genomic DNA sequences of H. contortus. A summary of the characteristics of the raw and assembled data is given in Table 1. For L3, 9,046 proteins were inferred from a total of 20,066 assembled nucleotide sequences (Table 2), of which 3,066 matched known proteins with 1,465 different domains (Figure 1a; Additional file 1). For xL3, 9,001 proteins were inferred from a total of 20,116 assembled nucleotide sequences (Table 2), of which 2,885 mapped to known proteins with 1,394 distinct domains (Figure 1a; Additional file 1). For both L3 and xL3, the 'transthyretin-like' (8%, IPR001534), 'calcium-binding EF-hand' (6.7%, IPR002048), 'nicotinamide adenine dinucleotide(phosphate) [NAD(P)]-binding domain' [NAD(P)-binding] (5%, IPR016040), and 'nucleotide-binding, alpha beta plait' (5%, IPR012677) domains were most commonly detected (see Figure 1a; Additional file 1).
The subsequent annotation of the proteins inferred for the L3 stage revealed 1,229 different Gene Ontology (GO) terms, of which 334 were linked to 'biological process', 116 'cellular component' and 779 'molecular function' (see Figure 1b; Additional file 2). For xL3, 1,213 GO terms were identified, of which 330 represented 'biological process', 93 'cellular component' and 790 'molecular function'. For both L3 and xL3, the most common GO terms were 'translation' (9.2%, GO:0006412) and 'metabolic process' (7.1%, GO:0008152) for 'biological process', 'intracellular' (18.7%, GO:0005622) and 'ribosome' (15.1%, GO:0005840) for 'cellular component' and 'structural constituent of ribosome' (9%, GO:0003735) and 'catalytic activity' (6.3%, GO:0003824) for 'molecular function' (see Figure 1b; Additional file 2). Pathway mapping using KEGG Orthology-Based Annotation System (KOBAS) predicted a total of 253 different pathways, of which the most represented were 'ribosome', 'oxidative phosphorylation', 'chaperones and folding catalysts' (see Additional file 3).
Subtractive bioinformatic analysis identified 560 and 685 nucleotide sequences as being uniquely transcribed in the L3 and xL3 stages, respectively [available via http://research.vet.unimelb.edu.au/gasserlab/index.html]. The uniqueness of selected (n = 30) transcripts was verified by independent reverse transcription (rt)PCR-coupled sequencing of amplicons. Of the sequences unique to L3, 336 had orthologues in Caenorhabditis elegans (Additional file 4), 15.7% of which encoded ribosomal proteins (n = 32), collagens and elongation factors (n = 21). Of the 555 functional domains predicted for all proteins, 'glutathione S-transferase/chloride channel' (4.7%, IPR017933), 'homeodomain-like' (3.7%, IPR009057) and 'MORN motif' (3%, IPR003409, Figure 1a; Additional file 1) were abundant. Among the 204 GO terms identified in the L3-'unique' dataset, the most common were 'isoprenoid biosynthetic process' (4.6%, GO:0008299) 'DNA recombination' (3%, GO:0006310) and 'DNA replication' for 'biological process'; 'alpha DNA polymerase, primase complex' (1.1%, GO:0005658) for 'cellular component'; and, 'cation binding' (6.8%, GO:0443169), 'metallo carboxypeptidase activity' (5%, GO:0004181) and '4-hydroxyphenylpyruvate dioxygenase activity' (1.8%, GO:0003868) for 'molecular function'(Figure 1b; Additional file 2). KOBAS mapping revealed links to 21 different pathways, of which 'biosynthesis and biodegradation of secondary metabolites' and 'cell motility and secretion' were the most commonly represented (Additional file 3).
Of the 685 sequences unique to xL3, 421 had known orthologues in C. elegans, the largest number (i.e. 5%) representing peptidases and other enzymes of the amino acid catabolism (n = 21) (Additional file 4). Among the 484 functional domains identified in the InterProScan analysis, 'globin' or 'globin-like' domains (19.7%, IPR000971, IPR012292 and IPR009050) had the greatest representation (Figure 1a; Additional file 1). The GO analysis identified a total of 188 terms, the most common being 'oxygen transport' (16.9%, GO:0015671), 'regulation of phosphorylation' (9.1%, GO:0042325) and 'nitrogen compound metabolic process' (4.7%, GO:0006807) for 'metabolic process'; 'transcription factor TFIIA complex' (10.3%, GO:0005672) for 'cellular component'; and, 'oxygen binding' (11%, GO:0019825), 'phosphoric ester hydrolase activity' (6.3%, GO:0042578) and 'dopamine beta-monooxygenase activity' (2.6%, GO:0004500) for 'molecular function' (Figure 1b; Additional file 2). Among the seven KOBAS pathways predicted for the xL3-specific data, 'cytokine-cytokine receptor interaction' and 'D-arginine and D-ornithine metabolism' were most commonly represented (Additional file 3).
Twenty-one C. elegans orthologues (i.e. genes elc-1, rps-22, sod-2, rps-14, rpl-24.1, rpl-32, pas-4, ced-10, T10B11.2, krs-1, rpl-31, pas-5, cul-3, F21D5.7, rps-27, smo-1, rps-11, Y37E3.8a, rpl-5, col-97 and F57A8.2b) unique to L3 (cf. Additional file 5) were predicted to interact directly with a total of 535 other genes (range: 1-227, see Figure 2a and Additional file 7). In contrast, seven C. elegans orthologues (i.e. ncbp-2, F57B10.3a, vab-9, ncs-1, cpr-6, Y52B11A.2a and nhr-80) unique to xL3 (see Additional file 6) were inferred to interact with 45 other genes (range 1-27; see Figure 2b and Additional file 7). The functional classification of the C. elegans orthologues of molecules which were unique to L3 revealed one predominant group linked to 'embryonic development' (within 'biological process'); this cluster contained molecules associated with 'regulation of development' and other pathways, including 'positive regulation of growth', 'RNA metabolism' and 'biosynthesis' (Figure 3a; Additional file 7). Orthologues unique to xL3 also revealed one main cluster associated with 'embryonic development' (within 'biological process') (Figure 2b; Additional file 7).
In H. contortus, the transition to parasitism was proposed to be accompanied by differences in the levels of transcription of particular genes that may facilitate the invasion of the host and the evasion of the immune response . Next-generation sequencing technology was used here to investigate the transcriptome of the H. contortus L3 before and after the exsheathment process in vitro. Interestingly, ~75% of the sequences determined did not have orthologues/homologues in any other organism for which sequence data are available. However, annotation was enhanced through the mining of genomic data available for H. contortus (http://www.sanger.ac.uk; August 2008), ultimately increasing the percentage of annotated sequences to ~50% for the L3 and xL3 datasets (see Results section). The likely explanation for this result is technical and would appear to relate to a 3'-bias in sequence reads , thus affecting the prediction of open reading frames (ORFs) and, thus, the identification of orthologues/homologues. In the absence of the complete genome sequence for H. contortus, genomic mining provided the sole option for maximizing the mapping and subsequent annotation of sequences.
The integrated approach of using both transcriptomic and available genomic data allowed the annotation of a total of 11,302 sequences, of which 5,117 matched unique C. elegans orthologues. Of the 3,010 unique InterPro domains identified, the most represented in both L3 and xL3 were the 'transthyretin-like' (InterPro: IPR001534), 'calcium-binding EF-hand' (IPR002048), 'NAD(P)-binding' (IPR016040) and 'sperm-coating protein (SCP)-like extracellular' (IPR014044) motifs (Additional file 1). Most 'transthyretin' proteins identified to date are conserved (across large evolutionary distances)  and have been recognized as enzymes of the purine catabolism that catalyse the conversion of 5-hydroxyisourate (HIU) to OHCU [29, 30]. Of the 4,000 'nematode-specific' protein families encoded by 'nematode-restricted' genes, the 'transthyretin-like' proteins (TTLs)  represent one of the largest groups . Indeed, members of the TTL family have been detected in plant parasitic nematodes, such as Xiphinema index, Heterodera glycines, Meloidogyne incognita and Radophilis similis [31, 33–35], the filarial nematode of humans, Brugia malayi , and Ostertagia ostertagi of ruminants (which is related to H. contortus) [37, 38]. In O. ostertagi, at least 18 ttl genes have been identified by database mining; most of these genes are constitutively transcribed from the free-living L3 stage through to adult males and females . In parasitic nematodes, TTLs might play a role as carriers of lipophilic substances or hormones [32, 39]. More recently, Jacob and coworkers  localized the expression of ttl genes to the ventral nerve of R. similis (i.e. Rs-ttl-2) and to the nervous system and hypodermis of C. elegans (gene code R13A5.6; http://www.wormbase.org). These findings, together with the sequence similarity between TTLs and some neuropeptides, suggest that these proteins play a yet uncharacterized role in the nematode nervous system . In H. contortus, a TTL has been identified previously using a proteomic approach and shown to be the most immunogenic protein in excretory/secretory products (E/S) of the adult stage . TTLs are also highly abundant in E/S of the canine hookworm, A. caninum . These latter two studies propose a central role for these proteins in the parasite-host interplay [40, 41], which warrants detailed investigations.
Some host-parasite interactions have been reported to rely on the essential role played by calcium-mediated signalling pathways in protein secretion, motility, cell invasion and/or differentiation [42–44]. These functions are controlled by various, specialized subcellular structures (such as the Golgi apparatus, and some channels and transporters) for the uptake and release of calcium, which acts as a secondary messenger for the activation of calcium-dependent proteins , particularly those containing 'EF-hand' domains. In the present study, the large number of sequences (equating to 6.7% overall) with such domains in the L3 and xL3 stages of H. contortus might reflect a substantial need for calcium ions in the cascade of developmental events occurring during larval growth, particularly those linked to the maturation of the nervous system, as proposed to function in the C. elegans dauer stage . Similarly, the redox reactions catalysed by oxidoreductases with a NAD(P)-binding domain are essential for all developmental processes . The NAD+, a coenzyme found in all living cells, is involved in redox reactions and carries electrons from one reaction to another. In contrast, the main function of NADP+ is as a reducing agent in anabolism, with this coenzyme being involved in key pathways, such as fatty acid synthesis and glucose metabolism .
Another relatively abundant group of sequences in both L3 and xL3 encoded key SCP/Tpx-1/Ag5/PR-1/Sc7 proteins (designated SCP/TAPS; Pfam accession number no. PF00188), characterized by one or more SCP-like extracellular domains (InterPro: IPR014044, p < 1e-05) . In parasitic nematodes, such proteins are also called Ancylostoma secreted proteins or activation-associated proteins (ASPs), mainly because they were originally discovered in hookworms [18, 49]. Although the numbers of sequences encoding SCP/TAPS (187 and 207, assembled into three and seven EST clusters, respectively) were similar in L3 and xL3 (cf. Table 3), there were some qualitative differences in transcription between these two stages. Homologues of the molecules designated Hc24 and Hc40, previously identified in excretory/secretory (ES) products of adult H. contortus [50, 51], were encoded in both of these stages. However, other previously undescribed SCP/TAPS, such as an inferred protein homologue (encoded by C. elegans gene F09B9.5) with no known homologues in any other nematode, were encoded only by sequences in xL3 (cf. Table 3). The finding of entirely novel SCP/TAPS supports a previous hypothesis, formulated based on observations from transcriptomic and proteomic analyses, that a broader range of molecules of this group is expressed in H. contortus . This finding also suggests a diversified, active and specific involvement of SCP/TAPS in the exsheathment process, establishment of H. contortus in its host and/or the parasite-host interplay. Current evidence indicates that the transcriptional profile for H. contortus in the transition from L3 to xL3 is quite distinct from that of A. caninum, in that there was a massive expansion in SCP/TAPS in the serum-activated L3 (i.e. xL3) in the hookworm . Although there were methodological differences between studies, particularly in relation to the use of normalized cDNA libraries herein, an apparent difference in the diversity of transcripts encoding SCP/TAPS might reflect distinct roles for various types of these molecules in H. contortus compared with A. caninum during the invasion process in the host. Specifically, L3s of H. contortus are ingested by the ruminant host, exsheath in the reticulo-rumen and establish in the abomasum, in which the xL3s start feeding on blood to subsequently develop into dioecious adults . In contrast, infective L3s of A. caninum penetrate skin of the canid host, exsheath and then migrate as xL3s via the circulatory system and the lung to finally reside as adults in the small intestine; the dioecious hookworm adults attach via their buccal capsule to the intestinal mucosa, rupture capillaries and feed on blood . Some xL3s also undergo somatic migration to then encyst in tissues, such as muscle and the mammary gland, and become arrested as L3s [53–55]. Therefore, the added biological complexity in the hookworm life cycle, particularly in the initial phases of invasion and migration as an xL3, compared with the relatively simple gastric infection process for H. contortus, supports the current proposal of a very distinct arsenal of SCP/TAPS (both qualitatively and quantitatively) between the two nematodes during their transition to parasitism. Nonetheless, only a direct comparison of transcription between L3 and xL3 for the two parasites using the same or a similar approach as employed herein will allow this hypothesis to be tested conclusively.
The in silico subtraction approach employed in the present study identified 560 and 421 sequences that were specific to H. contortus L3 and xL3, respectively. The specificity of a subset of these transcripts was unequivocally verified by rtPCR-coupled sequencing (see Results section). Most of the subset of 'L3-specific' molecules with known orthologues in C. elegans encoded ribosomal proteins (n = 32), collagens (n = 10) and elongation factors (n = 11) (cf. Table 2; Additional file 4). In nematodes, the synthesis of collagens has been observed to increase significantly prior to a moult , whereas ribosomal proteins have been reported to participate in various cellular processes besides protein biosynthesis. For instance, they can act as components of the translation apparatus and also regulate cell proliferation and apoptosis . Indeed, genetic interaction networking for C. elegans orthologues of H. contortus L3-specific transcripts predicted a number of genes coding for ribosomal proteins (i.e. rps-11, rps-14, rps-22, rps-27, rpl-5 and rpl-24.1) to interact directly with a GTPase (i.e. ced-10), which is required for phagocytosis during programmed cell death and for migration of the distal tip cells of the somatic gonad  (see Figure 2a; Additional file 7). In addition, the protein CED-10 has been demonstrated to be essential in the cascade of events which leads to the development of the nervous system in the growing and developing larva, by guiding the migrating cells and axonal growth cones to their final position in the body . A direct interaction was also predicted between ced-10 and another C. elegans orthologue of an L3-specific transcript (i.e. pas-4) encoding a proteasome alpha-type seven subunit of the core 20S proteasome subcomplex (see Figure 2a; Additional file 7;(http://www.wormbase.org) involved in the 'ubiquitin proteasome system' (UPS) . Loss-of-function in this system is known to be the primary cause or secondary consequence of pathological conditions of the nervous system (e.g., ). Moreover, gene silencing of pas-4 in C. elegans has been shown to result in a wide variety of defects, including embryonic and larval lethality, sterility, abnormal locomotion, slow growth, and abnormal transgene expression and subcellular localization  http://http://www.wormbase.org. Here, both ced-10 and pas-4 were predicted herein to interact with a C. elegans orthologue of SUMO (i.e. smo-1, Figure 2a; Additional file 7). SUMO is a small ubiquitin-like moiety that, when attached to protein substrates, regulates subcellular localization and activity [62, 63]. The loss of function of smo-1 results in developmental defects of the nervous and reproductive systems as well as embryonic or larval lethality [62, 63]. Based on the genetic interaction networks predicted for orthologues of H. contortus L3-specific molecules, it could be proposed that axon guidance and synapse formation, which are primarily controlled by UPS, through the regulation of protein turnover at the growth cone and the synapse , are crucial in the biology of the L3 of this parasitic nematode. This statement is supported by the fact that the exsheathment process is triggered by gaseous CO2, detected by chemosensory neurons of amphids, located in the anterior end of the L3 [7, 8], ultimately leading to the secretion of the neurotransmitter noradrenaline [1, 3].
Of the H. contortus xL3-specific transcripts, those predicted to be involved in mechanisms regulating chemotaxis during larval development were also highly represented. For example, the C. elegans orthologue ncs-1 encodes a neuronal calcium sensor protein, whose expression in C. elegans is primarily associated with two sensory organs, the head amphids and tail phasmids . This orthologue was predicted to interact with 11 distinct genes involved in axon guidance, neuron projection, embryonic and/or larval development (see Figure 2b; Additional file 7). However, the largest number of genetic interactions was predicted for a C. elegans orthologue of a H. contortus xL3-specific transcript encoding a subunit of the nuclear cap-binding complex (i.e. ncbp-2, Figure 2b; Additional file 7). This complex includes RNA-binding proteins that bind to the 5'-cap within the nucleus; when RNA is exported to the cytoplasm, the nuclear cap-binding protein complex is replaced by cytoplasmic cap-binding complex . Although ncbp-2 has been demonstrated to play key roles in mRNA decay, embryonic/larval development, reproduction and vulval morphogenesis , its precise role in the biology of the developing larva is still unclear.
The largest number of C. elegans orthologues of H. contortus xL3-specific transcripts encoded peptidases and other enzymes of the amino acid catabolism (cf. Additional files 4 and 6). A similar spectrum of proteases and other molecules linked to catalytic activity have been shown also to be highly represented in the transcriptome of the serum-activated xL3 stage of A. caninum by comparison to its L3 . This finding, for two haematophagous bursate nematodes with differing biologies, is likely to reflect the key roles that these molecules play in host tissue invasion, degradation and/or digestion [1, 18]. In particular, a cysteine protease (i.e. cpr-6) was amongst the C. elegans orthologues of xL3-specific molecules in H. contortus (see Additional file 6). This finding is supported by previous evidence, showing that cysteine proteases play a crucial role in the catabolism of globin peptides by the cleavage of haemoglobin in blood-feeding nematodes (e.g., A. caninum and Necator americanus) [67–70]. Given this role, these proteases have been considered as promising candidates for developing recombinant vaccines against H. contortus as well as hookworms [71, 72]. Typically, these proteases include a 'globin domain' that protects the haeme iron from rapid oxidation and regulates oxygen homoeostasis  in the gut of the host, which is characterized by a reduced oxygen tension . Indeed, this domain (represented by InterPro codes IPR000971, IPR012292 and IPR009050, cf. Figure 1a; Additional file 1) was identified as the most abundant motif amongst the xL3-specific molecules. In H. contortus, transcripts encoding cathepsin B cysteine protease-like (CBL) proteins have been reported to be most abundant (~16%) in the intestine of the adult female . In addition, a comparison of H. contortus isolates from Europe and North America has revealed a minor to moderate nucleotide diversity in cbl genes, which has been proposed to reflect antigenic variation among CBLs . Clearly, C. elegans provides a useful surrogate system  to express H. contortus cysteine-proteases and to assess their function(s).
The present study has provided new and exciting insights into the molecular biology of the L3 stage of H. contortus and has elucidated transcriptional alterations taking place during the transition from the free-living to the parasitic stage of this nematode. Although approximately half of the sequences generated by 454 sequencing remains uncharacterised, as a consequence of the absence of complete genomic information for this parasite, future sequencing efforts using, for example, Illumina technology, should provide the efficiency and depth of coverage required to define the complete transcriptomes of all developmental stages and various tissues as well as the genome of H. contortus. The determination of the genome sequence of H. contortus has major potential to accelerate large-scale studies of genes and gene products involved in nematode development and reproduction, parasite-host interactions and the disease caused by the parasite. Importantly, comprehensive genomic and transcriptomic data will also underpin future proteomic and metabolomic studies of H. contortus. Such an integrated approach should lead to important conceptual advances in our understanding of various aspects of nematode biology and should have major implications, in the medium to long term, for the development of novel strategies for parasite intervention, resulting in biotechnological outcomes (such as drugs and diagnostic tests). For instance, future work could focus on defining a spectrum of key molecules involved in pathways linked to the development of the nervous system in different stages of H. contortus and assessing their potential as drug targets. Moreover, determining the structure and function of SCP/TAPS homologues/orthologues could establish their role(s) in the invasion of and establishment of this parasite in the host animal, providing the prospect of interfering with the host-parasite relationship. Although the present study focused on H. contortus, the research findings and the integrated technological approach employed should find broad applicability to other parasitic nematode of major global importance.
L3s of H. contortus ("Haecon5 strain") were cultured  from the faeces from an experimental sheep with a monospecific infection of this strain using a protocol described previously . Animal ethics approval (AEC no. 0707528) was given by The University of Melbourne, and the care and maintenance of sheep followed this institution's guidelines. The L3s were sieved (mesh size: 20 μm), washed extensively in sterile water and then stored at 10°C for 18 days. L3s (n = 500,000) were exsheathed , consistently achieving an exsheathment rate of 99%. Exsheathment using CO2 in vitro  is recognised to represent the process in vivo within the host animal [10, 11]. L3s and xL3s were each snap frozen in liquid nitrogen and then stored at -70°C until use.
Library construction and next-generation sequencing
For each L3 and xL3 of H. contortus, a normalized cDNA library was constructed (Eurofins MWG Operon, Ebersberg, Germany). In brief, total RNA was extracted from either L3 or xL3 using the mirVana isolation kit (Ambion). The integrity of each RNA sample was verified using a Bioanalyzer 2100 (Agilent Technologies), and the yield determined spectrophotometrically (ND-1000 UV-VIS v.3.2.1, NanoDrop Technologies). RNA was then treated with DNA se I (DNAfree, Ambion) and polyA+ RNA purified from ~120 ng of total RNA. For each library constructed, first-strand cDNA was synthesized using an oligo(dT)-adapter primer, and the second-strand was produced using a random (N)6-adapter primer. The resultant double-stranded cDNA was amplified (21 cycles) by Long and Accurate PCR (LA-PCR)  and an aliquot (~200 ng) examined electrophoretically. Normalization was conducted using one cycle of denaturation and reassociation of the cDNA, followed by purification of the reassociated double-stranded cDNA on a hydroxylapatite column. The normalized cDNA (500-700 bases) was then amplified using 9 cycles of LA-PCR and sequenced in a Genome Sequencer™ (GS) FLX Instrument (Roche Diagnostics) employing a standard protocol .
ESTs determined from the L3 and xL3 libraries were subjected to separate bioinformatic analyses. For each stage, all ESTs were pre-processed (SeqClean ; RepeatMasker ), aligned and then clustered using the Contig Assembly Program v.3 (CAP3) , employing a minimum sequence overlap of 30 bases and an identity threshold of 95% (for the removal of the flanking adapter sequences), and assembled. The small number of sequences (n = 55; 0.24% of 23,245 contigs) with a perfect match to those available for Ovis aries [GenBank: GI:3288836-GI:280977729; e-value < 1e -15] were eliminated from each dataset. Both unassembled and assembled EST datasets were compared, at the nucleotide level using BLASTn (e-value < 1e -05), with genomic sequence data publicly available for H. contortus (http://www.sanger.ac.uk; 21st August 2008). Contigs and singletons in each EST dataset, and the corresponding genomic sequence(s), were subjected to BLASTx (NCBI: http://www.ncbi.nlm.nih.gov) and BLASTn (EMBL-EBI Parasite Genome BLAST Server: http://www.ebi.ac.uk) searches to identify putative orthologues/homologues in C. elegans, other nematodes, and organisms other than nematodes (e-value of < 1e-05). WormBase (http://www.wormbase.org) was interrogated for relevant information on C. elegans orthologues/homologues, including RNAi phenotypes as well as transcriptomic, proteomic and interactomic data. Orthologues/homologues predicted from individual ESTs were always consistent with those inferred from genomic sequences (not shown). ESTs with no match to any sequences in the NCBI and/or EMBL-EBI databases were subjected to BLASTn analysis against the genomic data for H. contortus to increase the likelihood of identifying orthologues/homologues in currently available databases.
Following the conceptual translation of contigs and singletons into peptides using ESTScan , functional annotation was performed by GO using BLAST2GO . Domains/motifs within translated peptides were identified via InterProScan  and linked to pathways in C. elegans using KOBAS . The open reading frames (ORFs) inferred from sequences with orthologues in C. elegans were also subjected to 'secretome analysis' using the program SignalP v.2.0 (available at http://www.cbs.dtu.dk/services/SignalP/), employing both the neural network and hidden Markov models to predict signal peptides and/or anchors [87–89]. In addition, transmembrane domains were inferred using the program TMHMM (http://www.cbs.dtu.dk/services/TMHMM/; [90–92]).
Identification of transcripts unique to either L3 or xL3 by in silico subtraction, and verification of specificity by PCR-coupled sequencing
In order to identify transcripts unique to each stage, the L3 and xL3 datasets were subtracted from one another (in both directions) using in silico subtraction (at both the nucleotide and the amino acid levels) employing the BLASTn and tBLASTx algorithms, set at a stringent cut-off (e-value of < 1e-15). Subsets of molecules (sequences) that were present in either L3 (n = 10) or xL3 (n = 10) or in both stages (n = 10) were selected for subsequent experimental verification of specificity by rtPCR. Forward and reverse primers designed to selected contig sequences and corresponding genomic sequences (available at http://www.sanger.ac.uk) (Additional file 8), ensuring that at least one intronic region in the genomic sequence was included (per sequence) to enable the detection by rtPCR of any possible residual genomic DNA in the RNA used. The specificity of each primer designed was evaluated in silico by BLASTn analysis against all presently available databases, including that containing H. contortus genomic data. The rtPCR was carried out as described previously  with the following modifications: (i) the cycling conditions were: 95°C, 2 min (initial denaturation) followed by 35 cycles of 95°C for 1 min (denaturation), 55°C for 1 min (annealing) and 72°C for 1 min (extension), and a final extension of 72°C for 5 min; (ii) in addition to part of the β-tubulin (250 bp) gene being used as a positive control, a portion of the elongation factor-1α (216 bp) gene was included for the amplification (cf. Additional file 8). Negative control reactions using template RNA from L3 and xL3 were also included. All amplicons produced were sequenced  to demonstrate unequivocally that they represented the correct sequence. The Fisher's exact test was used to confirm that the number of molecules selected from each dataset (obtained following in silico subtraction and containing sequences unique to either L3 or xL3) for the verification of differential transcription by rtPCR was representative statistically (p < 0.001).
Probabilistic functional gene networking
The method developed by Zhong and Sternberg  was used to predict the interaction networks among C. elegans orthologues of molecules transcribed in either L3 or xL3. Data regarding interactions, phenotypes, expression and GO for selected C. elegans gene orthologues/homologues, also incorporating data from Drosophila melanogaster (vinegar fly), Saccharomyces cerevisiae (yeast), Mus musculus (mouse) and Homo sapiens (human), were integrated using a naïve Bayesian model to predict genetic interactions among C. elegans genes using the recommended, stringent cut-off value of 4.6 [94, 95]. The predicted networks resulting from the analyses were saved in a graphic display file (gdf) format, examined using the graph exploration system available at http://graphexploration.cond.org/ and drawn using Adobe Illustrator CS2 (Adobe Systems Inc.). The genes predicted to interact with selected transcripts unique to L3 and xL3 were also classified according to the 'Biological process Gene Ontology annotations of their interacting partners' using the PROSTIDIN website (http://crfb.univ-mrs.fr/webdistin/; ).
Nikolaou S, Gasser RB: Prospects for exploring molecular developmental processes in Haemonchus contortus. Int J Parasitol. 2006, 36: 859-868. 10.1016/j.ijpara.2006.04.007.
Veglia F: The anatomy and life history of the Haemonchus contortus (Rud). Rep Dir Vet Res. 1915, 3-4: 347-500.
Davey KG, Sommerville RI, Rogers WP: The effect of ethoxyzolamide, an analogue of insect juvenile hormone, nor-adrenaline and iodine on changes in the optical path difference in the excretory cells and oesophagus during exsheathment in Haemonchus contortus. Int J Parasitol. 1982, 12: 509-513. 10.1016/0020-7519(82)90045-5.
Rogers WP, Petronijevic T: The infective stage and development of nematodes. Biology and Control of Endoparasites. Edited by: Symons LEA, Donald AD, Dineen JK. 1982, Academic Press, Australia, 3-28.
Petronijevic T, Rogers WP: Gene activity and the development of early parasitic stages of nematodes. Int J Parasitol. 1983, 13: 197-199. 10.1016/0020-7519(83)90012-7.
Hartfelder K: Insect juvenile hormone: from "status quo" to high society. Braz J Med Biol Res. 2000, 33: 157-177. 10.1590/S0100-879X2000000200003.
Ashton FT, Li J, Schad GA: Chemo- and thermosensory neurons: structure and function in animal parasitic nematodes. Vet Parasitol. 1999, 84: 297-316. 10.1016/S0304-4017(99)00037-0.
Sciacca J, Forbes WM, Ashton FT, Lombardini E, Gamble HR, Schad GA: Response to carbon dioxide by the infective larvae of three species of parasitic nematodes. Parasitol Int. 2002, 51: 53-62. 10.1016/S1383-5769(01)00105-2.
Li J, Zhu X, Ashton FT, Gamble HR, Schad GA: Sensory neuroanatomy of a passively ingested nematode parasite, Haemonchus contortus: amphidial neurons of the third-stage larva. J Parasitol. 2001, 87: 65-72.
Conder GA, Johnson SS: Viability of infective larvae of Haemonchus contortus, Ostertagia ostertagi, and Trichostrongylus colubriformis following exsheathment by various techniques. J Parasitol. 1996, 82: 100-102. 10.2307/3284123.
Sommerville RI: Effect of carbon dioxide on the development of third-stage larvae of Haemonchus contortus in vitro. Nature. 1964, 202: 316-317. 10.1038/202316a0.
Hoekstra R, Visser A, Otsen M, Tibben J, Lenstra JA, Roos MH: EST sequencing of the parasitic nematode Haemonchus contortus suggests a shift in gene expression during transition to the parasitic stages. Mol Biochem Parasitol. 1964, 110: 53-68. 10.1016/S0166-6851(00)00255-3.
Hartman D, Donald DR, Nikolaou S, Savin KW, Hasse D, Presidente PJ, Newton SE: Analysis of developmentally regulated genes of the parasite Haemonchus contortus. Int J Parasitol. 2001, 31: 1236-1245. 10.1016/S0020-7519(01)00248-X.
Geldhof P, Whitton C, Gregory WF, Blaxter M, Knox DP: Characterisation of the two most abundant genes in the Haemonchus contortus expressed sequence tag dataset. Int J Parasitol. 2005, 35: 513-522. 10.1016/j.ijpara.2005.02.009.
Hawdon JM, Hotez PJ: Hookworm: developmental biology of the infectious process. Curr Opin Genet Dev. 1996, 6: 618-623. 10.1016/S0959-437X(96)80092-X.
Moser JM, Freitas T, Arasu P, Gibson G: Gene expression profiles associated with the transition to parasitism in Ancylostoma caninum larvae. Mol Biochem Parasitol. 2005, 143: 39-48. 10.1016/j.molbiopara.2005.04.012.
Datu BJ, Gasser RB, Nagaraj SH, Ong EK, O'Donoghue P, McInnes R, Ranganathan S, Loukas A: Transcriptional changes in the hookworm, Ancylostoma caninum, during the transition from a free-living to a parasitic larva. PLoS Negl Trop Dis. 2008, 2: e130-10.1371/journal.pntd.0000130.
Datu BJ, Loukas A, Cantacessi C, O'Donoghue P, Gasser RB: Investigation of the regulation of transcriptional changes in Ancylostoma caninum larvae following serum activation, with a focus on the insulin-like signalling pathway. Vet Parasitol. 2009, 159: 139-148. 10.1016/j.vetpar.2008.10.026.
Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA: Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005, 437: 376-380.
Nagaraj SH, Gasser RB, Ranganathan S: A hitchhiker's guide to expressed sequence tag (EST) analysis. Brief Bioinform. 2007, 8: 6-21. 10.1093/bib/bbl015.
Nagaraj SH, Deshpande N, Gasser RB, Ranganathan S: ESTExplorer: an expressed sequence tag (EST) assembly and annotation platform. Nucleic Acids Res. 2007, 35: W143-W147. 10.1093/nar/gkm378.
Mardis ER: The impact of next-generation sequencing technology on genetics. Trends Genet. 2008, 24: 133-141.
Pandey V, Nutter RC, Prediger E: Applied Biosystems SOLiD™ System: Ligation-Based Sequencing. Next Generation Genome Sequencing: Towards Personalized Medicine. Edited by: Jantz M. 2008, Milton, Australia, Wiley, 29-41. full_text.
Bentley DR, Balasubramanian S, Swerdlow HP, Smith GP, Milton J, Brown CG: Accurate whole human genome sequencing using reversible terminator chemistry. Nature. 2008, 456: 53-59. 10.1038/nature07517.
Harris TD, Buzby PR, Babcock H, Beer E, Bowers J, Braslavsky I: Single-molecule DNA sequencing of a viral genome. Science. 2008, 320: 106-109. 10.1126/science.1150427.
Ranganathan S, Menon R, Gasser RB: Advanced in silico analysis of expressed sequence tag (EST) data for parasitic nematodes of major socio-economic importance--fundamental insights toward biotechnological outcomes. Biotechnol Adv. 2009, 27: 439-448. 10.1016/j.biotechadv.2009.03.005.
Shin H, Hirst M, Bainbridge MN, Magrini V, Mardis E, Moerman DG, Marra MA, Baillie DL, Jones SJ: Transcriptome analysis for Caenorhabditis elegans based on novel expressed sequence tags. BMC Biol. 2008, 6: 30-10.1186/1741-7007-6-30.
Rahat O, Yitzhaky A, Schreiber G: Cluster conservation as a novel tool for studying protein-protein interactions evolution. Proteins. 2008, 71: 621-630. 10.1002/prot.21749.
Lee Y, Lee DH, Kho CW, Lee AY, Jang M, Cho S, Lee CH, Lee JS, Myung PK, Park BC, Park SG: Transthyretin-related proteins function to facilitate the hydrolysis of 5-hydroxyisourate, the end product of the uricase reaction. FEBS Lett. 2005, 579: 4769-4774. 10.1016/j.febslet.2005.07.056.
Ramazzina I, Folli C, Secchi A, Berni R, Percudani R: Completing the uric acid degradation pathway through phylogenetic comparison of whole genomes. Nat Chem Biol. 2006, 2: 144-148. 10.1038/nchembio768.
Jacob J, Vanholme B, Haegeman A, Gheysen G: Four transthyretin-like genes of the migratory plant-parasitic nematode Radopholus similis: members of an extensive nematode-specific family. Gene. 2007, 402: 9-19. 10.1016/j.gene.2007.07.015.
Parkinson J, Mitreva M, Whitton C, Thomson M, Daub J, Martin J, Schmid R, Hall N, Barrell B, Waterston RH, McCarter JP, Blaxter ML: A transcriptomic analysis of the phylum Nematoda. Nat Genet. 2004, 36: 1259-1267. 10.1038/ng1472.
Gao B, Allen R, Maier T, Davis EL, Baum TJ, Hussey RS: The parasitome of the phytonematode Heterodera glycines. Mol Plant Microbe Interact. 2003, 16: 720-726. 10.1094/MPMI.2003.16.8.720.
McCarter JP, Mitreva MD, Martin J, Dante M, Wylie T, Rao U, Pape D, Bowers Y, Theising B, Murphy CV, Kloek AP, Chiapelli BJ, Clifton SW, Bird DM, Waterston RH: Analysis and functional classification of transcripts from the nematode Meloidogyne incognita. Genome Biol. 2004, 4: R26-10.1186/gb-2003-4-4-r26.
Furlanetto C, Cardle L, Brown DJF, Jones JT: Analysis of expressed sequence tags from the ectoparasitic nematode Xiphinema index. Nematology. 2005, 7: 95-104. 10.1163/1568541054192180.
Hewitson JP, Harcus YM, Curwen RS, Dowle AA, Atmadja AK, Ashton PD, Wilson A, Maizels RM: The secretome of the filarial parasite, Brugia malayi: proteomic profile of adult excretory-secretory products. Mol Biochem Parasitol. 2008, 160: 8-21. 10.1016/j.molbiopara.2008.02.007.
Vercauteren I, Geldhof P, Peelaers I, Claerebout E, Berx G, Vercruysse J: Identification of excretory-secretory products of larval and adult Ostertagia ostertagi by immunoscreening of cDNA libraries. Mol Biochem Parasitol. 2003, 126: 201-208. 10.1016/S0166-6851(02)00274-8.
Saverwyns H, Visser A, Van Durme J, Power D, Morgado I, Kennedy MW, Knox DP, Schymkowitz J, Rousseau F, Gevaert K, Vercruysse J, Claerebout E, Geldhof P: Analysis of the transthyretin-like (TTL) gene family in Ostertagia ostertagi - Comparison with other strongylid nematodes and Caenorhabditis elegans. Int J Parasitol. 2008, 38: 1545-1556. 10.1016/j.ijpara.2008.04.004.
McElwee JJ, Schuster E, Blanc E, Thomas JH, Gems D: Shared transcriptional signature in Caenorhabditis elegans Dauer larvae and long-lived daf-2 mutants implicates detoxification system in longevity assurance. J Biol Chem. 2004, 279: 44533-44543. 10.1074/jbc.M406207200.
Yatsuda AP, Krijgsveld J, Cornelissen AW, Heck AJ, de Vries E: Comprehensive analysis of the secreted proteins of the parasite Haemonchus contortus reveals extensive sequence variation and differential immune recognition. J Biol Chem. 2003, 278: 16941-16951. 10.1074/jbc.M212453200.
Mulvenna J, Hamilton B, Nagaraj SH, Smyth D, Loukas A, Gorman JJ: Proteomics analysis of the excretory/secretory component of the blood-feeding stage of the hookworm, Ancylostoma caninum. Mol Cell Proteomics. 2009, 8: 109-121. 10.1074/mcp.M800206-MCP200.
Guillou F, Roger E, Moné Y, Rognon A, Grunau C, Théron A, Mitta G, Coustau C, Gourbal BE: Excretory-secretory proteome of larval Schistosoma mansoni and Echinostoma caproni, two parasites of Biomphalaria glabrata. Mol Biochem Parasitol. 2007, 155: 45-56. 10.1016/j.molbiopara.2007.05.009.
Vaid A, Thomas DC, Sharma P: Role of Ca2+/calmodulin-PfPKB signaling pathway in erythrocyte invasion by Plasmodium falciparum. J Biol Chem. 2008, 283: 5589-5597. 10.1074/jbc.M708465200.
Bagnall N, Gough J, Cadogan L, Burns B, Kongsuwan K: Expression of intracellular calcium signalling genes in cattle skin during tick infestation. Parasite Immunol. 2009, 31: 177-187. 10.1111/j.1365-3024.2008.01092.x.
Nagamune K, Moreno SN, Chini EN, Sibley LD: Calcium regulation and signaling in apicomplexan parasites. Subcell Biochem. 2008, 47: 70-81. full_text.
Estevez AO, Cowie RH, Gardner KL, Estevez M: Both insulin and calcium channel signaling are required for developmental regulation of serotonin synthesis in the chemosensory ADF neurons of Caenorhabditis elegans. Dev Biol. 2006, 298: 32-44. 10.1016/j.ydbio.2006.06.005.
Jaeger T, Flohé L: The thiol-based redox networks of pathogens: unexploited targets in the search for new drugs. Biofactors. 2006, 27: 109-120. 10.1002/biof.5520270110.
Cantacessi C, Campbell BE, Visser A, Geldhof P, Nolan MJ, Nisbet AJ, Matthews JB, Loukas A, Hofmann A, Otranto D, Sternberg PW, Gasser RB: A portrait of the "SCP/TAPS" proteins of eukaryotes--developing a framework for fundamental research and biotechnological outcomes. Biotechnol Adv. 2009, 27: 376-388. 10.1016/j.biotechadv.2009.02.005.
Hawdon JM, Jones BF, Hoffman DR, Hotez PJ: Cloning and characterization of Ancylostoma-secreted protein. A novel protein associated with the transition to parasitism by infective hookworm larvae. J Biol Chem. 1996, 271: 6672-6678. 10.1074/jbc.271.12.6672.
Schallig HD, van Leeuwen MA, Cornelissen AW: Protective immunity induced by vaccination with two Haemonchus contortus excretory secretory proteins in sheep. Parasite Immunol. 1997, 19: 447-453. 10.1046/j.1365-3024.1997.d01-148.x.
Rehman A, Jasmer DP: A tissue specific approach for analysis of membrane and secreted protein antigens from Haemonchus contortus gut and its application to diverse nematode species. Mol Biochem Parasitol. 1998, 97: 55-68. 10.1016/S0166-6851(98)00132-7.
Gasser RB, Cantacessi C, Loukas A: DNA technological progress toward advanced diagnostic tools to support human hookworm control. Biotechnol Adv. 2008, 26: 35-45. 10.1016/j.biotechadv.2007.09.003.
Schad GA, Chowdhury AB, Dean CG, Kochar VK, Nawalinski TA, Thomas J, Tonascia JA: Arrested development in human hookworm infections: an adaptation to a seasonally unfavorable external environment. Science. 1973, 180: 502-504. 10.1126/science.180.4085.502.
Loukas A, Prociv P: Immune responses in hookworm infections. Clin Microbiol Rev. 2001, 14: 689-703. 10.1128/CMR.14.4.689-703.2001.
Hotez PJ, Brooker S, Bethony JM, Bottazzi ME, Loukas A, Xiao S: Hookworm infection. N Engl J Med. 2004, 351: 799-807. 10.1056/NEJMra032492.
Fetterer RH: Growth and cuticular synthesis in Ascaris suum larvae during development from third to fourth stage in vitro. Vet Parasitol. 1996, 65: 275-282. 10.1016/S0304-4017(96)00956-9.
Rhoads RE, Dinkova TD, Korneeva NL: Mechanism and regulation of translation in C. elegans. WormBook. 2006, 28: 1-18.
Reddien PW, Horvitz HR: CED-2/CrkII and CED-10/Rac control phagocytosis and cell migration in Caenorhabditis elegans. Nat Cell Biol. 2000, 2: 131-136. 10.1038/35004000.
Quinn CC, Pfeil DS, Wadsworth WG: CED-10/Rac1 mediates axon guidance by regulating the asymmetric distribution of MIG-10/lamellipodin. Curr Biol. 2008, 18: 808-813. 10.1016/j.cub.2008.04.050.
Kahn NW, Rea SL, Moyle S, Kell A, Johnson TE: Proteasomal dysfunction activates the transcription factor SKN-1 and produces a selective oxidative-stress response in Caenorhabditis elegans. Biochem J. 2008, 409: 205-213. 10.1042/BJ20070521.
Ding M, Shen K: The role of the ubiquitin proteasome system in synapse remodeling and neurodegenerative diseases. Bioessays. 2008, 30: 1075-1083. 10.1002/bies.20843.
Broday L, Kolotuev I, Didier C, Bhoumik A, Gupta BP, Sternberg PW, Podbilewicz B, Ronai Z: The small ubiquitin-like modifier (SUMO) is required for gonadal and uterine-vulval morphogenesis in Caenorhabditis elegans. Genes Dev. 2004, 18: 2380-2391. 10.1101/gad.1227104.
Bhalla N, Wynne DJ, Jantsch V, Dernburg AF: ZHP-3 acts at crossovers to couple meiotic recombination with synaptonemal complex disassembly and bivalent formation in C. elegans. PLoS Genet. 2008, 4: e1000235-10.1371/journal.pgen.1000235.
Gomez M, De Castro E, Guarin E, Sasakura H, Kuhara A, Mori I, Bartfai T, Bargmann CI, Nef P: Ca2+ signaling via the neuronal calcium sensor-1 regulates associative learning and memory in C. elegans. Neuron. 2001, 30: 241-248. 10.1016/S0896-6273(01)00276-8.
Cook A, Bono F, Jinek M, Conti E: Structural biology of nucleocytoplasmic transport. Annu Rev Biochem. 2007, 76: 647-671. 10.1146/annurev.biochem.76.052705.161529.
Lall S, Piano F, Davis RE: Caenorhabditis elegans decapping proteins: localization and functional analysis of Dcp1, Dcp2, and DcpS during embryogenesis. Mol Biol Cell. 2005, 16: 5880-5890. 10.1091/mbc.E05-07-0622.
Pratt D, Cox GN, Milhausen MJ, Boisvenue RJ: A developmentally regulated cysteine protease gene family in Haemonchus contortus. Mol Biochem Parasitol. 1990, 43: 181-191. 10.1016/0166-6851(90)90143-A.
Williamson AL, Lecchi P, Turk BE, Choe Y, Hotez PJ, McKerrow JH, Cantley LC, Sajid M, Craik CS, Loukas A: A multi-enzyme cascade of hemoglobin proteolysis in the intestine of blood-feeding hookworms. J Biol Chem. 2004, 279: 35950-35957. 10.1074/jbc.M405842200.
Ranjit N, Zhan B, Stenzel DJ, Mulvenna J, Fujiwara R, Hotez PJ, Loukas A: A family of cathepsin B cysteine proteases expressed in the gut of the human hookworm, Necator americanus. Mol Biochem Parasitol. 2008, 160: 90-99. 10.1016/j.molbiopara.2008.04.008.
Ranjit N, Zhan B, Hamilton B, Stenzel D, Lowther J, Pearson M, Gorman J, Hotez P, Loukas A: Proteolytic degradation of hemoglobin in the intestine of the human hookworm Necator americanus. J Infect Dis. 2009, 199: 904-912. 10.1086/597048.
Skuce PJ, Redmond DL, Liddell S, Stewart EM, Newlands GF, Smith WD, Knox DP: Molecular cloning and characterization of gut-derived cysteine proteinases associated with a host protective extract from Haemonchus contortus. Parasitology. 1999, 119: 405-412. 10.1017/S0031182099004813.
Loukas A, Bethony JM, Williamson AL, Goud GN, Mendez S, Zhan B, Hawdon JM, Bottazzi ME, Brindley PJ, Hotez PJ: Vaccination of dogs with a recombinant cysteine protease from the intestine of canine hookworms diminishes the fecundity and growth of worms. J Infect Dis. 2004, 189: 1952-1961. 10.1086/386346.
Freitas TA, Saito JA, Hou S, Alam M: Globin-coupled sensors, protoglobins, and the last universal common ancestor. J Inorg Biochem. 2005, 99: 23-33. 10.1016/j.jinorgbio.2004.10.024.
Fetterer RH, Hill DE, Rhoads ML: Characterization of a hemoglobin-like protein from adult Haemonchus contortus. J Parasitol. 1999, 85: 295-300. 10.2307/3285637.
Jasmer DP, Mitreva MD, McCarter JP: mRNA sequences for Haemonchus contortus intestinal cathepsin B-like cysteine proteases display an extreme in abundance and diversity compared with other adult mammalian parasitic nematodes. Mol Biochem Parasitol. 2004, 137: 297-305. 10.1016/j.molbiopara.2004.06.010.
Murray L, Geldhof P, Clark D, Knox DP, Britton C: Expression and purification of an active cysteine protease of Haemonchus contortus using Caenorhabditis elegans. Int J Parasitol. 2007, 37: 1117-1125. 10.1016/j.ijpara.2007.02.012.
Sommerville RI: Development of Haemonchus contortus in vitro and the stimulus from the host. J Parasitol. 1977, 63: 344-347. 10.2307/3280076.
Nikolaou S, Hartman D, Presidente PJ, Newton SE, Gasser RB: HcSTK, a Caenorhabditis elegans PAR-1 homologue from the parasitic nematode, Haemonchus contortus. Int J Parasitol. 2002, 32: 749-758. 10.1016/S0020-7519(02)00008-5.
Barnes WM: PCR amplification of up to 35-kb DNA with high fidelity and high yield from lambda bacteriophage templates. Proc Natl Acad Sci USA. 1994, 91: 2216-2220. 10.1073/pnas.91.6.2216.
Chen YA, Lin CC, Wang CD, Wu HB, Hwang PI: An optimized procedure greatly improves EST vector contamination removal. BMC Genomics. 2007, 8: 416-10.1186/1471-2164-8-416.
Chen N: Using RepeatMasker to identify repetitive elements in genomic sequences. Curr Protoc Bioinformatics. 2004, 4: Unit 4.10-
Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9: 868-877. 10.1101/gr.9.9.868.
Nagaraj SH, Gasser RB, Nisbet AJ, Ranganathan S: In silico analysis of expressed sequence tags from Trichostrongylus vitrinus (Nematoda): comparison of the automated ESTExplorer workflow platform with conventional database searches. BMC Bioinf. 2008, 9: S10-10.1186/1471-2105-9-S1-S10.
Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676. 10.1093/bioinformatics/bti610.
Hunter S, Apweiler R, Attwood TK, Bairoch A, Bateman A, Binns D, Bork P: InterPro: the integrative protein signature database. Nucleic Acids Res. 2009, 37: D211-D215. 10.1093/nar/gkn785.
Wu J, Mao X, Cai T, Luo J, Wei L: KOBAS server: a web-based platform for automated annotation and pathway identification. Nucleic Acids Res. 2006, 34: W720-W724. 10.1093/nar/gkl167.
Nielsen H, Engelbrecht J, Brunak S, von Heijne G: Identification of prokaryotic and eukaryotic signal peptides and prediction of their cleavage sites. Prot Eng. 1997, 10: 1-6. 10.1093/protein/10.1.1.
Nielsen H, Krogh A: Prediction of signal peptides and signal anchors by a hidden Markov model. Proceedings of the Sixth International Conference on Intelligent Systems for Molecular Biology (ISMB 6). 1998, AAAI Press, Menlo Park, California, 122-130.
Bendtsen JD, Nielsen H, von Heijne G, Brunak S: Improved prediction of signal peptides: SignalP 3.0. J Mol Biol. 2004, 340: 783-795. 10.1016/j.jmb.2004.05.028.
Sonnhammer ELL, von Heijne G, Krogh A: A hidden Markov model for predicting transmembrane helices in protein sequences. Proceedings of the Sixth International Conference on Intelligent Systems for Molecular Biology. 1998, Menlo Park, CA, AAAI, 175-182.
Krogh A, Larsson B, von Heijne G, Sonnhammer ELL: Predicting transmembrane protein topology with a hidden Markov model: Application to complete genomes. J Mol Biol. 2001, 305: 567-580. 10.1006/jmbi.2000.4315.
Moller S, Croning MDR, Apweiler R: Evaluation of methods for the prediction of membrane spanning regions. Bioinformatics. 2001, 17: 646-653. 10.1093/bioinformatics/17.7.646.
Hu M, Abs EL-Osta YG, Campbell BE, Boag PR, Nisbet AJ, Beveridge I, Gasser RB: Trichostrongylus vitrinus (Nematoda: Strongylida): molecular characterization and transcriptional analysis of Tv-stp-1, a serine/threonine phosphatase gene. Exp Parasitol. 2007, 117: 22-34. 10.1016/j.exppara.2007.03.008.
Campbell BE, Nagaraj SH, Hu M, Zhong W, Sternberg PW, Ong EK, Loukas A, Ranganathan S, Beveridge I, McInnes RL, Hutchinson GW, Gasser RB: Gender-enriched transcripts in Haemonchus contortus--predicted functions and genetic interactions based on comparative analyses with Caenorhabditis elegans. Int J Parasitol. 2008, 38: 65-83. 10.1016/j.ijpara.2007.07.001.
Zhong W, Sternberg PW: Genome-wide prediction of C. elegans genetic interactions. Science. 2006, 311: 1481-1484. 10.1126/science.1123287.
Baudot A, Martin D, Mouren P, Chevenet F, Guénoche A, Jacq B, Brun C: PRODISTIN Web Site: a tool for the functional classification of proteins from interaction networks. Bioinformatics. 2006, 22: 248-250. 10.1093/bioinformatics/bti757.
This work was supported by grants from the Australian Research Council (RBG, AL, PWS) and the Australian-American Fulbright Commission (RBG). CC is the grateful recipient of International Postgraduate Research Scholarship (IPRS) from the Australian Government and a fee-remission scholarship through the University of Melbourne as well as the Clunies Ross (2008) and Sue Newton (2009) awards from the School of Veterinary Science of the same university. The staff of WormBase is acknowledged for their contributions.
This article is dedicated to the late Paul JA Presidente.
CC performed the bioinformatic analyses, analysed the results and drafted the manuscript, BEC, NDY, ARJ and RSH also participated in the bioinformatic analyses, PJAP and JLZ provided the parasite material, WZ and BAM assisted in the probabilistic genetic interaction network predictions, AL and PWS contributed to the drafting of the manuscript. RBG conceived and designed the study, coordinated and supervised the project and drafted the manuscript with CC. All authors read and approved the final manuscript.