Mining nematode genome data for novel drug targets

Post-genomic progress in helminth parasitology

Helminth parasitology is an important discipline, which poses often unique technical challenges. One challenge is that helminth parasites, particularly those in humans, are often difficult to obtain alive and in sufficient quantities for study; another is the challenge of studying these organisms in vitro – no helminth parasite life cycle has been fully recapitulated outside of a host. Arguably, the key issue retarding progress in helminth parasitology has been a lack of experimental tools and resources, certainly relative to the riches that have driven many parasitologists to adopt free-living model organisms as surrogate systems. In response to these needs, the past 10–12 years have seen the beginnings of helminth parasitology's journey into the ‘omics’ era, with the release of abundant sequencing resources, and the functional genomics tools with which to test biological hypotheses. To reflect this progress, the 2019 Autumn Symposium of the British Society for Parasitology was held in Queen's University Belfast on the topic of ‘post-genomic progress in helminth parasitology’. This issue presents examples of the current state of play in the field, while this editorial summarizes how genomic datasets and functional genomic tools have stimulated impressive recent progress in our understanding of parasite biology.

Meloidogyne incognita (Nematoda: Meloidogynidae) sterol-binding protein Mi-SBP-1 as a target for its management

Meloidogyne incognita is a polyphagous plant-parasitic nematode that causes considerable yield loss in agricultural and horticultural crops. The management options available for M. incognita are extremely limited. Here we identified and characterised a M. incognita homolog of Caenorhabditis elegans sterol-binding protein (Mi-SBP-1), a transcriptional regulator of several lipogenesis pathway genes, and used RNA interference-mediated gene silencing to establish its utility as a target for the management of M. incognita. Mi-sbp-1 is predicted to be a helix-loop-helix domain containing DNA binding transcription factor, and is present in the M. incognita genome in three copies. The RNA-Seq analysis of Mi-sbp-1 silenced second stage juveniles confirmed the key role of this gene in lipogenesis regulation in M. incognita. In vitro and host-induced gene silencing of Mi-sbp-1 in M. incognita second stage juveniles resulted in loss of nematodes' ability to utilise the stored fat reserves, slower nematode development, and reduced parasitism on adzuki bean and tobacco plants. The multiplication factor for the Mi-sbp-1 silenced nematodes on adzuki bean plants was reduced by 51% compared with the control nematodes in which Mi-sbp-1 was not silenced. Transgenic expression of the double-stranded RNA construct of the Mi-sbp-1 gene in tobacco plants caused 40-45% reduction in M. incognita multiplication, 30-43.8% reduction in the number of egg masses, and 33-54% reduction in the number of eggs per egg mass compared with the wild type control plants. Our results confirm that Mi-sbp-1 is a key regulator of lipogenesis in M. incognita and suggest that it can be used as an effective target for its management. The findings of this study can be extended to develop methods to manage other economically important parasitic nematodes.

Functional analysis of a novel parasitic nematode-specific protein of Setaria digitata larvae in Culex quinquefasciatus by siRNA mediated RNA interference

BACKGROUND

Functional analysis of animal parasitic nematode genes is often quite challenging due to the unavailability of standardised in vitro culture conditions and lack of adequate tools to manipulate these genes. Therefore, this study was undertaken to investigate the suitability of Culex quinquefasciatus, as an in vivo culture platform for Setaria digitata larvae and RNA interference (RNAi), as a post-transcriptional gene silencing tool to study the roles of a vital gene that encodes a novel parasitic nematode-specific protein (SDNP).

RESULTS

The red colour fluorescence detected following RNAi injection to the thorax of C. quinquefasciatus indicated the uptake of dsRNA by S. digitata larvae. The reduction of SDNP transcripts in siRNA treated larvae compared to non-treated larvae, as determined by qPCR, indicated that the siRNA pathway is operational in S. digitata larvae. The observation of motility reductions and deformities during the development indicated the association of SDNP in larvae locomotion and development processes, respectively. The irregularities in the migration of larvae in mosquitoes and elevated survival rates of mosquitoes compared to their untreated counterparts indicated reduced parasitism of S. digitata larvae in mosquitoes upon targeted downregulation of SDNP by siRNA treatment.

