Towards a more accurate annotation of tyrosine-based site-specific recombinases in bacterial genomes.

BACKGROUND: Tyrosine-based site-specific recombinases (TBSSRs) are DNA breaking-rejoining enzymes. In bacterial genomes, they play a major role in the comings and goings of mobile genetic elements (MGEs), such as temperate phage genomes, integrated conjugative elements (ICEs) or integron cassettes. TBSSRs are also involved in the segregation of plasmids and chromosomes, the resolution of plasmid dimers and of co-integrates resulting from the replicative transposition of transposons. With the aim of improving the annotation of TBSSR genes in genomic sequences and databases, which so far is far from robust, we built a set of over 1,300 TBSSR protein sequences tagged with their genome of origin. We organized them in families to investigate: i) whether TBSSRs tend to be more conserved within than between classes of MGE types and ii) whether the (sub)families may help in understanding more about the function of TBSSRs associated in tandem or trios on plasmids and chromosomes. RESULTS: A total of 67% of the TBSSRs in our set are MGE type specific. We define a new class of actinobacterial transposons, related to Tn554, containing one abnormally long TBSSR and one of typical size, and we further characterize numerous TBSSRs trios present in plasmids and chromosomes of α- and ß-proteobacteria. CONCLUSIONS: The simple in silico procedure described here, which uses a set of reference TBSSRs from defined MGE types, could contribute to greatly improve the annotation of tyrosine-based site-specific recombinases in plasmid, (pro)phage and other integrated MGE genomes. It also reveals TBSSRs families whose distribution among bacterial taxa suggests they mediate lateral gene transfer.

A modular view of the bacteriophage genomic space: identification of host and lifestyle marker modules.

Bacteriophage genomes can be regarded as an ensemble of modules which are accessible to the whole phage population via recombination. The time spent by prophages in the bacterial host provides them with the opportunity to exchange modules with other prophages or infecting phages. Here we analyze the modular structure of a set of 457 phages and 760 prophages extracted from completely sequenced bacterial genomes using the ACLAME database and its associated tools. We identified 91 modules of proteins with similar phylogenetic profiles. Of these, 25 and 6 are associated with temperate and virulent phages, respectively; 57 are restricted to a host or small group of hosts; and 55 could be annotated with a phage function. We use the transposable phages as a study case and show how the inclusion of prophages allows us to unveil new types of genome organization (i.e. novel module combinations) and obtain insight into the host range for this particular group, highlighting the utility of prophage prediction to better characterize phage diversity.

Diversity of bacterial type II toxin-antitoxin systems: a comprehensive search and functional analysis of novel families.

Type II toxin-antitoxin (TA) systems are generally composed of two genes organized in an operon, encoding a labile antitoxin and a stable toxin. They were first discovered on plasmids where they contribute to plasmid stability by a phenomenon denoted as 'addiction', and subsequently in bacterial chromosomes. To discover novel families of antitoxins and toxins, we developed a bioinformatics approach based on the 'guilt by association' principle. Extensive experimental validation in Escherichia coli of predicted antitoxins and toxins increased significantly the number of validated systems and defined novel toxin and antitoxin families. Our data suggest that toxin families as well as antitoxin families originate from distinct ancestors that were assembled multiple times during evolution. Toxin and antitoxin families found on plasmids tend to be promiscuous and widespread, indicating that TA systems move through horizontal gene transfer. We propose that due to their addictive properties, TA systems are likely to be maintained in chromosomes even though they do not necessarily confer an advantage to their bacterial hosts. Therefore, addiction might play a major role in the evolutionary success of TA systems both on mobile genetic elements and in bacterial chromosomes.

Towards Viral Genome Annotation Standards, Report from the 2010 NCBI Annotation Workshop.

Improvements in DNA sequencing technologies portend a new era in virology and could possibly lead to a giant leap in our understanding of viral evolution and ecology. Yet, as viral genome sequences begin to fill the world's biological databases, it is critically important to recognize that the scientific promise of this era is dependent on consistent and comprehensive genome annotation. With this in mind, the NCBI Genome Annotation Workshop recently hosted a study group tasked with developing sequence, function, and metadata annotation standards for viral genomes. This report describes the issues involved in viral genome annotation and reviews policy recommendations presented at the NCBI Annotation Workshop.

