<title>Introduction</title> <p>The intestine is one of the major organs in nematodes, creating a key surface at the intestinal apical membrane that interacts with the environment. While specific cellular characteristics of the intestine can be diverse among nematode species, they typically conform to polarized epithelial cells with the apical membrane composed of microvilli lining the digestive tube. In apparent contrast to other surfaces of nematodes, digestive and assimilative functions, as well as various metabolic pathways and cellular trafficking, are expected to be extremely active at the intestinal surface. For example, an adult <italic>Caenorhabditis elegans</italic> is capable of producing oocytes with about the same total biomass as its own body per day ##UREF##0##[1]##, but the average intestinal residence time for foods was estimated to be less than two minutes in <italic>C. elegans</italic>
##UREF##1##[2]##, suggesting that the microvillous membrane must have an enormous capacity for nutrient digestion and absorption. In addition, the intestine has to offer innate immunity against invasive pathogens, and adaptations at the apical intestinal membrane may be required to protect parasitic nematodes against host immune systems. Furthermore, the nematode intestine has been suggested to be involved in other biological processes such as stress responses, body size control, and aging ##UREF##0##[1]##.</p> <p>Three lines of evidence indicate that the intestine is an important target for the control of parasitic nematodes. First, intestinal antigens enriched for apical membrane-associated proteins have been successfully used to immunize against <italic>Haemonchus contortus</italic>, a hematophagous nematode of small ruminants ##REF##1937800##[3]##–##REF##11522403##[7]##. Surface-bound nematode proteases are a dominant, but not exclusive, group of proteins that have been implicated in inducing this protection. A prospective mechanism of the immunity involves perturbing nutrient digestion and acquisition at the intestinal surface by the ingested host-derived antibodies capable of neutralizing parasite digestive proteases ##REF##11522403##[7]##. Further investigations conducted with hematophagous hookworms also produced similar effects ##REF##17067922##[8]##. Second, adult <italic>H. contortus</italic> intestinal cells are hypersensitive to benzimidazole anthelmintics, apparently through the target protein beta-tubulin isotype 1 ##REF##10613701##[9]##,##REF##12594960##[10]##. It was suggested that the drug inhibited vesicle transport in the apical secretory pathway, causing the intracellular release of the digestive enzymes destined for secretion and subsequent cytotoxic effects ##REF##10613701##[9]##. Third, parasite control has been demonstrated by inhibition of an intestinal enzyme, cathepsin L cysteine protease, by either RNA interference or a chemical inhibitor in the plant parasitic nematode <italic>Meloidogyne incognita</italic>
##REF##16996059##[11]##. These observations generate great interests to uncover the basic characteristics of the intestinal cells that might be further exploited for the broad control of parasitic nematodes. However, the dearth of relevant experimental systems and molecular information such as gene repertoires for many parasitic species has impeded rapid progress.</p> <p>Five major clades (I–V) are currently recognized to comprise the phylum Nematoda ##REF##17084644##[12]##,##REF##9510248##[13]##. So far, almost all studies of the intestine at the gene level have focused on the clade V nematodes. A small-scale sampling of expressed sequence tags (ESTs) from the dissected intestine from adult <italic>H. contortus</italic> females identified 51 intestinal genes including cysteine proteases ##REF##11522349##[14]##, this list was later expanded via a proteomic approach to include a number of apical intestinal membrane proteases from <italic>H. contortus</italic> and hookworms ##REF##12957519##[15]##. Intestinal EST libraries generated from laser-dissected materials from <italic>Necator americanus</italic> and <italic>Ancylostoma caninum</italic> allowed the identification of 544 intestine-expressed genes ##REF##16545815##[16]##. Although a more comprehensive dataset with &gt;5,000 intestinal genes is available in <italic>C. elegans</italic>
##UREF##2##[17]##–##REF##16354718##[19]##, it is unclear, given the evolutionary diversity within Nematoda, to what extent the molecular and cellular functions of the intestine can be extrapolated across nematode species.</p> <p>In this study, we sampled the transcribed genomes from several tissues and developmental stages from two parasitic nematodes: the clade III nematode <italic>Ascaris suum</italic>, which presumably feeds on the semi-digested contents in the host intestine, and the clade V blood-feeding parasite <italic>H. contortus</italic>. Nearly 10,000 and 5,000 genes were identified from the two nematodes, respectively. More importantly, given the attention to the intestine, we produced the largest collection of intestinal genes in parasitic nematodes by dissecting adult intestine from each species, a procedure that is not practical for many other nematodes because of their small sizes and the lack of laboratory culturing systems. Extensive cross-species comparisons were made among the adult intestinal genes from the parasites and those expressed in the adult intestine of the free-living bacterivore <italic>C. elegans</italic>. Both diversification and conservation of intestinal gene repertories were evident among the species investigated. The diversities of intestinal transcriptomes by clade and species may reflect the substantial life style differences among these nematodes. A group of 241 protein families were found conserved in the intestine of all three nematodes, accounting for ∼20% of the intestinal gene repertoires from the three species. These genes may include core intestinal functions that are indispensable among many nematodes. Functional annotations were generated for the intestinal genes. Molecular characteristics of the intestinal genes were further explored to highlight various physiological aspects of the nematode intestine.</p>
<title>Materials and Methods</title> <title>EST generation, clustering, and translation</title> <p>Dissection of the adult intestine was carefully performed under microscopy as described previously ##REF##11522349##[14]##,##REF##10498185##[20]##,##REF##9879887##[21]##. The samples used in this study had also passed another round of visual inspection microscopically to ensure they did not contain other tissues such as muscle, esophagus, or hypodermis. Detailed information on genetic materials and cDNA library construction are available at <ext-link ext-link-type="uri" xlink:href="http://www.nematode.net">www.nematode.net</ext-link>. ESTs were processed and clustered as described before ##REF##16697384##[22]##–##REF##14762059##[25]##. EST contig sequences were translated individually by Prot4EST, a 6-tier translation pipeline combining both similarity-based methods and <italic>de novo</italic> predictions ##REF##15571632##[26]##, for downstream analysis.</p> <title>Identification of sequence similarities</title> <p>Databases used for sequence comparisons were<italic>:</italic> i) <italic>Caenorhabditis spp.</italic>, all amino acid sequences in the complete genomes of <italic>C. elegans</italic> (Wormbase Release v150), <italic>C. briggsae</italic> (June, 2006), and <italic>C. remanei</italic> (June, 2006), ii) Other Nematoda, all non-<italic>Caenorhabditis</italic> nematode nucleic acid sequences in GenBank excluding those from <italic>A. suum</italic> (when analyzing <italic>A. suum</italic> sequences) or <italic>H. contortus</italic> (when querying <italic>H. contortus</italic> sequences) (October 18, 2006), and iii) Non-Nematoda, all amino acid sequences in the non-redundant protein database NR excluding those from nematode species (September 20, 2006). WU-BLASTP (wordmask = seg postwe B = 1000 topcomboN = 1) was used to query the translated sequences against protein databases, and WU-TBLASTN (wordmask = seg lcmask B = 1000 topcomboN = 1) for searching against nucleotide databases ##REF##2231712##[27]##. The <italic>E-</italic>value cutoff of 1.0e<sup>−5</sup> was used to accept sequence similarities in all BLAST searches.</p> <title>Identification of intestine-enriched genes</title> <p>Each intestinal EST cluster was assigned two counts according to the numbers of times it was sampled from either the intestinal or non-intestinal cDNA libraries, respectively. Similarly, each <italic>C. elegans</italic> intestinal gene was assigned two counts for the numbers of times it was sampled by SAGE tags from either the <italic>glp-4</italic> dissected gut or the <italic>glp-4</italic> adult whole worm, respectively. The mutants lack the gonad when raised at 25°C, therefore contamination by other tissues is less likely ##UREF##2##[17]##. The SAGE data was downloaded with sequence quality filter = 0.99, no normalization, duplicate ditags and ambiguous or antisense tags removed (April 19, 2006; mapped to Wormbase Release v150) ##UREF##2##[17]##. A Poisson-based enrichment test, considering both the total sampling sizes and random variations ##REF##9331369##[28]##, was implemented to compute an <italic>P</italic>-value to represent the likelihood of intestinal enrichment for each EST cluster or <italic>C. elegans</italic> gene using these two counts. The <italic>P</italic>-value cutoff of 0.001 was chosen to define the putative intestine-enriched genes from the three nematodes.</p> <title>Prediction of signal peptide (SP) and transmembrane (TM) domain</title> <p>A hidden Markov modeling-based algorithm, Phobius ##REF##15111065##[29]##, was used with default setting. Each query sequence was further annotated as TM-only, TM with SP, SP-only, or intracellular based on raw Phobius outputs. For each EST cluster, Phobius annotation was predicted for each contig and summarized at the EST cluster level. A modified Wormbase Release v150 containing only the longest splicing isoform at each gene loci was used as the complete gene set of the <italic>C. elegans</italic> genome.</p> <title>Identification of orthogous gene pairs between <italic>A. suum</italic> and <italic>C. elegans</italic>
</title> <p>For tissue-level comparisons made between intestine and gonad, InParanoid ##REF##11743721##[30]## was used at default settings to identify a total of 1,764 putative orthologous groups between all the <italic>A. suum</italic> EST clusters and the complete gene set of <italic>C. elegans</italic> (the modified Wormbase Release v150 containing only the longest splicing isoform at each gene loci). InParanoid-generated main orthologous pairs, which are essentially the mutual-best matches between all the available genes from the two species, were further screened against the 447,546 <italic>A. suum</italic> Genome Survey Sequences (GSSs) that were generated recently (Mitreva, unpublished), resulting in the final group of 1,652 putative main orthologous pairs in which the <italic>C. elegans</italic> members do not have better matches in GSSs than the <italic>A. suum</italic> EST partners assigned by InParanoid. <italic>C. elegans</italic> gonad-expressed genes were extracted from SAGE data generated from dissected gonad (March 12, 2007) ##UREF##2##[17]##.</p> <title>Identification of intestinal protein families</title> <p>An all-against-all WU-BLASTP was performed on all the 9,918 translated intestinal genes from the three species (including sequences for EST contigs from the two parasites and 5,056 <italic>C. elegans</italic> genes). Raw BLAST results were fed to a C-language implementation of Markov Cluster (MCL) Algorithm (<ext-link ext-link-type="uri" xlink:href="http://www.micans.org/mcl">www.micans.org/mcl</ext-link>), a fast and scalable unsupervised cluster algorithm based on simulation of flow in graphs ##REF##11917018##[31]##. An Inflation Fact of 1.6 was chosen for the MCL clustering. The MCL output was then summarized at the EST cluster level, during which we applied an additional filtering step to remove an EST cluster from a MCL protein family if less than 10% of its total contigs were clustered into that family. These parameters were based on manual inspection of the results on a test set consisting of the putative intestine-enriched genes with 210 parasite EST clusters and 247 <italic>C. elegans</italic> genes (false positive rate of 3%; data not shown).</p> <title>Gene Ontology mappings and identification of statistically enriched ontologies</title> <p>Default parameters for InterProScan v13.1 ##REF##15980438##[32]## were used to search against the InterPro database ##UREF##4##[33]##. Raw InterProScan results for the translated EST contigs were summarized at the EST cluster level. Gene ontology (GO) terms were further assigned and displayed graphically by the AmiGO browser with default parameters and the ontology data released on March 15, 2007##REF##10802651##[34]##. Complete GO mappings for the three intestinal transcriptomes are available at <ext-link ext-link-type="uri" xlink:href="http://www.nematode.net">www.nematode.net</ext-link>. For each GO term, its enrichment in an IntFam group (such as the IntFam-241 group) was measured over the complete set of 9,918 translated intestinal genes using a hypergeometric test, the <italic>p</italic>-value cutoff of 1.0e<sup>−5</sup> was chosen for enrichment. The less informative ontologies, including those at level 4 or higher for Biological Process or Molecular Function, and those at level 2 or higher for Cellular Component, were removed from the enrichment list. Also removed were redundant ontologies by keeping only the lower level more informative ontology if the same group of genes was mapped to more than one GO term.</p> <title>KEGG pathway analysis</title> <p>An empirical mixed approach was used for mapping the novel genes to canonical pathways. The <italic>E</italic>-value cut-off of 1.0e<sup>−10</sup> reported by WU-BLASTP against the Genes Database Release 39.0 from Kyoto Encyclopedia of Genes and Genomes (KEGG) was first used for finding homologous matches. Then the top match and all the matches within a range of 30% of the top BLAST score, if meeting the cut-off, were accepted for valid KEGG associations ##REF##9521924##[35]##–##REF##14681412##[37]##. A hypergeometric test, measuring the relative coverage of the KEGG-annotated orthologous groups assigned to a pathway, was implemented to identify the enriched pathways for each intestine ##REF##15817693##[38]##.</p> <title>Accession numbers</title> <p>Nucleotide sequences data reported in this paper are available in the GenBank, EMBL and DDBJ databases. The accession numbers for ESTs from <italic>A. suum</italic> are: BI781215-BI784439, BM032617-BM034650, BM280443-BM285290, BM318846-BM319958, BM515079-BM518821, BM566483-BM567588, BM568416-BM569529, BM732977-BM734435, BM964439-BM965448, BQ094886-BQ096565, BQ380669-BQ383404, BQ835081-BQ835723, BU965907-BU966430, CA849193-CA850481, CA953713-CA955182, CB100077-CB102042, DV018957-DV019894, EB186562-EB187079. The accession numbers for ESTs from <italic>H. contortus</italic> are: CA033335-CA034379, CA868595-CA870175, CA956361-CA959150, CB018493-CB022024, CB063882-CB065260, CB099467-CB100076, CB190871-CB192419, CB331948-CB333475.</p>
<title>Results/Discussion</title> <title>Generation and Clustering of ESTs from Clade III and V Parasitic Nematodes</title> <p>We constructed 18 <italic>A. suum</italic> and 6 <italic>H. contortus</italic> stage- or tissue-specific cDNA libraries, and sequenced 31,416 and 14,014 5-prime ESTs from the two species, respectively. These ESTs totaled to 13.6 and 6.3 million bases for <italic>A. suum</italic> and <italic>H. contortus</italic>, accounting for 77% and 63% of the total nucleotides from the two species currently available in public databases (##SUPPL##2##Table S1##). Supplemented by 9,354 <italic>A. suum</italic> and 8,146 <italic>H. contortus</italic> ESTs previously deposited in GenBank (retrieved in January, 2006), all available ESTs were grouped into 17,989 <italic>A. suum</italic> and 9,842 <italic>H. contortus</italic> EST contigs, each containing ESTs derived from nearly identical transcripts according to overlapping sequences to reduce sequence redundancy ##UREF##3##[23]##,##REF##15338281##[24]##. The contigs were further assembled into 9,947 <italic>A. suum</italic> and 5,058 <italic>H. contortus</italic> EST clusters based on sequence similarities identified among contigs as well as in previously identified genes (##SUPPL##2##Table S1##). Each EST cluster likely represents transcripts derived from a single genomic locus and therefore is approximated as one gene ##REF##16697384##[22]##–##REF##15338281##[24]##. Given that <italic>C. elegans</italic> and <italic>C. briggsae</italic> each contains ∼19,000 protein-coding loci, and between 14,500 and 17,800 genes were inferred from the <italic>Brugia malayi</italic> draft genome ##REF##17885136##[39]##, we have consequently identified a substantial portion of the complete gene sets from the two parasites. These data will vastly facilitate the genome assembly and annotation in the related nematode genome sequencing projects currently underway. Initial investigation of the identities of these novel genes was performed by comparing the translated sequences with known proteins from other organisms (##SUPPL##10##Text S1##; ##SUPPL##0##Figure S1##).</p> <title>Intestinal Transcriptomes from Adult <italic>A. suum</italic>, <italic>H. contortus</italic>, and <italic>C. elegans</italic>
</title> <p>To study the intestinal transcriptomes, four cDNA libraries (out of the 18) from <italic>A. suum</italic> and three (out of the 6) from <italic>H. contortus</italic> were constructed from dissected adult intestine with methods based on either Poly-A ##REF##15383298##[40]## or spliced leader sequences ##REF##15338281##[24]##. Among all the ESTs we generated, a total of 9,586 <italic>A. suum</italic> and 7,068 <italic>H. contortus</italic> ESTs were derived from these intestinal libraries. These ESTs occurred in 3,121 <italic>A. suum</italic> and 1,755 <italic>H. contortus</italic> EST clusters, accounting for about 30% of the total genes sampled in each nematode. Since these EST clusters contained ESTs sampled from the adult intestine, they were considered to represent adult intestinal genes, making this the largest tissue-level gene discovery in parasitic nematodes thus far (##TAB##0##Table 1##).</p> <p>In contrast to the two gastrointestinal parasites, the free-living model nematode <italic>C. elegans</italic> is a bacterivore obtaining nutrients primarily or exclusively from the consumption of bacteria. Two previous studies reported identification of genes expressed in the adult <italic>C. elegans</italic> intestine: i) sequence tags generated by serial analysis of gene expression (SAGE) from the dissected adult intestine were mapped to over 4,000 <italic>C. elegans</italic> genes ##UREF##2##[17]##,##REF##17113066##[18]##; ii) a study using mRNA tagging and microarray gene expression profiling identified ∼1,900 intestine-expressed genes ##REF##16354718##[19]##. Consolidating the two efforts provided us with a non-redundant set of 5,065 intestinal genes from adult <italic>C. elegans</italic>, covering over 25% of all coding loci in its entire genome (##TAB##0##Table 1##).</p> <p>The phylum Nematoda is ancient and diverse. Even though the evolutionary distance between clade III <italic>A. suum</italic> and clade V <italic>C. elegans</italic> was estimated to be ∼350 million years ##REF##8075836##[41]##, the nematode intestine has maintained high similarity in both tissue morphology and presumably physiology (i.e. involvement in feeding). However, it is unknown how much the intestine is conserved, or diversified, at the molecular level across species. The tissue-level gene sampling in this study offered an opportunity to investigate this question.</p> <title>Diversification among Intestinal Transcriptomes</title> <p>Differences in the intestinal gene repertoires were obvious among the three nematodes. In total, 39% of <italic>A. suum</italic> and 19% of <italic>H. contortus</italic> intestinal genes were found to be novel compared to all known proteins in the public databases (##FIG##0##Figure 1##). Such novel intestine-expressed parasite genes contained no match in the complete genome of the free-living <italic>C. elegans</italic>, thus not in the <italic>C. elegans</italic> intestine, making them unique by comparison to <italic>C. elegans</italic>. In addition, for the sampled intestinal genes from both parasites, the non-<italic>Caenorhabditis</italic> nematodes offered the largest numbers of homologous matches than either the <italic>Caenorhabditis</italic> species or the non-nematode organisms (##FIG##0##Figure 1##). Such differences may suggest the existence of lineage- or species-specific diversification in the nematode intestine.</p> <p>Furthermore, we observed higher levels of diversification in the putative intestine-enriched genes from the three nematodes. Taking into consideration sample size and random sampling fluctuation ##REF##9331369##[28]##, we identified 150 <italic>A. suum</italic>, 60 <italic>H. contortus</italic>, and 247 <italic>C. elegans</italic> putative intestine-enriched genes based on the “digital” expression levels revealed in EST and SAGE data (at the Poisson distribution-based <italic>P</italic>-value cutoff of 0.001) (##SUPPL##3##Table S2##; ##SUPPL##4##Table S3##; ##SUPPL##5##Table S4##). Many of these predicted enrichments suggested unique intestinal functions for the individual species. For example, the group of 60 genes from the blood-feeding <italic>H. contortus</italic> includes 2 fibrinogen-related proteins that may function as thrombin inhibitors to prevent clotting of ingested blood. Also included are putative enzymes that may be involved in the digestion of hemoglobin, one of the major food sources of blood-feeding parasites, including a serine-type protease, a metallopeptidase, and 13 different cysteine-type proteases that were reported previously ##REF##15383300##[42]## (##SUPPL##4##Table S3##). Interestingly, a significantly higher percentages of these genes (e.g. 15%-31% higher than all the sampled intestinal genes) encode proteins predicted as secreted or trans-membrane ##REF##15111065##[29]## (##FIG##1##Figure 2##), suggesting that they interact with the extracellular environment. However, 64%, 54%, and 69% of them, from the three species respectively, were distinct from members of the protein families conserved in the intestine of all three nematodes (IntFam-241; see below), indicating that a large portion of these putative intestine-enriched genes are specific to the intestine of individual nematode lineages or species. This further underlines the diversification of intestinal transcriptomes in accommodating the different life styles and feeding patterns among nematodes.</p> <title>Molecular Conservation of the Nematode Intestine</title> <p>To evaluate common characteristics of the nematode intestine, we first sought evidence for the molecular conservation of the tissue in the context of phylogeny. We made comparisons among genes expressed in the intestine of <italic>A. suum</italic> and <italic>C. elegans</italic> and those expressed in another tissue, namely the gonad. These two species have the largest numbers of sequences available, and they also represent the most distant relationship among the three nematodes investigated. The gonad was chosen because the next largest group of genes was sampled from this tissue in <italic>A. suum</italic> after the intestine. <italic>H. contortus</italic> was excluded from this analysis because a gonad-expressed gene set was not available from this nematode. Genes expressed in the intestine and gonad were divided into four putative tissue-specific groups: i) 2,453 <italic>A. suum</italic> and ii) 2,557 <italic>C. elegans</italic> genes expressed in the intestine but not in the gonad (the two intestine groups), and iii) 2,690 <italic>A. suum</italic> and iv) 2,589 <italic>C. elegans</italic> genes that were found in the gonad but not in the intestine (the two gonad groups). The use of the similar numbers of genes in each group is expected to reduce false results caused by over-representation from any single category.</p> <p>Molecular conservation was first evaluated by comparing the numbers of putative homologous pairs identified among the intestine and gonad gene groups. The number of the putative homologs between the two intestine groups was significantly larger than that between the intestine and gonad groups (<italic>p</italic>-value = 2.5e<sup>−04</sup> at the bit-score cutoff of 100 in a permutation two-tailed Z-test; <italic>p</italic>-value = 4.2e<sup>−08</sup> at the bit-score cutoff of 50; ##SUPPL##1##Figure S2##), suggesting that for genes expressed in the intestine of one nematode, their homologous matches in another species are significantly more likely to be expressed in the intestine than in the gonad of the second nematode. These results provide evidence for the molecular conservation of the intestine across these distantly related nematodes. In contrast, the number of putative homologs between the two gonad groups was not statistically different from that between the gonad and intestine (##SUPPL##1##Figure S2##), indicating that the gonad genes appeared to be less conserved than those expressed in the intestine in this two-tissue comparison.</p> <p>To increase the confidence of analysis, we next focused on the putative orthologous pairs predicted among the intestine and gonad gene groups, which was a smaller data set than the homologous pairs used above but with higher stringency. Among the total of 1,652 putative orthologous pairs predicted from <italic>A. suum</italic> and <italic>C. elegans</italic> (see <xref ref-type="sec" rid="s2">Materials and Methods</xref>), 289 were paired among genes from the intestine and gonad groups. They were used in a Chi Square statistical test, with random distribution of orthologous pairs as the null hypothesis. Compared to the expected numbers, there was a 31% enrichment of orthologous pairs observed between the <italic>A. suum</italic> and <italic>C. elegans</italic> intestine groups (##FIG##2##Figure 3##), whereas the enrichment between the two gonad groups was only marginal (5%), and the observed numbers of orthologous pairs between the gonad and intestine groups were less than expected (##FIG##2##Figure 3##). Overall, a significant χ<sup>2</sup> value of 11.9 rejects the null hypothesis at a confidence level higher than 99% (<italic>p</italic>-value &lt;0.01) ##UREF##5##[43]##, and selective pressure is evident on molecular conservation of the intestinal gene repertories.</p> <p>Although the use of the incomplete transcriptomes and a bias towards relatively abundant transcripts in EST sampling can affect results, analyses of either homologous or orthologous pairs both provide direct support for the molecular conservation of the nematode intestine. With the obvious pattern of diversification in the nematode intestine (discussed earlier), our results indicate that a subset of the intestinal gene repertoires, which likely contribute to the intestinal characteristics conserved across diverse nematode species, remain conserved during the evolution of Nematoda.</p> <p>Interestingly, genes expressed in the gonad appear to be less well conserved based on both analyses. However, these results do not suggest the lack of evidence for the conservation of the gonad. Instead, the two-tissue comparisons indicate that the levels of conservation are lower in the gonad than in the intestine, suggesting that the levels of molecular conservation may differ in different nematode tissues. In fact, the conserved characteristics of the gonad may become more evident with larger sample sizes and/or by comparisons with another tissue with a lower level of conservation than the intestine, when new sequence data becomes available. Similarly, differences at the levels of molecular conservation were observed in different tissues between human and mouse, which diverged only about 25 million years ago ##REF##17301236##[44]##. Future comparisons with more complete expression data across multiple tissues in different nematode species should offer additional insights into this aspect of nematode evolution.</p> <title>Identification of “Core” and Other Groups of Intestinal Protein Families</title> <p>To compare the intestinal transcriptomes of <italic>A. suum</italic>, <italic>H. contortus</italic>, and <italic>C. elegans</italic> in a single analysis, we built protein families from the complete set of 9,918 translated intestinal genes combined from the three nematodes. A total of 5,587 intestinal protein families (IntFam) were identified conservatively based on sequence similarities by MCL clustering ##REF##11917018##[31]## (##FIG##3##Figure 4##). Proteins assigned into the same protein family contain putative homologous or orthologous matches among the three species. Both diversification and conservation of the intestinal transcriptomes was obvious at the protein family level in this 3-species comparison. A total of 59% of all the sampled intestinal genes were members of the protein families containing proteins from only one nematode (##FIG##3##Figure 4##). Although the assignments for many of these single-species families are likely to change when more complete intestinal gene repertories become available, this group includes the genes contributing to the unique intestinal features in each species. The remaining 41% of the intestinal genes formed 910 multi-species protein families; they are conserved in the intestine of at least two nematodes. Among these multi-species families, 241 had members from all three species, accounting for ∼20% of all the intestinal genes under investigation (##FIG##3##Figure 4##). Given the differences in life styles and feeding patterns among the three nematodes, we propose that these 241 intestinal protein families represent an ancestral intestinal transcriptome involved in core cellular and physiological intestinal functions common to the investigated species or even across the Nematoda. Therefore, we referred to them as the “core” IntFam-241 group.</p> <title>Functional Analysis of Intestinal Genes Based on Gene Ontology</title> <p>The 9,918 translated intestinal genes sampled from the three nematodes were annotated and classified using Gene Ontology ##REF##10802651##[34]##,##REF##12654719##[45]##. Ontologies were assigned at a higher ratio (58%) to the <italic>C. elegans</italic> intestinal genes than to those from <italic>A. suum</italic> (31%) or <italic>H. contortus</italic> (35%; ##TAB##0##Table 1##). In addition, genes in the multi-species IntFam groups, which contained members from at least two nematodes, were annotated at higher ratios (47%–74%), whereas only 8% of the genes were annotated from the two single-species IntFam groups containing members only from <italic>A. suum</italic> or <italic>H. contortus</italic> (data not shown). These data may indicate that novel intestinal genes have independently evolved in relation to the different lineages of parasitism. Complete GO mappings for the three intestinal transcriptomes are presented in the searchable AmiGO browser at <ext-link ext-link-type="uri" xlink:href="http://www.nematode.net">www.nematode.net</ext-link>
##REF##14681448##[46]##. Furthermore, A hypergeometric test was implemented to identify ontologies that are statistically enriched, thus indicating enriched features, in the core IntFam-241 (##TAB##1##Table 2##) as well as other IntFam groups (##SUPPL##10##Text S1##; ##SUPPL##7##Table S6##).</p> <p>Five of the 17 enriched Molecular Function ontologies in IntFam-241 are related to protein kinases (##TAB##1##Table 2##; ##TAB##2##Table 3##). Protein kinases are one of the largest and most influential protein families, accounting for about 2% of genes in a variety of eukaryotic genomes including <italic>C. elegans</italic> and <italic>B. malayi</italic>. They regulate almost every aspect of cellular activities and may phosphorylate up to 30% of entire proteomes ##REF##17885136##[39]##,##REF##10570119##[47]##. Based on GO annotations, protein kinases were enriched by ∼3.5 fold in IntFam-241 over the complete set of intestinal genes (5.3% vs. 1.5% of the total genes for each group). Both serine/threonine- and tyrosine-types of protein kinases were found to be enriched. Novel protein kinases from the parasites were further classified based on their <italic>C. elegans</italic> homologs (##SUPPL##6##Table S5##). Interestingly, molecular functions such as adenyl nucleotide binding, ATP binding, and GTP binding were also enriched. The involvement of these functions in protein kinase activities further suggested key roles of cellular signaling in the nematode intestine (##TAB##1##Table 2##).</p> <p>The other major Molecular Function terms enriched in IntFam-241 were the proteases (##TAB##1##Table 2##; ##TAB##2##Table 3##). All but one of the six subtypes of proteases (glutamic acid-type proteases as the exception) had been identified in IntFam-241 (##TAB##2##Table 3##), suggesting conservation of essential protease functions, such as nutrient digestion and acquisition, among the three species or even across many species of Nematoda. Because each species feeds on distinct food sources, it is possible that related digestive proteases have evolved within each species to adapt for digestion of the different food types. Given the success of parasite control achieved by immunization with <italic>H. contortus</italic> and hookworm intestinal protease-type antigens ##REF##1937800##[3]##,##REF##17067922##[8]##, these proteases may warrant further investigation in <italic>A. suum</italic> and other parasites.</p> <p>Analysis of the IntFam groups other than IntFam-241 was also conducted. However, in absence of deeper sampling of the intestinal transcriptomes, it is difficult to interpret the results in relation to broadly conserved or lineage- and species-specific characteristics (##SUPPL##10##Text S1##; ##SUPPL##7##Table S6##).</p> <title>KEGG-based Pathway Analysis of Three Intestinal Transcriptomes</title> <p>To identify the biological pathways that are active in the nematode intestine, we mapped the 9,918 intestinal translated sequences, and for comparison, the complete <italic>C. elegans</italic> genes (Wormbase Release v150), to the reference canonical pathways in Kyoto Encyclopedia of Genes and Genomes ##REF##9521924##[35]##–##REF##14681412##[37]## (##TAB##3##Table 4##). Complete listing of all KEGG mappings including graphical representation is available for navigation at <ext-link ext-link-type="uri" xlink:href="http://www.nematode.net">www.nematode.net</ext-link>
##REF##14681448##[46]##.</p> <p>The enrichment of specific major KEGG pathways was evident for each intestine by comparisons to the complete KEGG mappings for all <italic>C. elegans</italic> genes (##TAB##3##Table 4##) ##REF##15817693##[38]##. Carbohydrate metabolism, energy metabolism, and translation were identified as the statistically enriched pathways in all three intestinal transcriptomes (at the <italic>p</italic>-value cutoff of 0.05). Interestingly, immune system was an enriched KEGG cellular process in the <italic>A. suum</italic> intestine; this pathway barely missed the cutoff for enrichment in <italic>H. contortus</italic> (with a <italic>p</italic>-value of 9.9e<sup>−02</sup>), but no enrichment was indicated for the <italic>C. elegance</italic> intestine (##TAB##3##Table 4##). The KEGG immune system was built based on studies in mammalian systems. Many of those from the two parasites were mapped to intracellular proteins of immune cells involved in, for example, intracellular signaling or antigen processing (##SUPPL##8##Table S7##; ##SUPPL##9##Table S8##). Therefore, the potential for their involvement in interactions with the host are not a primary suggestion here, but it cannot be completely excluded either.</p> <title>
<italic>C. elegans</italic> RNA Interference and Intestinal Genes</title> <p>RNA interference (RNAi) has been developed and successfully applied to genome-wide gene silencing to inhibit gene functions in <italic>C. elegans</italic>
##REF##9486653##[48]##–##REF##15791247##[51]##. <italic>C. elegans</italic> RNAi information can be further extrapolated in understanding functions of orthologous genes in other nematodes, especially in parasitic nematodes where high-throughput screening is not yet practical ##REF##16099532##[52]##. For the 3,455 IntFam protein families containing <italic>C. elegans</italic> genes, observed RNAi phenotypes for their <italic>C. elegans</italic> members (Wormbase Release v150) were extracted and extrapolated to a total of 45% of these IntFams (##TAB##4##Table 5##). Protein families from the IntFam-241 were assigned at a higher ratio (73%) than those from other IntFam groups with <italic>C. elegans</italic> members (##TAB##4##Table 5##). Among the IntFams-241 families with RNAi phenotypes assigned, 74% (131/176) had severe phenotypes including embryonic, larval, or adult lethal, sterile, sterile progeny, and larval or adult growth arrest (data not shown). Since the IntFam-241 families represent proteins conserved in all the three species, these results further support our hypothesis that the core IntFam-241 protein families likely play key roles in the nematode intestine across many species.</p> <title>Summary</title> <p>We have performed large-scale sampling of the transcribed genomes in <italic>A. suum</italic> and <italic>H. contortus</italic> from various tissues or developmental stages, accounting for 77% and 63% of total available bases for the two nematodes, respectively. The identification of 9,947 <italic>A. suum</italic> and 5,058 <italic>H. contortus</italic> genes in this study will vastly facilitate the related genome sequencing projects currently underway. The research has produced the largest samplings of the adult intestinal transcriptomes thus far in parasitic nematodes by identifying 3,121 <italic>A. suum</italic> and 1,755 <italic>H. contortus</italic> intestinal genes, making possible the extensive comparative studies with the adult intestinal transcriptome of the free-living <italic>C. elegans</italic>. We found that, even with the evolutionary distance of an estimated 350 million years separating clades III and V nematodes ##REF##8075836##[41]##, both significant conservation and diversification of gene repertories were evident for the intestine. A group of 241 intestinal protein families, each containing members from all three intestines, were further identified. The IntFam-241 group, containing ∼20% of all intestinal genes sampled from the three species, was proposed to represent an ancient intestinal transcriptome responsible for core cellular and physiological intestinal functions that are conserved in the investigated species or many other nematodes. In addition, various aspects of nematode intestinal physiology were revealed by GO and KEGG classifications of the intestinal transcriptomes, and the examination and extrapolation of available RNAi phenotypes from <italic>C. elegans</italic>. Overall, this study has contributed to a better understanding of nematode biology, providing central information for the development of novel and more effective parasite control strategies. Finally, the use of the <italic>C. elegans</italic> model to dissect basic parasite biology has been slow to evolve. Results presented here identified numerous specific areas of research where <italic>C. elegans</italic> might contribute in this way.</p>
<title>Results/Discussion</title> <title>Generation and Clustering of ESTs from Clade III and V Parasitic Nematodes</title> <p>We constructed 18 <italic>A. suum</italic> and 6 <italic>H. contortus</italic> stage- or tissue-specific cDNA libraries, and sequenced 31,416 and 14,014 5-prime ESTs from the two species, respectively. These ESTs totaled to 13.6 and 6.3 million bases for <italic>A. suum</italic> and <italic>H. contortus</italic>, accounting for 77% and 63% of the total nucleotides from the two species currently available in public databases (##SUPPL##2##Table S1##). Supplemented by 9,354 <italic>A. suum</italic> and 8,146 <italic>H. contortus</italic> ESTs previously deposited in GenBank (retrieved in January, 2006), all available ESTs were grouped into 17,989 <italic>A. suum</italic> and 9,842 <italic>H. contortus</italic> EST contigs, each containing ESTs derived from nearly identical transcripts according to overlapping sequences to reduce sequence redundancy ##UREF##3##[23]##,##REF##15338281##[24]##. The contigs were further assembled into 9,947 <italic>A. suum</italic> and 5,058 <italic>H. contortus</italic> EST clusters based on sequence similarities identified among contigs as well as in previously identified genes (##SUPPL##2##Table S1##). Each EST cluster likely represents transcripts derived from a single genomic locus and therefore is approximated as one gene ##REF##16697384##[22]##–##REF##15338281##[24]##. Given that <italic>C. elegans</italic> and <italic>C. briggsae</italic> each contains ∼19,000 protein-coding loci, and between 14,500 and 17,800 genes were inferred from the <italic>Brugia malayi</italic> draft genome ##REF##17885136##[39]##, we have consequently identified a substantial portion of the complete gene sets from the two parasites. These data will vastly facilitate the genome assembly and annotation in the related nematode genome sequencing projects currently underway. Initial investigation of the identities of these novel genes was performed by comparing the translated sequences with known proteins from other organisms (##SUPPL##10##Text S1##; ##SUPPL##0##Figure S1##).</p> <title>Intestinal Transcriptomes from Adult <italic>A. suum</italic>, <italic>H. contortus</italic>, and <italic>C. elegans</italic>
</title> <p>To study the intestinal transcriptomes, four cDNA libraries (out of the 18) from <italic>A. suum</italic> and three (out of the 6) from <italic>H. contortus</italic> were constructed from dissected adult intestine with methods based on either Poly-A ##REF##15383298##[40]## or spliced leader sequences ##REF##15338281##[24]##. Among all the ESTs we generated, a total of 9,586 <italic>A. suum</italic> and 7,068 <italic>H. contortus</italic> ESTs were derived from these intestinal libraries. These ESTs occurred in 3,121 <italic>A. suum</italic> and 1,755 <italic>H. contortus</italic> EST clusters, accounting for about 30% of the total genes sampled in each nematode. Since these EST clusters contained ESTs sampled from the adult intestine, they were considered to represent adult intestinal genes, making this the largest tissue-level gene discovery in parasitic nematodes thus far (##TAB##0##Table 1##).</p> <p>In contrast to the two gastrointestinal parasites, the free-living model nematode <italic>C. elegans</italic> is a bacterivore obtaining nutrients primarily or exclusively from the consumption of bacteria. Two previous studies reported identification of genes expressed in the adult <italic>C. elegans</italic> intestine: i) sequence tags generated by serial analysis of gene expression (SAGE) from the dissected adult intestine were mapped to over 4,000 <italic>C. elegans</italic> genes ##UREF##2##[17]##,##REF##17113066##[18]##; ii) a study using mRNA tagging and microarray gene expression profiling identified ∼1,900 intestine-expressed genes ##REF##16354718##[19]##. Consolidating the two efforts provided us with a non-redundant set of 5,065 intestinal genes from adult <italic>C. elegans</italic>, covering over 25% of all coding loci in its entire genome (##TAB##0##Table 1##).</p> <p>The phylum Nematoda is ancient and diverse. Even though the evolutionary distance between clade III <italic>A. suum</italic> and clade V <italic>C. elegans</italic> was estimated to be ∼350 million years ##REF##8075836##[41]##, the nematode intestine has maintained high similarity in both tissue morphology and presumably physiology (i.e. involvement in feeding). However, it is unknown how much the intestine is conserved, or diversified, at the molecular level across species. The tissue-level gene sampling in this study offered an opportunity to investigate this question.</p> <title>Diversification among Intestinal Transcriptomes</title> <p>Differences in the intestinal gene repertoires were obvious among the three nematodes. In total, 39% of <italic>A. suum</italic> and 19% of <italic>H. contortus</italic> intestinal genes were found to be novel compared to all known proteins in the public databases (##FIG##0##Figure 1##). Such novel intestine-expressed parasite genes contained no match in the complete genome of the free-living <italic>C. elegans</italic>, thus not in the <italic>C. elegans</italic> intestine, making them unique by comparison to <italic>C. elegans</italic>. In addition, for the sampled intestinal genes from both parasites, the non-<italic>Caenorhabditis</italic> nematodes offered the largest numbers of homologous matches than either the <italic>Caenorhabditis</italic> species or the non-nematode organisms (##FIG##0##Figure 1##). Such differences may suggest the existence of lineage- or species-specific diversification in the nematode intestine.</p> <p>Furthermore, we observed higher levels of diversification in the putative intestine-enriched genes from the three nematodes. Taking into consideration sample size and random sampling fluctuation ##REF##9331369##[28]##, we identified 150 <italic>A. suum</italic>, 60 <italic>H. contortus</italic>, and 247 <italic>C. elegans</italic> putative intestine-enriched genes based on the “digital” expression levels revealed in EST and SAGE data (at the Poisson distribution-based <italic>P</italic>-value cutoff of 0.001) (##SUPPL##3##Table S2##; ##SUPPL##4##Table S3##; ##SUPPL##5##Table S4##). Many of these predicted enrichments suggested unique intestinal functions for the individual species. For example, the group of 60 genes from the blood-feeding <italic>H. contortus</italic> includes 2 fibrinogen-related proteins that may function as thrombin inhibitors to prevent clotting of ingested blood. Also included are putative enzymes that may be involved in the digestion of hemoglobin, one of the major food sources of blood-feeding parasites, including a serine-type protease, a metallopeptidase, and 13 different cysteine-type proteases that were reported previously ##REF##15383300##[42]## (##SUPPL##4##Table S3##). Interestingly, a significantly higher percentages of these genes (e.g. 15%-31% higher than all the sampled intestinal genes) encode proteins predicted as secreted or trans-membrane ##REF##15111065##[29]## (##FIG##1##Figure 2##), suggesting that they interact with the extracellular environment. However, 64%, 54%, and 69% of them, from the three species respectively, were distinct from members of the protein families conserved in the intestine of all three nematodes (IntFam-241; see below), indicating that a large portion of these putative intestine-enriched genes are specific to the intestine of individual nematode lineages or species. This further underlines the diversification of intestinal transcriptomes in accommodating the different life styles and feeding patterns among nematodes.</p> <title>Molecular Conservation of the Nematode Intestine</title> <p>To evaluate common characteristics of the nematode intestine, we first sought evidence for the molecular conservation of the tissue in the context of phylogeny. We made comparisons among genes expressed in the intestine of <italic>A. suum</italic> and <italic>C. elegans</italic> and those expressed in another tissue, namely the gonad. These two species have the largest numbers of sequences available, and they also represent the most distant relationship among the three nematodes investigated. The gonad was chosen because the next largest group of genes was sampled from this tissue in <italic>A. suum</italic> after the intestine. <italic>H. contortus</italic> was excluded from this analysis because a gonad-expressed gene set was not available from this nematode. Genes expressed in the intestine and gonad were divided into four putative tissue-specific groups: i) 2,453 <italic>A. suum</italic> and ii) 2,557 <italic>C. elegans</italic> genes expressed in the intestine but not in the gonad (the two intestine groups), and iii) 2,690 <italic>A. suum</italic> and iv) 2,589 <italic>C. elegans</italic> genes that were found in the gonad but not in the intestine (the two gonad groups). The use of the similar numbers of genes in each group is expected to reduce false results caused by over-representation from any single category.</p> <p>Molecular conservation was first evaluated by comparing the numbers of putative homologous pairs identified among the intestine and gonad gene groups. The number of the putative homologs between the two intestine groups was significantly larger than that between the intestine and gonad groups (<italic>p</italic>-value = 2.5e<sup>−04</sup> at the bit-score cutoff of 100 in a permutation two-tailed Z-test; <italic>p</italic>-value = 4.2e<sup>−08</sup> at the bit-score cutoff of 50; ##SUPPL##1##Figure S2##), suggesting that for genes expressed in the intestine of one nematode, their homologous matches in another species are significantly more likely to be expressed in the intestine than in the gonad of the second nematode. These results provide evidence for the molecular conservation of the intestine across these distantly related nematodes. In contrast, the number of putative homologs between the two gonad groups was not statistically different from that between the gonad and intestine (##SUPPL##1##Figure S2##), indicating that the gonad genes appeared to be less conserved than those expressed in the intestine in this two-tissue comparison.</p> <p>To increase the confidence of analysis, we next focused on the putative orthologous pairs predicted among the intestine and gonad gene groups, which was a smaller data set than the homologous pairs used above but with higher stringency. Among the total of 1,652 putative orthologous pairs predicted from <italic>A. suum</italic> and <italic>C. elegans</italic> (see <xref ref-type="sec" rid="s2">Materials and Methods</xref>), 289 were paired among genes from the intestine and gonad groups. They were used in a Chi Square statistical test, with random distribution of orthologous pairs as the null hypothesis. Compared to the expected numbers, there was a 31% enrichment of orthologous pairs observed between the <italic>A. suum</italic> and <italic>C. elegans</italic> intestine groups (##FIG##2##Figure 3##), whereas the enrichment between the two gonad groups was only marginal (5%), and the observed numbers of orthologous pairs between the gonad and intestine groups were less than expected (##FIG##2##Figure 3##). Overall, a significant χ<sup>2</sup> value of 11.9 rejects the null hypothesis at a confidence level higher than 99% (<italic>p</italic>-value &lt;0.01) ##UREF##5##[43]##, and selective pressure is evident on molecular conservation of the intestinal gene repertories.</p> <p>Although the use of the incomplete transcriptomes and a bias towards relatively abundant transcripts in EST sampling can affect results, analyses of either homologous or orthologous pairs both provide direct support for the molecular conservation of the nematode intestine. With the obvious pattern of diversification in the nematode intestine (discussed earlier), our results indicate that a subset of the intestinal gene repertoires, which likely contribute to the intestinal characteristics conserved across diverse nematode species, remain conserved during the evolution of Nematoda.</p> <p>Interestingly, genes expressed in the gonad appear to be less well conserved based on both analyses. However, these results do not suggest the lack of evidence for the conservation of the gonad. Instead, the two-tissue comparisons indicate that the levels of conservation are lower in the gonad than in the intestine, suggesting that the levels of molecular conservation may differ in different nematode tissues. In fact, the conserved characteristics of the gonad may become more evident with larger sample sizes and/or by comparisons with another tissue with a lower level of conservation than the intestine, when new sequence data becomes available. Similarly, differences at the levels of molecular conservation were observed in different tissues between human and mouse, which diverged only about 25 million years ago ##REF##17301236##[44]##. Future comparisons with more complete expression data across multiple tissues in different nematode species should offer additional insights into this aspect of nematode evolution.</p> <title>Identification of “Core” and Other Groups of Intestinal Protein Families</title> <p>To compare the intestinal transcriptomes of <italic>A. suum</italic>, <italic>H. contortus</italic>, and <italic>C. elegans</italic> in a single analysis, we built protein families from the complete set of 9,918 translated intestinal genes combined from the three nematodes. A total of 5,587 intestinal protein families (IntFam) were identified conservatively based on sequence similarities by MCL clustering ##REF##11917018##[31]## (##FIG##3##Figure 4##). Proteins assigned into the same protein family contain putative homologous or orthologous matches among the three species. Both diversification and conservation of the intestinal transcriptomes was obvious at the protein family level in this 3-species comparison. A total of 59% of all the sampled intestinal genes were members of the protein families containing proteins from only one nematode (##FIG##3##Figure 4##). Although the assignments for many of these single-species families are likely to change when more complete intestinal gene repertories become available, this group includes the genes contributing to the unique intestinal features in each species. The remaining 41% of the intestinal genes formed 910 multi-species protein families; they are conserved in the intestine of at least two nematodes. Among these multi-species families, 241 had members from all three species, accounting for ∼20% of all the intestinal genes under investigation (##FIG##3##Figure 4##). Given the differences in life styles and feeding patterns among the three nematodes, we propose that these 241 intestinal protein families represent an ancestral intestinal transcriptome involved in core cellular and physiological intestinal functions common to the investigated species or even across the Nematoda. Therefore, we referred to them as the “core” IntFam-241 group.</p> <title>Functional Analysis of Intestinal Genes Based on Gene Ontology</title> <p>The 9,918 translated intestinal genes sampled from the three nematodes were annotated and classified using Gene Ontology ##REF##10802651##[34]##,##REF##12654719##[45]##. Ontologies were assigned at a higher ratio (58%) to the <italic>C. elegans</italic> intestinal genes than to those from <italic>A. suum</italic> (31%) or <italic>H. contortus</italic> (35%; ##TAB##0##Table 1##). In addition, genes in the multi-species IntFam groups, which contained members from at least two nematodes, were annotated at higher ratios (47%–74%), whereas only 8% of the genes were annotated from the two single-species IntFam groups containing members only from <italic>A. suum</italic> or <italic>H. contortus</italic> (data not shown). These data may indicate that novel intestinal genes have independently evolved in relation to the different lineages of parasitism. Complete GO mappings for the three intestinal transcriptomes are presented in the searchable AmiGO browser at <ext-link ext-link-type="uri" xlink:href="http://www.nematode.net">www.nematode.net</ext-link>
##REF##14681448##[46]##. Furthermore, A hypergeometric test was implemented to identify ontologies that are statistically enriched, thus indicating enriched features, in the core IntFam-241 (##TAB##1##Table 2##) as well as other IntFam groups (##SUPPL##10##Text S1##; ##SUPPL##7##Table S6##).</p> <p>Five of the 17 enriched Molecular Function ontologies in IntFam-241 are related to protein kinases (##TAB##1##Table 2##; ##TAB##2##Table 3##). Protein kinases are one of the largest and most influential protein families, accounting for about 2% of genes in a variety of eukaryotic genomes including <italic>C. elegans</italic> and <italic>B. malayi</italic>. They regulate almost every aspect of cellular activities and may phosphorylate up to 30% of entire proteomes ##REF##17885136##[39]##,##REF##10570119##[47]##. Based on GO annotations, protein kinases were enriched by ∼3.5 fold in IntFam-241 over the complete set of intestinal genes (5.3% vs. 1.5% of the total genes for each group). Both serine/threonine- and tyrosine-types of protein kinases were found to be enriched. Novel protein kinases from the parasites were further classified based on their <italic>C. elegans</italic> homologs (##SUPPL##6##Table S5##). Interestingly, molecular functions such as adenyl nucleotide binding, ATP binding, and GTP binding were also enriched. The involvement of these functions in protein kinase activities further suggested key roles of cellular signaling in the nematode intestine (##TAB##1##Table 2##).</p> <p>The other major Molecular Function terms enriched in IntFam-241 were the proteases (##TAB##1##Table 2##; ##TAB##2##Table 3##). All but one of the six subtypes of proteases (glutamic acid-type proteases as the exception) had been identified in IntFam-241 (##TAB##2##Table 3##), suggesting conservation of essential protease functions, such as nutrient digestion and acquisition, among the three species or even across many species of Nematoda. Because each species feeds on distinct food sources, it is possible that related digestive proteases have evolved within each species to adapt for digestion of the different food types. Given the success of parasite control achieved by immunization with <italic>H. contortus</italic> and hookworm intestinal protease-type antigens ##REF##1937800##[3]##,##REF##17067922##[8]##, these proteases may warrant further investigation in <italic>A. suum</italic> and other parasites.</p> <p>Analysis of the IntFam groups other than IntFam-241 was also conducted. However, in absence of deeper sampling of the intestinal transcriptomes, it is difficult to interpret the results in relation to broadly conserved or lineage- and species-specific characteristics (##SUPPL##10##Text S1##; ##SUPPL##7##Table S6##).</p> <title>KEGG-based Pathway Analysis of Three Intestinal Transcriptomes</title> <p>To identify the biological pathways that are active in the nematode intestine, we mapped the 9,918 intestinal translated sequences, and for comparison, the complete <italic>C. elegans</italic> genes (Wormbase Release v150), to the reference canonical pathways in Kyoto Encyclopedia of Genes and Genomes ##REF##9521924##[35]##–##REF##14681412##[37]## (##TAB##3##Table 4##). Complete listing of all KEGG mappings including graphical representation is available for navigation at <ext-link ext-link-type="uri" xlink:href="http://www.nematode.net">www.nematode.net</ext-link>
##REF##14681448##[46]##.</p> <p>The enrichment of specific major KEGG pathways was evident for each intestine by comparisons to the complete KEGG mappings for all <italic>C. elegans</italic> genes (##TAB##3##Table 4##) ##REF##15817693##[38]##. Carbohydrate metabolism, energy metabolism, and translation were identified as the statistically enriched pathways in all three intestinal transcriptomes (at the <italic>p</italic>-value cutoff of 0.05). Interestingly, immune system was an enriched KEGG cellular process in the <italic>A. suum</italic> intestine; this pathway barely missed the cutoff for enrichment in <italic>H. contortus</italic> (with a <italic>p</italic>-value of 9.9e<sup>−02</sup>), but no enrichment was indicated for the <italic>C. elegance</italic> intestine (##TAB##3##Table 4##). The KEGG immune system was built based on studies in mammalian systems. Many of those from the two parasites were mapped to intracellular proteins of immune cells involved in, for example, intracellular signaling or antigen processing (##SUPPL##8##Table S7##; ##SUPPL##9##Table S8##). Therefore, the potential for their involvement in interactions with the host are not a primary suggestion here, but it cannot be completely excluded either.</p> <title>
<italic>C. elegans</italic> RNA Interference and Intestinal Genes</title> <p>RNA interference (RNAi) has been developed and successfully applied to genome-wide gene silencing to inhibit gene functions in <italic>C. elegans</italic>
##REF##9486653##[48]##–##REF##15791247##[51]##. <italic>C. elegans</italic> RNAi information can be further extrapolated in understanding functions of orthologous genes in other nematodes, especially in parasitic nematodes where high-throughput screening is not yet practical ##REF##16099532##[52]##. For the 3,455 IntFam protein families containing <italic>C. elegans</italic> genes, observed RNAi phenotypes for their <italic>C. elegans</italic> members (Wormbase Release v150) were extracted and extrapolated to a total of 45% of these IntFams (##TAB##4##Table 5##). Protein families from the IntFam-241 were assigned at a higher ratio (73%) than those from other IntFam groups with <italic>C. elegans</italic> members (##TAB##4##Table 5##). Among the IntFams-241 families with RNAi phenotypes assigned, 74% (131/176) had severe phenotypes including embryonic, larval, or adult lethal, sterile, sterile progeny, and larval or adult growth arrest (data not shown). Since the IntFam-241 families represent proteins conserved in all the three species, these results further support our hypothesis that the core IntFam-241 protein families likely play key roles in the nematode intestine across many species.</p> <title>Summary</title> <p>We have performed large-scale sampling of the transcribed genomes in <italic>A. suum</italic> and <italic>H. contortus</italic> from various tissues or developmental stages, accounting for 77% and 63% of total available bases for the two nematodes, respectively. The identification of 9,947 <italic>A. suum</italic> and 5,058 <italic>H. contortus</italic> genes in this study will vastly facilitate the related genome sequencing projects currently underway. The research has produced the largest samplings of the adult intestinal transcriptomes thus far in parasitic nematodes by identifying 3,121 <italic>A. suum</italic> and 1,755 <italic>H. contortus</italic> intestinal genes, making possible the extensive comparative studies with the adult intestinal transcriptome of the free-living <italic>C. elegans</italic>. We found that, even with the evolutionary distance of an estimated 350 million years separating clades III and V nematodes ##REF##8075836##[41]##, both significant conservation and diversification of gene repertories were evident for the intestine. A group of 241 intestinal protein families, each containing members from all three intestines, were further identified. The IntFam-241 group, containing ∼20% of all intestinal genes sampled from the three species, was proposed to represent an ancient intestinal transcriptome responsible for core cellular and physiological intestinal functions that are conserved in the investigated species or many other nematodes. In addition, various aspects of nematode intestinal physiology were revealed by GO and KEGG classifications of the intestinal transcriptomes, and the examination and extrapolation of available RNAi phenotypes from <italic>C. elegans</italic>. Overall, this study has contributed to a better understanding of nematode biology, providing central information for the development of novel and more effective parasite control strategies. Finally, the use of the <italic>C. elegans</italic> model to dissect basic parasite biology has been slow to evolve. Results presented here identified numerous specific areas of research where <italic>C. elegans</italic> might contribute in this way.</p>
<p>Conceived and designed the experiments: YY RKW MM. Analyzed the data: YY JM SA MM. Contributed reagents/materials/analysis tools: YY JM SA ALS JPM DJ. Wrote the paper: YY ALS JPM RKW DJ MM.</p> <title>Background</title> <p>The nematode intestine is a major organ responsible for nutrient digestion and absorption; it is also involved in many other processes, such as reproduction, innate immunity, stress responses, and aging. The importance of the intestine as a target for the control of parasitic nematodes has been demonstrated. However, the lack of detailed knowledge on the molecular and cellular functions of the intestine and the level of its conservation across nematodes has impeded breakthroughs in this application.</p> <title>Methods and Findings</title> <p>As part of an extensive effort to investigate various transcribed genomes from <italic>Ascaris suum</italic> and <italic>Haemonchus contortus</italic>, we generated a large collection of intestinal sequences from parasitic nematodes by identifying 3,121 <italic>A. suum</italic> and 1,755 <italic>H. contortus</italic> genes expressed in the adult intestine through the generation of expressed sequence tags. Cross-species comparisons to the intestine of the free-living <italic>C. elegans</italic> revealed substantial diversification in the adult intestinal transcriptomes among these species, suggesting lineage- or species-specific adaptations during nematode evolution. In contrast, significant conservation of the intestinal gene repertories was also evident, despite the evolutionary distance of ∼350 million years separating them. A group of 241 intestinal protein families (IntFam-241), each containing members from all three species, was identified based on sequence similarities. These conserved proteins accounted for ∼20% of the sampled intestinal transcriptomes from the three nematodes and are proposed to represent conserved core functions in the nematode intestine. Functional characterizations of the IntFam-241 suggested important roles in molecular functions such as protein kinases and proteases, and biological pathways of carbohydrate metabolism, energy metabolism, and translation. Conservation in the core protein families was further explored by extrapolating observable RNA interference phenotypes in <italic>C. elegans</italic> to their parasitic counterparts.</p> <title>Conclusions</title> <p>Our study has provided novel insights into the nematode intestine and lays foundations for further comparative studies on biology, parasitism, and evolution within the phylum Nematoda.</p> <title>Author Summary</title> <p>Biological properties of the nematode intestine warrant in-depth investigation, the results of which can be utilized in the control of parasitic nematodes that infect humans, livestock, and plants. Both the importance of intestinal antigens from <italic>Haemonchus contortus</italic> in immunity and the damage to <italic>H. contortus</italic> intestine by anthelmintic fenbendazole have highlighted the versatility of the intestine as an emerging target. However, biological information regarding fundamental intestinal cell functions and mechanisms is currently limited. Conserved intestinal genes across nematode pathogens could offer molecular targets for broad parasite control. Furthermore, qualitative and quantitative comparisons on intestinal gene expression among species and lineages can identify basic adaptations relative to a critical selective force, the nutrient acquisition. This study begins to identify intestinal cell characteristics that are conserved across representatives of two clades of nematodes (V and III) and further clarifies diversities that likely reflect species- or lineage-specific adaptations. Results consistent with functional data on digestive enzymes from <italic>H. contortus</italic> and RNAi in <italic>Caenorhabditis elegans</italic>, as examples, support the potential for the comparative genomics approach to produce practical applications. This study provides a platform on which extensive investigation of intestinal genes and a more comprehensive understanding of the Nematoda can be gained.</p>
<title>Supporting Information</title>
<p>We thank Dr. Dante Zarlenga, Dr. Timothy Geary, Dr. Richard Davis for providing some materials, Brandi Chiapelli, and Claire Murphy for constructing cDNA libraries, Mike Dante for performing EST processing and clustering, Todd Wylie, David Messina, Dr. Jarret Glasscock, and Dr. Ken Chen for technical assistance and scientific discussions.</p>
<fig id="pntd-0000269-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pntd.0000269.g001</object-id><label>Figure 1</label><caption><title>Sequence Similarities Identified in the <italic>A. suum</italic> and <italic>H. contortus</italic> Intestinal Transcriptomes.</title><p>The three phylogenetically specific sequence groups used to identify sequence similarities of the intestinal genes were: i) <italic>Caenorhabditis spp.</italic>, amino acid sequences from the complete genomes of <italic>C. elegans</italic>, <italic>C. briggsae</italic>, and <italic>C. remanei</italic>, ii) Other Nematoda, non-<italic>Caenorhabditis</italic> nematode nucleic acid sequences excluding those from either <italic>A. suum</italic> or <italic>H. contortus</italic>, when sequences from <italic>A. suum</italic> or <italic>H. contortus</italic> were queried, respectively, and iii) Non-Nematoda, non-nematode amino acid sequences from the non-redundant protein database NR. In total, 61% (1,893/3,121) <italic>A. suu</italic>m and 81% (1,423/1,755) <italic>H. contortus</italic> intestinal genes contained primary sequence similarities to known proteins from other organisms.</p></caption></fig> <fig id="pntd-0000269-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pntd.0000269.g002</object-id><label>Figure 2</label><caption><title>Putative Secreted or Trans-membrane Proteins in <italic>A. suum</italic>, <italic>H. contortus</italic>, and <italic>C. elegans</italic>.</title><p>Larger percentages of the putative intestine-enriched genes (Int. Enriched) were predicted to be secreted with signal peptide (SP) or trans-membrane (TM) than either the complete set of intestinal genes (Int. Genes) or the complete set of all available genes (All Genes) in <italic>A. suum</italic>, <italic>H. contortus</italic>, and <italic>C. elegans</italic>.</p></caption></fig> <fig id="pntd-0000269-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pntd.0000269.g003</object-id><label>Figure 3</label><caption><title>Orthologous Genes Tend to Maintain Their Intestinal Expression Patterns across <italic>A. suum</italic> and <italic>C. elegans</italic>.</title><p>A total of 289 putative orthologous pairs were identified among the intestine or gonad gene groups from <italic>A. suum</italic> and <italic>C. elegans</italic>. Ninety such pairs were found between the two intestine gene groups, representing an enrichment of 31% over the expectation from a random distribution of orthologous pairs, and an enrichment of only 5% was detected between genes in the two gonad groups. The null hypothesis of random orthologous pairing was rejected at a confidence level of at least 99% with a χ<sup>2</sup> value of 11.9 between the observations and expectations.</p></caption></fig> <fig id="pntd-0000269-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pntd.0000269.g004</object-id><label>Figure 4</label><caption><title>Protein Families in the Adult Intestine from <italic>A. suum</italic>, <italic>H. contortus</italic>, and <italic>C. elegans</italic>.</title><p>In total, 5,587 intestinal protein families (IntFam) were built from the complete set of 9,918 translated intestinal genes sampled from the three species. Forty-one percent of all the intestinal genes were grouped into 910 multiple-species IntFam groups; 2,024 genes, including 752 from <italic>A. suum</italic>, 455 from <italic>H. contortus</italic>, and 817 from <italic>C. elegans</italic>, were found in a group of 241 families (IntFam-241) containing members from all three intestines. The IntFam-241 families likely represent an ancestral intestinal transcriptome involved in core cellular and physiological intestinal functions common to the investigated species or even to Nematoda, they are thus referred as the “core” IntFam-241 group.</p></caption></fig>
<table-wrap id="pntd-0000269-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pntd.0000269.t001</object-id><label>Table 1</label><caption><title>Intestinal Transcriptomes and Intestine-Enriched Genes from Three Nematodes.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Summary</td><td align="left" rowspan="1" colspan="1">
<italic>A. suum</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>H. contortus</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>C. elegans</italic>
</td></tr></thead><tbody><tr><td colspan="4" align="left" rowspan="1">
<bold>cDNA libraries</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Non-intestinal</td><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Intestinal</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td colspan="4" align="left" rowspan="1">
<bold>ESTs/SAGE Tags </bold>
<xref ref-type="table-fn" rid="nt101">1</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Total</td><td align="left" rowspan="1" colspan="1">38,978</td><td align="left" rowspan="1" colspan="1">22,152</td><td align="left" rowspan="1" colspan="1">104,756</td></tr><tr><td align="left" rowspan="1" colspan="1">Non-intestinal</td><td align="left" rowspan="1" colspan="1">29,392</td><td align="left" rowspan="1" colspan="1">15,084</td><td align="left" rowspan="1" colspan="1">53,197</td></tr><tr><td align="left" rowspan="1" colspan="1">Intestinal</td><td align="left" rowspan="1" colspan="1">9,586</td><td align="left" rowspan="1" colspan="1">7,068</td><td align="left" rowspan="1" colspan="1">51,559</td></tr><tr><td colspan="4" align="left" rowspan="1">
<bold>EST Clusters/Genes</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Total</td><td align="left" rowspan="1" colspan="1">9,947</td><td align="left" rowspan="1" colspan="1">5,058</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Non-intestinal</td><td align="left" rowspan="1" colspan="1">6,826</td><td align="left" rowspan="1" colspan="1">3,303</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Intestinal</td><td align="left" rowspan="1" colspan="1">3,121</td><td align="left" rowspan="1" colspan="1">1,755</td><td align="left" rowspan="1" colspan="1">5,065</td></tr><tr><td align="left" rowspan="1" colspan="1"> Int. genes w/Gene Ontology</td><td align="left" rowspan="1" colspan="1">954</td><td align="left" rowspan="1" colspan="1">620</td><td align="left" rowspan="1" colspan="1">2,911</td></tr><tr><td align="left" rowspan="1" colspan="1"> % Int. genes w/Gene Ontology</td><td align="left" rowspan="1" colspan="1">31%</td><td align="left" rowspan="1" colspan="1">35%</td><td align="left" rowspan="1" colspan="1">57%</td></tr><tr><td colspan="4" align="left" rowspan="1">
<bold>Intestine-enriched EST Clusters/Genes </bold>
<xref ref-type="table-fn" rid="nt102">2</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic> p</italic>-value &lt; = 0.001</td><td align="left" rowspan="1" colspan="1">150 (10)</td><td align="left" rowspan="1" colspan="1">60 (5)</td><td align="left" rowspan="1" colspan="1">247 (7)</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pntd-0000269-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pntd.0000269.t002</object-id><label>Table 2</label><caption><title>Gene Ontology Terms Statistically Enriched among 2,024 Intestinal Genes in IntFam-241.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">GO id</td><td colspan="4" align="left" rowspan="1">IntFam-241 Genes Mapped <xref ref-type="table-fn" rid="nt103">1</xref>
</td><td align="left" rowspan="1" colspan="1">All Int. Genes Mapped</td><td align="left" rowspan="1" colspan="1">
<italic>P</italic>_value</td><td align="left" rowspan="1" colspan="1">Gene Ontology</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<italic>As</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>Hc</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>Ce</italic>
</td><td align="left" rowspan="1" colspan="1">Total</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">
<bold>GO:0008150</bold>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>biological_process</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0044267</td><td align="left" rowspan="1" colspan="1">222</td><td align="left" rowspan="1" colspan="1">165</td><td align="left" rowspan="1" colspan="1">246</td><td align="left" rowspan="1" colspan="1">633</td><td align="left" rowspan="1" colspan="1">1,081</td><td align="left" rowspan="1" colspan="1">0.0E+00</td><td align="left" rowspan="1" colspan="1">cellular protein metabolic process</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0006412</td><td align="left" rowspan="1" colspan="1">154</td><td align="left" rowspan="1" colspan="1">62</td><td align="left" rowspan="1" colspan="1">74</td><td align="left" rowspan="1" colspan="1">290</td><td align="left" rowspan="1" colspan="1">451</td><td align="left" rowspan="1" colspan="1">0.0E+00</td><td align="left" rowspan="1" colspan="1">translation</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0006464</td><td align="left" rowspan="1" colspan="1">36</td><td align="left" rowspan="1" colspan="1">24</td><td align="left" rowspan="1" colspan="1">104</td><td align="left" rowspan="1" colspan="1">164</td><td align="left" rowspan="1" colspan="1">297</td><td align="left" rowspan="1" colspan="1">0.0E+00</td><td align="left" rowspan="1" colspan="1">protein modification process</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0006796</td><td align="left" rowspan="1" colspan="1">34</td><td align="left" rowspan="1" colspan="1">24</td><td align="left" rowspan="1" colspan="1">103</td><td align="left" rowspan="1" colspan="1">161</td><td align="left" rowspan="1" colspan="1">284</td><td align="left" rowspan="1" colspan="1">0.0E+00</td><td align="left" rowspan="1" colspan="1">phosphate metabolic process</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0016310</td><td align="left" rowspan="1" colspan="1">29</td><td align="left" rowspan="1" colspan="1">23</td><td align="left" rowspan="1" colspan="1">99</td><td align="left" rowspan="1" colspan="1">151</td><td align="left" rowspan="1" colspan="1">240</td><td align="left" rowspan="1" colspan="1">0.0E+00</td><td align="left" rowspan="1" colspan="1">phosphorylation</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0043687</td><td align="left" rowspan="1" colspan="1">28</td><td align="left" rowspan="1" colspan="1">18</td><td align="left" rowspan="1" colspan="1">101</td><td align="left" rowspan="1" colspan="1">147</td><td align="left" rowspan="1" colspan="1">243</td><td align="left" rowspan="1" colspan="1">0.0E+00</td><td align="left" rowspan="1" colspan="1">post-translational protein modification</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0006468</td><td align="left" rowspan="1" colspan="1">16</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">82</td><td align="left" rowspan="1" colspan="1">107</td><td align="left" rowspan="1" colspan="1">149</td><td align="left" rowspan="1" colspan="1">0.0E+00</td><td align="left" rowspan="1" colspan="1">protein amino acid phosphorylation</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0007264</td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">28</td><td align="left" rowspan="1" colspan="1">46</td><td align="left" rowspan="1" colspan="1">63</td><td align="left" rowspan="1" colspan="1">9.92E-11</td><td align="left" rowspan="1" colspan="1">small GTPase mediated signal transduction</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0006457</td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">22</td><td align="left" rowspan="1" colspan="1">43</td><td align="left" rowspan="1" colspan="1">62</td><td align="left" rowspan="1" colspan="1">5.95E-09</td><td align="left" rowspan="1" colspan="1">protein folding</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0006508</td><td align="left" rowspan="1" colspan="1">19</td><td align="left" rowspan="1" colspan="1">69</td><td align="left" rowspan="1" colspan="1">49</td><td align="left" rowspan="1" colspan="1">137</td><td align="left" rowspan="1" colspan="1">281</td><td align="left" rowspan="1" colspan="1">2.50E-08</td><td align="left" rowspan="1" colspan="1">proteolysis</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0042254</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">21</td><td align="left" rowspan="1" colspan="1">25</td><td align="left" rowspan="1" colspan="1">2.46E-07</td><td align="left" rowspan="1" colspan="1">ribosome biogenesis and assembly</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0006414</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">23</td><td align="left" rowspan="1" colspan="1">29</td><td align="left" rowspan="1" colspan="1">4.69E-07</td><td align="left" rowspan="1" colspan="1">translational elongation</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0006512</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1">29</td><td align="left" rowspan="1" colspan="1">43</td><td align="left" rowspan="1" colspan="1">4.62E-06</td><td align="left" rowspan="1" colspan="1">ubiquitin cycle</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0046034</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">16</td><td align="left" rowspan="1" colspan="1">38</td><td align="left" rowspan="1" colspan="1">62</td><td align="left" rowspan="1" colspan="1">5.33E-06</td><td align="left" rowspan="1" colspan="1">ATP metabolic process</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>GO:0005575</bold>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>cellular_component</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0005737</td><td align="left" rowspan="1" colspan="1">172</td><td align="left" rowspan="1" colspan="1">72</td><td align="left" rowspan="1" colspan="1">123</td><td align="left" rowspan="1" colspan="1">367</td><td align="left" rowspan="1" colspan="1">697</td><td align="left" rowspan="1" colspan="1">0.0E+00</td><td align="left" rowspan="1" colspan="1">cytoplasm</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0044444</td><td align="left" rowspan="1" colspan="1">164</td><td align="left" rowspan="1" colspan="1">68</td><td align="left" rowspan="1" colspan="1">95</td><td align="left" rowspan="1" colspan="1">327</td><td align="left" rowspan="1" colspan="1">590</td><td align="left" rowspan="1" colspan="1">0.0E+00</td><td align="left" rowspan="1" colspan="1">cytoplasmic part</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0005840</td><td align="left" rowspan="1" colspan="1">145</td><td align="left" rowspan="1" colspan="1">52</td><td align="left" rowspan="1" colspan="1">58</td><td align="left" rowspan="1" colspan="1">255</td><td align="left" rowspan="1" colspan="1">371</td><td align="left" rowspan="1" colspan="1">1.52E-12</td><td align="left" rowspan="1" colspan="1">ribosome</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0030529</td><td align="left" rowspan="1" colspan="1">149</td><td align="left" rowspan="1" colspan="1">58</td><td align="left" rowspan="1" colspan="1">67</td><td align="left" rowspan="1" colspan="1">274</td><td align="left" rowspan="1" colspan="1">408</td><td align="left" rowspan="1" colspan="1">1.81E-12</td><td align="left" rowspan="1" colspan="1">ribonucleoprotein complex</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0043232</td><td align="left" rowspan="1" colspan="1">161</td><td align="left" rowspan="1" colspan="1">58</td><td align="left" rowspan="1" colspan="1">69</td><td align="left" rowspan="1" colspan="1">288</td><td align="left" rowspan="1" colspan="1">480</td><td align="left" rowspan="1" colspan="1">2.36E-12</td><td align="left" rowspan="1" colspan="1">intracellular non-membrane-bound organelle</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0015935</td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1">21</td><td align="left" rowspan="1" colspan="1">24</td><td align="left" rowspan="1" colspan="1">5.74E-08</td><td align="left" rowspan="1" colspan="1">small ribosomal subunit</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0005740</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">19</td><td align="left" rowspan="1" colspan="1">38</td><td align="left" rowspan="1" colspan="1">61</td><td align="left" rowspan="1" colspan="1">3.03E-06</td><td align="left" rowspan="1" colspan="1">mitochondrial envelope</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0005874</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">5.55E-06</td><td align="left" rowspan="1" colspan="1">microtubule</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0044429</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">21</td><td align="left" rowspan="1" colspan="1">40</td><td align="left" rowspan="1" colspan="1">67</td><td align="left" rowspan="1" colspan="1">7.61E-06</td><td align="left" rowspan="1" colspan="1">mitochondrial part</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0005739</td><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">26</td><td align="left" rowspan="1" colspan="1">51</td><td align="left" rowspan="1" colspan="1">92</td><td align="left" rowspan="1" colspan="1">9.02E-06</td><td align="left" rowspan="1" colspan="1">mitochondrion</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>GO:0003674</bold>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>molecular_function</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0030554</td><td align="left" rowspan="1" colspan="1">44</td><td align="left" rowspan="1" colspan="1">24</td><td align="left" rowspan="1" colspan="1">161</td><td align="left" rowspan="1" colspan="1">229</td><td align="left" rowspan="1" colspan="1">468</td><td align="left" rowspan="1" colspan="1">0.0E+00</td><td align="left" rowspan="1" colspan="1">adenyl nucleotide binding</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0016773</td><td align="left" rowspan="1" colspan="1">18</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">84</td><td align="left" rowspan="1" colspan="1">112</td><td align="left" rowspan="1" colspan="1">186</td><td align="left" rowspan="1" colspan="1">0.0E+00</td><td align="left" rowspan="1" colspan="1">phosphotransferase activity, alcohol group as acceptor</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0004672</td><td align="left" rowspan="1" colspan="1">16</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">82</td><td align="left" rowspan="1" colspan="1">107</td><td align="left" rowspan="1" colspan="1">151</td><td align="left" rowspan="1" colspan="1">0.0E+00</td><td align="left" rowspan="1" colspan="1">protein kinase activity</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0004674</td><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1">63</td><td align="left" rowspan="1" colspan="1">84</td><td align="left" rowspan="1" colspan="1">103</td><td align="left" rowspan="1" colspan="1">0.0E+00</td><td align="left" rowspan="1" colspan="1">protein serine/threonine kinase activity</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0005525</td><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1">17</td><td align="left" rowspan="1" colspan="1">40</td><td align="left" rowspan="1" colspan="1">72</td><td align="left" rowspan="1" colspan="1">103</td><td align="left" rowspan="1" colspan="1">0.0E+00</td><td align="left" rowspan="1" colspan="1">GTP binding</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0004713</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">26</td><td align="left" rowspan="1" colspan="1">39</td><td align="left" rowspan="1" colspan="1">46</td><td align="left" rowspan="1" colspan="1">0.0E+00</td><td align="left" rowspan="1" colspan="1">protein-tyrosine kinase activity</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0000062</td><td align="left" rowspan="1" colspan="1">29</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">36</td><td align="left" rowspan="1" colspan="1">36</td><td align="left" rowspan="1" colspan="1">0.0E+00</td><td align="left" rowspan="1" colspan="1">acyl-CoA binding</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0004194</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">31</td><td align="left" rowspan="1" colspan="1">32</td><td align="left" rowspan="1" colspan="1">0.0E+00</td><td align="left" rowspan="1" colspan="1">pepsin A activity</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0005524</td><td align="left" rowspan="1" colspan="1">44</td><td align="left" rowspan="1" colspan="1">25</td><td align="left" rowspan="1" colspan="1">156</td><td align="left" rowspan="1" colspan="1">225</td><td align="left" rowspan="1" colspan="1">436</td><td align="left" rowspan="1" colspan="1">1.98E-13</td><td align="left" rowspan="1" colspan="1">ATP binding</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0016301</td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">88</td><td align="left" rowspan="1" colspan="1">119</td><td align="left" rowspan="1" colspan="1">213</td><td align="left" rowspan="1" colspan="1">1.04E-12</td><td align="left" rowspan="1" colspan="1">kinase activity</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0004190</td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">32</td><td align="left" rowspan="1" colspan="1">35</td><td align="left" rowspan="1" colspan="1">3.79E-12</td><td align="left" rowspan="1" colspan="1">aspartic-type endopeptidase activity</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0004197</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">44</td><td align="left" rowspan="1" colspan="1">16</td><td align="left" rowspan="1" colspan="1">65</td><td align="left" rowspan="1" colspan="1">95</td><td align="left" rowspan="1" colspan="1">4.26E-12</td><td align="left" rowspan="1" colspan="1">cysteine-type endopeptidase activity</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0008234</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">46</td><td align="left" rowspan="1" colspan="1">18</td><td align="left" rowspan="1" colspan="1">69</td><td align="left" rowspan="1" colspan="1">105</td><td align="left" rowspan="1" colspan="1">1.02E-11</td><td align="left" rowspan="1" colspan="1">cysteine-type peptidase activity</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0004175</td><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1">59</td><td align="left" rowspan="1" colspan="1">33</td><td align="left" rowspan="1" colspan="1">107</td><td align="left" rowspan="1" colspan="1">190</td><td align="left" rowspan="1" colspan="1">2.75E-11</td><td align="left" rowspan="1" colspan="1">endopeptidase activity</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0015035</td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1">29</td><td align="left" rowspan="1" colspan="1">41</td><td align="left" rowspan="1" colspan="1">9.97E-07</td><td align="left" rowspan="1" colspan="1">protein disulfide oxidoreductase activity</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0005544</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">1.84E-06</td><td align="left" rowspan="1" colspan="1">calcium-dependent phospholipid binding</td></tr><tr><td align="left" rowspan="1" colspan="1">GO:0016620</td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1">2.10E-06</td><td align="left" rowspan="1" colspan="1">oxidoreductase activity, acting on the aldehyde or oxo group of donors, NAD or NADP as acceptor</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pntd-0000269-t003" position="float"><object-id pub-id-type="doi">10.1371/journal.pntd.0000269.t003</object-id><label>Table 3</label><caption><title>Selected Gene Ontology Annotations in IntFam-241.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">GO Annotation</td><td align="left" rowspan="1" colspan="1">GO id</td><td align="left" rowspan="1" colspan="1">Genes</td></tr></thead><tbody><tr><td colspan="3" align="left" rowspan="1">
<bold>Protein Kinases</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">protein kinase activity</td><td align="left" rowspan="1" colspan="1">GO:0004672</td><td align="left" rowspan="1" colspan="1">107</td></tr><tr><td colspan="3" align="left" rowspan="1">
<bold>Serine/threonine-specific protein kinases</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1"> protein serine/threonine kinase activity</td><td align="left" rowspan="1" colspan="1">GO:0004674</td><td align="left" rowspan="1" colspan="1">84</td></tr><tr><td align="left" rowspan="1" colspan="1">  receptor signaling protein serine/threonine kinase activity</td><td align="left" rowspan="1" colspan="1">GO:0004702</td><td align="left" rowspan="1" colspan="1">4</td></tr><tr><td align="left" rowspan="1" colspan="1">   MAP kinase activity</td><td align="left" rowspan="1" colspan="1">GO:0004707</td><td align="left" rowspan="1" colspan="1">4</td></tr><tr><td colspan="3" align="left" rowspan="1">
<bold>Tyrosine-specific protein kinases</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1"> protein-tyrosine kinase activity</td><td align="left" rowspan="1" colspan="1">GO:0004713</td><td align="left" rowspan="1" colspan="1">39</td></tr><tr><td align="left" rowspan="1" colspan="1">  transmembrane receptor protein tyrosine kinase activity</td><td align="left" rowspan="1" colspan="1">GO:0004714</td><td align="left" rowspan="1" colspan="1">2</td></tr><tr><td align="left" rowspan="1" colspan="1">   epidermal growth factor receptor activity</td><td align="left" rowspan="1" colspan="1">GO:0005006</td><td align="left" rowspan="1" colspan="1">1</td></tr><tr><td align="left" rowspan="1" colspan="1">   vascular endothelial growth factor receptor activity</td><td align="left" rowspan="1" colspan="1">GO:0005021</td><td align="left" rowspan="1" colspan="1">1</td></tr><tr><td colspan="3" align="left" rowspan="1">
<bold>Protein Phosphatases</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">phosphoprotein phosphatase activity</td><td align="left" rowspan="1" colspan="1">GO:0004721</td><td align="left" rowspan="1" colspan="1">10</td></tr><tr><td align="left" rowspan="1" colspan="1"> protein tyrosine phosphatase activity</td><td align="left" rowspan="1" colspan="1">GO:0004725</td><td align="left" rowspan="1" colspan="1">10</td></tr><tr><td colspan="3" align="left" rowspan="1">
<bold>Proteases</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">peptidase activity</td><td align="left" rowspan="1" colspan="1">GO:0008233</td><td align="left" rowspan="1" colspan="1">137</td></tr><tr><td colspan="3" align="left" rowspan="1">
<bold>Serine proteases</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1"> serine-type peptidase activity</td><td align="left" rowspan="1" colspan="1">GO:0008236</td><td align="left" rowspan="1" colspan="1">9</td></tr><tr><td colspan="3" align="left" rowspan="1">
<bold>Threonine proteases</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1"> threonine endopeptidase activity</td><td align="left" rowspan="1" colspan="1">GO:0004298</td><td align="left" rowspan="1" colspan="1">8</td></tr><tr><td colspan="3" align="left" rowspan="1">
<bold>Cystine proteases</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1"> cysteine-type peptidase activity</td><td align="left" rowspan="1" colspan="1">GO:0008234</td><td align="left" rowspan="1" colspan="1">69</td></tr><tr><td align="left" rowspan="1" colspan="1"> cysteine-type endopeptidase activity</td><td align="left" rowspan="1" colspan="1">GO:0004197</td><td align="left" rowspan="1" colspan="1">65</td></tr><tr><td align="left" rowspan="1" colspan="1">  legumain activity</td><td align="left" rowspan="1" colspan="1">GO:0001509</td><td align="left" rowspan="1" colspan="1">3</td></tr><tr><td colspan="3" align="left" rowspan="1">
<bold>Aspartic acid proteases</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1"> aspartic-type endopeptidase activity</td><td align="left" rowspan="1" colspan="1">GO:0004190</td><td align="left" rowspan="1" colspan="1">32</td></tr><tr><td align="left" rowspan="1" colspan="1">  pepsin A activity</td><td align="left" rowspan="1" colspan="1">GO:0004194</td><td align="left" rowspan="1" colspan="1">31</td></tr><tr><td colspan="3" align="left" rowspan="1">
<bold>Metalloproteases</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1"> metallopeptidase activity</td><td align="left" rowspan="1" colspan="1">GO:0008237</td><td align="left" rowspan="1" colspan="1">16</td></tr><tr><td align="left" rowspan="1" colspan="1"> metalloendopeptidase activity</td><td align="left" rowspan="1" colspan="1">GO:0004222</td><td align="left" rowspan="1" colspan="1">1</td></tr><tr><td align="left" rowspan="1" colspan="1"> metalloexopeptidase activity</td><td align="left" rowspan="1" colspan="1">GO:0008235</td><td align="left" rowspan="1" colspan="1">11</td></tr><tr><td align="left" rowspan="1" colspan="1">  leucyl aminopeptidase activity</td><td align="left" rowspan="1" colspan="1">GO:0004178</td><td align="left" rowspan="1" colspan="1">1</td></tr><tr><td align="left" rowspan="1" colspan="1">  membrane alanyl aminopeptidase activity</td><td align="left" rowspan="1" colspan="1">GO:0004179</td><td align="left" rowspan="1" colspan="1">7</td></tr><tr><td align="left" rowspan="1" colspan="1">  metallocarboxypeptidase activity</td><td align="left" rowspan="1" colspan="1">GO:0004181</td><td align="left" rowspan="1" colspan="1">3</td></tr><tr><td colspan="3" align="left" rowspan="1">
<bold>Protease Inhibitors</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">protease inhibitor activity</td><td align="left" rowspan="1" colspan="1">GO:0030414</td><td align="left" rowspan="1" colspan="1">17</td></tr><tr><td align="left" rowspan="1" colspan="1"> serine-type endopeptidase inhibitor activity</td><td align="left" rowspan="1" colspan="1">GO:0004867</td><td align="left" rowspan="1" colspan="1">10</td></tr><tr><td align="left" rowspan="1" colspan="1"> cysteine protease inhibitor activity</td><td align="left" rowspan="1" colspan="1">GO:0004869</td><td align="left" rowspan="1" colspan="1">7</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pntd-0000269-t004" position="float"><object-id pub-id-type="doi">10.1371/journal.pntd.0000269.t004</object-id><label>Table 4</label><caption><title>KEGG Pathway Mappings for Intestinal Genes from Three Nematodes.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">KEGG Pathway</td><td colspan="3" align="left" rowspan="1">
<italic>A. suum</italic>
</td><td colspan="3" align="left" rowspan="1">
<italic>H. contortus</italic>
</td><td colspan="3" align="left" rowspan="1">
<italic>C. elegans</italic>
</td><td colspan="2" align="left" rowspan="1">Wormbase v150</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Genes</td><td align="left" rowspan="1" colspan="1">KOs <xref ref-type="table-fn" rid="nt104">1</xref>
</td><td align="left" rowspan="1" colspan="1">
<italic>P-</italic>value <xref ref-type="table-fn" rid="nt105">2</xref>
</td><td align="left" rowspan="1" colspan="1">Genes</td><td align="left" rowspan="1" colspan="1">KOs</td><td align="left" rowspan="1" colspan="1">
<italic>P-</italic>value</td><td align="left" rowspan="1" colspan="1">Genes</td><td align="left" rowspan="1" colspan="1">KOs</td><td align="left" rowspan="1" colspan="1">
<italic>P-</italic>value</td><td align="left" rowspan="1" colspan="1">Genes</td><td align="left" rowspan="1" colspan="1">KOs</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">
<bold>1. Metabolism</bold>
</td><td align="left" rowspan="1" colspan="1">247</td><td align="left" rowspan="1" colspan="1">176</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">167</td><td align="left" rowspan="1" colspan="1">152</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">414</td><td align="left" rowspan="1" colspan="1">316</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">788</td><td align="left" rowspan="1" colspan="1">480</td></tr><tr><td align="left" rowspan="1" colspan="1">1.1 Carbohydrate Metabolism</td><td align="left" rowspan="1" colspan="1">79</td><td align="left" rowspan="1" colspan="1">57</td><td align="left" rowspan="1" colspan="1">6.9E-03</td><td align="left" rowspan="1" colspan="1">46</td><td align="left" rowspan="1" colspan="1">45</td><td align="left" rowspan="1" colspan="1">4.3E-02</td><td align="left" rowspan="1" colspan="1">123</td><td align="left" rowspan="1" colspan="1">88</td><td align="left" rowspan="1" colspan="1">4.7E-03</td><td align="left" rowspan="1" colspan="1">230</td><td align="left" rowspan="1" colspan="1">122</td></tr><tr><td align="left" rowspan="1" colspan="1">1.2 Energy Metabolism</td><td align="left" rowspan="1" colspan="1">105</td><td align="left" rowspan="1" colspan="1">70</td><td align="left" rowspan="1" colspan="1">4.8E-05</td><td align="left" rowspan="1" colspan="1">80</td><td align="left" rowspan="1" colspan="1">70</td><td align="left" rowspan="1" colspan="1">6.7E-09</td><td align="left" rowspan="1" colspan="1">117</td><td align="left" rowspan="1" colspan="1">96</td><td align="left" rowspan="1" colspan="1">6.4E-03</td><td align="left" rowspan="1" colspan="1">182</td><td align="left" rowspan="1" colspan="1">135</td></tr><tr><td align="left" rowspan="1" colspan="1">1.3 Lipid Metabolism</td><td align="left" rowspan="1" colspan="1">25</td><td align="left" rowspan="1" colspan="1">22</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">22</td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">71</td><td align="left" rowspan="1" colspan="1">41</td><td align="left" rowspan="1" colspan="1">2.3E-02</td><td align="left" rowspan="1" colspan="1">133</td><td align="left" rowspan="1" colspan="1">55</td></tr><tr><td align="left" rowspan="1" colspan="1">1.4 Nucleotide Metabolism</td><td align="left" rowspan="1" colspan="1">25</td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">37</td><td align="left" rowspan="1" colspan="1">36</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">92</td><td align="left" rowspan="1" colspan="1">57</td></tr><tr><td align="left" rowspan="1" colspan="1">1.5 Amino Acid Metabolism</td><td align="left" rowspan="1" colspan="1">42</td><td align="left" rowspan="1" colspan="1">35</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">35</td><td align="left" rowspan="1" colspan="1">29</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">106</td><td align="left" rowspan="1" colspan="1">84</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">209</td><td align="left" rowspan="1" colspan="1">129</td></tr><tr><td align="left" rowspan="1" colspan="1">1.6 Metabolism of Other Amino Acids</td><td align="left" rowspan="1" colspan="1">23</td><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">51</td><td align="left" rowspan="1" colspan="1">30</td><td align="left" rowspan="1" colspan="1">2.6E-02</td><td align="left" rowspan="1" colspan="1">91</td><td align="left" rowspan="1" colspan="1">39</td></tr><tr><td align="left" rowspan="1" colspan="1">1.7 Glycan Biosynthesis and Metabolism</td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">30</td><td align="left" rowspan="1" colspan="1">21</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">72</td><td align="left" rowspan="1" colspan="1">50</td></tr><tr><td align="left" rowspan="1" colspan="1">1.8 Biosynthesis of Polyketides and Nonribosomal Peptides</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">2</td></tr><tr><td align="left" rowspan="1" colspan="1">1.9 Metabolism of Cofactors and Vitamins</td><td align="left" rowspan="1" colspan="1">19</td><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">17</td><td align="left" rowspan="1" colspan="1">16</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">47</td><td align="left" rowspan="1" colspan="1">31</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">98</td><td align="left" rowspan="1" colspan="1">47</td></tr><tr><td align="left" rowspan="1" colspan="1">1.10 Biosynthesis of Secondary Metabolites</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">24</td><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">58</td><td align="left" rowspan="1" colspan="1">21</td></tr><tr><td align="left" rowspan="1" colspan="1">1.11 Xenobiotics Biodegradation and Metabolism</td><td align="left" rowspan="1" colspan="1">21</td><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">58</td><td align="left" rowspan="1" colspan="1">26</td><td align="left" rowspan="1" colspan="1">6.1E-04</td><td align="left" rowspan="1" colspan="1">120</td><td align="left" rowspan="1" colspan="1">29</td></tr><tr><td align="left" rowspan="1" colspan="1">1.12 Enzyme Families</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>2. Genetic Information Processing</bold>
</td><td align="left" rowspan="1" colspan="1">280</td><td align="left" rowspan="1" colspan="1">133</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">135</td><td align="left" rowspan="1" colspan="1">120</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">248</td><td align="left" rowspan="1" colspan="1">215</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">432</td><td align="left" rowspan="1" colspan="1">336</td></tr><tr><td align="left" rowspan="1" colspan="1">2.1 Transcription</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">18</td><td align="left" rowspan="1" colspan="1">16</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">22</td><td align="left" rowspan="1" colspan="1">21</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">55</td><td align="left" rowspan="1" colspan="1">50</td></tr><tr><td align="left" rowspan="1" colspan="1">2.2 Translation</td><td align="left" rowspan="1" colspan="1">218</td><td align="left" rowspan="1" colspan="1">91</td><td align="left" rowspan="1" colspan="1">5.3E-10</td><td align="left" rowspan="1" colspan="1">79</td><td align="left" rowspan="1" colspan="1">73</td><td align="left" rowspan="1" colspan="1">6.0E-07</td><td align="left" rowspan="1" colspan="1">130</td><td align="left" rowspan="1" colspan="1">115</td><td align="left" rowspan="1" colspan="1">1.6E-04</td><td align="left" rowspan="1" colspan="1">194</td><td align="left" rowspan="1" colspan="1">155</td></tr><tr><td align="left" rowspan="1" colspan="1">2.3 Folding, Sorting and Degradation</td><td align="left" rowspan="1" colspan="1">50</td><td align="left" rowspan="1" colspan="1">31</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">35</td><td align="left" rowspan="1" colspan="1">28</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">88</td><td align="left" rowspan="1" colspan="1">71</td><td align="left" rowspan="1" colspan="1">6.9E-02</td><td align="left" rowspan="1" colspan="1">149</td><td align="left" rowspan="1" colspan="1">104</td></tr><tr><td align="left" rowspan="1" colspan="1">2.4 Replication and Repair</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">34</td><td align="left" rowspan="1" colspan="1">27</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>3. Environmental Information Processing</bold>
</td><td align="left" rowspan="1" colspan="1">76</td><td align="left" rowspan="1" colspan="1">55</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">42</td><td align="left" rowspan="1" colspan="1">36</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">170</td><td align="left" rowspan="1" colspan="1">99</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">940</td><td align="left" rowspan="1" colspan="1">207</td></tr><tr><td align="left" rowspan="1" colspan="1">3.1 Membrane Transport</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">16</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">47</td><td align="left" rowspan="1" colspan="1">18</td></tr><tr><td align="left" rowspan="1" colspan="1">3.2 Signal Transduction</td><td align="left" rowspan="1" colspan="1">44</td><td align="left" rowspan="1" colspan="1">27</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">23</td><td align="left" rowspan="1" colspan="1">22</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">79</td><td align="left" rowspan="1" colspan="1">61</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">168</td><td align="left" rowspan="1" colspan="1">103</td></tr><tr><td align="left" rowspan="1" colspan="1">3.3 Signaling Molecules and Interaction</td><td align="left" rowspan="1" colspan="1">32</td><td align="left" rowspan="1" colspan="1">28</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">19</td><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">86</td><td align="left" rowspan="1" colspan="1">37</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">752</td><td align="left" rowspan="1" colspan="1">107</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>4. Cellular Processes</bold>
</td><td align="left" rowspan="1" colspan="1">75</td><td align="left" rowspan="1" colspan="1">53</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">86</td><td align="left" rowspan="1" colspan="1">39</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">105</td><td align="left" rowspan="1" colspan="1">85</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">254</td><td align="left" rowspan="1" colspan="1">153</td></tr><tr><td align="left" rowspan="1" colspan="1">4.1 Cell Motility</td><td align="left" rowspan="1" colspan="1">16</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">17</td><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">39</td><td align="left" rowspan="1" colspan="1">28</td></tr><tr><td align="left" rowspan="1" colspan="1">4.2 Cell Growth and Death</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">24</td><td align="left" rowspan="1" colspan="1">22</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">51</td><td align="left" rowspan="1" colspan="1">42</td></tr><tr><td align="left" rowspan="1" colspan="1">4.3 Cell Communication</td><td align="left" rowspan="1" colspan="1">24</td><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">16</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">37</td><td align="left" rowspan="1" colspan="1">31</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">114</td><td align="left" rowspan="1" colspan="1">57</td></tr><tr><td align="left" rowspan="1" colspan="1">4.4 Endocrine System</td><td align="left" rowspan="1" colspan="1">18</td><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">33</td><td align="left" rowspan="1" colspan="1">27</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">68</td><td align="left" rowspan="1" colspan="1">42</td></tr><tr><td align="left" rowspan="1" colspan="1">4.5 Immune System</td><td align="left" rowspan="1" colspan="1">32</td><td align="left" rowspan="1" colspan="1">18</td><td align="left" rowspan="1" colspan="1">3.1E-02</td><td align="left" rowspan="1" colspan="1">58</td><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1">9.9E-02</td><td align="left" rowspan="1" colspan="1">21</td><td align="left" rowspan="1" colspan="1">16</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">55</td><td align="left" rowspan="1" colspan="1">34</td></tr><tr><td align="left" rowspan="1" colspan="1">4.6 Nervous System</td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">56</td><td align="left" rowspan="1" colspan="1">18</td></tr><tr><td align="left" rowspan="1" colspan="1">4.7 Sensory System</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">5</td></tr><tr><td align="left" rowspan="1" colspan="1">4.8 Development</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">22</td><td align="left" rowspan="1" colspan="1">20</td></tr><tr><td align="left" rowspan="1" colspan="1">4.9 Behavior</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">2</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>5. Human Diseases</bold>
</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">17</td><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">42</td><td align="left" rowspan="1" colspan="1">30</td></tr><tr><td align="left" rowspan="1" colspan="1">5.1 Neurodegenerative Disorders</td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">32</td><td align="left" rowspan="1" colspan="1">21</td></tr><tr><td align="left" rowspan="1" colspan="1">5.2 Metabolic Disorders</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1">11</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>7. Unclassified</bold>
</td><td align="left" rowspan="1" colspan="1">44</td><td align="left" rowspan="1" colspan="1">26</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">19</td><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">72</td><td align="left" rowspan="1" colspan="1">49</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">221</td><td align="left" rowspan="1" colspan="1">74</td></tr><tr><td align="left" rowspan="1" colspan="1">7.1 Metabolism</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">17</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">49</td><td align="left" rowspan="1" colspan="1">35</td><td align="left" rowspan="1" colspan="1">3.6E-02</td><td align="left" rowspan="1" colspan="1">145</td><td align="left" rowspan="1" colspan="1">47</td></tr><tr><td align="left" rowspan="1" colspan="1">7.2 Genetic Information Processing</td><td align="left" rowspan="1" colspan="1">18</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">1</td></tr><tr><td align="left" rowspan="1" colspan="1">7.3 Cellular Processes and Signaling</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">23</td><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">75</td><td align="left" rowspan="1" colspan="1">26</td></tr><tr><td align="left" rowspan="1" colspan="1">7.4 Poorly Characterized</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>Total Mapped Genes</bold>
</td><td align="left" rowspan="1" colspan="1">660</td><td align="left" rowspan="1" colspan="1">406</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">419</td><td align="left" rowspan="1" colspan="1">334</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">905</td><td align="left" rowspan="1" colspan="1">687</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">2393</td><td align="left" rowspan="1" colspan="1">1125</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pntd-0000269-t005" position="float"><object-id pub-id-type="doi">10.1371/journal.pntd.0000269.t005</object-id><label>Table 5</label><caption><title>RNA interference (RNAi) Phenotypes Assigned to IntFams through <italic>C. elegans</italic> Members.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">IntFam Group</td><td align="left" rowspan="1" colspan="1">Families</td><td align="left" rowspan="1" colspan="1">Families w/RNAi <xref ref-type="table-fn" rid="nt106">1</xref>
</td><td align="left" rowspan="1" colspan="1">Families w/Pheno <xref ref-type="table-fn" rid="nt107">2</xref>
</td><td align="left" rowspan="1" colspan="1">% Families w/Pheno</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">IntFam-241</td><td align="left" rowspan="1" colspan="1">241</td><td align="left" rowspan="1" colspan="1">241</td><td align="left" rowspan="1" colspan="1">176</td><td align="left" rowspan="1" colspan="1">73%</td></tr><tr><td align="left" rowspan="1" colspan="1">IntFam_<italic>As</italic>_<italic>Ce</italic>
</td><td align="left" rowspan="1" colspan="1">332</td><td align="left" rowspan="1" colspan="1">323</td><td align="left" rowspan="1" colspan="1">182</td><td align="left" rowspan="1" colspan="1">55%</td></tr><tr><td align="left" rowspan="1" colspan="1">IntFam_<italic>Hc_Ce</italic>
</td><td align="left" rowspan="1" colspan="1">290</td><td align="left" rowspan="1" colspan="1">286</td><td align="left" rowspan="1" colspan="1">168</td><td align="left" rowspan="1" colspan="1">58%</td></tr><tr><td align="left" rowspan="1" colspan="1">IntFam_<italic>Ce</italic>_only</td><td align="left" rowspan="1" colspan="1">2,592</td><td align="left" rowspan="1" colspan="1">2,488</td><td align="left" rowspan="1" colspan="1">1,037</td><td align="left" rowspan="1" colspan="1">40%</td></tr><tr><td align="left" rowspan="1" colspan="1">All <italic>Ce</italic>-containing IntFams</td><td align="left" rowspan="1" colspan="1">3,455</td><td align="left" rowspan="1" colspan="1">3,338</td><td align="left" rowspan="1" colspan="1">1,563</td><td align="left" rowspan="1" colspan="1">45%</td></tr></tbody></table></alternatives></table-wrap>
<supplementary-material content-type="local-data" id="pntd.0000269.s001"><label>Figure S1</label><caption><p>Distribution of Sequence Similarities Identified in <italic>A. suum</italic> and <italic>H. contortus</italic> EST Clusters. The three phylogenetically specific sequence groups used to identify sequence similarities of the intestinal genes were: i) <italic>Caenorhabditis</italic> spp., amino acid sequences from the complete genomes of <italic>C. elegans, C. briggsae</italic>, and <italic>C. remanei</italic>, ii) Other Nematoda, non-Caenorhabditis nematode nucleic acid sequences excluding those from either <italic>A. suum</italic> or <italic>H. contortus</italic>, when sequences from <italic>A. suum</italic> or <italic>H. contortus</italic> were queried, respectively, and iii) Non-Nematoda, non-nematode amino acid sequences from the non-redundant protein database NR. In total, 53% (5,303/9,947) <italic>A. suum</italic> and 75% (3,792/5,058) <italic>H. contortus</italic> EST clusters contained primary sequence similarities to known genes from other species, but similar distributions of the identified matches to various species groups were observed in the two parasites.</p><p>(0.05 MB PPT)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pntd.0000269.s002"><label>Figure S2</label><caption><p>Homologous Pairs between the Intestine and Gonad Gene Groups from <italic>A. suum</italic> and <italic>C. elegans</italic>. Significant larger number of genes in the <italic>A. suum</italic> intestine group had homologous counterparts in the <italic>C. elegans</italic> intestine group than in the <italic>C. elegans</italic> gonad group at BLAST bit-score cutoff of either 50 or 100, indicating the intestinal expression of homologous genes tend to be maintained across nematodes. However, the number of homologous pairs detected between the two gonad groups was not different from that between the gonad and intestine groups.</p><p>(0.14 MB PPT)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pntd.0000269.s003"><label>Table S1</label><caption><p>EST Generation and Clustering.</p><p>(0.02 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pntd.0000269.s004"><label>Table S2</label><caption><p>Identification of 150 Intestine-enriched Genes in <italic>A. suum</italic>.</p><p>(0.06 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pntd.0000269.s005"><label>Table S3</label><caption><p>Identification of 60 Intestine-enriched Genes in <italic>H. contortus.</italic>
</p><p>(0.04 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pntd.0000269.s006"><label>Table S4</label><caption><p>Identification of 247 Intestine-enriched Genes in <italic>C. elegans.</italic>.</p><p>(0.07 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pntd.0000269.s007"><label>Table S5</label><caption><p>Classification of Putative Protein Kinases in <italic>A. suum</italic> and <italic>H. contortus.</italic>
</p><p>(0.03 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pntd.0000269.s008"><label>Table S6</label><caption><p>Gene Ontology Terms Statistically Enriched in IntFam Groups Other Than IntFam-241.</p><p>(0.03 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pntd.0000269.s009"><label>Table S7</label><caption><p>Genes Mapped to KEGG Immune System in <italic>A. suum.</italic>
</p><p>(0.03 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pntd.0000269.s010"><label>Table S8</label><caption><p>Genes Mapped to KEGG Immune System in <italic>H. contortus.</italic>
</p><p>(0.03 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pntd.0000269.s011"><label>Text S1</label><caption><p>Supplementary Materials.</p><p>(0.04 MB DOC)</p></caption></supplementary-material>
<table-wrap-foot><fn id="nt101"><label>1</label><p>SAGE tags from <italic>glp-4</italic> adults (6,903 genes) and dissected gut (4,071 genes) for <italic>C. elegans.</italic>
</p></fn><fn id="nt102"><label>2</label><p>Included in parentheses are putative false positives.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt103"><label>1</label><p>
<italic>As</italic>: <italic>A. suum</italic>; <italic>Hc</italic>: <italic>H. contortus</italic>; <italic>Ce</italic>: <italic>C. elegans</italic>; Total: IntFam-241 genes from all the three species.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt104"><label>1</label><p>Each KEGG Orthology (KO) represent an ortholog/paralog group defined by KEGG.</p></fn><fn id="nt105"><label>2</label><p>Only <italic>P</italic>-values meeting the cutoff of 0.1 are shown here.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt106"><label>1</label><p>Families whose <italic>C. elegans</italic> members have RNAi results available.</p></fn><fn id="nt107"><label>2</label><p>Families whose <italic>C. elegans</italic> members have RNAi phenotypes reported.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p>
<italic>Ascaris suum</italic> and <italic>Haemonchus contortus</italic> EST sequencing at Washington University was supported by NIH-NIAID research grant AI 46593 to RKW. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</p></fn></fn-group>
<graphic id="pntd-0000269-t001-1" xlink:href="pntd.0000269.t001"/> <graphic xlink:href="pntd.0000269.g001"/> <graphic xlink:href="pntd.0000269.g002"/> <graphic xlink:href="pntd.0000269.g003"/> <graphic xlink:href="pntd.0000269.g004"/> <graphic id="pntd-0000269-t002-2" xlink:href="pntd.0000269.t002"/> <graphic id="pntd-0000269-t003-3" xlink:href="pntd.0000269.t003"/> <graphic id="pntd-0000269-t004-4" xlink:href="pntd.0000269.t004"/> <graphic id="pntd-0000269-t005-5" xlink:href="pntd.0000269.t005"/>
<media xlink:href="pntd.0000269.s001.ppt"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pntd.0000269.s002.ppt"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pntd.0000269.s003.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pntd.0000269.s004.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pntd.0000269.s005.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pntd.0000269.s006.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pntd.0000269.s007.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pntd.0000269.s008.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pntd.0000269.s009.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pntd.0000269.s010.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pntd.0000269.s011.doc"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["1"], "element-citation": ["\n"], "surname": ["McGhee"], "given-names": ["JD"], "year": ["2007"], "article-title": ["The "], "italic": ["C. elegans", "WormBook", "C. elegans"], "comment": ["The "], "ext-link": ["http://www.wormbook.org"]}, {"label": ["2"], "element-citation": ["\n"], "surname": ["Ghafouri", "McGhee"], "given-names": ["S", "JD"], "year": ["2007"], "article-title": ["Bacterial residence time in the intestine of "], "italic": ["Caenorhabditis elegans"], "source": ["Nematology"], "volume": ["9"], "fpage": ["87"], "lpage": ["91"]}, {"label": ["17"], "element-citation": ["\n"], "surname": ["McKay", "Johnson", "Khattra", "Asano", "Baillie"], "given-names": ["SJ", "R", "J", "J", "DL"], "year": ["2004"], "article-title": ["Gene expression profiling of cells, tissues and developmental stages of the nematode "], "italic": ["C. elegans"], "source": ["Cold Spring Harb Symp Quant Biol"], "volume": ["68"], "fpage": ["159"], "lpage": ["170"]}, {"label": ["23"], "element-citation": ["\n"], "surname": ["McCarter", "Dautova Mitreva", "Martin", "Dante", "Wylie"], "given-names": ["J", "M", "J", "M", "T"], "year": ["2003"], "article-title": ["Analysis and functional classification of transcripts from the nematode "], "italic": ["Meloidogyne incognita"], "source": ["Genome Biology"], "volume": ["4: R26"], "fpage": ["21"], "lpage": ["19"]}, {"label": ["33"], "element-citation": ["\n"], "surname": ["Mulder", "Apweiler", "Attwood", "Bairoch", "Bateman"], "given-names": ["NJ", "R", "TK", "A", "A"], "year": ["2005"], "article-title": ["InterPro, progress and status in 2005."], "source": ["Nucleic Acids Res 33 Database Issue"], "fpage": ["D201"], "lpage": ["205"]}, {"label": ["43"], "element-citation": ["\n"], "surname": ["Steel", "Torrie"], "given-names": ["RGD", "JH"], "year": ["1960"], "article-title": ["Principles and procedures of statistics: McGraw-Hill Book Company, Inc. New York Toronto London."]}]
PLoS Negl Trop Dis. 2008 Aug 6; 2(8):e269