CONCLUSION

SDNP plays vital roles in muscle contraction, locomotion, development processes, larval development and parasitism of S. digitata. Its ubiquitous presence in parasitic nematodes and its absence in their hosts provide a tantalising prospect of the possibility of targeting SDNP for future development of anthelmintic drugs. The susceptibility of the larval stages of S. digitata for RNAi in Culex quinquefasciatus was also demonstrated for the first time in this study.

Potential involvement of Brugia malayi cysteine proteases in the maintenance of the endosymbiotic relationship with Wolbachia

•

Bm-cpl-3 and Bm-cpl-6 are involved during development and endosymbiosis.

•

In tetracycline-treated worms Bm-cpl-3 and -6 are regulated in a bimodal pattern.

•

Reduction in Bm-cpl-3 and -6 levels resulted in hindered microfilarial development.

•

Reduction in Bm-cpl-3 and -6 levels resulted in reduced Wolbachia DNA levels.

•

These enzymes might be strong drug target candidates.

Brugia malayi, a parasitic nematode that causes lymphatic filariasis, harbors endosymbiotic intracellular bacteria, Wolbachia, that are required for the development and reproduction of the worm. The essential nature of this endosymbiosis led to the development of anti-Wolbachia chemotherapeutic approaches for the treatment of human filarial infections. Our study is aimed at identifying specific proteins that play a critical role in this endosymbiotic relationship leading to the identification of potential targets in the adult worms. Filarial cysteine proteases are known to be involved in molting and embryogenesis, processes shown to also be Wolbachia dependent. Based on the observation that cysteine protease transcripts are differentially regulated in response to tetracycline treatment, we focused on defining their role in symbiosis. We observe a bimodal regulation pattern of transcripts encoding cysteine proteases when in vitro tetracycline treated worms were examined. Using tetracycline-treated infertile female worms and purified embryos we established that the first peak of the bimodal pattern corresponds to embryonic transcripts while the second takes place within the hypodermis of the adult worms. Localization studies of the native proteins corresponding to Bm-cpl-3 and Bm-cpl-6 indicate that they are present in the area surrounding Wolbachia, and, in some cases, the proteins appear localized within the bacteria. Both proteins were also found in the inner bodies of microfilariae. The possible role of these cysteine proteases during development and endosymbiosis was further characterized using RNAi. Reduction in Bm-cpl-3 and Bm-cpl-6 transcript levels was accompanied by hindered microfilarial development and release, and reduced Wolbachia DNA levels, making these enzymes strong drug target candidates.

An analysis of the transcriptome of Teladorsagia circumcincta: its biological and biotechnological implications

BACKGROUND

Teladorsagia circumcincta (order Strongylida) is an economically important parasitic nematode of small ruminants (including sheep and goats) in temperate climatic regions of the world. Improved insights into the molecular biology of this parasite could underpin alternative methods required to control this and related parasites, in order to circumvent major problems associated with anthelmintic resistance. The aims of the present study were to define the transcriptome of the adult stage of T. circumcincta and to infer the main pathways linked to molecules known to be expressed in this nematode. Since sheep develop acquired immunity against T. circumcincta, there is some potential for the development of a vaccine against this parasite. Hence, we infer excretory/secretory molecules for T. circumcincta as possible immunogens and vaccine candidates.

RESULTS

