ConsPred: a rule-based (re-)annotation framework for prokaryotic genomes.

MOTIVATION: The rapidly growing number of available prokaryotic genome sequences requires fully automated and high-quality software solutions for their initial and re-annotation. Here we present ConsPred, a prokaryotic genome annotation framework that performs intrinsic gene predictions, homology searches, predictions of non-coding genes as well as CRISPR repeats and integrates all evidence into a consensus annotation. ConsPred achieves comprehensive, high-quality annotations based on rules and priorities, similar to decision-making in manual curation and avoids conflicting predictions. Parameters controlling the annotation process are configurable by the user. ConsPred has been used in the institutions of the authors for longer than 5 years and can easily be extended and adapted to specific needs. SUMMARY: The ConsPred algorithm for producing a consensus from the varying scores of multiple gene prediction programs approaches manual curation in accuracy. Its rule-based approach for choosing final predictions avoids overriding previous manual curations. AVAILABILITY AND IMPLEMENTATION: ConsPred is implemented in Java, Perl and Shell and is freely available under the Creative Commons license as a stand-alone in-house pipeline or as an Amazon Machine Image for cloud computing, see https://sourceforge.net/projects/conspred/. CONTACT: thomas.rattei@univie.ac.atSupplementary information: Supplementary data are available at Bioinformatics online.

Unity in variety--the pan-genome of the Chlamydiae.

Chlamydiae are evolutionarily well-separated bacteria that live exclusively within eukaryotic host cells. They include important human pathogens such as Chlamydia trachomatis as well as symbionts of protozoa. As these bacteria are experimentally challenging and genetically intractable, our knowledge about them is still limited. In this study, we obtained the genome sequences of Simkania negevensis Z, Waddlia chondrophila 2032/99, and Parachlamydia acanthamoebae UV-7. This enabled us to perform the first comprehensive comparative and phylogenomic analysis of representative members of four major families of the Chlamydiae, including the Chlamydiaceae. We identified a surprisingly large core gene set present in all genomes and a high number of diverse accessory genes in those Chlamydiae that do not primarily infect humans or animals, including a chemosensory system in P. acanthamoebae and a type IV secretion system. In S. negevensis, the type IV secretion system is encoded on a large conjugative plasmid (pSn, 132 kb). Phylogenetic analyses suggested that a plasmid similar to the S. negevensis plasmid was originally acquired by the last common ancestor of all four families and that it was subsequently reduced, integrated into the chromosome, or lost during diversification, ultimately giving rise to the extant virulence-associated plasmid of pathogenic chlamydiae. Other virulence factors, including a type III secretion system, are conserved among the Chlamydiae to variable degrees and together with differences in the composition of the cell wall reflect adaptation to different host cells including convergent evolution among the four chlamydial families. Phylogenomic analysis focusing on chlamydial proteins with homology to plant proteins provided evidence for the acquisition of 53 chlamydial genes by a plant progenitor, lending further support for the hypothesis of an early interaction between a chlamydial ancestor and the primary photosynthetic eukaryote.

B2G-FAR, a species-centered GO annotation repository.

MOTIVATION: Functional genomics research has expanded enormously in the last decade thanks to the cost reduction in high-throughput technologies and the development of computational tools that generate, standardize and share information on gene and protein function such as the Gene Ontology (GO). Nevertheless, many biologists, especially working with non-model organisms, still suffer from non-existing or low-coverage functional annotation, or simply struggle retrieving, summarizing and querying these data. RESULTS: The Blast2GO Functional Annotation Repository (B2G-FAR) is a bioinformatics resource envisaged to provide functional information for otherwise uncharacterized sequence data and offers data mining tools to analyze a larger repertoire of species than currently available. This new annotation resource has been created by applying the Blast2GO functional annotation engine in a strongly high-throughput manner to the entire space of public available sequences. The resulting repository contains GO term predictions for over 13.2 million non-redundant protein sequences based on BLAST search alignments from the SIMAP database. We generated GO annotation for approximately 150 000 different taxa making available 2000 species with the highest coverage through B2G-FAR. A second section within B2G-FAR holds functional annotations for 17 non-model organism Affymetrix GeneChips. CONCLUSIONS: B2G-FAR provides easy access to exhaustive functional annotation for 2000 species offering a good balance between quality and quantity, thereby supporting functional genomics research especially in the case of non-model organisms. AVAILABILITY: The annotation resource is available at http://www.b2gfar.org.