ACLAME: a CLAssification of Mobile genetic Elements, update 2010.

The ACLAME database is dedicated to the collection, analysis and classification of sequenced mobile genetic elements (MGEs, in particular phages and plasmids). In addition to providing information on the MGEs content, classifications are available at various levels of organization. At the gene/protein level, families group similar sequences that are expected to share the same function. Families of four or more proteins are manually assigned with a functional annotation using the GeneOntology and the locally developed ontology MeGO dedicated to MGEs. At the genome level, evolutionary cohesive modules group sets of protein families shared among MGEs. At the population level, networks display the reticulate evolutionary relationships among MGEs. To increase the coverage of the phage sequence space, ACLAME version 0.4 incorporates 760 high-quality predicted prophages selected from the Prophinder database. Most of the data can be downloaded from the freely accessible ACLAME web site (http://aclame.ulb.ac.be). The BLAST interface for querying the database has been extended and numerous tools for in-depth analysis of the results have been added.

Genetic diversity of Group A Streptococcus M protein: implications for typing and vaccine development.

Group A Streptococci (GAS) are classified into 180 emm-types based on the sequence of the amino-terminal hyper-variable region of the M surface protein. The genetic relatedness of the whole surface-exposed part of M was investigated in well-characterized Belgian and Brazilian GAS isolates which belong to different epidemiological and clinical landscapes. Despite a small number of different emm-types and an apparent low diversity in the Belgian isolates (as revealed by the emm-typing method), our data showed that the overall genetic diversity of the M proteins was higher than expected. On the contrary, the M proteins from the Brazilian isolates were genetically highly related. Since M is a multi-functional protein, an analysis of the whole surface-exposed sequence that takes into account the different functional domains may provide tools for typing as well as for analyzing the molecular mechanisms of M virulence or defining vaccine strategies.

Prophinder: a computational tool for prophage prediction in prokaryotic genomes.

UNLABELLED: Prophinder is a prophage prediction tool coupled with a prediction database, a web server and web service. Predicted prophages will help to fill the gaps in the current sparse phage sequence space, which should cover an estimated 100 million species. Systematic and reliable predictions will enable further studies of prophages contribution to the bacteriophage gene pool and to better understand gene shuffling between prophages and phages infecting the same host. AVAILABILITY: Softare is available at http://aclame.ulb.ac.be/prophinder

Reticulate representation of evolutionary and functional relationships between phage genomes.

Bacteriophage genomes show pervasive mosaicism, indicating the importance of horizontal gene exchange in their evolution. Phage genomes represent unique combinations of modules, each of them with a different phylogenetic history. The traditional classification, based on a variety of criteria such as nucleic acid type (single/double-stranded DNA/RNA), morphology, and host range, appeared inconsistent with sequence analyses. With the genomic era, an ever increasing number of sequenced phages cannot be classified, in part due to a lack of morphological information and in part to the intrinsic incapability of tree-based methods to efficiently deal with mosaicism. This problem led some virologists to call for a moratorium on the creation of additional taxa in the order Caudovirales, in order to let virologists discuss classification schemes that might better suit phage evolution. In this context, we propose a framework for a reticulate classification of phages based on gene content. Starting from gene families, we built a weighted graph, where nodes represent phages and edges represent phage-phage similarities in terms of shared genes. We then apply various measures of graph topology to analyze the resulting graph. Most double-stranded DNA phages are found in a single component. The values of the clustering coefficient and closeness distinguish temperate from virulent phages, whereas chimeric phages are characterized by a high betweenness coefficient. We apply a 2-step clustering method to this graph to generate a reticulate classification of phages: Each phage is associated with a membership vector, which quantitatively characterizes its membership to the set of clusters. Furthermore, we cluster genes based on their "phylogenetic profiles" to define "evolutionary cohesive modules." In virulent phages, evolutionary modules span several functional categories, whereas in temperate phages they correspond better to functional modules. Moreover, despite the fact that modules only cover a fraction of all phage genes, phage groups can be distinguished by their different combination of modules, serving the bases for a higher level reticulate classification. These 2 classification schemes provide an automatic and dynamic way of representing the relationships within the phage population and can be extended to include newly sequenced phage genomes, as well as other types of genetic elements.

PhiGO, a phage ontology associated with the ACLAME database.

As is the case for other genomes and metagenomes, complete nucleotide sequences of bacteriophages and archaeviruses have become increasingly numerous and require robust annotation tools. We present here the first version of a phage ontology, PhiGO, which should contribute to more informational annotation of gene products in phage and prophage sequences. PhiGO uses the Gene Ontology schema, the de facto standard for describing knowledge about gene products across many databases.

Analysis of the phage sequence space: the benefit of structured information.

Phages are the most abundant biological entities on Earth and are central players in the evolution of their bacterial hosts and the emergence of new pathogens. In addition, they bear an enormous potential for the development of new drugs, therapies or nanotechnologies. As a result, interest in phages is reviving. In the genomic era, our perspective on the phage sequence space remains incredibly sparse. The modular and combinatorial structure of phage genomes is largely documented. It is confirmed by new sequence information and it fuels a recurrent debate on the need to revise phage taxonomy. The absence of structured, computer readable information on phages is a major bottleneck for an extensive global analysis of phage genomes and their relationships, but such information is essential to reassess phage classification. Based on the ACLAME database, which is dedicated to the organization and analysis of prokaryotic mobile genetic elements, we discuss here how structured information on phage-encoded proteins helps global in silico analysis and allows the prediction of prophages in bacterial genome sequences, providing access to additional phage sequence information.

A first global analysis of plasmid encoded proteins in the ACLAME database.

Many plasmids are mobile genetic elements (MGEs) and, as other members of that group of DNA entities, their genomes display a mosaic and combinatorial structure, making their classification extremely difficult. As other MGEs, plasmids play a major role in horizontal transfer of genetic materials and genome reorganization. Yet, the full impact of such phenomenon on major properties of the host cell, such as pathogenicity, the ability to use new carbon sources or resistance to antibiotics, remains to be fully assessed. More and more complete plasmid genome sequences are available. However, in the absence of standards for storing plasmid sequence data and annotating genes and gene products on sequenced plasmid genomes, the resulting information remains rather limited. Using 503 sequenced plasmids organized in the ACLAME database, we discuss how, by structuring information on the genomes, their host and the proteins they code for, one can gain access to either global or more detailed analysis of the plasmid sequence information, as illustrated by a network representation of the relationships between plasmids.

Mobile genetic elements: the agents of open source evolution.

Horizontal genomics is a new field in prokaryotic biology that is focused on the analysis of DNA sequences in prokaryotic chromosomes that seem to have originated from other prokaryotes or eukaryotes. However, it is equally important to understand the agents that effect DNA movement: plasmids, bacteriophages and transposons. Although these agents occur in all prokaryotes, comprehensive genomics of the prokaryotic mobile gene pool or 'mobilome' lags behind other genomics initiatives owing to challenges that are distinct from cellular chromosomal analysis. Recent work shows promise of improved mobile genetic element (MGE) genomics and consequent opportunities to take advantage - and avoid the dangers - of these 'natural genetic engineers'. This review describes MGEs, their properties that are important in horizontal gene transfer, and current opportunities to advance MGE genomics.

Variations in plant metallothioneins: the heavy metal hyperaccumulator Thlaspi caerulescens as a study case.

Plant metallothioneins (MTs) are extremely diverse and are thought to be involved in metal homeostasis or detoxification. Thlaspi caerulescens is a model Zn/Cd hyperaccumulator and thus constitutes an ideal system to study the variability of these MTs. Two T. caerulescens cDNAs (accession: 665511; accession: 665515), that are highly homologous to type 1 and type 2 Arabidopsis thaliana MTs, have been isolated using a functional screen for plant cDNAs that confer Cd tolerance to yeast. However, TcMT1 has a much shorter N-terminal domain than that of A. thaliana and so lacks Cys motifs conserved through all the plant MTs classified as type 1. A systematic search in plant databases allowed the detection of MT-related sequences. Sixty-four percent fulfil the criteria for MT classification described in Cobbett and Goldsbrough (2002) and further extend our knowledge about other conserved residues that might play an important role in plant MT structure. In addition, 34% of the total MT-related sequences cannot be classified strictly as they display modifications in the conserved residues according to the current plant MTs' classification. The significance of this variability in plant MT sequences is discussed. Functional complementation in yeast was used to assess whether these variations may alter the MTs' function in T. caerulescens. Regulation of the expression of MTs in T. caerulescens was also investigated. TcMT1 and TcMT2 display higher expression in T. caerulescens than in A. thaliana. Moreover, their differential expression patterns in organs and in response to metal exposure, suggest that the two types of MTs may have diverse roles and functions in T. caerulescens.

Assessment of CAPRI predictions in rounds 3-5 shows progress in docking procedures.

The current status of docking procedures for predicting protein-protein interactions starting from their three-dimensional (3D) structure is reassessed by evaluating blind predictions, performed during 2003-2004 as part of Rounds 3-5 of the community-wide experiment on Critical Assessment of PRedicted Interactions (CAPRI). Ten newly determined structures of protein-protein complexes were used as targets for these rounds. They comprised 2 enzyme-inhibitor complexes, 2 antigen-antibody complexes, 2 complexes involved in cellular signaling, 2 homo-oligomers, and a complex between 2 components of the bacterial cellulosome. For most targets, the predictors were given the experimental structures of 1 unbound and 1 bound component, with the latter in a random orientation. For some, the structure of the free component was derived from that of a related protein, requiring the use of homology modeling. In some of the targets, significant differences in conformation were displayed between the bound and unbound components, representing a major challenge for the docking procedures. For 1 target, predictions could not go to completion. In total, 1866 predictions submitted by 30 groups were evaluated. Over one-third of these groups applied completely novel docking algorithms and scoring functions, with several of them specifically addressing the challenge of dealing with side-chain and backbone flexibility. The quality of the predicted interactions was evaluated by comparison to the experimental structures of the targets, made available for the evaluation, using the well-agreed-upon criteria used previously. Twenty-four groups, which for the first time included an automatic Web server, produced predictions ranking from acceptable to highly accurate for all targets, including those where the structures of the bound and unbound forms differed substantially. These results and a brief survey of the methods used by participants of CAPRI Rounds 3-5 suggest that genuine progress in the performance of docking methods is being achieved, with CAPRI acting as the catalyst.

Evidence for copper homeostasis function of metallothionein (MT3) in the hyperaccumulator Thlaspi caerulescens.

Metallothioneins chelate metals and consequently may be a control point of metal homeostasis. Homologous to type 3 metallothioneins, TcMT3 cDNA was identified in the Cd/Zn hyperaccumulator, Thlaspi caerulescens. TcMT3 amino acid sequence showed modifications in the Cys positions when compared with its Arabidopsis orthologue. A structural model established that the MT3 carboxyterminal domain is similar to the beta domain of animal metallothioneins and predicts a smaller cavity to chelate metals for A. thaliana than for T. caerulescens. Functional testing in yeast and Northern blot analysis added further evidence for adaptative variations of MT3 for the maintenance of Cu homeostasis in a metal hyperaccumulator.

ACLAME: a CLAssification of Mobile genetic Elements.

The ACLAME database (http://aclame.ulb.ac.be) is a collection and classification of prokaryotic mobile genetic elements (MGEs) from various sources, comprising all known phage genomes, plasmids and transposons. In addition to providing information on the full genomes and genetic entities, it aims to build a comprehensive classification of the functional modules of MGEs at the protein, gene and higher levels. This first version contains a comprehensive classification of 5069 proteins from 119 DNA bacteriophages into over 400 functional families. This classification was produced automatically using TRIBE-MCL, a graph-theory-based Markov clustering algorithm that uses sequence measures as input, and then manually curated. Manual curation was aided by consulting annotations available in public databases retrieved through additional sequence similarity searches using Psi-Blast and Hidden Markov Models. The database is publicly accessible and open to expert volunteers willing to participate in its curation. Its web interface allows browsing as well as querying the classification. The main objectives are to collect and organize in a rational way the complexity inherent to MGEs, to extend and improve the inadequate annotation currently associated with MGEs and to screen known genomes for the validation and discovery of new MGEs.

The biphenyl- and 4-chlorobiphenyl-catabolic transposon Tn4371, a member of a new family of genomic islands related to IncP and Ti plasmids.

The nucleotide sequence of the biphenyl catabolic transposon Tn4371 has been completed and analyzed. It confirmed that the element has a mosaic structure made of several building blocks. In addition to previously identified genes coding for a tyrosine recombinase related to phage integrases and for biphenyl degradation enzymes very similar to those of Achromobacter georgiopolitanum KKS102, Tn4371 carries many plasmid-related genes involved in replication, partition, and other, as-yet-unknown, plasmid functions. One gene cluster contains most of the genes required to express a type IV secretion-mating pair formation apparatus coupled with a TraG ATPase, all of which are related to those found on IncP and Ti plasmids. Orthologues of all Tn4371 plasmid-related genes and of the tyrosine recombinase gene were found, with a very similar organization, in the chromosome of Ralstonia solanacearum and on the yet-to-be-determined genomic sequences of Erwinia chrysanthemi and Azotobacter vinelandii. In each of these chromosomal segments, conserved segments were separated by different groups of genes, which also differed from the Tn4371 bph genes. The conserved blocks of genes were also identified, in at least two copies, in the chromosome of Ralstonia metallidurans CH34. Tn4371 thus appears to represent a new family of potentially mobile genomic islands with a broad host range since they reside in a wide range of soil proteobacteria, including plant pathogens.

Assessment of blind predictions of protein-protein interactions: current status of docking methods.

The current status of docking procedures for predicting protein-protein interactions starting from their three-dimensional structure is assessed from a first major evaluation of blind predictions. This evaluation was performed as part of a communitywide experiment on Critical Assessment of PRedicted Interactions (CAPRI). Seven newly determined structures of protein-protein complexes were available as targets for this experiment. These were the complexes between a kinase and its protein substrate, between a T-cell receptor beta-chain and a superantigen, and five antigen-antibody complexes. For each target, the predictors were given the experimental structures of the free components, or of one free and one bound component in a random orientation. The structure of the complex was revealed only at the time of the evaluation. A total of 465 predictions submitted by 19 groups were evaluated. These groups used a wide range of algorithms and scoring functions, some of which were completely novel. The quality of the predicted interactions was evaluated by comparing residue-residue contacts and interface residues to those in the X-ray structures and by analyzing the fit of the ligand molecules (the smaller of the two proteins in the complex) or of interface residues only, in the predicted versus target complexes. A total of 14 groups produced predictions, ranking from acceptable to highly accurate for five of the seven targets. The use of available biochemical and biological information, and in one instance structural information, played a key role in achieving this result. It was essential for identifying the native binding modes for the five correctly predicted targets, including the kinase-substrate complex where the enzyme changes conformation on association. But it was also the cause for missing the correct solution for the two remaining unpredicted targets, which involve unexpected antigen-antibody binding modes. Overall, this analysis reveals genuine progress in docking procedures but also illustrates the remaining serious limitations and points out the need for better scoring functions and more effective ways for handling conformational flexibility.

MaxBench: evaluation of sequence and structure comparison methods.

SUMMARY: MaxBench is a web-based system available for evaluating the results of sequence and structure comparison methods, based on the SCOP protein domain classification. The system makes it easy for developers to both compare the overall performance of their methods to standard algorithms and investigate the results of individual comparisons. AVAILABILITY: http://www.sanger.ac.uk/Users/lp1/MaxBench/