Nematode Genomes Research Articles

The plant parasite Pratylenchus coffeae carries a minimal nematode genome

Here we report the genome sequence of the lesion nematode, Pratylenchus coffeae, a significant pest of banana and other staple crops in tropical and sub-tropical regions worldwide. Initial analysis of the 19.67 Mb genome reveals 6712 protein encoding genes, the smallest number found in a metazoan, although sufficient to make a nematode. Significantly, no developmental or physiological pathways are obviously missing when compared to the model free-living nematode Caenorhabditis elegans, which possesses approximately 21 000 genes. The highly streamlined P. coffeae genome may reveal a remarkable functional plasticity in nematode genomes and may also indicate evolutionary routes to increased specialisation in other nematode genera. In addition, the P. coffeae genome may begin to reveal the core set of genes necessary to make a multicellular animal. Nematodes exhibit striking diversity in the niches they occupy, and the sequence of P. coffeae is a tool to begin to unravel the mechanisms that enable the extraordinary success of this phylum as both free-living and parasitic forms. Unlike the sedentary endoparasitic root-knot nematodes (Meloidogyne spp.), P. coffeae is a root-lesion nematode that does not establish a feeding site within the root. Because the P. coffeae nematode genome encodes fewer than half the number of genes found in the genomes of root-knot nematodes, comparative analysis to determine genes P. coffeae does not carry may help to define development of more sophisticated forms of nematode-plant interactions. The P. coffeae genome sequence may help to define timelines related to evolution of parasitism amongst nematodes. The genome of P. coffeae is a significant new tool to understand not only nematode evolution but animal biology in general.

What helminth genomes have taught us about parasite evolution.

The genomes of more than 20 helminths have now been sequenced. Here we perform a meta-analysis of all sequenced genomes of nematodes and Platyhelminthes, and attempt to address the question of what are the defining characteristics of helminth genomes. We find that parasitic worms lack systems for surface antigenic variation, instead maintaining infections using their surfaces as the first line of defence against the host immune system, with several expanded gene families of genes associated with the surface and tegument. Parasite excretory/secretory products evolve rapidly, and proteases even more so, with each parasite exhibiting unique modifications of its protease repertoire. Endoparasitic flatworms show striking losses of metabolic capabilities, not matched by nematodes. All helminths do however exhibit an overall reduction in auxiliary metabolism (biogenesis of co-factors and vitamins). Overall, the prevailing pattern is that there are few commonalities between the genomes of independently evolved parasitic worms, with each parasite having undergone specific adaptations for their particular niche.

Extensive duplication of the Wolbachia DNA in chromosome four of Drosophila ananassae.

BackgroundLateral gene transfer (LGT) from bacterial Wolbachia endosymbionts has been detected in ~20% of arthropod and nematode genome sequencing projects. Many of these transfers are large and contain a substantial part of the Wolbachia genome.ResultsHere, we re-sequenced three D. ananassae genomes from Asia and the Pacific that contain large LGTs from Wolbachia. We find that multiple copies of the Wolbachia genome are transferred to the Drosophila nuclear genome in all three lines. In the D. ananassae line from Indonesia, the copies of Wolbachia DNA in the nuclear genome are nearly identical in size and sequence yielding an even coverage of mapped reads over the Wolbachia genome. In contrast, the D. ananassae lines from Hawaii and India show an uneven coverage of mapped reads over the Wolbachia genome suggesting that different parts of these LGTs are present in different copy numbers. In the Hawaii line, we find that this LGT is underrepresented in third instar larvae indicative of being heterochromatic. Fluorescence in situ hybridization of mitotic chromosomes confirms that the LGT in the Hawaii line is heterochromatic and represents ~20% of the sequence on chromosome 4 (dot chromosome, Muller element F).ConclusionsThis collection of related lines contain large lateral gene transfers composed of multiple Wolbachia genomes that constitute >2% of the D. ananassae genome (~5 Mbp) and partially explain the abnormally large size of chromosome 4 in D. ananassae.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-1097) contains supplementary material, which is available to authorized users.

New insights into the evolution of vertebrate CRH (corticotropin-releasing hormone) and invertebrate DH44 (diuretic hormone 44) receptors in metazoans