SIMAP--a comprehensive database of pre-calculated protein sequence similarities, domains, annotations and clusters.

The prediction of protein function as well as the reconstruction of evolutionary genesis employing sequence comparison at large is still the most powerful tool in sequence analysis. Due to the exponential growth of the number of known protein sequences and the subsequent quadratic growth of the similarity matrix, the computation of the Similarity Matrix of Proteins (SIMAP) becomes a computational intensive task. The SIMAP database provides a comprehensive and up-to-date pre-calculation of the protein sequence similarity matrix, sequence-based features and sequence clusters. As of September 2009, SIMAP covers 48 million proteins and more than 23 million non-redundant sequences. Novel features of SIMAP include the expansion of the sequence space by including databases such as ENSEMBL as well as the integration of metagenomes based on their consistent processing and annotation. Furthermore, protein function predictions by Blast2GO are pre-calculated for all sequences in SIMAP and the data access and query functions have been improved. SIMAP assists biologists to query the up-to-date sequence space systematically and facilitates large-scale downstream projects in computational biology. Access to SIMAP is freely provided through the web portal for individuals (http://mips.gsf.de/simap/) and for programmatic access through DAS (http://webclu.bio.wzw.tum.de/das/) and Web-Service (http://mips.gsf.de/webservices/services/SimapService2.0?wsdl).

Complete genome sequence of Cronobacter turicensis LMG 23827, a food-borne pathogen causing deaths in neonates.

Here, we report the complete and annotated genome sequence of Cronobacter turicensis, an opportunistic food-borne pathogen, which is known as a rare but important cause of life-threatening neonatal infections. Among all proteins of C. turicensis, 223 have been annotated as virulence- and disease-related proteins.

Genomic insights into the metabolic potential of the polycyclic aromatic hydrocarbon degrading sulfate-reducing Deltaproteobacterium N47.

Anaerobic degradation of polycyclic aromatic hydrocarbons (PAHs) is an important process during natural attenuation of aromatic hydrocarbon spills. However, knowledge about metabolic potential and physiology of organisms involved in anaerobic degradation of PAHs is scarce. Therefore, we introduce the first genome of the sulfate-reducing Deltaproteobacterium N47 able to catabolize naphthalene, 2-methylnaphthalene, or 2-naphthoic acid as sole carbon source. Based on proteomics, we analysed metabolic pathways during growth on PAHs to gain physiological insights on anaerobic PAH degradation. The genomic assembly and taxonomic binning resulted in 17 contigs covering most of the sulfate reducer N47 genome according to general cluster of orthologous groups (COGs) analyses. According to the genes present, the Deltaproteobacterium N47 can potentially grow with the following sugars including d-mannose, d-fructose, d-galactose, α-d-glucose-1P, starch, glycogen, peptidoglycan and possesses the prerequisites for butanoic acid fermentation. Despite the inability for culture N47 to utilize NO(3) (-) as terminal electron acceptor, genes for nitrate ammonification are present. Furthermore, it is the first sequenced genome containing a complete TCA cycle along with the carbon monoxide dehydrogenase pathway. The genome contained a significant percentage of repetitive sequences and transposase-related protein domains enhancing the ability of genome evolution. Likewise, the sulfate reducer N47 genome contained many unique putative genes with unknown function, which are candidates for yet-unknown metabolic pathways.

Sequence-based prediction of type III secreted proteins.

The type III secretion system (TTSS) is a key mechanism for host cell interaction used by a variety of bacterial pathogens and symbionts of plants and animals including humans. The TTSS represents a molecular syringe with which the bacteria deliver effector proteins directly into the host cell cytosol. Despite the importance of the TTSS for bacterial pathogenesis, recognition and targeting of type III secreted proteins has up until now been poorly understood. Several hypotheses are discussed, including an mRNA-based signal, a chaperon-mediated process, or an N-terminal signal peptide. In this study, we systematically analyzed the amino acid composition and secondary structure of N-termini of 100 experimentally verified effector proteins. Based on this, we developed a machine-learning approach for the prediction of TTSS effector proteins, taking into account N-terminal sequence features such as frequencies of amino acids, short peptides, or residues with certain physico-chemical properties. The resulting computational model revealed a strong type III secretion signal in the N-terminus that can be used to detect effectors with sensitivity of approximately 71% and selectivity of approximately 85%. This signal seems to be taxonomically universal and conserved among animal pathogens and plant symbionts, since we could successfully detect effector proteins if the respective group was excluded from training. The application of our prediction approach to 739 complete bacterial and archaeal genome sequences resulted in the identification of between 0% and 12% putative TTSS effector proteins. Comparison of effector proteins with orthologs that are not secreted by the TTSS showed no clear pattern of signal acquisition by fusion, suggesting convergent evolutionary processes shaping the type III secretion signal. The newly developed program EffectiveT3 (http://www.chlamydiaedb.org) is the first universal in silico prediction program for the identification of novel TTSS effectors. Our findings will facilitate further studies on and improve our understanding of type III secretion and its role in pathogen-host interactions.

PEDANT covers all complete RefSeq genomes.

The PEDANT genome database provides exhaustive annotation of nearly 3000 publicly available eukaryotic, eubacterial, archaeal and viral genomes with more than 4.5 million proteins by a broad set of bioinformatics algorithms. In particular, all completely sequenced genomes from the NCBI's Reference Sequence collection (RefSeq) are covered. The PEDANT processing pipeline has been sped up by an order of magnitude through the utilization of precalculated similarity information stored in the similarity matrix of proteins (SIMAP) database, making it possible to process newly sequenced genomes immediately as they become available. PEDANT is freely accessible to academic users at http://pedant.gsf.de. For programmatic access Web Services are available at http://pedant.gsf.de/webservices.jsp.

The genome of the amoeba symbiont "Candidatus Amoebophilus asiaticus" reveals common mechanisms for host cell interaction among amoeba-associated bacteria.

Protozoa play host for many intracellular bacteria and are important for the adaptation of pathogenic bacteria to eukaryotic cells. We analyzed the genome sequence of "Candidatus Amoebophilus asiaticus," an obligate intracellular amoeba symbiont belonging to the Bacteroidetes. The genome has a size of 1.89 Mbp, encodes 1,557 proteins, and shows massive proliferation of IS elements (24% of all genes), although the genome seems to be evolutionarily relatively stable. The genome does not encode pathways for de novo biosynthesis of cofactors, nucleotides, and almost all amino acids. "Ca. Amoebophilus asiaticus" encodes a variety of proteins with predicted importance for host cell interaction; in particular, an arsenal of proteins with eukaryotic domains, including ankyrin-, TPR/SEL1-, and leucine-rich repeats, which is hitherto unmatched among prokaryotes, is remarkable. Unexpectedly, 26 proteins that can interfere with the host ubiquitin system were identified in the genome. These proteins include F- and U-box domain proteins and two ubiquitin-specific proteases of the CA clan C19 family, representing the first prokaryotic members of this protein family. Consequently, interference with the host ubiquitin system is an important host cell interaction mechanism of "Ca. Amoebophilus asiaticus". More generally, we show that the eukaryotic domains identified in "Ca. Amoebophilus asiaticus" are also significantly enriched in the genomes of other amoeba-associated bacteria (including chlamydiae, Legionella pneumophila, Rickettsia bellii, Francisella tularensis, and Mycobacterium avium). This indicates that phylogenetically and ecologically diverse bacteria which thrive inside amoebae exploit common mechanisms for interaction with their hosts, and it provides further evidence for the role of amoebae as training grounds for bacterial pathogens of humans.

Combined genomic and proteomic approaches identify gene clusters involved in anaerobic 2-methylnaphthalene degradation in the sulfate-reducing enrichment culture N47.

The highly enriched deltaproteobacterial culture N47 anaerobically oxidizes the polycyclic aromatic hydrocarbons naphthalene and 2-methylnaphthalene, with sulfate as the electron acceptor. Combined genome sequencing and liquid chromatography-tandem mass spectrometry-based shotgun proteome analyses were performed to identify genes and proteins involved in anaerobic aromatic catabolism. Proteome analysis of 2-methylnaphthalene-grown N47 cells resulted in the identification of putative enzymes catalyzing the anaerobic conversion of 2-methylnaphthalene to 2-naphthoyl coenzyme A (2-naphthoyl-CoA), as well as the reductive ring cleavage of 2-naphthoyl-CoA, leading to the formation of acetyl-CoA and CO(2). The glycyl radical-catalyzed fumarate addition to the methyl group of 2-methylnaphthalene is catalyzed by naphthyl-2-methyl-succinate synthase (Nms), composed of alpha-, beta-, and gamma-subunits that are encoded by the genes nmsABC. Located upstream of nmsABC is nmsD, encoding the Nms-activating enzyme, which harbors the characteristic [Fe(4)S(4)] cluster sequence motifs of S-adenosylmethionine radical enzymes. The bns gene cluster, coding for enzymes involved in beta-oxidation reactions converting naphthyl-2-methyl-succinate to 2-naphthoyl-CoA, was found four intervening open reading frames further downstream. This cluster consists of eight genes (bnsABCDEFGH) corresponding to 8.1 kb, which are closely related to genes for enzymes involved in anaerobic toluene degradation within the denitrifiers "Aromatoleum aromaticum" EbN1, Azoarcus sp. strain T, and Thauera aromatica. Another contiguous DNA sequence harbors the gene for 2-naphthoyl-CoA reductase (ncr) and 16 additional genes that were found to be expressed in 2-methylnaphthalene-grown cells. These genes code for enzymes that were supposed to catalyze the dearomatization and ring cleavage reactions converting 2-naphthoyl-CoA to acetyl-CoA and CO(2). Comparative sequence analysis of the four encoding subunits (ncrABCD) showed the gene product to have the closest similarity to the Azoarcus type of benzoyl-CoA reductase. The present work provides the first insight into the genetic basis of anaerobic 2-methylnaphthalene metabolism and delivers implications for understanding contaminant degradation.

Identification of enzymes involved in anaerobic benzene degradation by a strictly anaerobic iron-reducing enrichment culture.

Anaerobic benzene degradation was studied with a highly enriched iron-reducing culture (BF) composed of mainly Peptococcaceae-related Gram-positive microorganisms. The proteomes of benzene-, phenol- and benzoate-grown cells of culture BF were compared by SDS-PAGE. A specific benzene-expressed protein band of 60 kDa, which could not be observed during growth on phenol or benzoate, was subjected to N-terminal sequence analysis. The first 31 amino acids revealed that the protein was encoded by ORF 138 in the shotgun sequenced metagenome of culture BF. ORF 138 showed 43% sequence identity to phenylphosphate carboxylase subunit PpcA of Aromatoleum aromaticum strain EbN1. A LC/ESI-MS/MS-based shotgun proteomic analysis revealed other specifically benzene-expressed proteins with encoding genes located adjacent to ORF 138 on the metagenome. The protein products of ORF 137, ORF 139 and ORF 140 showed sequence identities of 37% to phenylphosphate carboxylase PpcD of A. aromaticum strain EbN1, 56% to benzoate-CoA ligase (BamY) of Geobacter metallireducens and 67% to 3-octaprenyl-4-hydroxybenzoate carboxy-lyase (UbiD/UbiX) of A. aromaticum strain EbN1 respectively. These genes are proposed as constituents of a putative benzene degradation gene cluster (∼ 17 kb) composed of carboxylase-related genes. The identified gene sequences suggest that the initial activation reaction in anaerobic benzene degradation is probably a direct carboxylation of benzene to benzoate catalysed by putative anaerobic benzene carboxylase (Abc). The putative Abc probably consists of several subunits, two of which are encoded by ORFs 137 and 138, and belongs to a family of carboxylases including phenylphosphate carboxylase (Ppc) and 3-octaprenyl-4-hydroxybenzoate carboxy-lyase (UbiD/UbiX).

SIMAP--structuring the network of protein similarities.

Protein sequences are the most important source of evolutionary and functional information for new proteins. In order to facilitate the computationally intensive tasks of sequence analysis, the Similarity Matrix of Proteins (SIMAP) database aims to provide a comprehensive and up-to-date dataset of the pre-calculated sequence similarity matrix and sequence-based features like InterPro domains for all proteins contained in the major public sequence databases. As of September 2007, SIMAP covers approximately 17 million proteins and more than 6 million non-redundant sequences and provides a complete annotation based on InterPro 16. Novel features of SIMAP include a new, portlet-based web portal providing multiple, structured views on retrieved proteins and integration of protein clusters and a unique search method for similar domain architectures. Access to SIMAP is freely provided for academic use through the web portal for individuals at http://mips.gsf.de/simap/and through Web Services for programmatic access at http://mips.gsf.de/webservices/services/SimapService2.0?wsdl.

Comprehensive in silico prediction and analysis of chlamydial outer membrane proteins reflects evolution and life style of the Chlamydiae.

BACKGROUND: Chlamydiae are obligate intracellular bacteria comprising some of the most important bacterial pathogens of animals and humans. Although chlamydial outer membrane proteins play a key role for attachment to and entry into host cells, only few have been described so far. We developed a comprehensive, multiphasic in silico approach, including the calculation of clusters of orthologues, to predict outer membrane proteins using conservative criteria. We tested this approach using Escherichia coli (positive control) and Bacillus subtilis (negative control), and applied it to five chlamydial species; Chlamydia trachomatis, Chlamydia muridarum, Chlamydia (a.k.a. Chlamydophila) pneumoniae, Chlamydia (a.k.a. Chlamydophila) caviae, and Protochlamydia amoebophila. RESULTS: In total, 312 chlamydial outer membrane proteins and lipoproteins in 88 orthologous clusters were identified, including 238 proteins not previously recognized to be located in the outer membrane. Analysis of their taxonomic distribution revealed an evolutionary conservation among Chlamydiae, Verrucomicrobia, Lentisphaerae and Planctomycetes as well as lifestyle-dependent conservation of the chlamydial outer membrane protein composition. CONCLUSION: This analysis suggested a correlation between the outer membrane protein composition and the host range of chlamydiae and revealed a common set of outer membrane proteins shared by these intracellular bacteria. The collection of predicted chlamydial outer membrane proteins is available at the online database pCOMP http://www.microbial-ecology.net/pcomp and might provide future guidance in the quest for anti-chlamydial vaccines.

probeCheck--a central resource for evaluating oligonucleotide probe coverage and specificity.

The web server probeCheck, freely accessible at http://www.microbial-ecology.net/probecheck, provides a pivotal forum for rapid specificity and coverage evaluations of probes and primers against selected databases of phylogenetic and functional marker genes. Currently, 24 widely used sequence collections including the Ribosomal Database Project (RDP) II, Greengenes, SILVA and the Functional Gene Pipeline/Repository can be queried. For this purpose, probeCheck integrates a new online version of the popular ARB probe match tool with free energy (DeltaG) calculations for each perfectly matched and mismatched probe-target hybrid, allowing assessment of the theoretical binding stabilities of oligo-target and non-target hybrids. For each output sequence, the accession number, the GenBank taxonomy and a link to the respective entry at GenBank, EMBL and, if applicable, the query database are displayed. Filtering options allow customizing results on the output page. In addition, probeCheck is linked with probe match tools of RDP II and Greengenes, NCBI blast, the Oligonucleotide Properties Calculator, the two-state folding tool of the DINAMelt server and the rRNA-targeted probe database probeBase. Taken together, these features provide a multifunctional platform with maximal flexibility for the user in the choice of databases and options for the evaluation of published and newly developed probes and primers.

SIMAP: the similarity matrix of proteins.

Similarity Matrix of Proteins (SIMAP) (http://mips.gsf.de/simap) provides a database based on a pre-computed similarity matrix covering the similarity space formed by >4 million amino acid sequences from public databases and completely sequenced genomes. The database is capable of handling very large datasets and is updated incrementally. For sequence similarity searches and pairwise alignments, we implemented a grid-enabled software system, which is based on FASTA heuristics and the Smith-Waterman algorithm. Our ProtInfo system allows querying by protein sequences covered by the SIMAP dataset as well as by fragments of these sequences, highly similar sequences and title words. Each sequence in the database is supplemented with pre-calculated features generated by detailed sequence analyses. By providing WWW interfaces as well as web-services, we offer the SIMAP resource as an efficient and comprehensive tool for sequence similarity searches.

SIMAP--the similarity matrix of proteins.

MOTIVATION: Sequence similarity searches are of great importance in bioinformatics. Exhaustive searches for homologous proteins in databases are computationally expensive and can be replaced by a database of pre-calculated homologies in many cases. Retrieving similarities from an incrementally updated database instead of repeatedly recalculating them should provide homologs much faster and frees computational resources for other purposes. RESULTS: We have implemented SIMAP-a database containing the similarity space formed by almost all amino acid sequences from public databases and completely sequenced genomes. The database is capable of handling very large datasets and allows incremental updates. We have implemented a powerful backbone for similarity computation, which is based on FASTA heuristics. By providing WWW interfaces as well as web services, we make our data accessible to the worldwide community. We have also adapted procedures to detect putative orthologs as example applications. AVAILABILITY: The SIMAP portal page providing links to SIMAP services is publicly available: http://mips.gsf.de/services/analysis/simap/. The web services can be accessed under http://mips.gsf.de/proj/hobitws/services/RPCSimapService?wsdl and http://mips.gsf.de/proj/hobitws/services/DocSimapService?wsdl

Distinct gene set in two different lineages of ammonia-oxidizing archaea supports the phylum Thaumarchaeota.

Globally distributed archaea comprising ammonia oxidizers of moderate terrestrial and marine environments are considered the most abundant archaeal organisms on Earth. Based on 16S rRNA phylogeny, initial assignment of these archaea was to the Crenarchaeota. By contrast, features of the first genome sequence from a member of this group suggested that they belong to a novel phylum, the Thaumarchaeota. Here, we re-investigate the Thaumarchaeota hypothesis by including two newly available genomes, that of the marine ammonia oxidizer Nitrosopumilus maritimus and that of Nitrososphaera gargensis, a representative of another evolutionary lineage within this group predominantly detected in terrestrial environments. Phylogenetic studies based on r-proteins and other core genes, as well as comparative genomics, confirm the assignment of these organisms to a separate phylum and reveal a Thaumarchaeota-specific set of core informational processing genes, as well as potentially ancestral features of the archaea.