Pathogenomic Inference of Virulence-Associated Genes in Leptospira interrogans

Leptospirosis is a globally important, neglected zoonotic infection caused by spirochetes of the genus Leptospira. Since genetic transformation remains technically limited for pathogenic Leptospira, a systems biology pathogenomic approach was used to infer leptospiral virulence genes by whole genome comparison of culture-attenuated Leptospira interrogans serovar Lai with its virulent, isogenic parent. Among the 11 pathogen-specific protein-coding genes in which non-synonymous mutations were found, a putative soluble adenylate cyclase with host cell cAMP-elevating activity, and two members of a previously unstudied ∼15 member paralogous gene family of unknown function were identified. This gene family was also uniquely found in the alpha-proteobacteria Bartonella bacilliformis and Bartonella australis that are geographically restricted to the Andes and Australia, respectively. How the pathogenic Leptospira and these two Bartonella species came to share this expanded gene family remains an evolutionary mystery. In vivo expression analyses demonstrated up-regulation of 10/11 Leptospira genes identified in the attenuation screen, and profound in vivo, tissue-specific up-regulation by members of the paralogous gene family, suggesting a direct role in virulence and host-pathogen interactions. The pathogenomic experimental design here is generalizable as a functional systems biology approach to studying bacterial pathogenesis and virulence and should encourage similar experimental studies of other pathogens.     

What's in the genome of a filamentous fungus? Analysis of the Neurospora genome sequence

The German Neurospora Genome Project has assembled sequences from ordered cosmid and BAC clones of linkage groups II and V of the genome of Neurospora crassa in 13 and 12 contigs, respectively. Including additional sequences located on other linkage groups a total of 12 Mb were subjected to a manual gene extraction and annotation process. The genome comprises a small number of repetitive elements, a low degree of segmental duplications and very few paralogous genes. The analysis of the 3218 identified open reading frames provides a first overview of the protein equipment of a filamentous fungus. Significantly, N.crassa possesses a large variety of metabolic enzymes including a substantial number of enzymes involved in the degradation of complex substrates as well as secondary metabolism. While several of these enzymes are specific for filamentous fungi many are shared exclusively with prokaryotes.     

Trichomonas Transmembrane Cyclases Result from Massive Gene Duplication and Concomitant Development of Pseudogenes

Background

Trichomonas vaginalis has an unusually large genome (∼160 Mb) encoding ∼60,000 proteins. With the goal of beginning to understand why some Trichomonas genes are present in so many copies, we characterized here a family of ∼123 Trichomonas genes that encode transmembrane adenylyl cyclases (TMACs).

Methodology/Principal Findings

The large family of TMACs genes is the result of recent duplications of a small set of ancestral genes that appear to be unique to trichomonads. Duplicated TMAC genes are not closely associated with repetitive elements, and duplications of flanking sequences are rare. However, there is evidence for TMAC gene replacements by homologous recombination. A high percentage of TMAC genes (∼46%) are pseudogenes, as they contain stop codons and/or frame shifts, or the genes are truncated. Numerous stop codons present in the genome project G3 strain are not present in orthologous genes of two other Trichomonas strains (S1 and B7RC2). Each TMAC is composed of a series of N-terminal transmembrane helices and a single C-terminal cyclase domain that has adenylyl cyclase activity. Multiple TMAC genes are transcribed by Trichomonas cloned by limiting dilution.

Conclusions/Significance

We conclude that one reason for the unusually large genome of Trichomonas is the presence of unstable families of genes such as those encoding TMACs that are undergoing massive gene duplication and concomitant development of pseudogenes.     

Natural variation of the Y chromosome suppresses sex ratio distortion and modulates testis-specific gene expression in Drosophila simulans

