ACLAME: a CLAssification of Mobile genetic Elements, update 2010

An automated approach for the identification of horizontal gene transfers from complete genomes reveals the rhizome of Rickettsiales

BACKGROUND

Horizontal gene transfer (HGT) is considered to be a major force driving the evolutionary history of prokaryotes. HGT is widespread in prokaryotes, contributing to the genomic repertoire of prokaryotic organisms, and is particularly apparent in Rickettsiales genomes. Gene gains from both distantly and closely related organisms play crucial roles in the evolution of bacterial genomes. In this work, we focus on genes transferred from distantly related species into Rickettsiales species.

RESULTS

We developed an automated approach for the detection of HGT from other organisms (excluding alphaproteobacteria) into Rickettsiales genomes. Our systematic approach consisted of several specialized features including the application of a parsimony method for inferring phyletic patterns followed by blast filter, automated phylogenetic reconstruction and the application of patterns for HGT detection. We identified 42 instances of HGT in 31 complete Rickettsiales genomes, of which 38 were previously unidentified instances of HGT from Anaplasma, Wolbachia, Candidatus Pelagibacter ubique and Rickettsia genomes. Additionally, putative cases with no phylogenetic support were assigned gene ontology terms. Overall, these transfers could be characterized as "rhizome-like".

CONCLUSIONS

Our analysis provides a comprehensive, systematic approach for the automated detection of HGTs from several complete proteome sequences that can be applied to detect instances of HGT within other genomes of interest.

Contribution of lateral gene transfers to the genome composition and parasitic ability of root-knot nematodes

Lateral gene transfers (LGT), species to species transmission of genes by means other than direct inheritance from a common ancestor, have played significant role in shaping prokaryotic genomes and are involved in gain or transfer of important biological processes. Whether LGT significantly contributed to the composition of an animal genome is currently unclear. In nematodes, multiple LGT are suspected to have favored emergence of plant-parasitism. With the availability of whole genome sequences it is now possible to assess whether LGT have significantly contributed to the composition of an animal genome and to establish a comprehensive list of these events. We generated clusters of homologous genes and automated phylogenetic inference, to detect LGT in the genomes of root-knot nematodes and found that up to 3.34% of the genes originate from LGT of non-metazoan origin. After their acquisition, the majority of genes underwent series of duplications. Compared to the rest of the genes in these species, several predicted functional categories showed a skewed distribution in the set of genes acquired via LGT. Interestingly, functions related to metabolism, degradation or modification of carbohydrates or proteins were substantially more frequent. This suggests that genes involved in these processes, related to a parasitic lifestyle, have been more frequently fixed in these parasites after their acquisition. Genes from soil bacteria, including plant-pathogens were the most frequent closest relatives, suggesting donors were preferentially bacteria from the rhizosphere. Several of these bacterial genes are plasmid-borne, pointing to a possible role of these mobile genetic elements in the transfer mechanism. Our analysis provides the first comprehensive description of the ensemble of genes of non-metazoan origin in an animal genome. Besides being involved in important processes regarding plant-parasitism, genes acquired via LGT now constitute a substantial proportion of protein-coding genes in these nematode genomes.

TnpPred: A Web Service for the Robust Prediction of Prokaryotic Transposases

Transposases (Tnps) are enzymes that participate in the movement of insertion sequences (ISs) within and between genomes. Genes that encode Tnps are amongst the most abundant and widely distributed genes in nature. However, they are difficult to predict bioinformatically and given the increasing availability of prokaryotic genomes and metagenomes, it is incumbent to develop rapid, high quality automatic annotation of ISs. This need prompted us to develop a web service, termed TnpPred for Tnp discovery. It provides better sensitivity and specificity for Tnp predictions than given by currently available programs as determined by ROC analysis. TnpPred should be useful for improving genome annotation. The TnpPred web service is freely available for noncommercial use.

Comparative sequence analysis of a multidrug-resistant plasmid from Aeromonas hydrophila