The corticotropin releasing hormone receptors (CRHR) and the arthropod diuretic hormone 44 receptors (DH44R) are structurally and functionally related members of the G protein-coupled receptors (GPCR) of the secretin-like receptor superfamily. We show here that they derive from a bilaterian predecessor. In protostomes, the receptor became DH44R that has been identified and functionally characterised in several arthropods but the gene seems to be absent from nematode genomes. Duplicate DH44R genes (DH44 R1 and DH44R2) have been described in some arthropods resulting from lineage-specific duplications. Recently, CRHR-DH44R-like receptors have been identified in the genomes of some lophotrochozoans (molluscs, which have a lineage-specific gene duplication, and annelids) as well as representatives of early diverging deuterostomes. Vertebrates have previously been reported to have two CRHR receptors that were named CRHR1 and CRHR2. To resolve their origin we have analysed recently assembled genomes from representatives of early vertebrate divergencies including elephant shark, spotted gar and coelacanth. We show here by analysis of synteny conservation that the two CRHR genes arose from a common ancestral gene in the early vertebrate tetraploidizations (2R) approximately 500 million years ago. Subsequently, the teleost-specific tetraploidization (3R) resulted in a duplicate of CRHR1 that has been lost in some teleost lineages. These results help distinguish orthology and paralogy relationships and will allow studies of functional conservation and changes during evolution of the individual members of the receptor family and their multiple native peptide agonists.

Operons are a conserved feature of nematode genomes.

The organization of genes into operons, clusters of genes that are co-transcribed to produce polycistronic pre-mRNAs, is a trait found in a wide range of eukaryotic groups, including multiple animal phyla. Operons are present in the class Chromadorea, one of the two main nematode classes, but their distribution in the other class, the Enoplea, is not known. We have surveyed the genomes of Trichinella spiralis, Trichuris muris, and Romanomermis culicivorax and identified the first putative operons in members of the Enoplea. Consistent with the mechanism of polycistronic RNA resolution in other nematodes, the mRNAs produced by genes downstream of the first gene in the T. spiralis and T. muris operons are trans-spliced to spliced leader RNAs, and we are able to detect polycistronic RNAs derived from these operons. Importantly, a putative intercistronic region from one of these potential enoplean operons confers polycistronic processing activity when expressed as part of a chimeric operon in Caenorhabditis elegans. We find that T. spiralis genes located in operons have an increased likelihood of having operonic C. elegans homologs. However, operon structure in terms of synteny and gene content is not tightly conserved between the two taxa, consistent with models of operon evolution. We have nevertheless identified putative operons conserved between Enoplea and Chromadorea. Our data suggest that operons and “spliced leader” (SL) trans-splicing predate the radiation of the nematode phylum, an inference which is supported by the phylogenetic profile of proteins known to be involved in nematode SL trans-splicing.

Bioinformatic exploration of RIO protein kinases of parasitic and free-living nematodes

Despite right open reading frame kinases (RIOKs) being essential for life, their functions, substrates and cellular pathways remain enigmatic. In the present study, gene structures were characterised for 26 RIOKs from draft genomes of parasitic and free-living nematodes. RNA-seq transcription profiles of riok genes were investigated for selected parasitic nematodes and showed that these kinases are transcribed in developmental stages that infect their mammalian host. Three-dimensional structural models of Caenorhabditis elegans RIOKs were predicted, and elucidated functional domains and conserved regions in nematode homologs. These findings provide prospects for functional studies of riok genes in C. elegans, and an opportunity for the design and validation of nematode-specific inhibitors of these enzymes in socioeconomic parasitic worms.

The evolution of parasitism in Nematoda.

Nematodes are abundant and diverse, and include many parasitic species. Molecular phylogenetic analyses have shown that parasitism of plants and animals has arisen at least 15 times independently. Extant nematode species also display lifestyles that are proposed to be on the evolutionary trajectory to parasitism. Recent advances have permitted the determination of the genomes and transcriptomes of many nematode species. These new data can be used to further resolve the phylogeny of Nematoda, and identify possible genetic patterns associated with parasitism. Plant-parasitic nematode genomes show evidence of horizontal gene transfer from other members of the rhizosphere, and these genes play important roles in the parasite-host interface. Similar horizontal transfer is not evident in animal parasitic groups. Many nematodes have bacterial symbionts that can be essential for survival. Horizontal transfer from symbionts to the nematode is also common, but its biological importance is unclear. Over 100 nematode species are currently targeted for sequencing, and these data will yield important insights into the biology and evolutionary history of parasitism. It is important that these new technologies are also applied to free-living taxa, so that the pre-parasitic ground state can be inferred, and the novelties associated with parasitism isolated.

Adjustments, extinction, and remains of selenocysteine incorporation machinery in the nematode lineage.