X-linked sex-ratio distorters that disrupt spermatogenesis can cause a deficiency in functional Y-bearing sperm and a female-biased sex ratio. Y-linked modifiers that restore a normal sex ratio might be abundant and favored when a X-linked distorter is present. Here we investigated natural variation of Y-linked suppressors of sex-ratio in the Winters systems and the ability of these chromosomes to modulate gene expression in Drosophila simulans. Seventy-eight Y chromosomes of worldwide origin were assayed for their resistance to the X-linked sex-ratio distorter gene Dox. Y chromosome diversity caused males to sire ∼63% to ∼98% female progeny. Genome-wide gene expression analysis revealed hundreds of genes differentially expressed between isogenic males with sensitive (high sex ratio) and resistant (low sex ratio) Y chromosomes from the same population. Although the expression of about 75% of all testis-specific genes remained unchanged across Y chromosomes, a subset of post-meiotic genes was upregulated by resistant Y chromosomes. Conversely, a set of accessory gland-specific genes and mitochondrial genes were downregulated in males with resistant Y chromosomes. The D. simulans Y chromosome also modulated gene expression in XXY females in which the Y-linked protein-coding genes are not transcribed. The data suggest that the Y chromosome might exert its regulatory functions through epigenetic mechanisms that do not require the expression of protein-coding genes. The gene network that modulates sex ratio distortion by the Y chromosome is poorly understood, other than that it might include interactions with mitochondria and enriched for genes expressed in post-meiotic stages of spermatogenesis.     

The Solanum commersonii Genome Sequence Provides Insights into Adaptation to Stress Conditions and Genome Evolution of Wild Potato Relatives

Here, we report the draft genome sequence of Solanum commersonii, which consists of ∼830 megabases with an N50 of 44,303 bp anchored to 12 chromosomes, using the potato (Solanum tuberosum) genome sequence as a reference. Compared with potato, S. commersonii shows a striking reduction in heterozygosity (1.5% versus 53 to 59%), and differences in genome sizes were mainly due to variations in intergenic sequence length. Gene annotation by ab initio prediction supported by RNA-seq data produced a catalog of 1703 predicted microRNAs, 18,882 long noncoding RNAs of which 20% are shown to target cold-responsive genes, and 39,290 protein-coding genes with a significant repertoire of nonredundant nucleotide binding site-encoding genes and 126 cold-related genes that are lacking in S. tuberosum. Phylogenetic analyses indicate that domesticated potato and S. commersonii lineages diverged ∼2.3 million years ago. Three duplication periods corresponding to genome enrichment for particular gene families related to response to salt stress, water transport, growth, and defense response were discovered. The draft genome sequence of S. commersonii substantially increases our understanding of the domesticated germplasm, facilitating translation of acquired knowledge into advances in crop stability in light of global climate and environmental changes.     

Detailed Analysis of a Contiguous 22-Mb Region of the Maize Genome

Most of our understanding of plant genome structure and evolution has come from the careful annotation of small (e.g., 100 kb) sequenced genomic regions or from automated annotation of complete genome sequences. Here, we sequenced and carefully annotated a contiguous 22 Mb region of maize chromosome 4 using an improved pseudomolecule for annotation. The sequence segment was comprehensively ordered, oriented, and confirmed using the maize optical map. Nearly 84% of the sequence is composed of transposable elements (TEs) that are mostly nested within each other, of which most families are low-copy. We identified 544 gene models using multiple levels of evidence, as well as five miRNA genes. Gene fragments, many captured by TEs, are prevalent within this region. Elimination of gene redundancy from a tetraploid maize ancestor that originated a few million years ago is responsible in this region for most disruptions of synteny with sorghum and rice. Consistent with other sub-genomic analyses in maize, small RNA mapping showed that many small RNAs match TEs and that most TEs match small RNAs. These results, performed on ∼1% of the maize genome, demonstrate the feasibility of refining the B73 RefGen_v1 genome assembly by incorporating optical map, high-resolution genetic map, and comparative genomic data sets. Such improvements, along with those of gene and repeat annotation, will serve to promote future functional genomic and phylogenomic research in maize and other grasses.     

Diversity of 23S rRNA Genes within Individual Prokaryotic Genomes

Background

The concept of ribosomal constraints on rRNA genes is deduced primarily based on the comparison of consensus rRNA sequences between closely related species, but recent advances in whole-genome sequencing allow evaluation of this concept within organisms with multiple rRNA operons.

Methodology/Principal Findings

Using the 23S rRNA gene as an example, we analyzed the diversity among individual rRNA genes within a genome. Of 184 prokaryotic species containing multiple 23S rRNA genes, diversity was observed in 113 (61.4%) genomes (mean 0.40%, range 0.01%–4.04%). Significant (1.17%–4.04%) intragenomic variation was found in 8 species. In 5 of the 8 species, the diversity in the primary structure had only minimal effect on the secondary structure (stem versus loop transition). In the remaining 3 species, the diversity significantly altered local secondary structure, but the alteration appears minimized through complex rearrangement. Intervening sequences (IVS), ranging between 9 and 1471 nt in size, were found in 7 species. IVS in Deinococcus radiodurans and Nostoc sp. encode transposases. T. tengcongensis was the only species in which intragenomic diversity >3% was observed among 4 paralogous 23S rRNA genes.