Aeromonas hydrophila is a pathogenic bacterium that has been implicated in fish, animal, and human disease. Recently, a multidrug resistance (MDR) plasmid, pR148, was isolated from A. hydrophila obtained from a tilapia (Oreochromis niloticus) farm in Thailand. pR148 is a 165,906-bp circular plasmid containing 147 coding regions showing highest similarity to pNDM-1_Dok1, an MDR plasmid isolated from a human pathogen. pR148 was also very similar to other IncA/C plasmids isolated from humans, animals, food, and fish. pR148 contains a mercuric resistance operon and encodes the complete set of genes for the type 4 secretion system. pR148 encodes a Tn21 type transposon. This transposon contains the drug resistance genes qacH, bla(OXA-10), aadA1, and sul1 in a class 1 integron; tetA and tetR in transposon Tn1721; and catA2 and a duplicate sul1 in a locus showing 100% similarity to IncU plasmids isolated from fish. The bla(OXA-10) and aadA1 genes showed 100% similarity to those from the Acinetobacter baumannii AYE genome. The similarity of pR148 to a human pathogen-derived plasmid indicates that the plasmids were either transferred between different genera or that they are derived from a common origin. Previous studies have shown that IncA/C plasmids retain a conserved backbone, while the accessory region points to lateral gene transfer. These observations point out the dangers of indiscriminate use of antibiotics in humans and in animals and the necessity of understanding how drug resistance determinants are disseminated and transferred.

Contribution of phage-derived genomic islands to the virulence of facultative bacterial pathogens

Facultative pathogens have extremely dynamic pan-genomes, to a large extent derived from bacteriophages and other mobile elements. We developed a simple approach to identify phage-derived genomic islands and apply it to show that pathogens from diverse bacterial genera are significantly enriched in clustered phage-derived genes compared with related benign strains. These findings show that genome expansion by integration of prophages containing virulence factors is a major route of evolution of facultative bacterial pathogens.

Large plasmids of Escherichia coli and Salmonella encode highly diverse arrays of accessory genes on common replicon families

Plasmids are important in evolution and adaptation of host bacteria, yet we lack a comprehensive picture of their own natural variation. We used replicon typing and RFLP analysis to assess diversity and distribution of plasmids in the ECOR, SARA, SARB and SARC reference collections of Escherichia coli and Salmonella. Plasmids, especially large (≥30 kb) plasmids, are abundant in these collections. Host species and genotype clearly impact plasmid prevalence; plasmids are more abundant in ECOR than SAR, but, within ECOR, subgroup B2 strains have the fewest large plasmids. The majority of large plasmids have unique RFLP patterns, suggesting high variation, even within dominant replicon families IncF and IncI1. We found only four conserved plasmid types within ECOR, none of which are widely distributed. Within SAR, conserved plasmid types are primarily serovar-specific, including a pSLT-like plasmid in 13 Typhimurium strains. Conservation of pSLT contrasts with variability of other plasmids, suggesting evolution of serovar-specific virulence plasmids is distinct from that of most enterobacterial plasmids. We sequenced a conserved serovar Heidelberg plasmid but did not detect virulence or antibiotic resistance genes. Our data illustrate the high degree of natural variation in large plasmids of E. coli and Salmonella, even among plasmids sharing backbone genes.

Bacteriophages and their genomes

Bacteriophages occupy a unique position in biology, representing an absolute majority of all organisms in the biosphere. Because their genomes are relatively small, elucidating the genetic diversity of the phage population, deciphering their origins, and identifying the evolutionary mechanisms that shape the population would seem readily feasible. And yet the pace of phage genome characterization has slowed over the past three years, reflecting in part a need to transition from sequencing known and well-characterized bacteriophages to the isolation and comparative analysis of new isolates. The current state of bacteriophage genomics shows that the genetic diversity of the population is very high, that phages have been actively evolving for billions of years with active engagement of horizontal genetic exchange, and that their genomes are consequently pervasively mosaic in their architectures. But we have barely scratched the surface and the next years of phage genome exploration promise to be especially revealing.