Selenocysteine (Sec) is encoded by an UGA codon with the help of a SECIS element present in selenoprotein mRNAs. SECIS-binding protein (SBP2/SCBP-2) mediates Sec insertion, but the roles of its domains and the impact of its deficiency on Sec insertion are not fully understood. We used Caenorhabditis elegans to examine SBP2 function since it possesses a single selenoprotein, thioredoxin reductase-1 (TRXR-1). All SBP2 described so far have an RNA-binding domain (RBD) and a Sec-incorporation domain (SID). Surprisingly, C. elegans SBP2 lacks SID and consists only of an RBD. An sbp2 deletion mutant strain ablated Sec incorporation demonstrating SBP2 essentiality for Sec incorporation. Further in silico analyses of nematode genomes revealed conservation of SBP2 lacking SID and maintenance of Sec incorporation linked to TRXR-1. Remarkably, parasitic plant nematodes lost the ability to incorporate Sec, but retained SecP43, a gene associated with Sec incorporation. Interestingly, both selenophosphate synthetase (SPS) genes are absent in plant parasitic nematodes, while only Cys-containing SPS2 is present in Sec-incorporating nematodes. Our results indicate that C. elegans and the nematode lineage provide key insights into Sec incorporation and the evolution of Sec utilization trait, selenoproteomes, selenoproteins, and Sec residues. Finally, our study provides evidence of noncanonical translation initiation in C. elegans, not previously known for this well-established animal model.

Characterization of the reniform nematode genome by shotgun sequencing.

The reniform nematode (RN), a major agricultural pest particularly on cotton in the United States, is among the major plant-parasitic nematodes for which limited genomic information exists. In this study, over 380 Mb of sequence data were generated from pooled DNA of four adult female RNs and assembled into 67,317 contigs, including 25,904 (38.5%) predicted coding contigs and 41,413 (61.5%) noncoding contigs. Most of the characterized repeats were of low complexity (88.9%), and 0.9% of the contigs matched with 53.2% of GenBank ESTs. The most frequent Gene Ontology (GO) terms for molecular function and biological process were protein binding (32%) and embryonic development (20%). Further analysis showed that 741 (1.1%), 94 (0.1%), and 169 (0.25%) RN genomic contigs matched with 1328 (13.9%), 1480 (5.4%), and 1330 (7.4%) supercontigs of Meloidogyne incognita, Brugia malayi, and Pristionchus pacificus, respectively. Chromosome 5 of Caenorhabditis elegans had the highest number of hits to the RN contigs. Seven putative detoxification genes and three carbohydrate-active enzymes (CAZymes) involved in cell wall degradation were studied in more detail. Additionally, kinases, G protein-coupled receptors, and neuropeptides functioning in physiological, developmental, and regulatory processes were identified in the RN genome.

Nematode and Arthropod Genomes Provide New Insights into the Evolution of Class 2 B1 GPCRs

Nematodes and arthropods are the most speciose animal groups and possess Class 2 B1 G-protein coupled receptors (GPCRs). Existing models of invertebrate Class 2 B1 GPCR evolution are mainly centered on Caenorhabditis elegans and Drosophila melanogaster and a few other nematode and arthropod representatives. The present study reevaluates the evolution of metazoan Class 2 B1 GPCRs and orthologues by exploring the receptors in several nematode and arthropod genomes and comparing them to the human receptors. Three novel receptor phylogenetic clusters were identified and designated cluster A, cluster B and PDF-R-related cluster. Clusters A and B were identified in several nematode and arthropod genomes but were absent from D. melanogaster and Culicidae genomes, whereas the majority of the members of the PDF-R-related cluster were from nematodes. Cluster A receptors were nematode and arthropod-specific but shared a conserved gene environment with human receptor loci. Cluster B members were orthologous to human GCGR, PTHR and Secretin members with which they probably shared a common origin. PDF-R and PDF-R related clusters were present in representatives of both nematodes and arthropods. The results of comparative analysis of GPCR evolution and diversity in protostomes confirm previous notions that C. elegans and D. melanogaster genomes are not good representatives of nematode and arthropod phyla. We hypothesize that at least four ancestral Class 2 B1 genes emerged early in the metazoan radiation, which after the protostome-deuterostome split underwent distinct selective pressures that resulted in duplication and deletion events that originated the current Class 2 B1 GPCRs in nematode and arthropod genomes.

Mitochondrial genomes of Meloidogyne chitwoodi and M. incognita (Nematoda: Tylenchina): Comparative analysis, gene order and phylogenetic relationships with other nematodes