Conclusions/Significance

These findings indicate tight ribosomal constraints on individual 23S rRNA genes within a genome. Although classification using primary 23S rRNA sequences could be erroneous, significant diversity among paralogous 23S rRNA genes was observed only once in the 184 species analyzed, indicating little overall impact on the mainstream of 23S rRNA gene-based prokaryotic taxonomy.     

Evidence for Transceptor Function of Cellodextrin Transporters in Neurospora crassa

Neurospora crassa colonizes burnt grasslands and metabolizes both cellulose and hemicellulose from plant cell walls. When switched from a favored carbon source to cellulose, N. crassa dramatically up-regulates expression and secretion of genes encoding lignocellulolytic enzymes. However, the means by which N. crassa and other filamentous fungi sense the presence of cellulose in the environment remains unclear. Previously, we have shown that a N. crassa mutant carrying deletions of three β-glucosidase enzymes (Δ3βG) lacks β-glucosidase activity, but efficiently induces cellulase gene expression and cellulolytic activity in the presence of cellobiose as the sole carbon source. These observations indicate that cellobiose, or a modified version of cellobiose, functions as an inducer of lignocellulolytic gene expression and activity in N. crassa. Here, we show that in N. crassa, two cellodextrin transporters, CDT-1 and CDT-2, contribute to cellulose sensing. A N. crassa mutant carrying deletions for both transporters is unable to induce cellulase gene expression in response to crystalline cellulose. Furthermore, a mutant lacking genes encoding both the β-glucosidase enzymes and cellodextrin transporters (Δ3βGΔ2T) does not induce cellulase gene expression in response to cellobiose. Point mutations that severely reduce cellobiose transport by either CDT-1 or CDT-2 when expressed individually do not greatly impact cellobiose induction of cellulase gene expression. These data suggest that the N. crassa cellodextrin transporters act as “transceptors” with dual functions - cellodextrin transport and receptor signaling that results in downstream activation of cellulolytic gene expression. Similar mechanisms of transceptor activity likely occur in related ascomycetes used for industrial cellulase production.     

Functional and Probiotic Attributes of an Indigenous Isolate of Lactobacillus plantarum

Background

Probiotic microorganisms favorably alter the intestinal microflora balance, promote intestinal integrity and mobility, inhibit the growth of harmful bacteria and increase resistance to infection. Probiotics are increasingly used in nutraceuticals, functional foods or in microbial interference treatment. However, the effectiveness of probiotic organism is considered to be population-specific due to variation in gut microflora, food habits and specific host-microbial interactions. Most of the probiotic strains available in the market are of western or European origin, and a strong need for exploring new indigenous probiotic organisms is felt.

Methods and Findings

An indigenous isolate Lp9 identified as Lactobacillus plantarum by molecular-typing methods was studied extensively for its functional and probiotic attributes, viz., acid and bile salt tolerance, cell surface hydrophobicity, autoaggregation and Caco-2 cell-binding as well as antibacterial and antioxidative activities. Lp9 isolate could survive 2 h incubation at pH 1.5–2.0 and toxicity of 1.5–2.0% oxgall bile. Lp9 could deconjugate major bile salts like glycocholate and deoxytaurocholate, indicating its potential to cause hypocholesterolemia. The isolate exhibited cell-surface hydrophobicity of ∼37% and autoaggregation of ∼31%. Presence of putative probiotic marker genes like mucus-binding protein (mub), fibronectin-binding protein (fbp) and bile salt hydrolase (bsh) were confirmed by PCR. Presence of these genes suggested the possibility of specific interaction and colonization potential of Lp9 isolate in the gut, which was also suggested by a good adhesion ratio of 7.4±1.3% with Caco-2 cell line. The isolate demonstrated higher free radical scavenging activity than standard probiotics L. johnsonii LA1 and L. acidophilus LA7. Lp9 also exhibited antibacterial activity against E. coli, L. monocytogenes, S. typhi, S. aureus and B. cereus.

Conclusion