Going viral: next-generation sequencing applied to phage populations in the human gut

Over the past decade, researchers have begun to characterize viral diversity using metagenomic methods. These studies have shown that viruses, the majority of which infect bacteria, are probably the most genetically diverse components of the biosphere. Here, we briefly review the incipient rise of a phage biology renaissance, which has been catalysed by advances in next-generation sequencing. We explore how work characterizing phage diversity and lifestyles in the human gut is changing our view of ourselves as supra-organisms. Finally, we discuss how a renewed appreciation of phage dynamics may yield new applications for phage therapies designed to manipulate the structure and functions of our gut microbiomes.

Stability of a Pseudomonas putida KT2440 bacteriophage-carried genomic island and its impact on rhizosphere fitness

The stability of seven genomic islands of Pseudomonas putida KT2440 with predicted potential for mobilization was studied in bacterial populations associated with the rhizosphere of corn plants by multiplex PCR. DNA rearrangements were detected for only one of them (GI28), which was lost at high frequency. This genomic island of 39.4 kb, with 53 open reading frames, shows the characteristic organization of genes belonging to tailed phages. We present evidence indicating that it corresponds to the lysogenic state of a functional bacteriophage that we have designated Pspu28. Integrated and rarely excised forms of Pspu28 coexist in KT2440 populations. Pspu28 is self-transmissible, and an excisionase is essential for its removal from the bacterial chromosome. The excised Pspu28 forms a circular element that can integrate into the chromosome at a specific location, att sites containing a 17-bp direct repeat sequence. Excision/insertion of Pspu28 alters the promoter sequence and changes the expression level of PP_1531, which encodes a predicted arsenate reductase. Finally, we show that the presence of Pspu28 in the lysogenic state has a negative effect on bacterial fitness in the rhizosphere under conditions of intraspecific competition, thus explaining why clones having lost this mobile element are recovered from that environment.

Dynamic viral populations in hypersaline systems as revealed by metagenomic assembly

Viruses of the Bacteria and Archaea play important roles in microbial evolution and ecology, and yet viral dynamics in natural systems remain poorly understood. Here, we created de novo assemblies from 6.4 Gbp of metagenomic sequence from eight community viral concentrate samples, collected from 12 h to 3 years apart from hypersaline Lake Tyrrell (LT), Victoria, Australia. Through extensive manual assembly curation, we reconstructed 7 complete and 28 partial novel genomes of viruses and virus-like entities (VLEs, which could be viruses or plasmids). We tracked these 35 populations across the eight samples and found that they are generally stable on the timescale of days and transient on the timescale of years, with some exceptions. Cross-detection of the 35 LT populations in three previously described haloviral metagenomes was limited to a few genes, and most previously sequenced haloviruses were not detected in our samples, though 3 were detected upon reducing our detection threshold from 90% to 75% nucleotide identity. Similar results were obtained when we applied our methods to haloviral metagenomic data previously reported from San Diego, CA: 10 contigs that we assembled from that system exhibited a variety of detection patterns on a timescale of weeks to 1 month but were generally not detected in LT. Our results suggest that most haloviral populations have a limited or, possibly, a temporally variable global distribution. This study provides high-resolution insight into viral biogeography and dynamics and it places "snapshot" viral metagenomes, collected at a single time and location, in context.

Characterization of the Salmonella bacteriophage vB_SenS-Ent1

The bacteriophage vB_SenS-Ent1 (Ent1) is a member of the family Siphoviridae of tailed bacteriophages and infects a broad range of serovars of the enteric pathogen Salmonella enterica. The virion particle is composed of an icosahedral head 64 nm in diameter and a flexible, non-contractile tail of 116 × 8.5 nm possessing terminal fibres. The adsorption rate constant at 37 °C is 6.73 × 10(-9) ml min(-1). Latent and eclipse periods are 25 and 20 min, respectively, and the burst size is 35 progeny particles per cell after 35 min at 37 °C. Sequencing revealed a circularly permuted, 42 391 bp dsDNA genome containing 58 ORFs organized into four major transcriptional units. Comparisons with the genome sequences of other bacteriophages revealed a high level of nucleotide sequence identity and shared orthologous proteins with the Salmonella phages SETP3, SE2 and KS7 (SS3e) and the Escherichia phages K1G, K1H, K1ind1 and K1ind3.