Root-knot nematodes (Meloidogyne spp.) are among the most important plant pathogens. In this study, the mitochondrial (mt) genomes of the root-knot nematodes, M. chitwoodi and M. incognita were sequenced. PCR analyses suggest that both mt genomes are circular, with an estimated size of 19.7 and 18.6–19.1kb, respectively. The mt genomes each contain a large non-coding region with tandem repeats and the control region. The mt gene arrangement of M. chitwoodi and M. incognita is unlike that of other nematodes. Sequence alignments of the two Meloidogyne mt genomes showed three translocations; two in transfer RNAs and one in cox2. Compared with other nematode mt genomes, the gene arrangement of M. chitwoodi and M. incognita was most similar to Pratylenchus vulnus. Phylogenetic analyses (Maximum Likelihood and Bayesian inference) were conducted using 78 complete mt genomes of diverse nematode species. Analyses based on nucleotides and amino acids of the 12 protein-coding mt genes showed strong support for the monophyly of class Chromadorea, but only amino acid-based analyses supported the monophyly of class Enoplea. The suborder Spirurina was not monophyletic in any of the phylogenetic analyses, contradicting the Clade III model, which groups Ascaridomorpha, Spiruromorpha and Oxyuridomorpha based on the small subunit ribosomal RNA gene. Importantly, comparisons of mt gene arrangement and tree-based methods placed Meloidogyne as sister taxa of Pratylenchus, a migratory plant endoparasitic nematode, and not with the sedentary endoparasitic Heterodera. Thus, comparative analyses of mt genomes suggest that sedentary endoparasitism in Meloidogyne and Heterodera is based on convergent evolution.

Signatures of adaptation to plant parasitism in nematode genomes.

Plant-parasitic nematodes cause considerable damage to global agriculture. The ability to parasitize plants is a derived character that appears to have independently emerged several times in the phylum Nematoda. Morphological convergence to feeding style has been observed, but whether this is emergent from molecular convergence is less obvious. To address this, we assess whether genomic signatures can be associated with plant parasitism by nematodes. In this review, we report genomic features and characteristics that appear to be common in plant-parasitic nematodes while absent or rare in animal parasites, predators or free-living species. Candidate horizontal acquisitions of parasitism genes have systematically been found in all plant-parasitic species investigated at the sequence level. Presence of peptides that mimic plant hormones also appears to be a trait of plant-parasitic species. Annotations of the few genomes of plant-parasitic nematodes available to date have revealed a set of apparently species-specific genes on every occasion. Effector genes, important for parasitism are frequently found among those species-specific genes, indicating poor overlap. Overall, nematodes appear to have developed convergent genomic solutions to adapt to plant parasitism.

Current Tropical Medicine Reports: A Path Forward to Highlight Research and Clinical Advances, New Trends and Innovations

Welcome to the first edition of Current Tropical Medicine Reports! This is an exciting time in tropical medicine; the ongoing advancements in basic and applied science to clinical research have now positioned the discipline beyond what the giants in our field may have expected in their wildest dreams. Manson, Gorgas, and Kean would be impressed with the audacious milestones that the scientific community has completed: from having complete sequences of nematode genomes to new innovations for the development of vaccines and drugs against pathogens that have plagued mankind for eons. The goal of this journal is to enhance the understanding of tropical medicine by delving into the latest advances in the wide expanse of the world of biosciences. This journal seeks to access a varied audience of readers, from the established principal investigators wanting to gain new ideas for basic or applied research to the first-year medicine resident or fellow wanting to learn more about these global diseases and apply that knowledge to the care of their patients. New issues will be released on a quarterly basis, and will focus on current highlights, hot topics, and a review of recent advances in the study of pathogens that affect morbidity and mortality worldwide. The topics will include virology, bacteriology, parasitology, mycology, entomology, immunology, cellular and molecular biology, epidemiology, and more; the journal will address the latest in diagnostics, treatment, and prevention of these global diseases and pathogens. To add variety, the quarterly selections will include one or more organisms from a pathogenesis group and describe the latest preventive and therapeutic strategies, evolving trends, and important developments in their burden of disease, the geographical distribution, and advances in research & development, all the way to clinical trials and clinical management. The articles will focus on work from the past three years and the advancements that are shaping the way we study and practice medicine. The section editors and authors have been selected from a cadre of scientists and physicians who are leaders in their field. They bring years of accumulated knowledge and research, but also their expert opinion and personal experience in their field. In addition the journal has compiled a diverse group of global leaders to participate as the editorial board. This inaugural edition covers viral tropical medicine with a focus onDengue virus, from emergence in the Americas to the regulation of its replication. Also included is a review of the helminth Toxocara, from development of molecular assays to clinical presentations of toxocariasis. As recently assessed in the Global Burden of Disease Study 2010, tropical diseases continue to have huge contributions towards the world’s morbidity (1). Furthermore, the notion that tropical medicine is limited to tropical climates or resource-limited areas in developing countries is now changing. Active transmission of these R. Mejia :M. E. Bottazzi (*) Sabin Vaccine Institute and Texas Children’s Hospital Center for Vaccine Development, National School of TropicalMedicine, Baylor College of Medicine, Houston, TX, USA e-mail: Bottazzi@bcm.edu