The indigenous Lactobacillus plantarum Lp9 exhibited high resistance against low pH and bile and possessed antibacterial, antioxidative and cholesterol lowering properties with a potential for exploitation in the development of indigenous functional food or nutraceuticals.     

Complete Genome Sequence of Halalkalicoccus jeotgali B3T,an Extremely Halophilic Archaeon

Halalkalicoccus jeotgali B3^T, isolated from salt-fermented seafood from South Korea, is an extremely halophilic archaeon belonging to the family Halobacteriaceae. Here, we present the complete genome sequence of the type strain H. jeotgali B3^T (3,698,650 bp, with a G+C content of 62.5%), which consists of one chromosome and six plasmids. This is the first complete genome sequence of the Halalkalicoccus species.Extremely halophilic archaea (haloarchaea) are adapted to hypersaline environments and grow optimally in NaCl solutions of 2.6 M or higher (12). These haloarchaea are classified within the family Halobacteriaceae in the order Halobacteriales; currently, this family comprises 28 genera (3), and only 11 complete genome sequences in Halobacteriaceae have been reported. In a study of archaeal diversity in salt-fermented small shrimp or shellfish from South Korea, our laboratory isolated and characterized 5 novel, extremely halophilic archaeal strains of Halobacteriaceae. These strains included Natronococcus jeotgali (9), Halalkalicoccus jeotgali (11), Halorubrum cibi (7), Haloterrigena jeotgali (10) and Haladaptatus cibarius (8). We have now sequenced the genome of Halalkalicoccus jeotgali B3^T; genome sequencing had not been completed or initiated for any strain in this genus when our sequencing project was begun. The genus Halalkalicoccus currently contains only two species, Halalkalicoccus tibetensis (13) and H. jeotgali, and these species exhibit 98.6% gene sequence similarity in their 16S rRNA. The genome of H. jeotgali B3^T is the first of this genus to be sequenced.The complete genome sequence of H. jeotgali B3^T was determined by a whole-genome shotgun strategy using Roche 454 GS (FLX Titanium) pyrosequencing (898,168 reads totaling ∼348 Mb; ∼94-fold coverage of the genome) and a fosmid library (514 reads totaling ∼680 kb) at the Genome Resource Center, KRIBB (Korea Research Institute of Bioscience and Biotechnology). Genome sequences from pyrosequencing were processed by Roche''s software according to the manufacturer''s instructions, and sequences from the fosmid library were processed by PESTAS (6). A total of 898,196 reads were assembled using Newbler Assembler 2.3 (454 Life Science), which generated 54 large contigs (>100 bp in size) with bases having quality scores of 40 and above. The gaps between contigs were closed by primer walking and sequencing of PCR products across the gaps. The annotation was done by merging results obtained from the RAST (Rapid Annotation using Subsystem Technology) pipeline (1), Glimmer 3.02 (2), tRNAscan-SE 1.21 (5), and RNAmmer 1.2 (4).The H. jeotgali B3^T genome is 3,698,650 bases long with a 62.5% G+C content. The chromosome consists of a single circular chromosome (2,809,118 bp, with a G+C content of 65.0%) and six plasmids (406,285 bp, 55.3%; 363,534 bp, 54.2%; 44,576 bp, 58.9%; 44,459 bp, 54.9%; 23,727 bp, 47.6%; 6,951 bp, 60.6%). The genome contains 3,860 predicted coding sequences and 52 RNA genes (determined using RAST). The chromosome is predicted to contain 3,101 coding sequences with a coding intensity of 90.0%, including 47 tRNA genes, 1 5S rRNA gene, 1 16S rRNA gene, and 1 23S rRNA gene. The largest plasmid contains 466 coding sequences with a coding intensity of 81.2% and 2 tRNA genes, while the other five plasmids contain 425, 44, 48, 29, and 5 coding sequences with coding intensities of 80.2%, 84.2%, 83.0%, 69.6%, and 22.8%, respectively (determined using Glimmer3). More detailed analysis of this genome and comparative analysis with other haloarchaea will provide further insight into the genomic differences and metabolism of the extremely halophilic archaea.     

Increased Expression and Protein Divergence in Duplicate Genes Is Associated with Morphological Diversification