Differences in DNA curvature-related sequence periodicity between prokaryotic chromosomes and phages, and relationship to chromosomal prophage content

Background

Periodic spacing of A-tracts (short runs of A or T) with the DNA helical period of ~10–11 bp is characteristic of intrinsically bent DNA. In eukaryotes, the DNA bending is related to chromatin structure and nucleosome positioning. However, the physiological role of strong sequence periodicity detected in many prokaryotic genomes is not clear.

Results

We developed measures of intensity and persistency of DNA curvature-related sequence periodicity and applied them to prokaryotic chromosomes and phages. The results indicate that strong periodic signals present in chromosomes are generally absent in phage genomes. Moreover, chromosomes containing prophages are less likely to possess a persistent periodic signal than chromosomes with no prophages.

Conclusions

Absence of DNA curvature-related sequence periodicity in phages could arise from constraints associated with DNA packaging in the viral capsid. Lack of prophages in chromosomes with persistent periodic signal suggests that the sequence periodicity and concomitant DNA curvature could play a role in protecting the chromosomes from integration of phage DNA.

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.

Prokaryote genome fluidity: toward a system approach of the mobilome

The importance of horizontal/lateral gene transfer (LGT) in shaping the genomes of prokaryotic organisms has been recognized in recent years as a result of analysis of the increasing number of available genome sequences. LGT is largely due to the transfer and recombination activities of mobile genetic elements (MGEs). Bacterial and archaeal genomes are mosaics of vertically and horizontally transmitted DNA segments. This generates reticulate relationships between members of the prokaryotic world that are better represented by networks than by "classical" phylogenetic trees. In this review we summarize the nature and activities of MGEs, and the problems that presently limit their analysis on a large scale. We propose routes to improve their annotation in the flow of genomic and metagenomic sequences that currently exist and those that become available. We describe network analysis of evolutionary relationships among some MGE categories and sketch out possible developments of this type of approach to get more insight into the role of the mobilome in bacterial adaptation and evolution.

Sequence Analysis of Inducible Prophage phIS3501 Integrated into the Haemolysin II Gene of Bacillus thuringiensis var israelensis ATCC35646

Diarrheic food poisoning by bacteria of the Bacillus cereus group is mostly due to several toxins encoded in the genomes. One of them, cytotoxin K, was recently identified as responsible for severe necrotic syndromes. Cytotoxin K is similar to a class of proteins encoded by genes usually annotated as haemolysin II (hlyII) in the majority of genomes of the B. cereus group. The partially sequenced genome of Bacillus thuringiensis var israelensis ATCC35646 contains several potentially induced prophages, one of them integrated into the hlyII gene. We determined the complete sequence and established the genomic organization of this prophage-designated phIS3501. During induction of excision of this prophage with mitomycin C, intact hlyII gene is formed, thus providing to cells a genetic ability to synthesize the active toxin. Therefore, this prophage, upon its excision, can be implicated in the regulation of synthesis of the active toxin and thus in the virulence of bacterial host. A generality of selection for such systems in bacterial pathogens is indicated by the similarity of this genetic arrangement to that of Staphylococcus aureus  β-haemolysin.

Reticulate classification of mosaic microbial genomes using NeAT website

