On the Extent and Origins of Genic Novelty in the Phylum Nematoda
[['The vast majority of species are unlikely to be selected for whole genome sequencing , whatever their importance in terms of evolution , health and ecology .', 'The few eukaryote species selected for such projects , despite their utility in laboratory investigation , are unlikely to be representative of the genomic diversity of speciose phyla .', 'For example , Arthropoda and Nematoda have over one million species each [1] , [2] and the ∼20 genomes completed [3]–[7] or in sequencing will illuminate only small parts of their diversity .', 'Expressed sequence tags ( ESTs ) have proved to be a cost-effective and rapid method for identification of the genes from a target species [8] .', "Although the largest EST collections have been generated primarily for the annotation of complete genome sequences ( e . g . human and mouse ) , more than half the sequences in GenBank's EST depository ( dbEST ) [9] are from otherwise neglected genomes .", 'One phylum that has benefited from an EST sequencing approach is the Nematoda [10]–[13] .', 'Nematodes ( or round worms ) are abundant and diverse in terms of biology and ecology [14] .', 'They are ubiquitous members of the meiofauna and play a core role in nutrient recycling .', 'Parasitic species of this phylum are the causative agents of six of the thirteen neglected tropical diseases which afflict around 2 . 7 billion people [15]–[19] .', 'The diseases caused by nematodes are extremely varied , and include anaemia and malnutrition ( caused by hookworms such Ancylostoma ceylanicum ) , African river blindness ( caused by the filarial nematode Onchocerca volvulus ) and elephantiasis ( caused by the filarial nematode Brugia malayi ) .', 'In terms of disability adjusted life years ( DALYS ) , the burden of lymphatic filariasis ( 5 . 8 million DALYs ) , onchocerciasis ( 0 . 5 million DALYs ) and intestinal nematode infections ( 3 million DALYs ) is significant .', 'Among school aged children ( 5–14 ) , the impact of intestinal nematodes is even greater than malaria [20] .', 'Parasites are also responsible for substantial losses in agriculture .', 'Plant-parasitic nematodes , such as the root-knot nematodes ( Meloidogyne spp . ) , are major crop pathogens throughout the world , impacting both the quantity and quality of marketable yields , causing an estimated US$80bn in damage annually [21] , and parasites of livestock are the cause for severe economic losses .', 'The fully sequenced genomes of the free-living nematodes C . elegans and C . briggsae makes the analysis of EST datasets from parasitic nematode species particularly informative , in that both elements of core biology and particular adaptations specific to parasitism can be investigated .', 'Already more than a dozen species- or family-specific analyses of nematode EST datasets have been published , considering parasites of humans [22]–[24] , animals [25]–[27] and plants [28] , [29] .', 'The first whole-phylum meta-analysis was based on 265 , 000 sequences from 30 species , defining 93 , 645 putative genes [11] .', "Surprisingly , 30–70% of each species' dataset was found to have no significant similarity ( as defined by BLAST searches ) with any other sequence either within or outwith the sampled nematodes .", 'Do these sequences define new genes , with new functions in nematodes ?', 'Or are they transcriptional noise derived from non-coding sequence with no functional significance ?', "The majority of functional annotations have been assigned through sequence similarity to other proteins [30] , and thus a large number of nematode proteins lack clues as to their importance to the organism's survival .", 'In the absence of annotation , these data are limited in their practical use , for example , in identifying the lead novel targets for anthelminthic drugs .', "One indication of a gene's significance , in worm survival , is its presence in a number of nematode species .", 'Proteins with essential , conserved functions will tend to be conserved between species , and thus will be members of protein families .', 'Protein families restricted to the Nematoda , but found in a number of species , invite further study to reveal their function .', 'Proteins often share local regions or similarity despite being non-orthologous [31] , with the interplay between these domains underpinning their function .', 'There are a number of widely used protein domain databases [32]–[35] which provide domain models to search .', 'In addition , it is possible to identify new domains through similarity searches [36] , and nematode-restricted novel domains may yield novel insights into avenues for control of parasites .', 'EST datasets have been considered less than ideal for such analyses , due to the occurrence of frame-shifts , ambiguous base calls and untranslated regions [37]–[39] .', 'However , coding regions can be accurately predicted from EST cluster consensuses using a hierarchical approach such as that employed by prot4EST [39] .', 'A great deal of care must be taken when translating sequences that do not have sequence similarity to known proteins .', 'ESTScan , incorporated in the prot4EST pipeline , locates ( and corrects ) coding regions through the identification of frames that have oligonucleotide frequencies resembling those of the training dataset .', 'However , by definition few sequence data are available in the public repositories for neglected species such as parasitic nematodes .', 'Here we have inferred protein translations for over 120 , 000 putative genes from EST data from 37 species of nematodes using both high quality codon usage tables for each species [40] and synthetic training sets .', 'This protein dataset , NemPep3 , is employed here to investigate protein family ( NemFam3 ) and protein domain ( NemDom3 ) composition of nematodes , and presented in an online database NEMBASE3 .', 'Our key findings are:'], ['Sequence data were sourced from EMBL/GenBank/DDBJ and from WormBase ( http://www . wormbase . org ) as follows: We used TRIBE-MCL to generate protein families from NemPep3 [50] .', 'In TRIBE-MCL , the Inflation parameter defines the tightness of the clusters .', 'No single Inflation parameter value will correctly return all protein families , just as no single molecular clock exists to describe the evolution of all genes .', 'Therefore we repeated the clustering procedure over a range of values and recorded all the clusters , following a previous study of prokaryote proteins [51] .', 'The input to TRIBE-MCL was an all-against-all BLAST report .', 'The number of families generated varied from 42 , 865 to 71 , 867 .', 'All five sets of protein families are stored in NemBase3 .', "We used NemFam3 to investigate how sampling from additional species affected the discovery of protein families , generating a “collector's curve” of discovery of novelty .", 'First we took those families for which the only nematode species present was C . elegans .', 'We then added new families identified in each species in turn , adding them in the approximate order of their phylogenetic distance from C . elegans: Rhabditoidea ( CBG; see Figure 1 for three letter species codes ) ; Strongyloidea ( ACP , AYP , HCP , NAP , NBP , OOP , TDP ) ; Diplogasteromorpha ( PPP ) ; Panagrolaimomorpha ( PTP , SRP , SSP ) ; Tylenchomorpha ( GPP , GRP , HGP , HSP , MAP , MCP , MHP , MIP , MPP , PEP , PVP , RSP ) ; Cephalobomorpha ( ZPP ) ; Ascaridomorpha ( ALP , ASP , TCP ) ; Spiruromorpha ( BMP , DIP , LSP , OVP , WBP ) ; Trichinellida ( TMP , TVP , TSP ) ; Dorylaimida ( XIP ) .', 'All EST derived proteins were annotated with matches to the KEGG database [52] with a script developed in house which makes use of BLAST comparisons .', 'We wanted to identify metabolic processes absent in C . elegans but present in other nematodes .', 'To do this we compiled two separate lists of metabolites that are substrates of enzymes in C . elegans and in the other nematodes .', 'This step was important to reduce redundancy , as more than one enzyme ( EC number ) can be assigned to the same step of a pathway .', 'Next we compared the two lists and extracted those substrates missing from C . elegans , highlighting the enzymes that catalyse transformation of these molecules .', 'The Enzyme Commission ( EC ) identifiers of these proteins were obtained through the KEGG database .', 'Assignment of signal peptides was done using the SignalP3 . 0 web-interface [53] with the following parameters: organism group - eukaryotes; method - both neural networks and hidden Markov models; truncation - first 70 residues .', 'We used three Boolean tests provided by SignalP3 . 0 to determine if a signal peptide was present: first ‘D’ must be true; secondly , we considered ‘Cmax’ and ‘Ymax’ , if both were true then we deemed this strong evidence and weaker evidence if only one category was true .', 'Analyses of the secreted proteomes have been carried out previously for Nippostrongylus brasiliensis [54] and H . schachtii [55] .', 'Compared with these studies , and despite using more conservative parameters , we identified a larger number of signal peptide-containing proteins in N . brasiliensis ( 96 were identified , compared with 87 from Harcus et al . [54] ) and H . schachtii ( 105 identified compared with 65 from Vanholme et al . [55] ) .', 'This increase is likely to derive from more robust coding region predictions producing proteins that were more likely to contain the correct N-terminus .', 'NemPep3 proteins were annotated with protein domains using existing domain databases ( PfamA and ProDom ) and by de novo identification of domains in unannotated sequence .', 'NEMBASE3 is a relational database built using the PostgreSQL database manager ( http://postgresql . org ) .', 'It holds all the data types described above , including sequences , clustering information , consensuses derived from EST clusters , peptide predictions , protein families and protein domains .', 'All peptides have been annotated with extensive BLAST-based similarity data , as well as quality scored functional annotation ( GO , EC and KEGG identifiers ) derived from GOtcha [61] and annot8r [62] analyses .', 'The database is available through the www using custom php scripts from http://www . nematodes . org/ .'], ['Coding regions for EST cluster consensuses derived from NEMBASE [12] from 37 species from the phylum Nematoda were predicted using prot4EST , yielding a total of 121 , 694 polypeptide sequences ( Figure 1 ) .', 'For each species , specific codon usage tables [40] were used to reverse translate the C . elegans proteome , providing synthetic training-set transcriptomes ( see Methods ) .', 'To assess the accuracy of synthetic transcriptomes , partial datasets built for C . elegans [39] were translated in a similar fashion .', 'Comparison with a complete collection of coding sequences showed only a slight reduction in prediction using synthetic transcriptomes ( data not shown ) .', 'Importantly , for most species the simulated training sets were more accurate than simply using the complete C . elegans or C . briggsae transcriptomes .', 'The mean length of translation for the EST datasets ( excluding the caenorhabditids ) was 137 amino acids ( aa ) ( standard deviation 65 aa ) , and 84% of the bases in the EST cluster consensuses contributed to translations .', 'The regions not covered are likely to be predominantly untranslated regions , as well as regions of low-complexity sequence .', 'Previously , we have shown that the most accurate translations are obtained using similarity to a known protein or the prot4EST implementation of the ESTScan algorithm [38] , [39] .', 'For most nematode species , over 90% of EST cluster consensuses were translated using these two methods ( Figure 1 ) .', 'However , three Spiruromorph species had much lower rates of translation by these methods: Brugia malayi ( 71% translated using similarity or ESTScan methods ) , Onchocerca volvulus ( 78% ) and Wuchereria bancrofti ( 68% ) ( Figure 1; ‘percentage accepted’ ) .', 'These low rates appear to arise from two features of these data .', "Firstly , a relatively low proportion ( ∼40% ) of these species' EST cluster consensuses had significant similarity to protein sequences in UniRef100 [46] .", 'Secondly , only ∼54% of the novel sequences had compositions that matched models derived from known coding regions , simulated transcriptomes , or , in the case of B . malayi where a first pass annotation of the whole genome sequence is available [63] , an extensive transcriptome dataset .', 'Our inability to derive high quality translations for a significant number of clusters from these taxa could be due to a major biological difference and to the quality of the training set used or to the quality of the sequence data .', 'Other species that had similarly low proportions of sequence similarity matches , had higher rates of compositionally-identified coding regions ( e . g . Trichuris vulpis with 80% of the novel sequences translated by ESTScan and Meloidogyne javanica with 97% ) .', 'The addition of a 12 , 000-transcript , orthologous training set [63] did not improve the proportion of B . malayi cluster consensuses that yielded a translation .', 'For these three problem species , we noted that singleton cluster consensuses were much less likely to be robustly translated , but these species did not have an excess of singletons compared to the other nematodes .', 'The proportions of ESTs lacking detectable coding regions were compared between the source cDNA libraries .', 'Of 25 B . malayi libraries , five were significantly enriched for ESTs not translated ( G-statistic\u200a=\u200a682; p≪0 . 001 ) .', 'Two libraries from the eight available for O . volvulus and two for W . bancrofti were also shown to contain an excess of ESTs without a coding region .', 'Strikingly , 93% of the untranslatable sequences from B . malayi came from the highlighted five libraries , while the O . volvulus and W . bancrofti libraries accounted for around 30% of each species suspect contigs .', "We conclude that some of the unique features of the three species' data derive from the relative quality of some cDNA libraries sampled .", 'To ensure that subsequent analyses were performed on the most accurate collection of polypeptides , we excluded EST cluster consensuses that could not be translated with either the sequence similarity or ESTScan components of prot4EST .', 'Addition of the proteomes from the fully-sequenced C . elegans and C . briggsae yielded a high quality dataset ( NemPep3 ) .', 'The current release of NemPep , version 3 , includes 154 , 501 polypeptide sequences ( Figure 1 ) , with a mean length of 220 amino acids .', 'NemPep3 is available for download from NEMBASE3 ( http://www . nematodes . org/nembase3/ ) .', 'We used TRIBE-MCL [50] to derive putative protein families ( NemFam3 ) from NemPep3 .', 'These families were compared to proteins from the UniProt database [46] to identify overlap with previously defined protein families .', 'The results of the clustering algorithm , MCL , can be tuned with an Inflation parameter .', 'In the context of protein clusters , this value determines how tight , or strict , the clustering is ( see Methods ) .', 'No single parameter set for TRIBE-MCL can be used to accurately identify all ( or even most ) families and so we generated independent estimates at five different Inflation values .', 'To simplify analyses presented here , we have examined in detail the 65 , 179 protein families generated using an Inflation value of 3 . 0 , the default used for the TRIBE-MCL database [64] .', 'Despite having a large sample ( 37 species and over 150 , 000 individual sequences ) we found no evidence of having exhausted the diversity of nematode ‘protein space’ .', 'There was a near-linear increase in the number of protein families identified with addition of sequences and species ( Figure 2 ) .', 'This finding is congruent with that of Parkinson et al . ( 2004b ) but here we have used a rigorous protein family definition schema rather than simply BLAST matches .', 'Analyses of complete prokaryote proteomes also show an increase in the number of novel proteins as further species are sequenced [65] , although as a proportion of all prokaryote proteins the number of novel proteins is decreasing [66] .', 'This trend is not apparent in the nematode dataset ( Figure 2 ) .', 'The distribution of size of the NemFam3 protein families can be described by a power law , matching that of many protein family databases ( Figure 3a ) [67] .', 'We identified protein families that were restricted to all levels of nematode taxonomy , from species-specific to phylum-specific ( Figure 3b ) .', 'By comparing NemFam3 families to proteins from non-nematode species , we divided them into three classes: NemFam3 families that were unique to the Nematoda ( region A of Figure 2 ) ; NemFam3 families that were not found in C . elegans but did have homologues in other phyla ( region B ) ; and NemFam3 families that included C . elegans members and had homologues in other phyla ( region C ) .', 'Region C presumably encompasses proteins with core metabolic functions shared with other phyla .', 'Gene loss is a common feature of genome evolution [68]–[70] .', 'Gene gain by horizontal gene transfer is common in non-eukaryotes , but its role in eukaryotes , and particularly in metazoans , is still controversial [71]–[73] .', 'Gene loss in C . elegans has been reported previously [6] , [74]–[76] .', 'For example , orthologues of the Hox genes Antennapedia and Hox3 are absent in C . elegans but present in B . malayi and other invertebrates [76] .', 'Comparison of C . elegans and C . briggsae [6] identified a large number of proteins in each species that did not have an orthologue in the other .', 'Using NemPep3 and the UniProt database ( release 5 ) reduces the number of orphan proteins in C . elegans from 2 , 108 to 1 , 846 and from 2 , 141 to 1 , 961 in C . briggsae .', 'Comparison with proteomes from additional Caenorhabditis sp .', 'genomes currently being sequenced will clarify the patterns of gene gain and loss in this lineage .', 'We identified 4 , 864 protein families ( containing 6 , 903 proteins ) that had significant sequence similarity to proteins from outside the Nematoda but that contained no C . elegans representatives ( ‘loss/gain candidates’ ) .', 'To investigate the effect of using partial sequences , we compared loss/gain candidate EST cluster consensuses to the C . elegans genome .', 'Thirty-nine loss/gain candidate families ( 92 sequences ) could be aligned to the genome ( using the program BLAT [77] ) and overlapped an annotated coding sequence: the failure of TRIBE-MCL to group the C . elegans proteins with their loss/gain candidate matches was because their BLAST alignments had a low similarity score .', 'Three loss/gain candidate families ( eight proteins ) matched regions of the C . elegans genome that were not part of a coding region: these may correspond to valid but unannotated genes in C . elegans .', 'Thus the majority of loss/gain candidate families are absent from C . elegans .', 'Gene Ontology ( GO ) annotation of the loss/gain candidate families showed that a large number are involved in metabolism .', 'One hundred and fourteen individual Enzyme Commission ( EC ) classifications could be assigned to 240 families ( a full list of these annotations is available in Table S1 ) .', 'Some of these putative functions complemented gaps in the metabolic map of C . elegans .', 'For example , C . elegans lacks a canonical DNA methylation pathway enzyme , cytosine-5′-methyltransferase [78] .', 'Homologues of cytosine-5′-methyltransferase were identified in Ostertagia ostertagi , Teladorsagia .', 'circumcincta and Xiphinema index ( Figure 4 ) , and a homologue has also been identified in Pristionchus pacificus [79] .', 'It will be informative to examine additional nematode genomes for the features of DNA methylation and thus identify when , and perhaps why , this core regulatory mechanism was lost .', 'While the above examples reveal the process of gene loss , we also identified putative gain of genes by horizontal transfer from other organisms ( Table 1 ) .', "Plant-parasitic nematodes modulate their host's metabolism and induce development of feeding sites ( for example induction of syncytia by cyst nematodes , and of giant cells by root-knot nematodes ) .", 'These modifications involve the secretion by the nematode of exoenzymes such as pectinases , proteinases and cellulases ( reviewed by Vanholme and colleagues [80] ) .', 'Putative effectors have been identified using directed cloning of nematode secretory gland products , including beta-1 , 4-endoglucanases from Globodera rostochiensis [81] , Heterodera schachtii [82] and Meloidogyne incognita [83] .', 'Analyses of plant-parasitic EST data also identified beta-1 , 4-endoglucanase , beta-1 , 4-xylanases [84] and pectate lyases [85] .', 'We identified two Meloidogyne orthologues ( M . javanica and M . hapla ) of a polygalacturonase previously reported from M . incognita [83] .', 'Beta-1 , 4-endoglucanases were identified in seven species , including Pratylenchus vulnus .', "The enzyme's presence across most Tylenchid genera studied ( missing in the small Rhadopholus similis dataset ) suggests that the acquisition of this endoglucanase gene occurred in an ancient tylench ancestor .", 'We identified seven additional protein families from plant parasitic nematodes that are similar to enzymes found in plants but not previously identified in non-nematode metazoans .', 'The activities that may be carried out by these genes fall into two classes .', 'Four genes , all from Tylenchomorpha , are enzymes that catabolise plant cell wall or starch carbohydrates ( polygalacturonase , beta-amylase and cellulase ) , and may mediate parasite modification or digestion of the root cell walls .', 'Three genes , from the dorylaim X . index and the tylenchomorph M . incognita , encode activities that could modify plant signaling or second metabolites ( flavonol synthase , scopoletin glucosyltransferase and polyneuridine-aldehyde esterase ) , and may represent ‘anti-immunity’ mediators secreted by the parasite in order to subvert the necrotic or other responses of the host .', 'Another mechanism of “gene gain” is de novo evolution of functional proteins .', 'While it is clear that this mechanism has been active on the scale of phyla and kingdoms , its ongoing role in genome evolution is unclear [86] .', 'We identified 56 , 407 protein families ( including 94 , 343 proteins ) restricted to nematodes ( NR families ) .', 'Analyses of novel proteins in other species have shown that they are characterized by a significant reduction in average length compared to proteins with homologues in other taxa [65] .', 'However , the average length of the NR family proteins ( 200 aa ) is only slightly shorter than those with homologues elsewhere ( 220 aa ) .', 'It might be expected that novel genes would be expressed at low levels , and that they might thus be indistinguishable from aberrant transcripts from non-coding regions of the genome .', 'Over 80% of the NR families contained an EST-derived sequence; not restricted to the caenorhabditids .', 'Of these 69% were derived from a single EST ( data not shown ) .', 'For loss-gain candidate protein families , 68% were derived from a single EST , while of families with matches in C . elegans and elsewhere , only 35% were derived from single ESTs .', 'Thus , while the NR family sequences are expressed at low levels compared to core nematode genes , their expression levels are comparable to those of genes with wide phylogenetic distribution .', 'We analyzed further the 2 , 098 NR families with at least five members .', 'The number of NR families restricted to each taxonomic family or species correlated well with the depth of sequencing for each taxon ( Table S2 ) .', 'We note that despite cogent evidence for gene loss in the caenorhabditids [74]–[76] , many NR families with a disjoint distribution in Nematoda are likely to be present in additional species , but as yet unsampled by ESTs .', 'For example , 388 protein families ( 2 , 985 proteins ) were restricted to the complete proteomes of the caenorhabditids ( Family Rhabditoidea ) .', 'The lack of homologues in other nematodes is likely to result in part from the depth of EST sampling , as only 1 , 385 ( 46% ) of these proteins had corresponding C . elegans ESTs ( out of 346 , 064 EST sequences ) .', 'All nine nematode taxonomic families in this study had taxon-restricted protein families .', 'For example , of 35 protein families that were restricted to Spiruromorpha , only three were species-specific ( one restricted to B . malayi and two to O . volvulus; data not shown ) .', 'Fourteen of the spiruromorph protein families occurred in four species and one ( NemFam3 family 3 . 0_3062 ) contained all five species ( whose multiple sequence alignment is shown in Figure S1 ) .', 'Many ( 630 ) NR protein families with at least five members did not contain a protein from the complete proteome of C . elegans .', 'The processes of ‘gene invention’ ( and high rate of protein evolution ) are ongoing in Nematoda .', 'Indeed , the preponderance of apparently species-specific proteins is just what we would predict from this process , given the pull towards new functions , and thus may not be simply due to lack of representation in EST data .', 'However , compared with our previous analysis [11] , many sequences once thought to be species-specific now have inferred nodes of origin deeper in the nematode phylogeny , and we would expect this trend to continue as additional data are collected .', 'It has been hypothesized that the secreted subset of parasitic nematode proteomes may be especially enriched in novel proteins , through rapid evolution to perform novel functions such as interactions with the host and other environmental challenges [54] , [55] .', 'The protein families restricted to the nematodes were significantly enriched for signal peptides ( 19% ) compared to those that had homologues in other phyla ( 12% ) ( Figure 5 ) .', 'Within the class of nematode protein families that did have homologues in other phyla ( non-NR ) , 2 , 490 proteins ( 28% ) were predicted to have signal peptides .', 'Surprisingly , aligning these signal peptide-containing nematode proteins to homologues from other phyla revealed that 1 , 883 nematode proteins ( from 856 NemFam3 families , both NR and non-NR ) appear to have gained an N-terminal signal peptide .', 'For two thirds of these protein families , C . elegans and C . briggsae proteins do not contain a signal peptide , suggesting that the acquisition of a signal peptide did not occur in the caenorhabditid lineage .', 'The T . circumcincta proteome was the most enriched with signal peptides in both nematode-restricted and shared proteins .', 'Mapping these T . circumcincta proteins onto NR families identified 48 strongylomorph-restricted families where signal peptide-containing proteins predominated .', 'Despite the incomplete sampling of nematode protein space it is likely that many of these protein families are involved in specializations of the parasitic mode of life in strongylids .', 'Domains are the basic functional and structural units of proteins and , while primary sequence diversity is expected to be huge , the diversity of domains has been predicted to be rather small [87] , [88] .', 'As novel genes are being evolved in nematodes , we predicted that there might be de novo or accelerated evolution of protein domains .', 'Identification of protein domains typically involves comparing sequences to a library of protein domain alignments [32] , [33] , [35] .', 'These alignments are characterized either as hidden Markov models ( HMM ) or position-specific scoring matrices ( PSSM ) .', 'Such an approach is well suited for full-length sequences , where a match , global ( i . e . full-length ) with respect to the domain , is usually considered necessary .', 'However , proteomes derived from EST projects contain incomplete sequences , where only part of the domain is present making these global searches problematic .', 'In particular it is difficult to robustly recognize domains that extend over the termini of partial translations .', 'We devised a heuristic approach to assigning domain presence , based on different scoring thresholds available for domain models , in order to return a high coverage of domain annotation while keeping number of false positives to a minimum ( see Materials and Methods ) .', 'The resulting nematode domain classification ( NemDom3 ) contained 39 , 944 unique domains ( Table 2 ) of which 2 , 593 were from PfamA and 10 , 684 from ProDom .', 'The majority of these domains were derived from the complete caenorhabditid genomes , but more than half were found in the EST-derived proteome .', 'Previously , 348 PfamA domains had been identified in non-caenorhabditid nematodes .', 'We found 2 , 300 PfamA domain matches in the EST-derived proteomes of which 214 domains ( increased from thirteen ) were absent in C . elegans and C . briggsae .', 'All but eight of these domains were exclusive to protein sequences that we had already identified as loss/gain candidates ( described above ) , including those restricted to plant-parasites: cellulase ( PF00150 ) and pectate lyase ( PF03211 ) .', 'Of the eight domains identified in protein families that include Caenorhabditis sp .', 'members , two of these , domains associated with the ribosomal large subunit protein 6 ( PF03868 ) and NADH:ubiqunione oxidoreductase ( PF08122 ) , have been reported in C . elegans [89] .', 'However their sequences have been so diverged from the domain model as not to be recognized .', 'Seventy-seven PfamA domains were found only in nematodes , with six found in species other than exclusive to C . elegans or C . briggsae ( Table S3 ) .', 'With the exception of the abundant larval transcript ( ALT ) domain ( PF05535 ) , all nematode-restricted ( NR ) domains were first identified in C . elegans [90]–[93] .', 'Surprisingly , we were able to expand the species-distribution in only 24 of the 77 domains .', 'It is possible that the remaining NR domains are restricted to the caenorhabditid lineage .', 'However , it is more likely that many , if not most , are present in other nematode species , but were not yet represented in EST data , or were not recognized by domain models that were too constrained .', 'Inspection of the multiple sequence alignments of caenorhabditid-specific NR domains revealed often extremely high levels of identity .', 'These alignments may generate hidden Markov models ( HMMs ) that cannot identify more divergent members .', 'To illustrate this , we returned to the ALT domain ( PF05535 ) , which was , expectedly , identified in proteins from filarial species , but the searches did not find the known instance in C . elegans [90] , [91] .', 'Using the Pfam alignment for this domain ( based on five filarial sequences ) , we constructed a PSSM and performed a RPS-BLAST search .', 'This identified ALT domains in C . elegans as well as predicted proteins from Ascaris suum and A . lumbricoides .', 'We defined over 23 , 000 protein domains seemingly unique to nematodes .', 'Nearly half of these are found in non-caenorhabditid species .', 'Many of these new domains are found as part of multi-domain architectures , with 15 , 152 ( 65% ) present with at least one different domain ( all classes ) and 6 , 625 associated with a PfamA domain .', 'Profile searches with these novel domains ( see Methods ) identified 3 , 694 domains that matched non-nematode UniProt proteins .', 'The most common distribution of these domains was the 270 domains found throughout the Ecdysozoa .', 'However many domains had disjoint distributions , such as the 56 novel domains apparently exclusive to the nematodes and Viridplantae .', 'Ten of these domains were found in plant-parasite nematode species ( Table 3 ) .', 'The presence of putative homologues for three of these domains in C . elegans confuses of the issue of their origin .', 'The absence of these domains in other metazoans suggests that they were either acquired through horizontal gene transfer or diverged from an ancestral nematode domain .', 'Convergent evolution has been reported previously in nematodes [94] , [95] .', 'Are these domains real , conserved units ?', 'Of the 1 , 652 novel domains that were exclusive to the Spiruromorpha , 824 were found in at least two species of t
['Introduction', 'Materials and Methods', 'Results/Discussion']
['The phylum Nematoda is biologically diverse , including parasites of plants and animals as well as free-living taxa .', 'Underpinning this diversity will be commensurate diversity in expressed genes , including gene sets associated specifically with evolution of parasitism .', 'Here we have analyzed the extensive expressed sequence tag data ( available for 37 nematode species , most of which are parasites ) and define over 120 , 000 distinct putative genes from which we have derived robust protein translations .', 'Combined with the complete proteomes of Caenorhabditis elegans and Caenorhabditis briggsae , these proteins have been grouped into 65 , 000 protein families that in turn contain 40 , 000 distinct protein domains .', 'We have mapped the occurrence of domains and families across the Nematoda and compared the nematode data to that available for other phyla .', 'Gene loss is common , and in particular we identify nearly 5 , 000 genes that may have been lost from the lineage leading to the model nematode C . elegans .', 'We find a preponderance of novelty , including 56 , 000 nematode-restricted protein families and 26 , 000 nematode-restricted domains .', 'Mapping of the latest time-of-origin of these new families and domains across the nematode phylogeny revealed ongoing evolution of novelty .', 'A number of genes from parasitic species had signatures of horizontal transfer from their host organisms , and parasitic species had a greater proportion of novel , secreted proteins than did free-living ones .', 'These classes of genes may underpin parasitic phenotypes , and thus may be targets for development of effective control measures .']
['The high-throughput sequencing of messenger RNA from parasitic organisms has permitted large-scale sequence analyses typically reserved for complete genome studies .', 'Such expressed sequence tags ( ESTs ) have previously been generated for 37 species from the phylum Nematoda , of which 35 were from parasitic species .', 'These datasets were combined with the complete genomes of Caenorhabditis elegans and C . briggsae .', 'The sequences were assembled into 65 , 000 protein families , and decorated with 40 , 000 distinct protein domains .', 'These annotations were analysed in the context of the nematode phylogeny .', 'We identified massive gene loss in the model nematode , C . elegans , as well as plant-like proteins in nematodes that cause crop damage .', 'Furthermore , many protein families were found in small groups of closely related species and may represent innovations necessary to sustain their parasitic ecologies .', 'All of these data are presented at NemBase ( www . nematodes . org ) and will aid researchers working on this important group of parasites .']
['genetics and genomics/gene discovery', 'evolutionary biology/evolutionary and comparative genetics', 'infectious diseases/helminth infections', 'evolutionary biology/bioinformatics', 'computational biology/genomics', 'genetics and genomics/bioinformatics']