The differentiation of both gene expression and protein function is thought to be important as a mechanism of the functionalization of duplicate genes. However, it has not been addressed whether expression or protein divergence of duplicate genes is greater in those genes that have undergone functionalization compared with those that have not. We examined a total of 492 paralogous gene pairs associated with morphological diversification in a plant model organism (Arabidopsis thaliana). Classifying these paralogous gene pairs into high, low, and no morphological diversification groups, based on knock-out data, we found that the divergence rate of both gene expression and protein sequences were significantly higher in either high or low morphological diversification groups compared with those in the no morphological diversification group. These results strongly suggest that the divergence of both expression and protein sequence are important sources for morphological diversification of duplicate genes. Although both mechanisms are not mutually exclusive, our analysis suggested that changes of expression pattern play the minor role (33%–41%) and that changes of protein sequence play the major role (59%–67%) in morphological diversification. Finally, we examined to what extent duplicate genes are associated with expression or protein divergence exerting morphological diversification at the whole-genome level. Interestingly, duplicate genes randomly chosen from A. thaliana had not experienced expression or protein divergence that resulted in morphological diversification. These results indicate that most duplicate genes have experienced minor functionalization.     

New Assembly,Reannotation and Analysis of the Entamoeba histolytica Genome Reveal New Genomic Features and Protein Content Information

Background

In order to maintain genome information accurately and relevantly, original genome annotations need to be updated and evaluated regularly. Manual reannotation of genomes is important as it can significantly reduce the propagation of errors and consequently diminishes the time spent on mistaken research. For this reason, after five years from the initial submission of the Entamoeba histolytica draft genome publication, we have re-examined the original 23 Mb assembly and the annotation of the predicted genes.

Principal Findings

The evaluation of the genomic sequence led to the identification of more than one hundred artifactual tandem duplications that were eliminated by re-assembling the genome. The reannotation was done using a combination of manual and automated genome analysis. The new 20 Mb assembly contains 1,496 scaffolds and 8,201 predicted genes, of which 60% are identical to the initial annotation and the remaining 40% underwent structural changes. Functional classification of 60% of the genes was modified based on recent sequence comparisons and new experimental data. We have assigned putative function to 3,788 proteins (46% of the predicted proteome) based on the annotation of predicted gene families, and have identified 58 protein families of five or more members that share no homology with known proteins and thus could be entamoeba specific. Genome analysis also revealed new features such as the presence of segmental duplications of up to 16 kb flanked by inverted repeats, and the tight association of some gene families with transposable elements.

Significance

This new genome annotation and analysis represents a more refined and accurate blueprint of the pathogen genome, and provides an upgraded tool as reference for the study of many important aspects of E. histolytica biology, such as genome evolution and pathogenesis.     

Metabolic potential of a single cell belonging to one of the most abundant lineages in freshwater bacterioplankton

Actinobacteria within the acI lineage are often numerically dominating in freshwater ecosystems, where they can account for >50% of total bacteria in the surface water. However, they remain uncultured to date. We thus set out to use single-cell genomics to gain insights into their genetic make-up, with the aim of learning about their physiology and ecological niche. A representative from the highly abundant acI-B1 group was selected for shotgun genomic sequencing. We obtained a draft genomic sequence in 75 larger contigs (sum=1.16 Mb), with an unusually low genomic G+C mol% (∼42%). Actinobacteria core gene analysis suggests an almost complete genome recovery. We found that the acI-B1 cell had a small genome, with a rather low percentage of genes having no predicted functions (∼15%) as compared with other cultured and genome-sequenced microbial species. Our metabolic reconstruction hints at a facultative aerobe microorganism with many transporters and enzymes for pentoses utilization (for example, xylose). We also found an actinorhodopsin gene that may contribute to energy conservation under unfavorable conditions. This project reveals the metabolic potential of a member of the global abundant freshwater Actinobacteria.     

Rapid Evolution of Virulence and Drug Resistance in the Emerging Zoonotic Pathogen Streptococcus suis

Background

Streptococcus suis is a zoonotic pathogen that infects pigs and can occasionally cause serious infections in humans. S. suis infections occur sporadically in human Europe and North America, but a recent major outbreak has been described in China with high levels of mortality. The mechanisms of S. suis pathogenesis in humans and pigs are poorly understood.

Methodology/Principal Findings