A total of 407,357 ESTs were assembled yielding 39,852 putative gene sequences. Conceptual translation predicted 24,013 proteins, which were then subjected to detailed annotation which included pathway mapping of predicted proteins (including 112 excreted/secreted [ES] and 226 transmembrane peptides), domain analysis and GO annotation was carried out using InterProScan along with BLAST2GO. Further analysis was carried out for secretory signal peptides using SignalP and non-classical sec pathway using SecretomeP tools. For ES proteins, key pathways, including Fc epsilon RI, T cell receptor, and chemokine signalling as well as leukocyte transendothelial migration were inferred to be linked to immune responses, along with other pathways related to neurodegenerative diseases and infectious diseases, which warrant detailed future studies. KAAS could identify new and updated pathways like phagosome and protein processing in endoplasmic reticulum. Domain analysis for the assembled dataset revealed families of serine, cysteine and proteinase inhibitors which might represent targets for parasite intervention. InterProScan could identify GO terms pertaining to the extracellular region. Some of the important domain families identified included the SCP-like extracellular proteins which belong to the pathogenesis-related proteins (PRPs) superfamily along with C-type lectin, saposin-like proteins. The 'extracellular region' that corresponds to allergen V5/Tpx-1 related, considered important in parasite-host interactions, was also identified. Six cysteine motif (SXC1) proteins, transthyretin proteins, C-type lectins, activation-associated secreted proteins (ASPs), which could represent potential candidates for developing novel anthelmintics or vaccines were few other important findings. Of these, SXC1, protein kinase domain-containing protein, trypsin family protein, trypsin-like protease family member (TRY-1), putative major allergen and putative lipid binding protein were identified which have not been reported in the published T. circumcincta proteomics analysis. Detailed analysis of 6,058 raw EST sequences from dbEST revealed 315 putatively secreted proteins. Amongst them, C-type single domain activation associated secreted protein ASP3 precursor, activation-associated secreted proteins (ASP-like protein), cathepsin B-like cysteine protease, cathepsin L cysteine protease, cysteine protease, TransThyretin-Related and Venom-Allergen-like proteins were the key findings.

CONCLUSIONS

We have annotated a large dataset ESTs of T. circumcincta and undertaken detailed comparative bioinformatics analyses. The results provide a comprehensive insight into the molecular biology of this parasite and disease manifestation which provides potential focal point for future research. We identified a number of pathways responsible for immune response. This type of large-scale computational scanning could be coupled with proteomic and metabolomic studies of this parasite leading to novel therapeutic intervention and disease control strategies. We have also successfully affirmed the use of bioinformatics tools, for the study of ESTs, which could now serve as a benchmark for the development of new computational EST analysis pipelines.

Targeting the Wolbachia cell division protein FtsZ as a new approach for antifilarial therapy

The use of antibiotics targeting the obligate bacterial endosymbiont Wolbachia of filarial parasites has been validated as an approach for controlling filarial infection in animals and humans. Availability of genomic sequences for the Wolbachia (wBm) present in the human filarial parasite Brugia malayi has enabled genome-wide searching for new potential drug targets. In the present study, we investigated the cell division machinery of wBm and determined that it possesses the essential cell division gene ftsZ which was expressed in all developmental stages of B. malayi examined. FtsZ is a GTPase thereby making the protein an attractive Wolbachia drug target. We described the molecular characterization and catalytic properties of Wolbachia FtsZ. We also demonstrated that the GTPase activity was inhibited by the natural product, berberine, and small molecule inhibitors identified from a high-throughput screen. Furthermore, berberine was also effective in reducing motility and reproduction in B. malayi parasites in vitro. Our results should facilitate the discovery of selective inhibitors of FtsZ as a novel anti-symbiotic approach for controlling filarial infection. NOTE: The nucleotide sequences reported in this paper are available in GenBank™ Data Bank under the accession number wAlB-FtsZ (JN616286).

Computational prediction of essential genes in an unculturable endosymbiotic bacterium, Wolbachia of Brugia malayi

BACKGROUND

Wolbachia (wBm) is an obligate endosymbiotic bacterium of Brugia malayi, a parasitic filarial nematode of humans and one of the causative agents of lymphatic filariasis. There is a pressing need for new drugs against filarial parasites, such as B. malayi. As wBm is required for B. malayi development and fertility, targeting wBm is a promising approach. However, the lifecycle of neither B. malayi nor wBm can be maintained in vitro. To facilitate selection of potential drug targets we computationally ranked the wBm genome based on confidence that a particular gene is essential for the survival of the bacterium.