Evolutionarily dynamic roles of a PUF RNA‐binding protein in the somatic development of Caenorhabditis briggsae

Gene duplication and divergence has emerged as an important aspect of developmental evolution. The genomes of Caenorhabditis nematodes encode an ancient family of PUF RNA-binding proteins. Most have been implicated in germline development, and are often redundant with paralogs of the same sub-family. An exception is Cbr-puf-2 (one of three Caenorhabditis briggsae PUF-2 sub-family paralogs), which is required for development past the second larval stage. Here, we provide a detailed functional characterization of Cbr-puf-2. The larval arrest of Cbr-puf-2 mutant animals is caused by inefficient breakdown of bacterial food, which leads to starvation. Cbr-puf-2 is required for the normal grinding cycle of the muscular terminal bulb during early larval stages, and is transiently expressed in this tissue. In addition, rescue of larval arrest reveals that Cbr-puf-2 also promotes normal vulval development. It is expressed in the anchor cell (which induces vulval fate) and vulval muscles, but not in the vulva precursor cells (VPCs) themselves. This contrasts with the VPC-autonomous repression of vulval development described for the Caenorhabditis elegans homologs fbf-1/2. These different roles for PUF proteins occur even as the vulva and pharynx maintain highly conserved anatomies across Caenorhabditis, indicating pervasive developmental system drift (DSD). Because Cbr-PUF-2 shares RNA-binding specificity with its paralogs and with C. elegans FBF, we suggest that functional novelty of RNA-binding proteins evolves through changes in the site of their expression, perhaps in concert with cis-regulatory evolution in target mRNAs.

Identification of Novel Target Genes for Safer and More Specific Control of Root-Knot Nematodes from a Pan-Genome Mining

Root-knot nematodes are globally the most aggressive and damaging plant-parasitic nematodes. Chemical nematicides have so far constituted the most efficient control measures against these agricultural pests. Because of their toxicity for the environment and danger for human health, these nematicides have now been banned from use. Consequently, new and more specific control means, safe for the environment and human health, are urgently needed to avoid worldwide proliferation of these devastating plant-parasites. Mining the genomes of root-knot nematodes through an evolutionary and comparative genomics approach, we identified and analyzed 15,952 nematode genes conserved in genomes of plant-damaging species but absent from non target genomes of chordates, plants, annelids, insect pollinators and mollusks. Functional annotation of the corresponding proteins revealed a relative abundance of putative transcription factors in this parasite-specific set compared to whole proteomes of root-knot nematodes. This may point to important and specific regulators of genes involved in parasitism. Because these nematodes are known to secrete effector proteins in planta, essential for parasitism, we searched and identified 993 such effector-like proteins absent from non-target species. Aiming at identifying novel targets for the development of future control methods, we biologically tested the effect of inactivation of the corresponding genes through RNA interference. A total of 15 novel effector-like proteins and one putative transcription factor compatible with the design of siRNAs were present as non-redundant genes and had transcriptional support in the model root-knot nematode Meloidogyne incognita. Infestation assays with siRNA-treated M. incognita on tomato plants showed significant and reproducible reduction of the infestation for 12 of the 16 tested genes compared to control nematodes. These 12 novel genes, showing efficient reduction of parasitism when silenced, constitute promising targets for the development of more specific and safer control means.

Intra-genomic variation in the ribosomal repeats of nematodes.