The tree of life is the classical representation of the evolutionary relationships between existent species. A tree is appropriate to display the divergence of species through mutation, i.e., by vertical descent. However, lateral gene transfer (LGT) is excluded from such representations. When LGT contribution to genome evolution cannot be neglected (e.g., for prokaryotes and mobile genetic elements), the tree becomes misleading. Networks appear as an intuitive way to represent both vertical and horizontal relationships, while overlapping groups within such graphs are more suitable for their classification. Here, we describe a method to represent both vertical and horizontal relationships. We start with a set of genomes whose coded proteins have been grouped into families based on sequence similarity. Next, all pairs of genomes are compared, counting the number of proteins classified into the same family. From this comparison, we derive a weighted graph where genomes with a significant number of similar proteins are linked. Finally, we apply a two-step clustering of this graph to produce a classification where nodes can be assigned to multiple clusters. The procedure can be performed using the Network Analysis Tools (NeAT) website.

ICEberg: a web-based resource for integrative and conjugative elements found in Bacteria

ICEberg (http://db-mml.sjtu.edu.cn/ICEberg/) is an integrated database that provides comprehensive information about integrative and conjugative elements (ICEs) found in bacteria. ICEs are conjugative self-transmissible elements that can integrate into and excise from a host chromosome. An ICE contains three typical modules, integration and excision, conjugation, and regulation modules, that collectively promote vertical inheritance and periodic lateral gene flow. Many ICEs carry likely virulence determinants, antibiotic-resistant factors and/or genes coding for other beneficial traits. ICEberg offers a unique, highly organized, readily explorable archive of both predicted and experimentally supported ICE-relevant data. It currently contains details of 428 ICEs found in representatives of 124 bacterial species, and a collection of >400 directly related references. A broad range of similarity search, sequence alignment, genome context browser, phylogenetic and other functional analysis tools are readily accessible via ICEberg. We propose that ICEberg will facilitate efficient, multi-disciplinary and innovative exploration of bacterial ICEs and be of particular interest to researchers in the broad fields of prokaryotic evolution, pathogenesis, biotechnology and metabolism. The ICEberg database will be maintained, updated and improved regularly to ensure its ongoing maximum utility to the research community.

Phamerator: a bioinformatic tool for comparative bacteriophage genomics

BACKGROUND

Bacteriophage genomes have mosaic architectures and are replete with small open reading frames of unknown function, presenting challenges in their annotation, comparative analysis, and representation.

RESULTS

We describe here a bioinformatic tool, Phamerator, that assorts protein-coding genes into phamilies of related sequences using pairwise comparisons to generate a database of gene relationships. This database is used to generate genome maps of multiple phages that incorporate nucleotide and amino acid sequence relationships, as well as genes containing conserved domains. Phamerator also generates phamily circle representations of gene phamilies, facilitating analysis of the different evolutionary histories of individual genes that migrate through phage populations by horizontal genetic exchange.

CONCLUSIONS

Phamerator represents a useful tool for comparative genomic analysis and comparative representations of bacteriophage genomes.

Solving the Problem: Genome Annotation Standards before the Data Deluge

The promise of genome sequencing was that the vast undiscovered country would be mapped out by comparison of the multitude of sequences available and would aid researchers in deciphering the role of each gene in every organism. Researchers recognize that there is a need for high quality data. However, different annotation procedures, numerous databases, and a diminishing percentage of experimentally determined gene functions have resulted in a spectrum of annotation quality. NCBI in collaboration with sequencing centers, archival databases, and researchers, has developed the first international annotation standards, a fundamental step in ensuring that high quality complete prokaryotic genomes are available as gold standard references. Highlights include the development of annotation assessment tools, community acceptance of protein naming standards, comparison of annotation resources to provide consistent annotation, and improved tracking of the evidence used to generate a particular annotation. The development of a set of minimal standards, including the requirement for annotated complete prokaryotic genomes to contain a full set of ribosomal RNAs, transfer RNAs, and proteins encoding core conserved functions, is an historic milestone. The use of these standards in existing genomes and future submissions will increase the quality of databases, enabling researchers to make accurate biological discoveries.

Do phages efficiently shuttle transposable elements among prokaryotes?

The distribution, dynamics, and evolution of insertion sequences (IS), the most frequent class of prokaryotic transposable elements, are conditioned by their ability to horizontally transfer between cells. IS horizontal transfer (HT) requires shuttling by other mobile genetic elements. It is widely assumed in the literature that these vectors are phages and plasmids. By examining the relative abundance of IS in 454 plasmid and 446 phage genomes, we found that IS are very frequent in plasmids but, surprisingly, very rare in phages. Our results indicate that IS rarity in phages reflects very strong and efficient postinsertional purifying selection, mainly caused by a higher density of deleterious insertion sites in phages compared to plasmids. As they do not tolerate IS insertions, we conclude that phages may be rather poor vectors of IS HT in prokaryotes, in sharp contrast with the conventional view.

The human gut virome: inter-individual variation and dynamic response to diet

Immense populations of viruses are present in the human gut and other body sites. Understanding the role of these populations (the human "virome") in health and disease requires a much deeper understanding of their composition and dynamics in the face of environmental perturbation. Here, we investigate viromes from human subjects on a controlled feeding regimen. Longitudinal fecal samples were analyzed by metagenomic sequencing of DNA from virus-like particles (VLP) and total microbial communities. Assembly of 336 Mb of VLP sequence yielded 7175 contigs, many identifiable as complete or partial bacteriophage genomes. Contigs were rich in viral functions required in lytic and lysogenic growth, as well as unexpected functions such as viral CRISPR arrays and genes for antibiotic resistance. The largest source of variance among virome samples was interpersonal variation. Parallel deep-sequencing analysis of bacterial populations showed covaration of the virome with the larger microbiome. The dietary intervention was associated with a change in the virome community to a new state, in which individuals on the same diet converged. Thus these data provide an overview of the composition of the human gut virome and associate virome structure with diet.

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.

Analysis of antibiotic resistance regions in Gram-negative bacteria

Antibiotic resistance in Gram-negative bacteria is often due to the acquisition of resistance genes from a shared pool. In multiresistant isolates these genes, together with associated mobile elements, may be found in complex conglomerations on plasmids or on the chromosome. Analysis of available sequences reveals that these multiresistance regions (MRR) are modular, mosaic structures composed of different combinations of components from a limited set arranged in a limited number of ways. Components common to different MRR provide targets for homologous recombination, allowing these regions to evolve by combinatorial evolution, but our understanding of this process is far from complete. Advances in technology are leading to increasing amounts of sequence data, but currently available automated annotation methods usually focus on identifying ORFs and predicting protein function by homology. In MRR, where the genes are often well characterized, the challenge is to identify precisely which genes are present and to define the boundaries of complete and fragmented mobile elements. This review aims to summarize the types of mobile elements involved in multiresistance in Gram-negative bacteria and their associations with particular resistance genes, to describe common components of MRR and to illustrate methods for detailed analysis of these regions.

Genome-wide protein localization prediction strategies for gram negative bacteria

BACKGROUND

Genome-wide prediction of protein subcellular localization is an important type of evidence used for inferring protein function. While a variety of computational tools have been developed for this purpose, errors in the gene models and use of protein sorting signals that are not recognized by the more commonly accepted tools can diminish the accuracy of their output.

RESULTS

As part of an effort to manually curate the annotations of 19 strains of Shewanella, numerous insights were gained regarding the use of computational tools and proteomics data to predict protein localization. Identification of the suite of secretion systems present in each strain at the start of the process made it possible to tailor-fit the subsequent localization prediction strategies to each strain for improved accuracy. Comparisons of the computational predictions among orthologous proteins revealed inconsistencies in the computational outputs, which could often be resolved by adjusting the gene models or ortholog group memberships. While proteomic data was useful for verifying start site predictions and post-translational proteolytic cleavage, care was needed to distinguish cellular versus sample processing-mediated cleavage events. Searches for lipoprotein signal peptides revealed that neither TatP nor LipoP are designed for identification of lipoprotein substrates of the twin arginine translocation system and that the +2 rule for lipoprotein sorting does not apply to this Genus. Analysis of the relationships between domain occurrence and protein localization prediction enabled identification of numerous location-informative domains which could then be used to refine or increase confidence in location predictions. This collective knowledge was used to develop a general strategy for predicting protein localization that could be adapted to other organisms.

CONCLUSION

Improved localization prediction accuracy is not simply a matter of developing better computational algorithms. It also entails gathering key knowledge regarding the host architecture and translocation machinery and associated substrate recognition via experimentation and integration of diverse computational analyses from many proteins and, where possible, that are derived from different species within the same genus.

Gene flow, mobile genetic elements and the recruitment of antibiotic resistance genes into Gram-negative pathogens

Antibiotics were one of the great discoveries of the 20th century. However, resistance appeared even in the earliest years of the antibiotic era. Antibiotic resistance continues to become worse, despite the ever-increasing resources devoted to combat the problem. One of the most important factors in the development of resistance to antibiotics is the remarkable ability of bacteria to share genetic resources via Lateral Gene Transfer (LGT). LGT occurs on a global scale, such that in theory, any gene in any organism anywhere in the microbial biosphere might be mobilized and spread. With sufficiently strong selection, any gene may spread to a point where it establishes a global presence. From an antibiotic resistance perspective, this means that a resistance phenotype can appear in a diverse range of infections around the globe nearly simultaneously. We discuss the forces and agents that make this LGT possible and argue that the problem of resistance can ultimately only be managed by understanding the problem from a broad ecological and evolutionary perspective. We also argue that human activities are exacerbating the problem by increasing the tempo of LGT and bacterial evolution for many traits that are important to humans.

Comparative genomic analysis of the Streptococcus dysgalactiae species group: gene content, molecular adaptation, and promoter evolution

Comparative genomics of closely related bacterial species with different pathogenesis and host preference can provide a means of identifying the specifics of adaptive differences. Streptococcus dysgalactiae (SD) is comprised of two subspecies: S. dysgalactiae subsp. equisimilis is both a human commensal organism and a human pathogen, and S. dysgalactiae subsp. dysgalactiae is strictly an animal pathogen. Here, we present complete genome sequences for both taxa, with analyses involving other species of Streptococcus but focusing on adaptation in the SD species group. We found little evidence for enrichment in biochemical categories of genes carried by each SD strain, however, differences in the virulence gene repertoire were apparent. Some of the differences could be ascribed to prophage and integrative conjugative elements. We identified approximately 9% of the nonrecombinant core genome to be under positive selection, some of which involved known virulence factors in other bacteria. Analyses of proteomes by pooling data across genes, by biochemical category, clade, or branch, provided evidence for increased rates of evolution in several gene categories, as well as external branches of the tree. Promoters were primarily evolving under purifying selection but with certain categories of genes evolving faster. Many of these fast-evolving categories were the same as those associated with rapid evolution in proteins. Overall, these results suggest that adaptation to changing environments and new hosts in the SD species group has involved the acquisition of key virulence genes along with selection of orthologous protein-coding loci and operon promoters.

Genome annotation and intraviral interactome for the Streptococcus pneumoniae virulent phage Dp-1

Streptococcus pneumoniae causes several diseases, including pneumonia, septicemia, and meningitis. Phage Dp-1 is one of the very few isolated virulent S. pneumoniae bacteriophages, but only a partial characterization is currently available. Here, we confirmed that Dp-1 belongs to the family Siphoviridae. Then, we determined its complete genomic sequence of 56,506 bp. It encodes 72 open reading frames, of which 44 have been assigned a function. We have identified putative promoters, Rho-independent terminators, and several genomic clusters. We provide evidence that Dp-1 may be using a novel DNA replication system as well as redirecting host protein synthesis through queuosine-containing tRNAs. Liquid chromatography-mass spectrometry analysis of purified phage Dp-1 particles identified at least eight structural proteins. Finally, using comprehensive yeast two-hybrid screens, we identified 156 phage protein interactions, and this intraviral interactome was used to propose a structural model of Dp-1.

TADB: a web-based resource for Type 2 toxin-antitoxin loci in bacteria and archaea

TADB (http://bioinfo-mml.sjtu.edu.cn/TADB/) is an integrated database that provides comprehensive information about Type 2 toxin-antitoxin (TA) loci, genetic features that are richly distributed throughout bacterial and archaeal genomes. Two-gene and much less frequently three-gene Type 2 TA loci code for cognate partners that have been hypothesized or demonstrated to play key roles in stress response, bacterial physiology and stabilization of horizontally acquired genetic elements. TADB offers a unique compilation of both predicted and experimentally supported Type 2 TA loci-relevant data and currently contains 10,753 Type 2 TA gene pairs identified within 1240 prokaryotic genomes, and details of over 240 directly relevant scientific publications. A broad range of similarity search, sequence alignment, genome context browser and phylogenetic tools are readily accessible via TADB. We propose that TADB will facilitate efficient, multi-disciplinary and innovative exploration of the bacteria and archaea Type 2 TA space, better defining presently recognized TA-related phenomena and potentially even leading to yet-to-be envisaged frontiers. The TADB database, envisaged as a one-stop shop for Type 2 TA-related research, will be maintained, updated and improved regularly to ensure its ongoing maximum utility to the research community.

mGenomeSubtractor: a web-based tool for parallel in silico subtractive hybridization analysis of multiple bacterial genomes

mGenomeSubtractor performs an mpiBLAST-based comparison of reference bacterial genomes against multiple user-selected genomes for investigation of strain variable accessory regions. With parallel computing architecture, mGenomeSubtractor is able to run rapid BLAST searches of the segmented reference genome against multiple subject genomes at the DNA or amino acid level within a minute. In addition to comparison of protein coding sequences, the highly flexible sliding window-based genome fragmentation approach offered can be used to identify short unique sequences within or between genes. mGenomeSubtractor provides powerful schematic outputs for exploration of identified core and accessory regions, including searches against databases of mobile genetic elements, virulence factors or bacterial essential genes, examination of G+C content and binucleotide distribution bias, and integrated primer design tools. mGenomeSubtractor also allows for the ready definition of species-specific gene pools based on available genomes. Pan-genomic arrays can be easily developed using the efficient oligonucleotide design tool. This simple high-throughput in silico ‘subtractive hybridization’ analytical tool will support the rapidly escalating number of comparative bacterial genomics studies aimed at defining genomic biomarkers of evolutionary lineage, phenotype, pathotype, environmental adaptation and/or disease-association of diverse bacterial species. mGenomeSubtractor is freely available to all users without any login requirement at: http://bioinfo-mml.sjtu.edu.cn/mGS/.

Detection of novel recombinases in bacteriophage genomes unveils Rad52, Rad51 and Gp2.5 remote homologs

Homologous recombination is a key in contributing to bacteriophages genome repair, circularization and replication. No less than six kinds of recombinase genes have been reported so far in bacteriophage genomes, two (UvsX and Gp2.5) from virulent, and four (Sak, Redβ, Erf and Sak4) from temperate phages. Using profile–profile comparisons, structure-based modelling and gene-context analyses, we provide new views on the global landscape of recombinases in 465 bacteriophages. We show that Sak, Redβ and Erf belong to a common large superfamily adopting a shortcut Rad52-like fold. Remote homologs of Sak4 are predicted to adopt a shortcut Rad51/RecA fold and are discovered widespread among phage genomes. Unexpectedly, within temperate phages, gene-context analyses also pinpointed the presence of distant Gp2.5 homologs, believed to be restricted to virulent phages. All in all, three major superfamilies of phage recombinases emerged either related to Rad52-like, Rad51-like or Gp2.5-like proteins. For two newly detected recombinases belonging to the Sak4 and Gp2.5 families, we provide experimental evidence of their recombination activity in vivo. Temperate versus virulent lifestyle together with the importance of genome mosaicism is discussed in the light of these novel recombinases. Screening for these recombinases in genomes can be performed at http://biodev.extra.cea.fr/virfam.