RESULTS

wBm protein sequences were aligned using BLAST to the Database of Essential Genes (DEG) version 5.2, a collection of 5,260 experimentally identified essential genes in 15 bacterial strains. A confidence score, the Multiple Hit Score (MHS), was developed to predict each wBm gene's essentiality based on the top alignments to essential genes in each bacterial strain. This method was validated using a jackknife methodology to test the ability to recover known essential genes in a control genome. A second estimation of essentiality, the Gene Conservation Score (GCS), was calculated on the basis of phyletic conservation of genes across Wolbachia's parent order Rickettsiales. Clusters of orthologous genes were predicted within the 27 currently available complete genomes. Druggability of wBm proteins was predicted by alignment to a database of protein targets of known compounds.

CONCLUSION

Ranking wBm genes by either MHS or GCS predicts and prioritizes potentially essential genes. Comparison of the MHS to GCS produces quadrants representing four types of predictions: those with high confidence of essentiality by both methods (245 genes), those highly conserved across Rickettsiales (299 genes), those similar to distant essential genes (8 genes), and those with low confidence of essentiality (253 genes). These data facilitate selection of wBm genes for entry into drug design pipelines.

Bioinformatic analysis of abundant, gender-enriched transcripts of adult Ascaris suum (Nematoda) using a semi-automated workflow platform

Expressed sequence tag (EST) data representing transcripts with a high level of differential hybridization in suppressive-subtractive hybridization (SSH)-based microarray analysis between adult female and male Ascaris suum were subjected to detailed bioinformatic analysis. A total of 361 ESTs clustered into 209 sequences, of which 52 and 157 represented transcripts that were enriched in female and male A. suum, respectively. Thirty (57.7%) of the 'female' subset of 52 sequences had orthologues/homologues in other parasitic nematodes and/or Caenorhabditis elegans, 13 (25%) exclusively in other parasitic nematodes and nine (17.3%) had no match in any other organism for which sequence data are currently available; the C. elegans orthologues encoded molecules involved in reproduction as well as embryonic and gamete development, such as vitellogenins and chitin-binding proteins. Of the 'male' subset of 157 sequences, 73 (46.5%) had orthologues/homologues in other parasitic nematodes and/or C. elegans, 57 (37.5%) in other parasitic nematodes only, and 22 (14.5%) had no significant similarity match in any other organism; the C. elegans orthologues encoded predominantly major sperm proteins (MSPs), kinases and phosphatases, actins, myosins and an Ancylostoma secreted protein-like molecule. The findings of the present study should support further genomic investigations of A. suum.

Functional analysis of the cathepsin-like cysteine protease genes in adult Brugia malayi using RNA interference

Background

Cathepsin-like enzymes have been identified as potential targets for drug or vaccine development in many parasites, as their functions appear to be essential in a variety of important biological processes within the host, such as molting, cuticle remodeling, embryogenesis, feeding and immune evasion. Functional analysis of Caenorhabditis elegans cathepsin L (Ce-cpl-1) and cathepsin Z (Ce-cpz-1) has established that both genes are required for early embryogenesis, with Ce-cpl-1 having a role in regulating in part the processing of yolk proteins. Ce-cpz-1 also has an important role during molting.

Methods and Findings

RNA interference assays have allowed us to verify whether the functions of the orthologous filarial genes in Brugia malayi adult female worms are similar. Treatment of B. malayi adult female worms with Bm-cpl-1, Bm-cpl-5, which belong to group Ia of the filarial cpl gene family, or Bm-cpz-1 dsRNA resulted in decreased numbers of secreted microfilariae in vitro. In addition, analysis of the intrauterine progeny of the Bm-cpl-5 or Bm-cpl Pro dsRNA- and siRNA-treated worms revealed a clear disruption in the process of embryogenesis resulting in structural abnormalities in embryos and a varied differential development of embryonic stages.

