The genome of Romanomermis culicivorax: revealing fundamental changes in the core developmental genetic toolkit in Nematoda
- Philipp H Schiffer1Email author,
- Michael Kroiher1,
- Christopher Kraus1,
- Georgios D Koutsovoulos2,
- Sujai Kumar2,
- Julia I R Camps1,
- Ndifon A Nsah1,
- Dominik Stappert3,
- Krystalynne Morris4,
- Peter Heger1,
- Janine Altmüller5,
- Peter Frommolt5,
- Peter Nürnberg5,
- W Kelley Thomas4,
- Mark L Blaxter2 and
- Einhard Schierenberg1
© Schiffer et al.; licensee BioMed Central Ltd. 2013
Received: 7 May 2013
Accepted: 17 December 2013
Published: 27 December 2013
The genetics of development in the nematode Caenorhabditis elegans has been described in exquisite detail. The phylum Nematoda has two classes: Chromadorea (which includes C. elegans) and the Enoplea. While the development of many chromadorean species resembles closely that of C. elegans, enoplean nematodes show markedly different patterns of early cell division and cell fate assignment. Embryogenesis of the enoplean Romanomermis culicivorax has been studied in detail, but the genetic circuitry underpinning development in this species has not been explored.
We generated a draft genome for R. culicivorax and compared its gene content with that of C. elegans, a second enoplean, the vertebrate parasite Trichinella spiralis, and a representative arthropod, Tribolium castaneum. This comparison revealed that R. culicivorax has retained components of the conserved ecdysozoan developmental gene toolkit lost in C. elegans. T. spiralis has independently lost even more of this toolkit than has C. elegans. However, the C. elegans toolkit is not simply depauperate, as many novel genes essential for embryogenesis in C. elegans are not found in, or have only extremely divergent homologues in R. culicivorax and T. spiralis. Our data imply fundamental differences in the genetic programmes not only for early cell specification but also others such as vulva formation and sex determination.
Despite the apparent morphological conservatism, major differences in the molecular logic of development have evolved within the phylum Nematoda. R. culicivorax serves as a tractable system to contrast C. elegans and understand how divergent genomic and thus regulatory backgrounds nevertheless generate a conserved phenotype. The R. culicivorax draft genome will promote use of this species as a research model.
KeywordsNematode Genome Evolution Development Caenorhabditis Mermithida Romanomermis
Nematodes have a generally conserved body plan. Their typical form is dictated by the presence of a single-chamber hydroskeleton, where longitudinal muscles act against an inextensible extracellular cuticle. The conservation of organ systems between nematode species is even more striking, with, for example, the nervous system, the somatic gonad and the vulva having very similar general organisations and cellular morphologies across the phylum. It might be thought that these similarities arise from highly stereotypical developmental programmes, but comparative studies challenge this “all nematodes are equal” view.
Regulative development is common among Metazoa, and is also observed in other Ecdysozoa, including Arthropoda. Indeed, in several enoplean species, early embryos have been found to not display polarised early divisions, arguing against a strongly determined mode of development in this group [13, 14]. The determined mode found in C. elegans is thus likely to be derived even within Nematoda , implying that the core developmental system in Nematoda has changed, while maintaining a very similar organismal output. This phenomenon, termed “developmental system drift” , reveals independent selection on the developmental mechanism and the final form produced.
To explore the genetics of development of enoplean and other non-rhabditid nematodes requires tractable experimental systems with a suitable set of methodological tools and extensive genomic data. While C. elegans and its embryos are relatively easily manipulated and observed, and the C. elegans genome has been fully sequenced , embryos from the Enoplia and Dorylaimia are much harder to culture and manipulate. Few viable laboratory cultures exist and obtaining large numbers of embryos from wild material is difficult. Functional molecular analysis of most nematodes, in particular Enoplea, is further hindered by the lack of genetic tools such as mutant analysis or gene-knockdown via RNAi. Performing detailed comparative experimental embryology on a phylogenetically representative set of species across the phylum Nematoda thus remains a distant goal.
The genetic toolkit utilised by a species is represented in its genome, and direct assessment of the genetic capabilities of an organism can thus be assessed through analysis of genome data. Using the background knowledge of pathways and modules used in other taxa, the underpinning logic of a species’ developmental system can be inferred from its genome, and the developmental toolkits of different species can be compared. These comparisons can reveal changes in developmental logic between taxa by identifying gene losses during evolution that must result in changed pathway functioning, and similarly identify genes recruited to developmental regulatory roles in particular lineages.
Efficient generation of genomic resources for non-model species, and the inference of developmental regulatory pathways from the encoded gene sets, is now possible. The majority of the fifteen nematode genomes published to date have been from Rhabditida (Figure 1) [18–26]. The single enoplean genome sequences is from the mammalian parasite Trichinella spiralis (Dorylaimia; order Trichocephalida) . T. spiralis is ovoviviparous, proper development requires intrauterine environment, and early blastomeres are extremely transparent  such that individual nuclei are hard to identify (E.S., unpublished observations). Hence, this species is of very limited value for light microscopical image analysis and experimental investigation correlating cell dynamics with the molecular circuitry regulating early development.
Although the genomes of many additional nematode species are being sequenced [29, 30], even in this wider sampling of the phylum, Enoplea remains neglected. The enoplean Romanomermis culicivorax (order Mermithida within Dorylaimia) has been established in culture for decades. It infects and kills the larvae of many different mosquito species , and is being investigated for its potential as a biocontrol agent of malaria and other disease vectors [31, 32]. R. culicivorax and T. spiralis differ fundamentally in many life-cycle and phenotypic characters. R. culicivorax reproduces sexually. A single female can produce more than a thousand eggs, and embryos are easily studied under laboratory conditions. They display a developmental pattern that differs markedly from C. elegans. As in other Enoplea [14, 33] the first division is equal, and not asymmetric as in C. elegans. R. culicivorax also shows an inversion of dorso-ventral axis polarity compared to C. elegans, while a predominantly monoclonal fate distribution indicates fewer modifying inductions between blastomeres [33, 34]. Generation of the hypodermis involves repetitive cell elements extending from posterior to anterior over the remainder of the embryo, a process distinct from that observed in C. elegans.
We here catalogue the R. culicivorax developmental toolkit derived from annotation of a draft genome sequence. We contrast genes and proteins identified in R. culicivorax and T. spiralis with those of C. elegans, and other Ecdysozoa, represented by the arthropod Tribolium castaneum. We conclude that major changes in the regulatory logic of development have taken place during nematode evolution, possibly as a consequence of developmental system drift, and that the model species C. elegans is considerably derived compared to an ecdysozoan (and possibly metazoan) ground system. However, we are still able to define conserved gene sets that may act in “phylotypic” developmental stages.
Results and discussion
Romanomermis culicivoraxhas a large and repetitive genome
A draft genome assembly for R. culicivorax was generated from 26.9 gigabases (Gb) of raw data (filtered from a total of 41 Gb sequenced; Additional file 1: Table S1). The assembly has a contig span of 267 million base pairs (Mb) and a scaffold span of 323 Mb. The 52 Mb of spanned gaps are likely inflated estimates derived from use of the SSPACE scaffolder. We do not currently have a validated independent estimate of genome size for R. culicivorax, but preliminary measurements with Feulgen densitometry suggest a size greater than 320 Mb (Elizabeth Martínez Salazar pers. comm.). The R. culicivorax genome is thus three times bigger than that of C. elegans, and five times that of T. spiralis (Table 1). The assembly is currently in 62,537 scaffolds and contigs larger than 500 bp, with an N50 of 17.6 kilobases (kb). The N50 for scaffolds larger than 10 kb is 29.9 kb, and the largest scaffold is over 200 kb. The GC content is 36%, comparable to 38% of C. elegans and 34% in T. spiralis. We identified 47% of the R. culicivorax genome as repetitive. To validate this estimate we applied our repeat-finding approach to previously published genomes and achieved good accordance with these data (Table 1). The non-repetitive content of the R. culicivorax genome is thus approximately twice that of C. elegans and three times that of T. spiralis. T. spiralis thus stands out as having the least complex nematode genome sequenced so far, and the contrast with R. culicivorax indicates that small genomes are not characteristic of Dorylaimia.
Assembly and annotation statistics
Contigs >100bp span
Scaffolds >500bp span
N50 contigs/scaffolds >500bp
N50 scaffolds >500bp
Max contig length
Max scaffold length
Mean transcript length
Mean protein length
MAKER AUGUSTUS predictions
MAKER SNAP predictions
Num. ESTs (isogroups)
Mean EST length
80% BLAT sequence coverage
CEGMA compl. completeness
CEGMA Group 1 part. compl.
CEGMA Group 2 part. compl.
CEGMA Group 3 part. compl.
CEGMA Group 4 part. compl.
We annotated 1,443 tRNAs in the R. culicivorax genome using INFERNAL  and tRNAscan-SE , of which 382 were pseudogenes (see Additional file 1: Table S2 for details). In comparison, T. spiralis has 134 tRNAs of which 7 are pseudogenes, while C. elegans has 606 tRNAs with 36 pseudogenes . Threonine (Thr) tRNAs were particularly overrepresented (676 copies), a finding echoed in the genomes of Meloidogyne incognita and Meloidogyne floridensis[24, 43] and in P. pacificus. The latter has also an overrepresentation of Arginine tRNAs .
We have made available the annotated R. culicivorax genome, with functional categorisations of predicted genes and proteins and annotation features, in a dedicated genome browser at http://romanomermis.nematod.es.
The R.culicivorax gene set is more representative of Dorylaimia than T. spiralis
The phylogenetic placement of R. culicivorax makes its genome attractive for exploring the likely genetic complexity of an ancestral nematode. With T. spiralis, it can be used to reveal the idiosyncrasies of the several genomes available for Rhabditida. To polarise this comparison, we used the arthropod Tribolium castaneum, for which a high quality genome sequence is available . T. castaneum development is considered less derived than that of the major arthropod model Drosophila melanogaster. The OrthoMCL pipeline accurately clusters orthologous proteins, facilitating the complex task of grouping proteins that are likely to share biological function in divergent organisms , and performs better than approaches that simply use domain presence information or aggregative approaches such as psiBLAST . We used the OrthoMCL pipeline to generate a set of protein clusters for the four species (R. culicivorax, T. spiralis, C. elegans and T. castaneum). While the large divergence between these species may obscure relationships between protein sequences, making inference of orthology problematic [48–50], the parameters used were most inclusive [50–52]. Additionally, as the R. culicivorax genome assembly may not be complete we based inference of absence on shared loss in both R. culicivorax and T. spiralis. Additionally, we validated inferences of absence from the OrthoMCL analyses by performing detailed sequence comparisons using BLAST+  (Additional file 2).
Clusters containing only R. culicivorax and T. spiralis proteins might identify functions distinct to these dorylaim nematodes. In the 461 T. spiralis and 806 R. culicivorax proteins in these clusters, a total of 65 GO terms were found to be overrepresented (single test p <0.05 by Fisher’s exact test). While C. elegans has a reduced ability to methylate DNA , we found four methylation-associated GO terms among the 64 overrepresented. We also detected significant enrichment (single test p <0.05) for GO terms describing chromatin and DNA methylation functions in the set of R. culicivorax proteins that lacked homologues in C. elegans (see Additional file 3). Important roles for methylation and changes in methylation patterns during T. spiralis development have been inferred from transcriptional profiling . Methylation is important for the silencing of transposable elements [57, 58] and could play a crucial role in the highly repetitive R. culicivorax genome.
The clusters that contained R. culicivorax, T. spiralis and T. castaneum proteins but no C. elegans orthologues might contain proteins involved in core ecdysozoan processes lost in C. elegans. In these clusters we identified 40 GO terms overrepresented (single test p <0.05) compared to the C. elegans proteome (see Additional file 3). Some of these GO terms were linked to chromatin remodelling and methylation (e.g. Ino80 complex, histone arginine methylation). Other overrepresented GO terms were related to cell signalling (the Wnt receptor pathway; the C. elegans Wnt signalling system is distinct from other metazoa ), and ecdysone receptor holocomplex (potentially a basic ecdysozoan function ).
The genetic background of development in R. culicivorax and T. spiralis differs markedly from that of C. elegans
In a recent multi-species developmental time course expression analysis within several Caenorhabditis species, conserved sets of genes were found to have conserved patterns of differential expression in discrete phases in the timeline from zygote to the hatching larva .
Core developmental pathways differ between nematodes
There are important differences in cell behaviour during early embryogenesis between R. culicivorax and C. elegans[33, 34]. We used the genomic data to follow up on some of the striking contrasts between the dorylaim and the rhabditid patterns of development: establishment of primary axis polarity, segregation of maternal message within the early embryo, hypodermis formation, the vulval specification pathways, epigenetic pathways (especially DNA methylation), sex determination and light sensing (see Additional file 1).
Presence and absence of selected ∗ C. elegans proteins in Dorylaimia
The hypodermis in C. elegans is derived from specific descendants of the anterior and posterior founder cells . In contrast, in R. culicivorax hypodermis is derived from descendants of a single cell . Several C. elegans genes expressed in the hypodermis or associated with hypodermal development were absent from R. culicivorax and T. spiralis (see Table 3 and Additional file 3). For example the GATA-like transcription factors ELT-1 and ELT-3 act redundantly in C. elegans. ELT-3 was absent from the dorylaim species, but ELT-1 was conserved in R. culicivorax, T. spiralis and T. castaneum. Thus, ELT-3 appears to be an innovation in the rhabditid lineage, suggesting changes of interaction complexity during nematode evolution.
In C. elegans, vulva formation is highly dependent on initial inductive signals from the anchor cell that activate a complex gene regulatory network, which drives tissue specific cell division and differentiation. The evolutionary plasticity of this system has been explored in rhabditid nematodes, revealing the changing relative importance of cell-cell interactions, inductions, and lineage-autonomous specifications [90, 91]. The signal transduction pathways include a RTK/RAS/MAPK cascade, activated by EGF- and wnt-signalling . Among the downstream targets in C. elegans are for example LIN-1 and the β-catenin BAR-1, which in turn regulates the HOX-5 orthologue LIN-39 [93–95]. These important regulators of vulva development are completely absent from the genomes of R. culicivorax and T. spiralis (Table 3 and Additional file 2). We identified a R. culicivorax protein with low similarity to C. elegans BAR-1 (24% sequence identity). However, this protein is not clustered with other dorylaim proteins, and appears to be either a duplication of the β-catenin ortholog HMP-2 or another armadillo repeat-containing protein rather than an orthologue of BAR-1 (see Additional file 5). These shared patterns of absence again suggest that similar morphological structures can be generated with very different genetic underpinnings. Vulva formation in the dorylaims may be regulated without the BAR-1 – LIN-39 interaction, as observed in P. pacificus. In C. elegans Hox gene expression is cell-lineage dependent [97, 98], organised so that the cells that express specific Hox genes are clustered along the anterior-posterior axis (see e.g. ). It will be informative to test whether in R. culicivorax and other non-rhabditid nematodes Hox genes act in an axis position-dependent, but cell lineage-independent manner, as observed in many other animals, notably arthropods [100, 101]. Epigenetic regulation is key to developmental processes in many animals, but its roles in C. elegans are more muted (see above). Notably C. elegans is depleted for chromatin re-modelling genes of the Polycomb and Trithorax groups . It is intriguing that we found orthologues of T. castaneum pleiohomeotic in R. culicivorax and T. spiralis, and orthologs of T. castaneum trithorax and Sex comb on midleg (Scm) in R. culicivorax. This suggests that dorylaim chromatin restructuring mechanisms may be more arthropod-like than in C. elegans. The presence of an intact methylation machinery and conserved chromatin re-modelling factors opens the prospects for a role for epigenetic modification in developmental regulation of dorylaim nematodes.
Defining a set of potential phylotypic stage genes
While the examples above demonstrate considerable developmental system drift in Nematoda, we also identified many sets of orthologous proteins conserved between Dorylaimia and C. elegans. We asked if these could be correlated with functions in distinct developmental phases with a conserved phenotype. Shortly before the start of morphogenesis, at the point of ventral enclosure, nematode embryos from Chromadorea and Enoplea share a similar morphology . Levin et al.  found that in five Caenorhabditis species a distinct set of genes had elevated expression around the ventral enclosure stage (their stage 7) (Figure 3) and proposed that this constitutes a “phylotypic stage" for nematodes. We used T. spiralis and R. culicivorax gene sets to refine and restrict this set of phylotypic stage genes. Of the 834 C. elegans genes with elevated expression between stages 6 to 8 , 355 had no orthologue in R. culicivorax, T. spiralis or T. castaneum. The remaining 479 phylotypic stage candidates from C. elegans were present in 279 of our OrthoMCL clusters. Of these clusters 93 were nematode-restricted containing 186 C. elegans proteins grouped with 129 R. culicivorax and 113 T. spiralis homologues. The remaining 186 clusters were part of the conserved ecdysozoan core proteome (see above) and contained 330 C. elegans proteins together with 248 R. culicivorax, 248 T. spiralis and 621 T. castaneum proteins (Figure 3; The total number of C. elegans candidates is larger than 479 due to the inclusion of co-orthologues in this species). In the set of phylotypic stage genes identified by Levin et al.  are proteins functioning in processes such as muscle and neuron formation, signalling between cells, and morphogenesis. This pattern was retained in the conserved clusters (see Additional file 5). Although time-resolved expression data will be needed to confirm the activity of these genes in developmental stages of R. culicivorax, their retention in the Dorylaimia supports their general importance. We can now sub-classify the set of conserved proteins expressed at the potential nematode phylotypic stage. A first, nematode-restricted set includes many proteins that are important for cuticle formation (e.g. collagen proteins) and some hedgehog-like proteins, expressed in the C. elegans hypodermis . As cuticle formation follows ventral enclosure in nematodes, these proteins may be involved in this nematode-specific function. The second set, comprising clusters conserved between the nematodes and T. castaneum, contains many important developmental transcription factors, such as the Hox gene mab-5, other homeobox genes, and helix-loop-helix and C2H2-type zinc finger transcription factors. This second set may represent a genetic backbone driving formation of phylotypic stage in diverse animal taxa, in accordance with the recent extension of the concept to Metazoa [104–106].
To be useful as a contrasting system to the canonical C. elegans model, any nematode species must be accessible to both descriptive and manipulative investigation. The reference genome for R. culicivorax lays bare the core machinery available for developmental regulation, and we have demonstrated that in situ hybridisation approaches are feasible for this species. Along with the long established, robust laboratory cultures, this makes R. culicivorax an attractive and tractable alternative model for understanding the evolutionary dynamics of nematode development. By combining the R. culicivorax genome with that of T. spiralis, we have been able to explore the molecular diversity of Dorylaimia, and provide robust contrasts with the intensively studied Rhabditida. Particularly surprising are the differences between R. culicivorax and T. spiralis. The R. culicivorax genome is much larger than that of T. spiralis, and contained a high proportion of repetitive sequence, including many transposable elements. Despite the phylogenetic and lifestyle affinities between the two dorylaims compared to C. elegans, the R. culicivorax genome retained many more genes in common with C. elegans than did T. spiralis. We suggest that T. spiralis may be an atypical representative of dorylaim nematodes, perhaps due to its highly derived life cycle.
Our analyses identified many genes apparently absent from the dorylaim genomes, despite relaxed analysis parameters. In particular, for genes identified as critical to C. elegans development but apparently absent from the dorylaims, we were unable to identify credible orthologues using sensitive search strategies. In this phylum-spanning comparison, inferences of gene orthology can be obscured by levels of divergence. In addition, the gene family birth rate in the chromadorean lineage leading to C. elegans is high [25, 27], and therefore C. elegans was expected to have many genes absent from the dorylaim species. Thus, we might not have found a R. culicivorax orthologue for a specific gene for three reasons: it may have arisen in the branch leading to C. elegans; its sequence divergence may be too great to permit clustering with potential homologs; or it was not assembled in the draft dorylaim genomes. The case of C. elegans PAR-3 and D. melanogaster bazooka illustrate some of these difficulties: the possible R. culicivorax orthologue was highly divergent. Whether or not we have been able to identify all the orthologues of the key C. elegans genes present in the R. culicivorax and T. spiralis genomes, the absence of an identified orthologue maximally implies loss from the genome, and minimally implies significant sequence divergence. In the latter case this would most likely cause changes in the networks and pathways in which genes interact to deliver biological function.
Between the model nematode C. elegans and arthropod models such as T. castaneum many key mechanisms governing early cell patterning are divergent . Our data strongly support the view that major variation also exists within Nematoda. T. spiralis and R. culicivorax both lack orthologues of genes involved in core developmental processes in C. elegans, and many of these C. elegans genes appear to be restricted to the Rhabditida. It is thus doubtful that these processes are regulated by same molecular interactions across the phylum. We suggest that developmental system drift has played a major role in nematode evolution. The phenotypic conservatism associated with the vermiform morphology of nematodes  has fostered unjustified expectations concerning the conservation of genetic programmes that determine these morphologies. Despite this divergence in developmental systems, we were able to define two sets of conserved genes possibly active in a taxon-specific phase of ventral enclosure and cuticle formation in Nematoda, and in a potential phylotypic stage of Ecdysozoa. The advent of robust, affordable and rapid genome sequencing also opens the vista of large-scale comparative genomics of development across the phylum Nematoda  to better understand the diversity of the phylum and also place the remarkable C. elegans model in context of its peers. It will next be necessary to extend these studies to a broader sampling of developmental pathway genes from a wider and representative sampling of nematode genomes across the full diversity of the phylum. We have highlighted a few of the possible avenues a research programme could follow: early axis formation and polarisation, the specification of hypodermis, sex determination, vulva formation, the roles of epigenetic processes in developmental regulation and the confirmation of potential “phylotypic stage genes” with expression analysis in R. culicivorax.
Sequencing and genome assembly
Genomic DNA was extracted from several hundred, mixed-sex, adult R. culicivorax specimens from a culture first established in Ed Platzer’s laboratory in Riverside, California. Illumina paired end and mate pair sequencing with libraries of varying insert sizes, and Roche 454 single end sequencing, was performed at the Cologne Center for Genomics (CCG: http://www.ccg.uni-koeln.de). A Roche 454 dataset of transcriptome reads from cDNA synthesised from mixed developmental stages and sexes was also generated (see Additional file 1: Table S1 for details of data generated).
The quality of the raw data was assessed with FastQC (v.0.9; http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Adapter sequences and low quality data were trimmed from the Illumina paired end data with custom scripts (see http://github.com/sujaikumar/assemblage) and from the mate pair libraries with Cutadapt (v.1.0) . We constructed a preliminary genome assembly, with relaxed insert size parameters, from the paired end Illumina libraries with the de-novo-assemble option of the clcAssemblyCell (v.4.03b) . We validated the actual insert sizes of our libraries by mapping back the reads to this preliminary assembly using clcAssemblyCell. The preliminary assembly was also used to screen out bacterial and other contaminant data . The transcriptome data were assembled with Roche GSAssembler (Newbler; version 2.5). For the production assembly, we explored assembly parameters using different mixes of our data, evaluating each for total span, maximal contig lengths, N50, number of contigs, representation of the transcriptome, and conserved eukaryotic gene content (using the CEGMA pipeline v.2.1 ). The most promising assembly was scaffolded with the filtered Illumina mate pair read sets using SSPACE (v.1.2) . As our genomic DNA derived from a population of nematodes of unknown genetic diversity, we removed short contigs that mapped entirely within larger ones using Cd-hit (v.4.5.7)  at a 95% cutoff. A final round of superscaffolding was performed, linking scaffolds that had logically consistent matches to the transcriptome data based on BLAT  hits and processed with SCUBAT (B. Elsworth, pers. comm.; http://github.com/elswob/SCUBAT). The final genome assembly was again assessed for completeness by assessing the mapping of the transcriptome contigs and with the CEGMA pipeline .
RepeatMasker (v.3.3.0) [112, 113], RepeatFinder  and RepeatModeler (v.1.0.5; http://www.repeatmasker.org/RepeatModeler.html; combining RECON (v.1.07)  and RepeatScout (v.1.05) ), were used to identify known and novel repetitive elements in the R. culicivorax genome. We employed the MAKER pipeline to find genes in the R. culicivorax genome assembly. In a first pass, the SNAP gene predictor included in MAKER was trained with a CEGMA  derived output of predicted highly conserved genes. As additional evidence we included the transcriptome assembly and a set of approximately 15,000 conserved nematode proteins derived from the NEMBASE4 database (recalculated by J. Parkinson; pers. comm.). In the second, definitive, pass we used the gene set derived from this first MAKER iteration to train AUGUSTUS  inside the MAKER pipeline for a second run, also including evidence from transcriptome to genome mapping obtained with GenomeThreader . Codon usage in R. culicivorax, T. spiralis, and C. elegans was calculated using INCA (v2.1) . Results were then compared to data from  (see Additional files 1 and 6).
We used Blast2GO (Blast2GO4Pipe, v.2.5, January 2012 database issue)  to annotate the gene set with Gene Ontology terms , based on BLAST matches with expect values less than 1e -5 to the UniProt/SwissProt database (March 2012 snapshot), and domain annotations derived from the InterPro database . Comparison of annotations between three nematode species (R. culicivorax, C. elegans, and T. spiralis) and, as a reference outgroup, the holometabolous coleopteran arthropod Tribolium castaneum, was based on GO Slim data retrieved with Blast2GO. RNA genes were predicted using INFERNAL (v.1.0.2) and the Rfam database , and tRNAscan-SE (v.1.3.1) .
We inferred clusters of orthologous proteins between R. culicivorax, T. spiralis, and C. elegans, and the beetle T. castaneum using OrthoMCL (v.2.0.3) . T. spiralis, C. elegans and T. castaneum protein sets were downloaded from NCBI and WormBase (see Additional file 1: Table S3) and redundancy screened with Cd-hit at the 99% threshold. We selected an inflation parameter of 1.5 for MCL clustering (based on [126, 127]) within OrthoMCL to generate an inclusive clusterings in our analysis likely to contain even highly diverged representatives from the four species. In analyses of selected developmental genes, clusters were manually validated using NCBI-BLAST+ . We affirmed the uniqueness of C. elegans proteins identified as lacking homologues in the enoplean nematodes by comparing them to the R. culicivorax proteome using BLAST. Those with no significant matches at all (all matches with E-values > 1e -5) were classified as confirmed absent. Those having matches with E-values < 1e -5 were investigated further by surveying the cluster memberships of the R. culicivorax matches. If the R. culicivorax protein was found to cluster with a different C. elegans protein, the uniqueness to C. elegans was again confirmed. If the R. culicivorax protein did not cluster with an alternative C. elegans protein, we reviewed the BLAST statistics (E-value, identity and sequence coverage) of the match and searched the GenBank non redundant protein database for additional evidence of possible orthology. Only if these tests yielded no indication of direct orthology was the C. elegans protein designated absent from the enoplean set. Further details of the process are given in Additional file 5.
We identified the protein sequences of 1,725 genes differentially expressed in C. elegans developmental stages  and selected, using our OrthoMCL clustering, those apparently lacking orthologues in R. culicivorax and T. spiralis (verified as above). Using Wormbase (http://www.wormbase.org, release WS233) we surveyed the C. elegans-restricted genes for their experimentally-defined roles in development.
Custom Perl scripts were used to group orthoMCL clusters on the basis of species membership patterns. The sets of clusters that contained (i) both T. spiralis and R. culicivorax members but no C. elegans members and (ii) T. spiralis and R. culicivorax and T. castaneum members but no C. elegans members were surveyed for GO annotations enriched in comparison to the whole C. elegans proteome (sets i and ii) and the T. castaneum proteome (set i), conducting Fisher’s exact test as implemented in Blast2GO. Due to the small size of both sets compared to the large reference set, p-values could not be corrected for multiple testing. To improve annotation reliability, these proteins were recompared (using BLAST) to the UniProt/SwissProt database and run through the Blast2GO pipeline as described above.
Whole-mount in situ hybridization
For in situ hybridisation we modified the freeze-crack procedure described previously for C. elegans and revised by Maduro et al. (2007; http://www.faculty.ucr.edu/~mmaduro/resources.htm). In particular, to achieve reliable penetration of the durable R. culicivorax egg envelopes we initially partly removed the protective layer by incubation in alkaline bleach solution (see ). Digoxygenine-labeled sense and antisense RNA probes were generated from linearized pBs vectors (Stratagene, La Jolla, USA) containing a 400 bp fragment of R. culicivorax mex-3 via run off in vitro transcription with T7 or T3 RNA-polymerase according to the manufacturer’s protocol (Roche, Mannheim, Germany). The concentration of the labelled probes was about 300 n g×m l-1.
We are indebted to E. Platzer, Riverside, for the continuous supply with R. culcivorax nematodes. We thank J. Schulze, Cologne, for advice on nematode cultivation and C. Becker and K. Konrad for expert technical assistance in the genome sequencing experiments. We are also grateful to H. Oezden, Cologne for assistance with In-situ hybridisations. We thank J. Parkinson, Toronto, for providing a conserved NEMBASE4 protein set, Elizabeth Martínez Salazar, Zacatecas, Mexico, for Feulgen C-value data and A. H. Jay Burr, Vancouver, Canada for sharing preliminary results on phototaxis in R. culicivorax.
Assemblies and other computations were conducted on the HPC cluster “CHEOPS” at the University of Cologne (http://rrzk.uni-koeln.de/cheops.html).
This work was partly funded through the SFB 680: “Molecular Basis of Evolutionary Innovations”. Philipp H. Schiffer is funded by the VolkswagenStiftung in the “Förderinitiative Evolutionsbiologie”. Gerogios D. Koutsovoulos is funded by a UK BBSRC Research Studentship and an Overseas Reasearch Studentship from the University of Edinburgh. Additional Funding came through the BMBF-Projekt “NGSgoesHPC”.
Raw genome and transcriptome sequence data reported in this manuscript have been deposited in the ERA under accession ERP002111 (http://www.ebi.ac.uk/ena/data/view/ERP002111), and assembled genomic contigs deposited in the ENA INSDC database under accession numbers CAQS01000001-CAQS01062537. Annotation information and additional data are available through (http://romanomermis.nematod.es).
- Maduro MF: Cell fate specification in the C. elegans embryo. Dev Dyn. 2010, 239: 1315-1329.PubMedGoogle Scholar
- Sulston JE, Schierenberg E, White JG, Thomson JN: The embryonic cell lineage of the nematode Caenorhabditis elegans. Dev Biol. 1983, 100: 64-119. 10.1016/0012-1606(83)90201-4.PubMedGoogle Scholar
- Blaxter M, de Ley P, Garey J, Liu L, Scheldeman P, Vierstraete A, Vanfleteren J, Mackey L, Dorris M, Frisse L: A molecular evolutionary framework for the phylum Nematoda. Nature. 1998, 392 (6671): 71-75. 10.1038/32160.PubMedGoogle Scholar
- Meldal B, Debenham N, de Ley P, de Ley I, Vanfleteren J, Vierstraete A, Bert W, Borgonie G, Moens T, Tyler P, Austen M, Blaxter M, Rodgers A, Lambshead P: An improved molecular phylogeny of the Nematoda with special emphasis on marine taxa. Mol Phylogenet Evol. 2007, 42 (3): 622-636. 10.1016/j.ympev.2006.08.025.PubMedGoogle Scholar
- Boveri T: Die Entwicklung von Ascaris megalocephala mit besonderer Ruecksicht auf die Kernverhaeltnisse. 1899, Fischer: Festschrift fuer Carl von Kupffer; JenaGoogle Scholar
- Müller H: Beitrag zur Embryonalentwicklung von Ascaris megalocephala. Zoologica. 1903, 41: 60-Google Scholar
- Vangestel S, Houthoofd W, Bert W, Borgonie G: The early embryonic development of the satellite organism Pristionchus pacificus: differences and similarities with Caenorhabditis elegans. Nematology. 2008, 10: 301-312. 10.1163/156854108783900267.Google Scholar
- Skiba F, Schierenberg E: Cell lineages, developmental timing, and spatial pattern formation in embryos of free-living soil nematodes. Dev Biol. 1992, 151 (2): 597-610. 10.1016/0012-1606(92)90197-O.PubMedGoogle Scholar
- Lahl V, Schulze J, Schierenberg E: Differences in embryonic pattern formation between Caenorhabditis elegans and its close parthenogenetic relative Diploscapter coronatus. Int J Dev Biol. 2009, 53 (4): 507-515. 10.1387/ijdb.082718vl.PubMedGoogle Scholar
- Brauchle M, Kiontke K, Macmenamin P, Fitch DHA, Piano F: Evolution of early embryogenesis in rhabditid nematodes. Dev Biol. 2009, 335: 253-262. 10.1016/j.ydbio.2009.07.033.PubMed CentralPubMedGoogle Scholar
- Wiegner O, Schierenberg E: Specification of gut cell fate differs significantly between the NematodesAcrobeloides nanusandCaenorhabditis elegans. Dev Biol. 1998, 204: 3-14. 10.1006/dbio.1998.9054.PubMedGoogle Scholar
- Wiegner O, Schierenberg E: Regulative development in a nematode embryo: a hierarchy of cell fate transformations. Dev Biol. 1999, 215: 1-12. 10.1006/dbio.1999.9423.PubMedGoogle Scholar
- Voronov DA, Panchin YV: Cell lineage in marine nematode Enoplus brevis. Dev. 1998, 125: 143-150.Google Scholar
- Schulze J, Schierenberg E: Evolution of embryonic development in nematodes. EvoDevo. 2011, 2: 18-10.1186/2041-9139-2-18.PubMed CentralPubMedGoogle Scholar
- Schulze J, Houthoofd W, Uenk J, Vangestel S, Schierenberg E: Plectus - a stepping stone in embryonic cell lineage evolution of nematodes. EvoDevo. 2012, 3: 13-10.1186/2041-9139-3-13.PubMed CentralPubMedGoogle Scholar
- True JR, Haag ES: Developmental system drift and flexibility in evolutionary trajectories. Evol Dev. 2001, 3 (2): 109-119. 10.1046/j.1525-142x.2001.003002109.x.PubMedGoogle Scholar
- C elegans Sequencing Consortium: Genome sequence of the nematode C. elegans: a platform for investigating biology. Science. 1998, 282 (5396): 2012-2018.Google Scholar
- Stein LD, Bao Z, Blasiar D, Blumenthal T, Brent MR, Chen N, Chinwalla A, Clarke L, Clee C, Coghlan A, Coulson A, D’eustachio P, Fitch DHA, Fulton LA, Fulton RE, Griffiths-Jones S, Harris TW, Hillier LW, Kamath R, Kuwabara PE, Mardis ER, Marra MA, Miner TL, Minx P, Mullikin JC, Plumb RW, Rogers J, Schein JE, Sohrmann M, Spieth J, et al: The genome sequence of, Caenorhabditis briggsae: a platform for comparative genomics. PLoS Biol. 2003, 1 (2): E45-PubMed CentralPubMedGoogle Scholar
- Mortazavi A, Schwarz EM, Williams B, Schaeffer L, Antoshechkin I, Wold BJ, Sternberg PW: Scaffolding a Caenorhabditis nematode genome with RNA-seq. Genome Res. 2010, 20 (12): 1740-1747. 10.1101/gr.111021.110.PubMed CentralPubMedGoogle Scholar
- Dieterich C, Clifton SW, Schuster LN, Chinwalla A, Delehaunty K, Dinkelacker I, Fulton L, Fulton R, Godfrey J, Minx P, Mitreva M, Roeseler W, Tian H, Witte H, Yang SP, Wilson RK, Sommer RJ: The Pristionchus pacificus genome provides a unique perspective on nematode lifestyle and parasitism. Nature Genet. 2008, 40 (10): 1193-1198. 10.1038/ng.227.PubMedGoogle Scholar
- Jex AR, Liu S, Li B, Young ND, Hall RS, Li Y, Yang L, Zeng N, Xu X, Xiong Z, Chen F, Wu X, Zhang G, Fang X, Kang Y, Anderson GA, Harris TW, Campbell BE, Vlaminck J, Wang T, Cantacessi C, Schwarz EM, Ranganathan S, Geldhof P, Nejsum P, Sternberg PW, Yang H, Wang J, Wang J, Gasser RB: Ascaris suum draft genome. Nature. 2011, 479 (7374): 529-533. 10.1038/nature10553.PubMedGoogle Scholar
- Ghedin E, Wang S, Spiro D, Caler E, Zhao Q, Crabtree J, Allen JE, Delcher AL, Guiliano DB, Miranda-Saavedra D, Angiuoli SV, Creasy T, Amedeo P, Haas B, El-Sayed NM, Wortman JR, Feldblyum T, Tallon L, Schatz M, Shumway M, Koo H, Salzberg SL, Schobel S, Pertea M, Pop M, White O, Barton GJ, Carlow CKS, Crawford MJ, Daub J, et al: Draft genome of the filarial nematode parasite Brugia malayi. Science. 2007, 317 (5845): 1756-1760. 10.1126/science.1145406.PubMed CentralPubMedGoogle Scholar
- Godel C, Kumar S, Koutsovoulos G, Ludin P, Nilsson D, Comandatore F, Wrobel N, Thompson M, Schmid CD, Goto S, Bringaud F, Wolstenholme A, Bandi C, Epe C, Kaminsky R, Blaxter M, Mäser P: The genome of the heartworm, Dirofilaria immitis, reveals drug and vaccine targets. FASEB J. 2012, 26 (11): 4650-4661. 10.1096/fj.12-205096.PubMed CentralPubMedGoogle Scholar
- Abad P, Gouzy J, Aury JM, Castagnone-Sereno P, Danchin EGJ, Deleury E, Perfus-Barbeoch L, Anthouard V, Artiguenave F, Blok VC, Caillaud MC, Coutinho PM, Dasilva C, De Luca F, Deau F, Esquibet M, Flutre T, Goldstone JV, Hamamouch N, Hewezi T, Jaillon O, Jubin C, Leonetti P, Magliano M, Maier TR, Markov GV, Mcveigh P, Pesole G, Poulain J, Robinson-Rechavi M, et al: Genome sequence of the metazoan plant-parasitic nematode, Meloidogyne incognita. Nature Biotechnol. 2008, 26 (8): 909-915. 10.1038/nbt.1482.Google Scholar
- Kikuchi T, Cotton JA, Dalzell JJ, Hasegawa K, Kanzaki N, Mcveigh P, Takanashi T, Tsai IJ, Assefa SA, Cock PJA, Otto TD, Hunt M, Reid AJ, Sanchez-Flores A, Tsuchihara K, Yokoi T, Larsson MC, Miwa J, Maule AG, Sahashi N, Jones JT, Berriman M: Genomic insights into the origin of parasitism in the emerging plant pathogenBursaphelenchus xylophilus. PLoS Pathogens. 2011, 7 (9): e1002219-10.1371/journal.ppat.1002219.PubMed CentralPubMedGoogle Scholar
- Srinivasan J, Dillman AR, Macchietto MG, Heikkinen L, Lakso M, Fracchia KM, Antoshechkin I, Mortazavi A, Wong G, Sternberg PW: The draft genome and transcriptome of Panagrellus redivivus are shaped by the harsh demands of a free-living lifestyle. Genetics. 2013, 193 (4): 1279-1295. 10.1534/genetics.112.148809.PubMed CentralPubMedGoogle Scholar
- Mitreva M, Jasmer DP, Zarlenga DS, Wang Z, Abubucker S, Martin J, Taylor CM, Yin Y, Fulton LA, Minx P, Yang SP, Warren WC, Fulton RS, Bhonagiri V, Zhang X, Hallsworth-Pepin K, Clifton SW, Mccarter JP, Appleton J, Mardis ER, Wilson RK: The draft genome of the parasitic nematode Trichinella spiralis. Nature Genet. 2011, 43 (3): 228-235. 10.1038/ng.769.PubMed CentralPubMedGoogle Scholar
- Hope IA: Embryology, Developmental Biology and the Genome. The Biology of Nematodes. Edited by: Lee DL. 2002, New York: Tayler & Francis, 121-145.Google Scholar
- Kumar S, Schiffer PH, Blaxter M: 959 Nematode Genomes: a semantic wiki for coordinating sequencing projects. Nucl Acids Res. 2012, 40 (D1): D1295-D1300. 10.1093/nar/gkr826.PubMed CentralPubMedGoogle Scholar
- Kumar S, Koutsovoulos G, Kaur G, Blaxter M: Toward 959 nematode genomes. Worm. 2012, 1: 0-8.Google Scholar
- Petersen JJ: Nematodes as biological control agents: Part I. Mermithidae. Advances in Parasitology. Edited by: Baker JR, Muller R. 1985, London: Academic Press, 307-346.Google Scholar
- Petersen JJ, Chapman HC, Willis OR, Fukuda T: Release of Romanomermis culicivorax for the control of Anopheles albimanus in El Salvador II. Application of the nematode. Ame J Trop Med Hyg. 1978, 27 (6): 1268-1273.Google Scholar
- Schulze J, Schierenberg E: Cellular pattern formation, establishment of polarity and segregation of colored cytoplasm in embryos of the nematode Romanomermis culicivorax. Dev Biol. 2008, 315 (2): 426-436. 10.1016/j.ydbio.2007.12.043.PubMedGoogle Scholar
- Schulze J, Schierenberg E: Embryogenesis of Romanomermis culicivorax: an alternative way to construct a nematode. Dev Biol. 2009, 334: 10-21. 10.1016/j.ydbio.2009.06.009.PubMedGoogle Scholar
- Kent W: BLAT—the BLAST-like alignment tool. Genome Res. 2002, 12 (4): 656-PubMed CentralPubMedGoogle Scholar
- Parra G, Bradnam K, Korf I: CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics. 2007, 23 (9): 1061-1067. 10.1093/bioinformatics/btm071.PubMedGoogle Scholar
- Cantarel BL, Korf I, Robb SMC, Parra G, Ross E, Moore B, Holt C, Sánchez Alvarado A, Yandell M: MAKER: an easy-to-use annotation pipeline designed for emerging model organism genomes. Genome Res. 2008, 18: 188-196.PubMed CentralPubMedGoogle Scholar
- Li W, Godzik A: Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006, 22 (13): 1658-1659. 10.1093/bioinformatics/btl158.PubMedGoogle Scholar
- Stanke M, Waack S: Gene prediction with a hidden Markov model and a new intron submodel. Bioinformatics. 2003, 19 (2): ii215-ii225.PubMedGoogle Scholar
- Wang J, Mitreva M, Berriman M, Thorne A, Magrini V, Koutsovoulos G, Kumar S, Blaxter ML, Davis RE: Silencing of germline-expressed genes by DNA elimination in somatic cells. Dev Cell. 2012, 23 (5): 1072-1080. 10.1016/j.devcel.2012.09.020.PubMed CentralPubMedGoogle Scholar
- Nawrocki EP, Kolbe DL, Eddy SR: Infernal 1.0: inference of RNA alignments. Bioinformatics. 2009, 25 (10): 1335-1337. 10.1093/bioinformatics/btp157.PubMed CentralPubMedGoogle Scholar
- Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucl Acids Res. 1997, 25 (5): 0955-0964.Google Scholar
- Kumar S: Data for PhD thesis on next generation nematode genomes. 2012, http://dx.doi.org/10.6084/m9.figshare.96089,Google Scholar
- Richards S, The T castaneum Genome Consortium: The genome of the model beetle and pest Tribolium castaneum. Nature. 2008, 452 (7190): 949-955. 10.1038/nature06784.PubMedGoogle Scholar
- Schröder R, Beermann A, Wittkopp N, Lutz R: From development to biodiversity—Tribolium castaneum, an insect model organism for short germband development. Dev Genes Evol. 2008, 218 (3-4): 119-126. 10.1007/s00427-008-0214-3.PubMedGoogle Scholar
- Chen F, Mackey AJ, Vermunt JK, Roos DS: Assessing performance of orthology detection strategies applied to eukaryotic genomes. PLoS ONE. 2007, 2 (4): e383-10.1371/journal.pone.0000383.PubMed CentralPubMedGoogle Scholar
- Altschul S, Madden T, Schäffer A: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucl Acids …. 1997, [http://nar.oxfordjournals.org/content/25/17/3389.short],Google Scholar
- Jensen RA: Orthologs and paralogs - we need to get it right. Genome Biol. 2001, 2 (8): interactions1002.1-1002.3.Google Scholar
- Koonin E: Orthologs, paralogs, and evolutionary genomics. Ann Rev Genet. 2005, 39: 309-338. 10.1146/annurev.genet.39.073003.114725.PubMedGoogle Scholar
- Moreno-Hagelsieb G, Latimer K: Choosing BLAST options for better detection of orthologs as reciprocal best hits. Bioinformatics. 2008, 24 (3): 319-324. 10.1093/bioinformatics/btm585.PubMedGoogle Scholar
- Shaye DS, Greenwald I: OrthoList: A compendium of C. elegans Genes with human Orthologs. PLoS ONE. 2011, 6 (5): e20085-10.1371/journal.pone.0020085.PubMed CentralPubMedGoogle Scholar
- Tautz D, Domazet-Los̆o T: The evolutionary origin of orphan genes. Nat Rev Genet. 2011, 12 (10): 692-702. 10.1038/nrg3053.PubMedGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.PubMedGoogle Scholar
- Keeling PJ, Corradi N, Morrison HG, Haag KL, Ebert D, Weiss LM, Akiyoshi DE, Tzipori S: The reduced genome of the parasitic microsporidian enterocytozoon bieneusi lacks genes for core carbon metabolism. Genome Biol Evol. 2010, 2 (0): 304-309. 10.1093/gbe/evq022.PubMed CentralPubMedGoogle Scholar
- Bird A: DNA methylation patterns and epigenetic memory. Genes & Dev. 2002, 16: 6-21. 10.1101/gad.947102.Google Scholar
- Gao F, Liu X, Wu XP, Wang XL, Gong D, Lu H, Xia Y, Song Y, Wang J, Du J, Liu S, Han X, Tang Y, Yang H, Jin Q, Zhang X, Liu M: Differential DNA methylation in discrete developmental stages of the parasitic nematode Trichinella spiralis. Genome Biol. 2012, 13 (10): R100-10.1186/gb-2012-13-10-r100.PubMed CentralPubMedGoogle Scholar
- Tran RK, Zilberman D, de Bustos C, Ditt RF, Henikoff JG, Lindroth AM, Delrow J, Boyle T, Kwong S, Bryson TD, Jacobsen SE, Henikoff S: Chromatin and siRNA pathways cooperate to maintain DNA methylation of small transposable elements in Arabidopsis. Genome Biol. 2005, 6 (11): R90-10.1186/gb-2005-6-11-r90.PubMed CentralPubMedGoogle Scholar
- Martienssen RA, Colot V: DNA methylation and epigenetic inheritance in plants and filamentous fungi. Science. 2001, 293 (5532): 1070-1074. 10.1126/science.293.5532.1070.PubMedGoogle Scholar
- Eisenmann DM: Wnt signaling. WormBook. 2005, [http://www.wormbook.org/chapters/www_wntsignaling/wntsignaling.html],Google Scholar
- Graham LD, Kotze AC, Fernley RT, Hill RJ: Molecular & biochemical parasitology. Mol Biochem Parasitol. 2010, 171 (2): 104-107. 10.1016/j.molbiopara.2010.03.003. [http://dx.doi.org/10.1016/j.molbiopara.2010.03.003],PubMedGoogle Scholar
- Levin M, Hashimshony T, Wagner F, Yanai I: Developmental milestones punctuate gene expression in the caenorhabditis embryo. Dev Cell. 2012, 22 (5): 1101-1108. 10.1016/j.devcel.2012.04.004.PubMedGoogle Scholar
- Phelan P: Innexins: members of an evolutionarily conserved family of gap-junction proteins. Biochimica et Biophysica Acta (BBA) - Biomembranes. 2005, 1711 (2): 225-245. 10.1016/j.bbamem.2004.10.004.Google Scholar
- Jones AK, Sattelle DB: Functional genomics of the nicotinic acetylcholine receptor gene family of the nematode, Caenorhabditis elegans. BioEssays. 2003, 26: 39-49.Google Scholar
- Antebi A: Nuclear hormone receptors in C. elegans. Wormbook. Edited by: The C. elegans Research Community. 2006, WormBook, http://www.wormbook.org,Google Scholar
- Altun ZF, Chen B, Wang ZW, Hall DH: High resolution map of Caenorhabditis elegansgap junction proteins. Dev Dyn. 2009, 238 (8): 1936-1950. 10.1002/dvdy.22025.PubMed CentralPubMedGoogle Scholar
- Haag E: The evolution of nematode sex determination: C. elegans as a reference point for comparative biology. WormBook. Edited by: The C. elegans Research Community. 2005, WormBook, http://www.wormbook.org,Google Scholar
- Tingley GA, Anderson RM: Environmental sex determination and density-dependent population regulation in the entomogenous nematode Romanomermis culicivorax. Parasitol. 1986, 92: 431-449. 10.1017/S0031182000064192.Google Scholar
- Powell JR, Jow MM, Meyer BJ: The T-box transcription factor SEA-1 is an autosomal element of the X:A signal that determines C. elegans sex. Dev Cell. 2005, 9 (3): 339-349. 10.1016/j.devcel.2005.06.009.PubMed CentralPubMedGoogle Scholar
- Chu DS, Dawes HE, Lieb JD, Chan RC, Kuo AF, Meyer BJ: A molecular link between gene-specific and chromosome-wide transcriptional repression. Genes & Dev. 2002, 16 (7): 796-805. 10.1101/gad.972702.Google Scholar
- Meyer B: X-Chromosome dosage compensation. WormBook. Edited by: The C. elegans Research Community. 2005, WormBook, http://www.wormbook.org,Google Scholar
- Zarkower D: Somatic sex determination. WormBook. Edited by: The C. elegans Research Community. 2006, WormBook, http://www.wormbook.org,Google Scholar
- Kuwabara PE, Okkema PG, Kimble J: tra-2 encodes a membrane protein and may mediate cell communication in the Caenorhabditis elegans sex determination pathway. Mol Biol Cell. 1992, 3 (4): 461-473. 10.1091/mbc.3.4.461.PubMed CentralPubMedGoogle Scholar
- Goodwin EB, Ellis RE: Turning clustering loops: sex determination in Caenorhabditis elegans. Curr Biol. 2002, 12 (3): R111-R120. 10.1016/S0960-9822(02)00675-9.PubMedGoogle Scholar
- Baldi C, Cho S, Ellis RE: Mutations in two independent pathways are sufficient to create hermaphroditic nematodes. Science. 2009, 326 (5955): 1002-1005. 10.1126/science.1176013.PubMedGoogle Scholar
- Goldstein B, Macara IG: The PAR proteins: fundamental players in animal cell polarization. Dev Cell. 2007, 13 (5): 609-622. 10.1016/j.devcel.2007.10.007.PubMed CentralPubMedGoogle Scholar
- Bowerman B: Embryonic polarity: Protein stability in asymmetric cell division. Curr Biol. 2000, 10 (17): R637-R641. 10.1016/S0960-9822(00)00660-6.PubMedGoogle Scholar
- Severson AF, Bowerman B: J Cell Biol. 2003, 161: 21-26. 10.1083/jcb.200210171.PubMed CentralPubMedGoogle Scholar
- Gönczy P, Rose LS: Asymmetric cell division and axis formation in the embryo. WormBook. Edited by: The C. elegans Research Community. WormBook, 2005-2005. http://www.wormbook.org,
- Cheng NN, Kirby CM, Kemphues KJ: Control of cleavage spindle orientation in Caenorhabditis elegans: the role of the genes par-2 and par-3. Genetics. 1995, 139 (2): 549-559.PubMed CentralPubMedGoogle Scholar
- Wu JC, Rose LS: PAR-3 and PAR-1 inhibit LET-99 localization to generate a cortical band important for spindle positioning in Caenorhabditis elegans embryos. Mol Biol Cell. 2007, 18 (11): 4470-4482. 10.1091/mbc.E07-02-0105.PubMed CentralPubMedGoogle Scholar
- Doerflinger H, Vogt N, Torres IL, Mirouse V, Koch I, Nusslein-Volhard C, St Johnston D: Bazooka is required for polarisation of the Drosophila anterior-posterior axis. Development. 2010, 137 (10): 1765-1773. 10.1242/dev.045807.PubMed CentralPubMedGoogle Scholar
- Goldstein B, Frisse L, Thomas W: Embryonic axis specification in nematodes: evolution of the first step in development. Curr Biol. 1998, 8 (3): 157-160. 10.1016/S0960-9822(98)70062-4.PubMedGoogle Scholar
- Draper BW, Mello CC, Bowerman B, Hardin J, Priess JR: MEX-3 is a KH domain protein that regulates blastomere identity in early C. elegans embryos. Cell. 1996, 87 (2): 205-216. 10.1016/S0092-8674(00)81339-2.PubMedGoogle Scholar
- Huang N, Mootz D, Walhout A, Vidal M, Hunter CP: MEX-3 interacting proteins link cell polarity to asymmetric gene expression in Caenorhabditis elegans. Development. 2002, 129 (3): 747-759.PubMedGoogle Scholar
- Evans TC, Hunter CP: Translational control of maternal RNAs. WormBook. Edited by: The C. elegans Research Community. WormBook, 2005-2005. http://www.wormbook.org,
- Gomes JE, Encalada SE, Swan KA, Shelton CA, Carter JC, Bowerman B: The maternal gene spn-4 encodes a predicted RRM protein required for mitotic spindle orientation and cell fate patterning in early C. elegans embryos. Dev. 2001, 128 (21): 4301-4314.Google Scholar
- Labbé JC, Goldstein B: Embryonic development: A New SPN on cell fate specification. Curr Biol. 2002, 12 (11): R396-R398. 10.1016/S0960-9822(02)00884-9.PubMedGoogle Scholar
- Simske JS, Hardin J: Getting into shape: epidermal morphogenesis in Caenorhabditis elegans embryos. BioEssays. 2001, 23: 12-23. 10.1002/1521-1878(200101)23:1<12::AID-BIES1003>3.3.CO;2-I.PubMedGoogle Scholar
- Gilleard JS, McGhee JD: Activation of hypodermal differentiation in the Caenorhabditis elegans embryo by GATA transcription factors ELT-1 and ELT-3. Mol Cell Biol. 2001, 21 (7): 2533-2544. 10.1128/MCB.21.7.2533-2544.2001.PubMed CentralPubMedGoogle Scholar
- Sommer R: As good as they get: cells in nematode vulva development and evolution. Curr Opin Cell Biol. 2001, 13 (6): 715-720. 10.1016/S0955-0674(00)00275-1.PubMedGoogle Scholar
- Kiontke K, Barriere A, Kolotuev I, Podbilewicz B, Sommer R, Fitch DHA, Félix MA: Trends, stasis, and drift in the evolution of nematode vulva development. Curr Biol : CB. 2007, 17 (22): 1925-1937. 10.1016/j.cub.2007.10.061.PubMedGoogle Scholar
- Sternberg PW: Vulval development. WormBook. Edited by: The C. elegans Research Community. 2005, WormBook, http://www.wormbook.org,Google Scholar
- Salser SJ, Loer CM, Kenyon C: Multiple HOM-C gene interactions specify cell fates in the nematode central nervous system. Genes & Dev. 1993, 7 (9): 1714-1724. 10.1101/gad.7.9.1714.Google Scholar
- Eisenmann DM, Kim SK: Protruding vulva mutants identify novel loci and Wnt signaling factors that function during Caenorhabditis elegans vulva development. Genetics. 2000, 156 (3): 1097-1116.PubMed CentralPubMedGoogle Scholar
- Shemer G, Podbilewicz B: LIN-39/Hox triggers cell division and represses EFF-1/fusogen-dependent vulval cell fusion. Genes & Dev. 2002, 16 (24): 3136-3141. 10.1101/gad.251202.Google Scholar
- Tian H, Schlager B, Xiao H, Sommer RJ: Wnt signaling induces vulva development in the nematode Pristionchus pacificus. Curr Biol. 2008, 18 (2): 142-146. 10.1016/j.cub.2007.12.048.PubMedGoogle Scholar
- Streit A, Kohler R, Marty T, Belfiore M, Takacs-Vellai K, Vigano MA, Schnabel R, Affolter M, Müller F: Conserved Regulation of the Caenorhabditis elegans labial/Hox1 Gene ceh-13. Dev Biol. 2002, 242 (2): 96-108. 10.1006/dbio.2001.0544.PubMedGoogle Scholar
- Aboobaker AA, Blaxter ML: Hox Gene Loss during dynamic evolution of the nematode cluster. Curr Biol. 2003, 13: 37-40. 10.1016/S0960-9822(02)01399-4.PubMedGoogle Scholar
- Chisholm A: Control of cell fate in the tail region of C. elegans by the gene egl-5. Deve. 1991, 111 (4): 921-932.Google Scholar
- Aboobaker A, Blaxter M: Hox gene evolution in nematodes: novelty conserved. Current Opinion in Genet & Dev. 2003, 13: 593-598. 10.1016/j.gde.2003.10.009.Google Scholar
- Lemons D, McGinnis W: Genomic evolution of Hox gene clusters. Science. 2006, 313 (5795): 1918-1922. 10.1126/science.1132040.PubMedGoogle Scholar
- Chamberlin HM, Thomas JH: The bromodomain protein LIN-49 and trithorax-related protein LIN-59 affect development and gene expression in Caenorhabditis elegans. Development. 2000, 127 (4): 713-723.PubMedGoogle Scholar
- Aspöck G, Kagoshima H, Niklaus G, Bürglin TR: Caenorhabditis elegans Has scores of hedgehog related genes: sequence and expression analysis. Genome Res. 1999, 9: 909-923. 10.1101/gr.9.10.909. http://genome.cshlp.org,PubMedGoogle Scholar
- Domazet-Los̆o T, Tautz D: A phylogenetically based transcriptome age index mirrors ontogenetic divergence patterns. Nature. 2010, 468 (7325): 815-818. 10.1038/nature09632.PubMedGoogle Scholar
- Richardson MK: A phylotypic stage for all animals?. Dev Cell. 2012, 22 (5): 903-904. 10.1016/j.devcel.2012.05.001.PubMedGoogle Scholar
- Kalinka AT, Tomancak P: The evolution of early animal embryos: conservation or divergence?. Trends In Ecology & Evolution. 2012, 27 (7): 385-393. 10.1016/j.tree.2012.03.007.Google Scholar
- De Ley P: A quick tour of nematode diversity and the backbone of nematode phylogeny. WormBook. Edited by: The C. elegans Research Community. 2006, WormBook, http://www.wormbook.org,Google Scholar
- Martin M: Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.Journal. 2011, 17 (1):Google Scholar
- CLCbio: White Paper on de Novo Assembly in CLC Assembly Cell. 2010, CLC Bio: WhitepaperGoogle Scholar
- Kumar S, Blaxter ML: Simultaneous genome sequencing of symbionts and their hosts. Symbiosis. 2012, 55 (3): 119-126.PubMed CentralGoogle Scholar
- Boetzer M, Henkel CV, Jansen HJ, Butler D, Pirovano W: Scaffolding pre-assembled contigs using SSPACE. Bioinformatics. 2011, 27 (4): 578-579. 10.1093/bioinformatics/btq683.PubMedGoogle Scholar
- Smit AFA, Hubley R, Green P: RepeatMasker Open-3.0. 1996-2010. http://www.repeatmasker.org,
- Jurka J, Kapitonov VV, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J: Repbase Update, a database of eukaryotic repetitive elements. Cytogenetic and Genome Res. 2005, 110 (1-4): 462-467. 10.1159/000084979.Google Scholar
- Volfovsky N, Haas BJ, Salzberg SL: A clustering method for repeat analysis in DNA sequences. Genome Biol. 2001, 2 (8): research0027.1-research0027.11.Google Scholar
- Bao Z, Eddy SR: Automated de novo identification of repeat sequence families in sequenced genomes. Genome Res. 2002, 12 (8): 1269-1276. 10.1101/gr.88502.PubMed CentralPubMedGoogle Scholar
- Price AL, Jones NC, Pevzner PA: De novo identification of repeat families in large genomes. Bioinformatics. 2005, 21 (Suppl 1): i351-i358. 10.1093/bioinformatics/bti1018.PubMedGoogle Scholar
- Parkinson J, Mitreva M, Whitton C, Thomson M: A transcriptomic analysis of the phylum Nematoda. Nature Genet. 2004, 36: 1259-1267. 10.1038/ng1472.PubMedGoogle Scholar
- Gremme G, Brendel V, Sparks ME, Kurtz S: Engineering a software tool for gene structure prediction in higher organisms. Information and Software Technol. 2005, 47 (15): 965-978. 10.1016/j.infsof.2005.09.005.Google Scholar
- Supek F, Vlahovicek K: INCA: synonymous codon usage analysis and clustering by means of self-organizing map. Bioinformatics. 2004, 20 (14): 2329-2330. 10.1093/bioinformatics/bth238.PubMedGoogle Scholar
- Cutter AD, Wasmuth JD, Blaxter ML: The evolution of biased codon and amino acid usage in nematode genomes. Mol Biol Evol. 2006, 23 (12): 2303-2315. 10.1093/molbev/msl097.PubMedGoogle Scholar
- 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 (18): 3674-3676. 10.1093/bioinformatics/bti610.PubMedGoogle Scholar
- The Gene Ontology Consortium: Gene Ontology: tool for the unification of biology. Nature Genetics. 2000, 25: 25-29. 10.1038/75556.PubMed CentralGoogle Scholar
- Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler R, Lopez R: InterProScan: protein domains identifier. Nucl Acids Res. 2005, 33 (Web Server): W116-W120. 10.1093/nar/gki442.PubMed CentralPubMedGoogle Scholar
- Griffiths-Jones S, Moxon S, Marshall M, Khanna A, Eddy S, Bateman A: Rfam: annotating non-coding RNAs in complete genomes. Nucl Acids Res. 2004, 33 (Database issue): D121-D124. 10.1093/nar/gki081.PubMed CentralGoogle Scholar
- Li L, Stoeckert CJ, Roos DS: OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Research. 2003, 13 (9): 2178-2189. 10.1101/gr.1224503.PubMed CentralPubMedGoogle Scholar
- van Dongen S: A Cluster algorithm for graphs. Report - Information Syst. 2000, 10: 1-40.Google Scholar
- van Dongen S: Graph Clustering by Flow Simulation. PhD thesis,. University of Utrecht 2000. http://www.wormbook.org,
- Seydoux G, Fire A: Whole-mount in situ hybridization for the detection of RNA in Caenorhabditis elegans embryos. Methods in cell biology. 1995, 48: 323-337.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.