The sequencing of whole genomes of S. suis isolates provides opportunities to investigate the genetic basis of infection. Here we describe whole genome sequences of three S. suis strains from the same lineage: one from European pigs, and two from human cases from China and Vietnam. Comparative genomic analysis was used to investigate the variability of these strains. S. suis is phylogenetically distinct from other Streptococcus species for which genome sequences are currently available. Accordingly, ∼40% of the ∼2 Mb genome is unique in comparison to other Streptococcus species. Finer genomic comparisons within the species showed a high level of sequence conservation; virtually all of the genome is common to the S. suis strains. The only exceptions are three ∼90 kb regions, present in the two isolates from humans, composed of integrative conjugative elements and transposons. Carried in these regions are coding sequences associated with drug resistance. In addition, small-scale sequence variation has generated pseudogenes in putative virulence and colonization factors.

Conclusions/Significance

The genomic inventories of genetically related S. suis strains, isolated from distinct hosts and diseases, exhibit high levels of conservation. However, the genomes provide evidence that horizontal gene transfer has contributed to the evolution of drug resistance.     

Genome Erosion in a Nitrogen-Fixing Vertically Transmitted Endosymbiotic Multicellular Cyanobacterium

Background

An ancient cyanobacterial incorporation into a eukaryotic organism led to the evolution of plastids (chloroplasts) and subsequently to the origin of the plant kingdom. The underlying mechanism and the identities of the partners in this monophyletic event remain elusive.

Methodology/Principal Findings

To shed light on this evolutionary process, we sequenced the genome of a cyanobacterium residing extracellularly in an endosymbiosis with a plant, the water-fern Azolla filiculoides Lam. This symbiosis was selected as it has characters which make it unique among extant cyanobacterial plant symbioses: the cyanobacterium lacks autonomous growth and is vertically transmitted between plant generations. Our results reveal features of evolutionary significance. The genome is in an eroding state, evidenced by a large proportion of pseudogenes (31.2%) and a high frequency of transposable elements (∼600) scattered throughout the genome. Pseudogenization is found in genes such as the replication initiator dnaA and DNA repair genes, considered essential to free-living cyanobacteria. For some functional categories of genes pseudogenes are more prevalent than functional genes. Loss of function is apparent even within the ‘core’ gene categories of bacteria, such as genes involved in glycolysis and nutrient uptake. In contrast, serving as a critical source of nitrogen for the host, genes related to metabolic processes such as cell differentiation and nitrogen-fixation are well preserved.

Conclusions/Significance

This is the first finding of genome degradation in a plant symbiont and phenotypically complex cyanobacterium and one of only a few extracellular endosymbionts described showing signs of reductive genome evolution. Our findings suggest an ongoing selective streamlining of this cyanobacterial genome which has resulted in an organism devoted to nitrogen fixation and devoid of autonomous growth. The cyanobacterial symbiont of Azolla can thus be considered at the initial phase of a transition from free-living organism to a nitrogen-fixing plant entity, a transition process which may mimic what drove the evolution of chloroplasts from a cyanobacterial ancestor.     

Bioinformatics Analysis Identify Novel OB Fold Protein Coding Genes in C. elegans

Background

The C. elegans genome has been extensively annotated by the WormBase consortium that uses state of the art bioinformatics pipelines, functional genomics and manual curation approaches. As a result, the identification of novel genes in silico in this model organism is becoming more challenging requiring new approaches. The Oligonucleotide-oligosaccharide binding (OB) fold is a highly divergent protein family, in which protein sequences, in spite of having the same fold, share very little sequence identity (5–25%). Therefore, evidence from sequence-based annotation may not be sufficient to identify all the members of this family. In C. elegans, the number of OB-fold proteins reported is remarkably low (n = 46) compared to other evolutionary-related eukaryotes, such as yeast S. cerevisiae (n = 344) or fruit fly D. melanogaster (n = 84). Gene loss during evolution or differences in the level of annotation for this protein family, may explain these discrepancies.

Methodology/Principal Findings

This study examines the possibility that novel OB-fold coding genes exist in the worm. We developed a bioinformatics approach that uses the most sensitive sequence-sequence, sequence-profile and profile-profile similarity search methods followed by 3D-structure prediction as a filtering step to eliminate false positive candidate sequences. We have predicted 18 coding genes containing the OB-fold that have remarkably partially been characterized in C. elegans.

Conclusions/Significance

This study raises the possibility that the annotation of highly divergent protein fold families can be improved in C. elegans. Similar strategies could be implemented for large scale analysis by the WormBase consortium when novel versions of the genome sequence of C. elegans, or other evolutionary related species are being released. This approach is of general interest to the scientific community since it can be used to annotate any genome.