Conclusions

Our studies suggest that these filarial cathepsin-like cysteine proteases are likely to be functional orthologs of the C. elegans genes. This functional conservation may thus allow for a more thorough investigation of their distinct functions and their development as potential drug targets.

Filarial nematodes are an important group of human pathogens, causing lymphatic filariasis and onchocerciasis, and infecting around 150 million people throughout the tropics with more than 1.5 billion at risk of infection. Control of filariasis currently relies on mass drug administration (MDA) programs using drugs which principally target the microfilarial life-cycle stage. These control programs are facing major challenges, including the absence of a drug with macrofilaricidal or permanent sterilizing activity, and the possibility of the development of drug-resistance against the drugs available. Cysteine proteases are essential enzymes which play important roles in a wide range of cellular processes, and the cathepsin-like cysteine proteases have been identified as potential targets for drug or vaccine development in many parasites. Here we have studied the function of several of the cathepsin-like enzymes in the filarial nematode, B. malayi, and demonstrate that these cysteine proteases are involved in the development of embryos, show similar functions to their counterparts in C. elegans, and therefore, provide an important target for future drug development targeted to eliminate filariasis.

Genomic-bioinformatic analysis of transcripts enriched in the third-stage larva of the parasitic nematode Ascaris suum

Differential transcription in Ascaris suum was investigated using a genomic-bioinformatic approach. A cDNA archive enriched for molecules in the infective third-stage larva (L3) of A. suum was constructed by suppressive-subtractive hybridization (SSH), and a subset of cDNAs from 3075 clones subjected to microarray analysis using cDNA probes derived from RNA from different developmental stages of A. suum. The cDNAs (n = 498) shown by microarray analysis to be enriched in the L3 were sequenced and subjected to bioinformatic analyses using a semi-automated pipeline (ESTExplorer). Using gene ontology (GO), 235 of these molecules were assigned to ‘biological process’ (n = 68), ‘cellular component’ (n = 50), or ‘molecular function’ (n = 117). Of the 91 clusters assembled, 56 molecules (61.5%) had homologues/orthologues in the free-living nematodes Caenorhabditis elegans and C. briggsae and/or other organisms, whereas 35 (38.5%) had no significant similarity to any sequences available in current gene databases. Transcripts encoding protein kinases, protein phosphatases (and their precursors), and enolases were abundantly represented in the L3 of A. suum, as were molecules involved in cellular processes, such as ubiquitination and proteasome function, gene transcription, protein–protein interactions, and function. In silico analyses inferred the C. elegans orthologues/homologues (n = 50) to be involved in apoptosis and insulin signaling (2%), ATP synthesis (2%), carbon metabolism (6%), fatty acid biosynthesis (2%), gap junction (2%), glucose metabolism (6%), or porphyrin metabolism (2%), although 34 (68%) of them could not be mapped to a specific metabolic pathway. Small numbers of these 50 molecules were predicted to be secreted (10%), anchored (2%), and/or transmembrane (12%) proteins. Functionally, 17 (34%) of them were predicted to be associated with (non-wild-type) RNAi phenotypes in C. elegans, the majority being embryonic lethality (Emb) (13 types; 58.8%), larval arrest (Lva) (23.5%) and larval lethality (Lvl) (47%). A genetic interaction network was predicted for these 17 C. elegans orthologues, revealing highly significant interactions for nine molecules associated with embryonic and larval development (66.9%), information storage and processing (5.1%), cellular processing and signaling (15.2%), metabolism (6.1%), and unknown function (6.7%). The potential roles of these molecules in development are discussed in relation to the known roles of their homologues/orthologues in C. elegans and some other nematodes. The results of the present study provide a basis for future functional genomic studies to elucidate molecular aspects governing larval developmental processes in A. suum and/or the transition to parasitism.