Ribosomal loci represent a major tool for investigating environmental diversity and community structure via high-throughput marker gene studies of eukaryotes (e.g. 18S rRNA). Since the estimation of species’ abundance is a major goal of environmental studies (by counting numbers of sequences), understanding the patterns of rRNA copy number across species will be critical for informing such high-throughput approaches. Such knowledge is critical, given that ribosomal RNA genes exist within multi-copy repeated arrays in a genome. Here we measured the repeat copy number for six nematode species by mapping the sequences from whole genome shotgun libraries against reference sequences for their rRNA repeat. This revealed a 6-fold variation in repeat copy number amongst taxa investigated, with levels of intragenomic variation ranging from 56 to 323 copies of the rRNA array. By applying the same approach to four C. elegans mutation accumulation lines propagated by repeated bottlenecking for an average of ~400 generations, we find on average a 2-fold increase in repeat copy number (rate of increase in rRNA estimated at 0.0285-0.3414 copies per generation), suggesting that rRNA repeat copy number is subject to selection. Within each Caenorhabditis species, the majority of intragenomic variation found across the rRNA repeat was observed within gene regions (18S, 28S, 5.8S), suggesting that such intragenomic variation is not a product of selection for rRNA coding function. We find that the dramatic variation in repeat copy number among these six nematode genomes would limit the use of rRNA in estimates of organismal abundance. In addition, the unique pattern of variation within a single genome was uncorrelated with patterns of divergence between species, reflecting a strong signature of natural selection for rRNA function. A better understanding of the factors that control or affect copy number in these arrays, as well as their rates and patterns of evolution, will be critical for informing estimates of global biodiversity.

Komplexität im Kleinen — Nematoden-Genome im Vergleich

The genome of Caenorhabditis eleganswas the first to be sequenced from a multicellular organism. One major finding was that despite a much smaller genome, C. eleganshas only marginally fewer genes than man. Sequencing of other nematode genomes demonstrated that this is a general feature of nematodes. With their high gene turnover and numerous examples of horizontal gene transfer, nematode genomes form a highly attractive system for comparative studies of genome structure and evolution.

Research progress of conserved non-coding elements in metazoan

In addition to protein coding sequence, the organism genome contains a significant amount of regulatory DNA. Comparative genomics reveals that the organism genomes of vertebrates, tunicate, cephalochordate, flies, and nematodes contain cis-regulatory elements with highly conserved non-coding elements (CNEs). CNEs that cluster around trans-dev genes are part of core gene regulatory networks (GRNs), and usually, they can act as transcriptional enhancers. In this review, we described the identification of CNEs and summarized their key properties across the metazoans, and then discussed the evolution of CNEs after large-scale genome duplication events and the role of CNEs in the evolution of ani-mal body plan.

Badger—an accessible genome exploration environment

Summary: High-quality draft genomes are now easy to generate, as sequencing and assembly costs have dropped dramatically. However, building a user-friendly searchable Web site and database for a newly annotated genome is not straightforward. Here we present Badger, a lightweight and easy-to-install genome exploration environment designed for next generation non-model organism genomes.Availability: Badger is released under the GPL and is available at http://badger.bio.ed.ac.uk/. We show two working examples: (i) a test dataset included with the source code, and (ii) a collection of four filarial nematode genomes.Contact: mark.blaxter@ed.ac.uk

A Lover and a Fighter: The Genome Sequence of an Entomopathogenic Nematode Heterorhabditis bacteriophora

Heterorhabditis bacteriophora are entomopathogenic nematodes that have evolved a mutualism with Photorhabdus luminescens bacteria to function as highly virulent insect pathogens. The nematode provides a safe harbor for intestinal symbionts in soil and delivers the symbiotic bacteria into the insect blood. The symbiont provides virulence and toxins, metabolites essential for nematode reproduction, and antibiotic preservation of the insect cadaver. Approximately half of the 21,250 putative protein coding genes identified in the 77 Mbp high quality draft H. bacteriophora genome sequence were novel proteins of unknown function lacking homologs in Caenorhabditis elegans or any other sequenced organisms. Similarly, 317 of the 603 predicted secreted proteins are novel with unknown function in addition to 19 putative peptidases, 9 peptidase inhibitors and 7 C-type lectins that may function in interactions with insect hosts or bacterial symbionts. The 134 proteins contained mariner transposase domains, of which there are none in C. elegans, suggesting an invasion and expansion of mariner transposons in H. bacteriophora. Fewer Kyoto Encyclopedia of Genes and Genomes Orthologies in almost all metabolic categories were detected in the genome compared with 9 other sequenced nematode genomes, which may reflect dependence on the symbiont or insect host for these functions. The H. bacteriophora genome sequence will greatly facilitate genetics, genomics and evolutionary studies to gain fundamental knowledge of nematode parasitism and mutualism. It also elevates the utility of H. bacteriophora as a bridge species between vertebrate parasitic nematodes and the C. elegans model.