In the present study, we constructed a cDNA library enriched for molecules of the infective third-stage larva (L3) of Ascaris suum, the common roundworm of pigs. Using the method of suppressive-subtractive hybridization (SSH), we explored transcription of a subset of molecules by microarray analysis and conducted bioinformatic analyses to characterize these molecules, map them to biochemical pathways, and predict genetic interactions based on comparisons with Caenorhabditis elegans and/or other organisms. The results provide interesting insights into early molecular processes in A. suum. Approximately 60% of the L3-enriched molecules discovered had homologues in C. elegans. Probabilistic analyses suggested that a complex genetic network regulates or controls larval growth and development in A. suum L3s, some of which might be involved in or regulate the switch from the free-living to the parasitic stage. Functional studies of these molecules to elucidate developmental processes in Ascaris could assist in identifying new targets for intervention.

Genes important in the parasitic life of the nematode Strongyloides ratti

Parasitic nematodes are important pathogens of humans and other animals. The genus Strongyloides has both a parasitic and a free-living adult generation. S. ratti infections of its rat host are negatively affected by the host immune response, such that a month after infection, worms are lost from the hosts. Here we have investigated the changes in parasite gene expression that occur as the anti-S. ratti immune pressure increases. Existing S. ratti expressed sequence tags were used to construct a microarray consisting of 2227 putative genes. This was probed with cDNA prepared from parasites subject to low or high immune pressures. There are significant changes in the gene expression of S. ratti when subject to different immune pressures. Most of the genes whose expression changes have no significant alignment to known genes. These data together with previous S. ratti EST data were then used to identify genes that we hypothesise are central to the parasitic life of S. ratti and, perhaps, other parasitic nematodes. These analyses have identified genes likely to play a key role in the parasitic life of S. ratti; these genes should be the priority for further investigation.

Mining predicted essential genes of Brugia malayi for nematode drug targets

We report results from the first genome-wide application of a rational drug target selection methodology to a metazoan pathogen genome, the completed draft sequence of Brugia malayi, a parasitic nematode responsible for human lymphatic filariasis. More than 1.5 billion people worldwide are at risk of contracting lymphatic filariasis and onchocerciasis, a related filarial disease. Drug treatments for filariasis have not changed significantly in over 20 years, and with the risk of resistance rising, there is an urgent need for the development of new anti-filarial drug therapies. The recent publication of the draft genomic sequence for B. malayi enables a genome-wide search for new drug targets. However, there is no functional genomics data in B. malayi to guide the selection of potential drug targets. To circumvent this problem, we have utilized the free-living model nematode Caenorhabditis elegans as a surrogate for B. malayi. Sequence comparisons between the two genomes allow us to map C. elegans orthologs to B. malayi genes. Using these orthology mappings and by incorporating the extensive genomic and functional genomic data, including genome-wide RNAi screens, that already exist for C. elegans, we identify potentially essential genes in B. malayi. Further incorporation of human host genome sequence data and a custom algorithm for prioritization enables us to collect and rank nearly 600 drug target candidates. Previously identified potential drug targets cluster near the top of our prioritized list, lending credibility to our methodology. Over-represented Gene Ontology terms, predicted InterPro domains, and RNAi phenotypes of C. elegans orthologs associated with the potential target pool are identified. By virtue of the selection procedure, the potential B. malayi drug targets highlight components of key processes in nematode biology such as central metabolism, molting and regulation of gene expression.

Molecular and biochemical characterization of nematode cofactor independent phosphoglycerate mutases

Phosphoglycerate mutase (PGM, EC 5.4.2.1) catalyzes the isomerization of 3-phosphoglycerate and 2-phosphoglycerate in glycolysis and gluconeogenesis. Two distinct types of PGM exist in nature, one that requires 2,3-bisphosphoglycerate as a cofactor (dPGM) and another that does not (iPGM). The two enzymes are structurally distinct and possess different mechanisms of action. In any particular organism, one form may exist or both. Nematodes possess the iPGM form whereas mammals have dPGM. In the present study, we have cloned and expressed iPGM from Onchocerca volvulus and described the catalytic properties of O. volvulus, Brugia malayi and Caenorhabditis elegans iPGM enzymes. Temperature and pH optima were determined for each enzyme. Like other iPGM enzymes, the activities of the nematode iPGM enzymes were dependent on the presence of divalent ions. Inactivation by EDTA could be restored most effectively by magnesium and manganese ions. Kinetic parameters and specific activities of the various recombinant enzymes were determined. The high similarity in catalytic properties among the enzymes indicates that a single enzyme inhibitor would likely be effective against all nematode enzymes. Inhibition of iPGM activity in vivo may lead to lethality as indicated by RNAi studies in C. elegans. Our results support the development of iPGM as a promising drug target in parasitic nematodes.

Oesophagostomum dentatum: potential as a model for genomic studies of strongylid nematodes, with biotechnological prospects

There are substantial gaps in the knowledge of the molecular processes of development and reproduction in parasitic nematodes, despite the fact that understanding such processes could lead to novel ways of treating and controlling parasitic diseases, through blocking or disrupting key biological pathways. Biotechnological advances through large-scale sequencing projects, approaches for the analysis of differential gene and protein expression and functional genomics (e.g., double-stranded RNA interference) now provide opportunities to investigate the molecular basis of developmental processes in some parasitic nematodes. The porcine nodule worm, Oesophagostomum dentatum (order Strongylida), may provide a platform for testing the function of genes from this and related nematodes, given that this species can be grown and maintained in culture in vitro for periods longer than other nematodes of the same order. In this article, we review relevant biological, biochemical and molecular biological and genomic information about O. dentatum and propose that the O. dentatum - pig system provides an attractive model for exploring molecular developmental and reproductive processes in strongylid nematodes, leading toward new intervention methods and biotechnological outcomes.

Construction of gender-enriched cDNA archives for adult Oesophagostomum dentatum by suppressive-subtractive hybridization and a microarray analysis of expressed sequence tags

In the present study, we constructed gender-enriched cDNA libraries for the adult stage of the parasitic nematode Oesophagostomum dentatum (order Strongylida) using suppressive-subtractive hybridization (SSH), sequenced clones from the female-library and male-library (480 from each) and conducted bioinformatic and microarray analyses of the expressed sequence tags (ESTs). In total, 873 ESTs (440 male and 433 female) were obtained, achieving a sequencing success of 91%The nucleotide sequences reported in this article (Tables 1-5) have been deposited in the EMBL, GenBank and DDJB databases under the Accession nos. AM157797-AM158083. Microarray analyses of 516 unique ESTs representing both gender-enriched libraries revealed differential hybridization for 391 of them (75.8%). Of these, 220 (56.3%) had significantly greater signal intensities in the female than in the male, and 154 (70%) of these were predicted to have homologues in C. elegans. These homologues were predicted to be involved in key biological processes, including embryonic nutrition, gametogenesis, molecular binding/transport or metabolism, nucleic acid synthesis and function, and signal transduction. Of the 171 ESTs with statistically higher signal intensities in male O. dentatum, 43.8% had homologues in C. elegans. These homologues included major sperm proteins (MSPs) or MSP-like molecules, keratin-like molecules, molecules involved in metabolism, PDZ domain-containing proteins, sugar binding proteins, protein kinases, serine proteases or protease inhibitors, molecules involved in proteolysis and other proteins, such as enzymes and various putative proteins. Of the 287 ESTs (from both gender-enriched cDNA libraries) with no known homologues in C. elegans, 50 (17.4%) had homologues in other nematodes, 8 had homologues in various other organisms and 104 (36.2%) had no homology to any sequence in current gene databases. The present study lays a foundation for the isolation and molecular, biochemical and functional characterization of selected genes from the gender-enriched cDNA archives established for O. dentatum.