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

Diverse Durham collection phages demonstrate complex BREX defense responses

Bacteriophages (phages) outnumber bacteria ten-to-one and cause infections at a rate of 1025 per second. The ability of phages to reduce bacterial populations makes them attractive alternative antibacterials for use in combating the rise in antimicrobial resistance. This effort may be hindered due to bacterial defenses such as Bacteriophage Exclusion (BREX) that have arisen from the constant evolutionary battle between bacteria and phages. For phages to be widely accepted as therapeutics in Western medicine, more must be understood about bacteria-phage interactions and the outcomes of bacterial phage defense. Here, we present the annotated genomes of 12 novel bacteriophage species isolated from water sources in Durham, UK, during undergraduate practical classes. The collection includes diverse species from across known phylogenetic groups. Comparative analyses of two novel phages from the collection suggest they may be founding members of a new genus. Using this Durham phage collection, we determined that particular BREX defense systems were likely to confer a varied degree of resistance against an invading phage. We concluded that the number of BREX target motifs encoded in the phage genome was not proportional to the degree of susceptibility. IMPORTANCE Bacteriophages have long been the source of tools for biotechnology that are in everyday use in molecular biology research laboratories worldwide. Phages make attractive new targets for the development of novel antimicrobials. While the number of phage genome depositions has increased in recent years, the expected bacteriophage diversity remains underrepresented. Here we demonstrate how undergraduates can contribute to the identification of novel phages and that a single City in England can provide ample phage diversity and the opportunity to find novel technologies. Moreover, we demonstrate that the interactions and intricacies of the interplay between bacterial phage defense systems such as Bacteriophage Exclusion (BREX) and phages are more complex than originally thought. Further work will be required in the field before the dynamic interactions between phages and bacterial defense systems are fully understood and integrated with novel phage therapies.

A widespread family of WYL-domain transcriptional regulators co-localizes with diverse phage defence systems and islands

Bacteria are under constant assault by bacteriophages and other mobile genetic elements. As a result, bacteria have evolved a multitude of systems that protect from attack. Genes encoding bacterial defence mechanisms can be clustered into 'defence islands', providing a potentially synergistic level of protection against a wider range of assailants. However, there is a comparative paucity of information on how expression of these defence systems is controlled. Here, we functionally characterize a transcriptional regulator, BrxR, encoded within a recently described phage defence island from a multidrug resistant plasmid of the emerging pathogen Escherichia fergusonii. Using a combination of reporters and electrophoretic mobility shift assays, we discovered that BrxR acts as a repressor. We present the structure of BrxR to 2.15 Å, the first structure of this family of transcription factors, and pinpoint a likely binding site for ligands within the WYL-domain. Bioinformatic analyses demonstrated that BrxR-family homologues are widespread amongst bacteria. About half (48%) of identified BrxR homologues were co-localized with a diverse array of known phage defence systems, either alone or clustered into defence islands. BrxR is a novel regulator that reveals a common mechanism for controlling the expression of the bacterial phage defence arsenal.

Bacteriophages in Biological Wastewater Treatment Systems: Occurrence, Characterization, and Function

Phage bacteria interactions can affect structure, dynamics, and function of microbial communities. In the context of biological wastewater treatment (BWT), the presence of phages can alter the efficiency of the treatment process and influence the quality of the treated effluent. The active role of phages in BWT has been demonstrated, but many questions remain unanswered regarding the diversity of phages in these engineered environments, the dynamics of infection, the determination of bacterial hosts, and the impact of their activity in full-scale processes. A deeper understanding of the phage ecology in BWT can lead the improvement of process monitoring and control, promote higher influent quality, and potentiate the use of phages as biocontrol agents. In this review, we highlight suitable methods for studying phages in wastewater adapted from other research fields, provide a critical overview on the current state of knowledge on the effect of phages on structure and function of BWT bacterial communities, and highlight gaps, opportunities, and priority questions to be addressed in future research.

The phage defence island of a multidrug resistant plasmid uses both BREX and type IV restriction for complementary protection from viruses

Bacteria have evolved a multitude of systems to prevent invasion by bacteriophages and other mobile genetic elements. Comparative genomics suggests that genes encoding bacterial defence mechanisms are often clustered in 'defence islands', providing a concerted level of protection against a wider range of attackers. However, there is a comparative paucity of information on functional interplay between multiple defence systems. Here, we have functionally characterised a defence island from a multidrug resistant plasmid of the emerging pathogen Escherichia fergusonii. Using a suite of thirty environmentally-isolated coliphages, we demonstrate multi-layered and robust phage protection provided by a plasmid-encoded defence island that expresses both a type I BREX system and the novel GmrSD-family type IV DNA modification-dependent restriction enzyme, BrxU. We present the structure of BrxU to 2.12 Å, the first structure of the GmrSD family of enzymes, and show that BrxU can utilise all common nucleotides and a wide selection of metals to cleave a range of modified DNAs. Additionally, BrxU undergoes a multi-step reaction cycle instigated by an unexpected ATP-dependent shift from an intertwined dimer to monomers. This direct evidence that bacterial defence islands can mediate complementary layers of phage protection enhances our understanding of the ever-expanding nature of phage-bacterial interactions.

Genome-Wide Identification and Analysis of Chromosomally Integrated Putative Prophages Associated with Clinical Klebsiella pneumoniae Strains

Klebsiella pneumoniae, an opportunistic pathogen found in the environment and human mucosal surfaces, is a leading cause of nosocomial infections. K. pneumoniae is now considered a global threat owing to the emergence of multidrug-resistant strains making its infections untreatable. In this study, 254 strains of K. pneumoniae were screened for the presence of prophages using the PHASTER tool. Very few strains lacked prophages (3.1%), while the remaining harboured both intact (811) and defective prophages (709). A subset of 42 unique strains of K. pneumoniae was chosen for further analysis. Our analysis revealed the presence of 110 complete prophages which were further classified as belonging to Myoviridae (67.3%), Siphoviridae (28.2%) and Podoviridae family (4.5%). An alignment of the 110 complete, prophage genome sequences clustered the prophages into 16 groups and 3 singletons. While none of the prophages encoded for virulence factors, 2 (1.8%) prophages were seen to encode for the antibiotic resistance-related genes. The CRISPR-Cas system was prevalent in 10 (23.8%) out of the 42 strains. Further analysis of the CRISPR spacers revealed 11.42% of the total spacers integrated in K. pneumoniae chromosome to match prophage protein sequences.

Biological species in the viral world

Due to their dependence on cellular organisms for metabolism and replication, viruses are typically named and assigned to species according to their genome structure and the original host that they infect. But because viruses often infect multiple hosts and the numbers of distinct lineages within a host can be vast, their delineation into species is often dictated by arbitrary sequence thresholds, which are highly inconsistent across lineages. Here we apply an approach to determine the boundaries of viral species based on the detection of gene flow within populations, thereby defining viral species according to the biological species concept (BSC). Despite the potential for gene transfer between highly divergent genomes, viruses, like the cellular organisms they infect, assort into reproductively isolated groups and can be organized into biological species. This approach revealed that BSC-defined viral species are often congruent with the taxonomic partitioning based on shared gene contents and host tropism, and that bacteriophages can similarly be classified in biological species. These results open the possibility to use a single, universal definition of species that is applicable across cellular and acellular lifeforms.

CRISPR-Cas-Mediated Phage Resistance Enhances Horizontal Gene Transfer by Transduction

A powerful contributor to prokaryotic evolution is horizontal gene transfer (HGT) through transformation, conjugation, and transduction, which can be advantageous, neutral, or detrimental to fitness. Bacteria and archaea control HGT and phage infection through CRISPR-Cas (clustered regularly interspaced short palindromic repeats-CRISPR-associated proteins) adaptive immunity. Although the benefits of resisting phage infection are evident, this can come at a cost of inhibiting the acquisition of other beneficial genes through HGT. Despite the ability of CRISPR-Cas to limit HGT through conjugation and transformation, its role in transduction is largely overlooked. Transduction is the phage-mediated transfer of bacterial DNA between cells and arguably has the greatest impact on HGT. We demonstrate that in Pectobacterium atrosepticum, CRISPR-Cas can inhibit the transduction of plasmids and chromosomal loci. In addition, we detected phage-mediated transfer of a large plant pathogenicity genomic island and show that CRISPR-Cas can inhibit its transduction. Despite these inhibitory effects of CRISPR-Cas on transduction, its more common role in phage resistance promotes rather than diminishes HGT via transduction by protecting bacteria from phage infection. This protective effect can also increase transduction of phage-sensitive members of mixed populations. CRISPR-Cas systems themselves display evidence of HGT, but little is known about their lateral dissemination between bacteria and whether transduction can contribute. We show that, through transduction, bacteria can acquire an entire chromosomal CRISPR-Cas system, including cas genes and phage-targeting spacers. We propose that the positive effect of CRISPR-Cas phage immunity on enhancing transduction surpasses the rarer cases where gene flow by transduction is restricted.IMPORTANCE The generation of genetic diversity through acquisition of DNA is a powerful contributor to microbial evolution and occurs through transformation, conjugation, and transduction. Of these, transduction, the phage-mediated transfer of bacterial DNA, is arguably the major route for genetic exchange. CRISPR-Cas adaptive immune systems control gene transfer by conjugation and transformation, but transduction has been mostly overlooked. Our results indicate that CRISPR-Cas can impede, but typically enhances the transduction of plasmids, chromosomal genes, and pathogenicity islands. By limiting wild-type phage replication, CRISPR-Cas immunity increases transduction in both phage-resistant and -sensitive members of mixed populations. Furthermore, we demonstrate mobilization of a chromosomal CRISPR-Cas system containing phage-targeting spacers by generalized transduction, which might partly account for the uneven distribution of these systems in nature. Overall, the ability of CRISPR-Cas to promote transduction reveals an unexpected impact of adaptive immunity on horizontal gene transfer, with broader implications for microbial evolution.

The Double-Stranded DNA Virosphere as a Modular Hierarchical Network of Gene Sharing

Virus genomes are prone to extensive gene loss, gain, and exchange and share no universal genes. Therefore, in a broad-scale study of virus evolution, gene and genome network analyses can complement traditional phylogenetics. We performed an exhaustive comparative analysis of the genomes of double-stranded DNA (dsDNA) viruses by using the bipartite network approach and found a robust hierarchical modularity in the dsDNA virosphere. Bipartite networks consist of two classes of nodes, with nodes in one class, in this case genomes, being connected via nodes of the second class, in this case genes. Such a network can be partitioned into modules that combine nodes from both classes. The bipartite network of dsDNA viruses includes 19 modules that form 5 major and 3 minor supermodules. Of these modules, 11 include tailed bacteriophages, reflecting the diversity of this largest group of viruses. The module analysis quantitatively validates and refines previously proposed nontrivial evolutionary relationships. An expansive supermodule combines the large and giant viruses of the putative order "Megavirales" with diverse moderate-sized viruses and related mobile elements. All viruses in this supermodule share a distinct morphogenetic tool kit with a double jelly roll major capsid protein. Herpesviruses and tailed bacteriophages comprise another supermodule, held together by a distinct set of morphogenetic proteins centered on the HK97-like major capsid protein. Together, these two supermodules cover the great majority of currently known dsDNA viruses. We formally identify a set of 14 viral hallmark genes that comprise the hubs of the network and account for most of the intermodule connections.

IMPORTANCE

Viruses and related mobile genetic elements are the dominant biological entities on earth, but their evolution is not sufficiently understood and their classification is not adequately developed. The key reason is the characteristic high rate of virus evolution that involves not only sequence change but also extensive gene loss, gain, and exchange. Therefore, in the study of virus evolution on a large scale, traditional phylogenetic approaches have limited applicability and have to be complemented by gene and genome network analyses. We applied state-of-the art methods of such analysis to reveal robust hierarchical modularity in the genomes of double-stranded DNA viruses. Some of the identified modules combine highly diverse viruses infecting bacteria, archaea, and eukaryotes, in support of previous hypotheses on direct evolutionary relationships between viruses from the three domains of cellular life. We formally identify a set of 14 viral hallmark genes that hold together the genomic network.

The complete genome sequence and analysis of vB_VorS-PVo5, a Vibrio phage infectious to the pathogenic bacterium Vibrio ordalii ATCC-33509

The bacterium Vibrio ordalii is best known as the causative agent of vibriosis outbreaks in fish and thus recognized for generating serious production losses in aquaculture systems. Here we report for the first time on the isolation and the genome sequencing of phage vB_VorS-PVo5, infectious to Vibrio ordalii ATCC 33509. The features as well as the complete genome sequence and annotation of the Vibrio phage are described; vB_VorS-PVo5 consists of a lineal double stranded DNA totaling ~ 80.6 Kb in length. Considering its ability to lyse Vibrio ordalii ATCC 33509, the phage is likely to gain importance in future aquaculture applications by controlling the pathogen and as such replacing antibiotics as the treatment of choice.

Structural plasticity and in vivo activity of Cas1 from the type I-F CRISPR-Cas system

CRISPR-Cas systems are adaptive immune systems in prokaryotes that provide protection against viruses and other foreign DNA. In the adaptation stage, foreign DNA is integrated into CRISPR (clustered regularly interspaced short palindromic repeat) arrays as new spacers. These spacers are used in the interference stage to guide effector CRISPR associated (Cas) protein(s) to target complementary foreign invading DNA. Cas1 is the integrase enzyme that is central to the catalysis of spacer integration. There are many diverse types of CRISPR-Cas systems, including type I-F systems, which are typified by a unique Cas1-Cas2-3 adaptation complex. In the present study we characterize the Cas1 protein of the potato phytopathogen Pectobacterium atrosepticum, an important model organism for understanding spacer acquisition in type I-F CRISPR-Cas systems. We demonstrate by mutagenesis that Cas1 is essential for adaptation in vivo and requires a conserved aspartic acid residue. By X-ray crystallography, we show that although P. atrosepticum Cas1 adopts a fold conserved among other Cas1 proteins, it possesses remarkable asymmetry as a result of structural plasticity. In particular, we resolve for the first time a flexible, asymmetric loop that may be unique to type I-F Cas1 proteins, and we discuss the implications of these structural features for DNA binding and enzymatic activity.

Computational clustering for viral reference proteomes

MOTIVATION

The enormous number of redundant sequenced genomes has hindered efforts to analyze and functionally annotate proteins. As the taxonomy of viruses is not uniformly defined, viral proteomes pose special challenges in this regard. Grouping viruses based on the similarity of their proteins at proteome scale can normalize against potential taxonomic nomenclature anomalies.

RESULTS

We present Viral Reference Proteomes (Viral RPs), which are computed from complete virus proteomes within UniProtKB. Viral RPs based on 95, 75, 55, 35 and 15% co-membership in proteome similarity based clusters are provided. Comparison of our computational Viral RPs with UniProt's curator-selected Reference Proteomes indicates that the two sets are consistent and complementary. Furthermore, each Viral RP represents a cluster of virus proteomes that was consistent with virus or host taxonomy. We provide BLASTP search and FTP download of Viral RP protein sequences, and a browser to facilitate the visualization of Viral RPs.

AVAILABILITY AND IMPLEMENTATION

http://proteininformationresource.org/rps/viruses/

CONTACT

chenc@udel.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Eco-Evolutionary Dynamics of Episomes among Ecologically Cohesive Bacterial Populations

Although plasmids and other episomes are recognized as key players in horizontal gene transfer among microbes, their diversity and dynamics among ecologically structured host populations in the wild remain poorly understood. Here, we show that natural populations of marine Vibrionaceae bacteria host large numbers of families of episomes, consisting of plasmids and a surprisingly high fraction of plasmid-like temperate phages. Episomes are unevenly distributed among host populations, and contrary to the notion that high-density communities in biofilms act as hot spots of gene transfer, we identified a strong bias for episomes to occur in free-living as opposed to particle-attached cells. Mapping of episomal families onto host phylogeny shows that, with the exception of all phage and a few plasmid families, most are of recent evolutionary origin and appear to have spread rapidly by horizontal transfer. Such high eco-evolutionary turnover is particularly surprising for plasmids that are, based on previously suggested categorization, putatively nontransmissible, indicating that this type of plasmid is indeed frequently transferred by currently unknown mechanisms. Finally, analysis of recent gene transfer among plasmids reveals a network of extensive exchange connecting nearly all episomes. Genes functioning in plasmid transfer and maintenance are frequently exchanged, suggesting that plasmids can be rapidly transformed from one category to another. The broad distribution of episomes among distantly related hosts and the observed promiscuous recombination patterns show how episomes can offer their hosts rapid assembly and dissemination of novel functions.

IMPORTANCE

Plasmids and other episomes are an integral part of bacterial biology in all environments, yet their study is heavily biased toward their role as vectors for antibiotic resistance genes. This study presents a comprehensive analysis of all episomes within several coexisting bacterial populations of Vibrionaceae from the coastal ocean and represents the largest-yet genomic survey of episomes from a single bacterial family. The host population framework allows analysis of the eco-evolutionary dynamics at unprecedented resolution, yielding several unexpected results. These include (i) discovery of novel, nonintegrative temperate phages, (ii) revision of a class of episomes, previously termed "nontransmissible," as highly transmissible, and (iii) surprisingly high evolutionary turnover of episomes, manifest as frequent birth, spread, and loss.

Whole genome comparison of a large collection of mycobacteriophages reveals a continuum of phage genetic diversity

The bacteriophage population is large, dynamic, ancient, and genetically diverse. Limited genomic information shows that phage genomes are mosaic, and the genetic architecture of phage populations remains ill-defined. To understand the population structure of phages infecting a single host strain, we isolated, sequenced, and compared 627 phages of Mycobacterium smegmatis. Their genetic diversity is considerable, and there are 28 distinct genomic types (clusters) with related nucleotide sequences. However, amino acid sequence comparisons show pervasive genomic mosaicism, and quantification of inter-cluster and intra-cluster relatedness reveals a continuum of genetic diversity, albeit with uneven representation of different phages. Furthermore, rarefaction analysis shows that the mycobacteriophage population is not closed, and there is a constant influx of genes from other sources. Phage isolation and analysis was performed by a large consortium of academic institutions, illustrating the substantial benefits of a disseminated, structured program involving large numbers of freshman undergraduates in scientific discovery.

DOI:http://dx.doi.org/10.7554/eLife.06416.001

Viruses are unable to replicate independently. To generate copies of itself, a virus must instead invade a target cell and commandeer that cell's replication machinery. Different viruses are able to invade different types of cell, and a group of viruses known as bacteriophages (or phages for short) replicate within bacteria. The enormous number and diversity of phages in the world means that they play an important role in virtually every ecosystem.

Despite their importance, relatively little is known about how different phage populations are related to each other and how they evolved. Many phages contain their genetic information in the form of strands of DNA. Using genetic sequencing to find out where and how different genes are encoded in the DNA can reveal information about how different viruses are related to each other. These relationships are particularly complicated in phages, as they can exchange genes with other viruses and microbes.

Previous studies comparing the genomes—the complete DNA sequence—of reasonably small numbers of phages that infect the Mycobacterium group of bacteria have found that the phages can be sorted into ‘clusters’ based on similarities in their genes and where these are encoded in their DNA. However, the number of phages investigated so far has been too small to conclude how different clusters are related. Are the clusters separate, or do they form a ‘continuum’ with different genes and DNA sequences shared between different clusters?

Here, Pope, Bowman, Russell et al. compare the individual genomes of 627 bacteriophages that infect the bacterial species Mycobacterium smegmatis. This is by far the largest number of phage genomes analyzed from a single host species. The large number of genomes analyzed allowed a much clearer understanding of the complexity and diversity of these phages to be obtained. The isolation, sequencing and analysis of the hundreds of M. smegmatis bacteriophage genomes was performed by an integrated research and education program, called the Science Education Alliance Phage Hunters Advancing Genomics and Evolutionary Science (SEA-PHAGES) program. This enabled thousands of undergraduate students from different institutions to contribute to the phage discovery and sequencing project, and co-author the report. SEA-PHAGES therefore shows that it is possible to successfully incorporate genuine scientific research into an undergraduate course, and that doing so can benefit both the students and researchers involved.

The results show that while the genomes could be categorized into 28 clusters, the genomes are not completely unrelated. Instead, a spread of diversity is seen, as genes and groups of genes are shared between different clusters. Pope, Bowman, Russell et al. further reveal that the phage population is in a constant state of change, and continuously acquires genes from other microorganisms and viruses.

DOI:http://dx.doi.org/10.7554/eLife.06416.002

Escherichia coli and Staphylococcus phages: effect of translation initiation efficiency on differential codon adaptation mediated by virulent and temperate lifestyles

Rapid biosynthesis is key to the success of bacteria and viruses. Highly expressed genes in bacteria exhibit a strong codon bias corresponding to the differential availability of tRNAs. However, a large clade of lambdoid coliphages exhibits relatively poor codon adaptation to the host translation machinery, in contrast to other coliphages that exhibit strong codon adaptation to the host. Three possible explanations were previously proposed but dismissed: (1) the phage-borne tRNA genes that reduce the dependence of phage translation on host tRNAs, (2) lack of time needed for evolving codon adaptation due to recent host switching, and (3) strong strand asymmetry with biased mutation disrupting codon adaptation. Here, we examined the possibility that phages with relatively poor codon adaptation have poor translation initiation which would weaken the selection on codon adaptation. We measured translation initiation by: (1) the strength and position of the Shine-Dalgarno (SD) sequence, and (2) the stability of the secondary structure of sequences flanking the SD and start codon known to affect accessibility of the SD sequence and start codon. Phage genes with strong codon adaptation had significantly stronger SD sequences than those with poor codon adaptation. The former also had significantly weaker secondary structure in sequences flanking the SD sequence and start codon than the latter. Thus, lambdoid phages do not exhibit strong codon adaptation because they have relatively inefficient translation initiation and would benefit little from increased elongation efficiency. We also provided evidence suggesting that phage lifestyle (virulent versus temperate) affected selection intensity on the efficiency of translation initiation and elongation.

Viral molecular mimicry circumvents abortive infection and suppresses bacterial suicide to make hosts permissive for replication

The global interplay between bacteria and bacteriophages has generated many macromolecules useful in biotechnology, through the co-evolutionary see-saw of bacterial defense and viral counter-attack measures. Bacteria can protect themselves using abortive infection systems, which induce altruistic suicide in an infected cell and therefore protect the clonal population at the expense of the infected individual. Our recent paper describes how bacteriophage ΦTE successfully subverted the activity of a plasmid-borne abortive infection system. ΦTE evolved mimics of the small RNA antitoxin that naturally inhibits the active toxin component of this anti-viral mechanism. These mutant phages further manipulated the behavior of the host population, through transduction of the plasmid encoding the abortive infection system. Transductants thereby became enslaved by the abortive infection system, committing suicide in response to infection by the original phage population. In effect, the new host was infected by an “addictive altruism,” to the advantage of the resistant bacteriophage.

Bacteriophage ϕMAM1, a viunalikevirus, is a broad-host-range, high-efficiency generalized transducer that infects environmental and clinical isolates of the enterobacterial genera Serratia and Kluyvera

Members of the enterobacterial genus Serratia are ecologically widespread, and some strains are opportunistic human pathogens. Bacteriophage ϕMAM1 was isolated on Serratia plymuthica A153, a biocontrol rhizosphere strain that produces the potently bioactive antifungal and anticancer haterumalide oocydin A. The ϕMAM1 phage is a generalized transducing phage that infects multiple environmental and clinical isolates of Serratia spp. and a rhizosphere strain of Kluyvera cryocrescens. Electron microscopy allowed classification of ϕMAM1 in the family Myoviridae. Bacteriophage ϕMAM1 is virulent, uses capsular polysaccharides as a receptor, and can transduce chromosomal markers at frequencies of up to 7 × 10(-6) transductants per PFU. We also demonstrated transduction of the complete 77-kb oocydin A gene cluster and heterogeneric transduction of a plasmid carrying a type III toxin-antitoxin system. These results support the notion of the potential ecological importance of transducing phages in the acquisition of genes by horizontal gene transfer. Phylogenetic analyses grouped ϕMAM1 within the ViI-like bacteriophages, and genomic analyses revealed that the major differences between ϕMAM1 and other ViI-like phages arise in a region encoding the host recognition determinants. Our results predict that the wider genus of ViI-like phages could be efficient transducing phages, and this possibility has obvious implications for the ecology of horizontal gene transfer, bacterial functional genomics, and synthetic biology.

sRNA antitoxins: more than one way to repress a toxin

Bacterial toxin-antitoxin loci consist of two genes: one encodes a potentially toxic protein, and the second, an antitoxin to repress its function or expression. The antitoxin can either be an RNA or a protein. For type I and type III loci, the antitoxins are RNAs; however, they have very different modes of action. Type I antitoxins repress toxin protein expression through interacting with the toxin mRNA, thereby targeting the mRNA for degradation or preventing its translation or both; type III antitoxins directly bind to the toxin protein, sequestering it. Along with these two very different modes of action for the antitoxin, there are differences in the functions of the toxin proteins and the mobility of these loci between species. Within this review, we discuss the major differences as to how the RNAs repress toxin activity, the potential consequences for utilizing different regulatory strategies, as well as the confirmed and potential biological roles for these loci across bacterial species.

Characterization and complete genome sequence analysis of novel bacteriophage IME-EFm1 infecting Enterococcus faecium

We isolated and characterized a novel virulent bacteriophage, IME-EFm1, specifically infecting multidrug-resistant Enterococcus faecium. IME-EFm1 is morphologically similar to members of the family Siphoviridae. It was found capable of lysing a wide range of our E. faecium collections, including two strains resistant to vancomycin. One-step growth tests revealed the host lysis activity of phage IME-EFm1, with a latent time of 30 min and a large burst size of 116 p.f.u. per cell. These biological characteristics suggested that IME-EFm1 has the potential to be used as a therapeutic agent. The complete genome of IME-EFm1 was 42 597 bp, and was linear, with terminally non-redundant dsDNA and a G+C content of 35.2 mol%. The termini of the phage genome were determined with next-generation sequencing and were further confirmed by nuclease digestion analysis. To our knowledge, this is the first report of a complete genome sequence of a bacteriophage infecting E. faecium. IME-EFm1 exhibited a low similarity to other phages in terms of genome organization and structural protein amino acid sequences. The coding region corresponded to 90.7 % of the genome; 70 putative ORFs were deduced and, of these, 29 could be functionally identified based on their homology to previously characterized proteins. A predicted metallo-β-lactamase gene was detected in the genome sequence. The identification of an antibiotic resistance gene emphasizes the necessity for complete genome sequencing of a phage to ensure it is free of any undesirable genes before use as a therapeutic agent against bacterial pathogens.

Differential codon adaptation between dsDNA and ssDNA phages in Escherichia coli

Because phages use their host translation machinery, their codon usage should evolve toward that of highly expressed host genes. We used two indices to measure codon adaptation of phages to their host, r_RSCU (the correlation in relative synonymous codon usage [RSCU] between phages and their host) and Codon Adaptation Index (CAI) computed with highly expressed host genes as the reference set (because phage translation depends on host translation machinery). These indices used for this purpose are appropriate only when hosts exhibit little mutation bias, so only phages parasitizing Escherichia coli were included in the analysis. For double-stranded DNA (dsDNA) phages, both r_RSCU and CAI decrease with increasing number of transfer RNA genes encoded by the phage genome. r_RSCU is greater for dsDNA phages than for single-stranded DNA (ssDNA) phages, and the low r_RSCU values are mainly due to poor concordance in RSCU values for Y-ending codons between ssDNA phages and the E. coli host, consistent with the predicted effect of C→T mutation bias in the ssDNA phages. Strong C→T mutation bias would improve codon adaptation in codon families (e.g., Gly) where U-ending codons are favored over C-ending codons (“U-friendly” codon families) by highly expressed host genes but decrease codon adaptation in other codon families where highly expressed host genes favor C-ending codons against U-ending codons (“U-hostile” codon families). It is remarkable that ssDNA phages with increasing C→T mutation bias also increased the usage of codons in the “U-friendly” codon families, thereby achieving CAI values almost as large as those of dsDNA phages. This represents a new type of codon adaptation.

SP10 infectivity is aborted after bacteriophage SP10 infection induces nonA transcription on the prophage SPβ region of the Bacillus subtilis genome

Bacteria have developed various strategies for phage resistance. Infection with phage induces the transcription of part of the phage resistance gene, but the regulatory mechanisms of such transcription remain largely unknown. The phage resistance gene nonA is located on the SPβ prophage region of the Bacillus subtilis Marburg strain genome. The nonA transcript was detected at the late stage of SP10 infection but is undetectable in noninfected cells. The nonA transcript was detected after the induction of the sigma factor Orf199-Orf200 (σ(Orf199-200)), when sigma factors encoded in the SP10 genome were expressed from a xylose-inducible plasmid. Thus, the SP10 sigma factor is an activator of a set of SP10 genes and nonA. The nonA gene encodes a 72-amino-acid protein with a transmembrane motif and has no significant homology with any protein in any database. NonA overexpression halted cell growth and reduced the efficiency of B. subtilis colony formation and respiration activity. In addition, SP10 virion protein synthesis was inhibited in the nonA(+) strain, and SP10 virion particles were scarce in it. These results indicate that NonA is a novel protein that can abort SP10 infection, and its transcription was regulated by SP10 sigma factor.

Phylogenomic network and comparative genomics reveal a diverged member of the ΦKZ-related group, marine vibrio phage ΦJM-2012

Bacteriophages are the largest reservoir of genetic diversity. Here we describe the novel phage ΦJM-2012. This natural isolate from marine Vibrio cyclitrophicus possesses very few gene contents relevant to other well-studied marine Vibrio phages. To better understand its evolutionary history, we built a mathematical model of pairwise relationships among 1,221 phage genomes, in which the genomes (nodes) are linked by edges representing the normalized number of shared orthologous protein families. This weighted network revealed that ΦJM-2012 was connected to only five members of the Pseudomonas ΦKZ-like phage family in an isolated network, strongly indicating that it belongs to this phage group. However, comparative genomic analyses highlighted an almost complete loss of colinearity with the ΦKZ-related genomes and little conservation of gene order, probably reflecting the action of distinct evolutionary forces on the genome of ΦJM-2012. In this phage, typical conserved core genes, including six RNA polymerase genes, were frequently displaced and the hyperplastic regions were rich in both unique genes and predicted unidirectional promoters with highly correlated orientations. Further, analysis of the ΦJM-2012 genome showed that segments of the conserved N-terminal parts of ΦKZ tail fiber paralogs exhibited evidence of combinatorial assortment, having switched transcriptional orientation, and there was recruitment and/or structural changes among phage endolysins and tail spike protein. Thus, this naturally occurring phage appears to have branched from a common ancestor of the ΦKZ-related groups, showing a distinct genomic architecture and unique genes that most likely reflect adaptation to its chosen host and environment.

The secret lives of mycobacteriophages

The study of mycobacteriophages provides insights into viral diversity and evolution, as well as the genetics and physiology of their pathogenic hosts. Genomic characterization of 80 mycobacteriophages reveals a high degree of genetic diversity and an especially rich reservoir of interesting genes. These include a vast number of genes of unknown function that do not match known database entries and many genes whose functions can be predicted but which are not typically found as components of phage genomes. Thus many mysteries surround these genomes, such as why the genes are there, what do they do, how are they expressed and regulated, how do they influence the physiology of the host bacterium, and what forces of evolution directed them to their genomic homes? Although the genetic diversity and novelty of these phages is full of intrigue, it is a godsend for the mycobacterial geneticist, presenting an abundantly rich toolbox that can be exploited to devise new and effective ways for understanding the genetics and physiology of human tuberculosis. As the number of sequenced genomes continues to grow, their mysteries continue to thicken, and the time has come to learn more about the secret lives of mycobacteriophages.

The role of prophage in plant-pathogenic bacteria

A diverse set of phage lineages is associated with the bacterial plant-pathogen genomes sequenced to date. Analysis of 37 genomes revealed 5,169 potential genes (approximately 4.3 Mbp) of phage origin, and at least 50% had no function assigned or are nonessential to phage biology. Some phytopathogens have transcriptionally active prophage genes under conditions that mimic plant infection, suggesting an association between plant disease and prophage transcriptional modulation. The role of prophages within genomes for cell biology varies. For pathogens such as Pectobacterium, Pseudomonas, Ralstonia, and Streptomyces, involvement of prophage in disease symptoms has been demonstrated. In Xylella and Xanthomonas, prophage activity is associated with genome rearrangements and strain differentiation. For other pathogens, prophage roles are yet to be established. This review integrates available information in a unique interface ( http://propnav.esalq.usp.br ) that may be assessed to improve research in prophage biology and its association with genome evolution and pathogenicity.

Genomic and proteomic analyses of the terminally redundant genome of the Pseudomonas aeruginosa phage PaP1: establishment of genus PaP1-like phages

We isolated and characterized a new Pseudomonas aeruginosa myovirus named PaP1. The morphology of this phage was visualized by electron microscopy and its genome sequence and ends were determined. Finally, genomic and proteomic analyses were performed. PaP1 has an icosahedral head with an apex diameter of 68–70 nm and a contractile tail with a length of 138–140 nm. The PaP1 genome is a linear dsDNA molecule containing 91,715 base pairs (bp) with a G+C content of 49.36% and 12 tRNA genes. A strategy to identify the genome ends of PaP1 was designed. The genome has a 1190 bp terminal redundancy. PaP1 has 157 open reading frames (ORFs). Of these, 143 proteins are homologs of known proteins, but only 38 could be functionally identified. Sodium dodecyl sulfate-polyacrylamide gel electrophoresis and high-performance liquid chromatography-mass spectrometry allowed identification of 12 ORFs as structural protein coding genes within the PaP1 genome. Comparative genomic analysis indicated that the Pseudomonas aeruginosa phage PaP1, JG004, PAK_P1 and vB_PaeM_C2-10_Ab1 share great similarity. Besides their similar biological characteristics, the phages contain 123 core genes and have very close phylogenetic relationships, which distinguish them from other known phage genera. We therefore propose that these four phages be classified as PaP1-like phages, a new phage genus of Myoviridae that infects Pseudomonas aeruginosa.

A novel virus genome discovered in an extreme environment suggests recombination between unrelated groups of RNA and DNA viruses

BACKGROUND

Viruses are known to be the most abundant organisms on earth, yet little is known about their collective origin and evolutionary history. With exceptionally high rates of genetic mutation and mosaicism, it is not currently possible to resolve deep evolutionary histories of the known major virus groups. Metagenomics offers a potential means of establishing a more comprehensive view of viral evolution as vast amounts of new sequence data becomes available for comparative analysis.

RESULTS

Bioinformatic analysis of viral metagenomic sequences derived from a hot, acidic lake revealed a circular, putatively single-stranded DNA virus encoding a major capsid protein similar to those found only in single-stranded RNA viruses. The presence and circular configuration of the complete virus genome was confirmed by inverse PCR amplification from native DNA extracted from lake sediment. The virus genome appears to be the result of a RNA-DNA recombination event between two ostensibly unrelated virus groups. Environmental sequence databases were examined for homologous genes arranged in similar configurations and three similar putative virus genomes from marine environments were identified. This result indicates the existence of a widespread but previously undetected group of viruses.

CONCLUSIONS

This unique viral genome carries implications for theories of virus emergence and evolution, as no mechanism for interviral RNA-DNA recombination has yet been identified, and only scant evidence exists that genetic exchange occurs between such distinct virus lineages.

REVIEWERS

This article was reviewed by EK, MK (nominated by PF) and AM. For the full reviews, please go to the Reviewers' comments section.

Genomics of bacterial and archaeal viruses: dynamics within the prokaryotic virosphere

Prokaryotes, bacteria and archaea, are the most abundant cellular organisms among those sharing the planet Earth with human beings (among others). However, numerous ecological studies have revealed that it is actually prokaryotic viruses that predominate on our planet and outnumber their hosts by at least an order of magnitude. An understanding of how this viral domain is organized and what are the mechanisms governing its evolution is therefore of great interest and importance. The vast majority of characterized prokaryotic viruses belong to the order Caudovirales, double-stranded DNA (dsDNA) bacteriophages with tails. Consequently, these viruses have been studied (and reviewed) extensively from both genomic and functional perspectives. However, albeit numerous, tailed phages represent only a minor fraction of the prokaryotic virus diversity. Therefore, the knowledge which has been generated for this viral system does not offer a comprehensive view of the prokaryotic virosphere. In this review, we discuss all families of bacterial and archaeal viruses that contain more than one characterized member and for which evolutionary conclusions can be attempted by use of comparative genomic analysis. We focus on the molecular mechanisms of their genome evolution as well as on the relationships between different viral groups and plasmids. It becomes clear that evolutionary mechanisms shaping the genomes of prokaryotic viruses vary between different families and depend on the type of the nucleic acid, characteristics of the virion structure, as well as the mode of the life cycle. We also point out that horizontal gene transfer is not equally prevalent in different virus families and is not uniformly unrestricted for diverse viral functions.

Large-scale analysis of plasmid relationships through gene-sharing networks

Plasmids are vessels of genetic exchange in microbial communities. They are known to transfer between different host organisms and acquire diverse genetic elements from chromosomes and/or other plasmids. Therefore, they constitute an important element in microbial evolution by rapidly disseminating various genetic properties among different communities. A paradigmatic example of this is the dissemination of antibiotic resistance (AR) genes that has resulted in the emergence of multiresistant pathogenic bacterial strains. To globally analyze the evolutionary dynamics of plasmids, we built a large graph in which 2,343 plasmids (nodes) are connected according to the proteins shared by each other. The analysis of this gene-sharing network revealed an overall coherence between network clustering and the phylogenetic classes of the corresponding microorganisms, likely resulting from genetic barriers to horizontal gene transfer between distant phylogenetic groups. Habitat was not a crucial factor in clustering as plasmids from organisms inhabiting different environments were often found embedded in the same cluster. Analyses of network metrics revealed a statistically significant correlation between plasmid mobility and their centrality within the network, providing support to the observation that mobile plasmids are particularly important in spreading genes in microbial communities. Finally, our study reveals an extensive (and previously undescribed) sharing of AR genes between Actinobacteria and Gammaproteobacteria, suggesting that the former might represent an important reservoir of AR genes for the latter.

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.

Identification of bacterial plasmids based on mobility and plasmid population biology

Plasmids contain a backbone of core genes that remains relatively stable for long evolutionary periods, making sense to speak about plasmid species. The identification and characterization of the core genes of a plasmid species has a special relevance in the study of its epidemiology and modes of transmission. Besides, this knowledge will help to unveil the main routes that genes, for example antibiotic resistance (AbR) genes, use to travel from environmental reservoirs to human pathogens. Global dissemination of multiple antibiotic resistances and virulence traits by plasmids is an increasing threat for the treatment of many bacterial infectious diseases. To follow the dissemination of virulence and AbR genes, we need to identify the causative plasmids and follow their path from reservoirs to pathogens. In this review, we discuss how the existing diversity in plasmid genetic structures gives rise to a large diversity in propagation strategies. We would like to propose that, using an identification methodology based on plasmid mobility types, we can follow the propagation routes of most plasmids in Gammaproteobacteria, as well as their cargo genes, in complex ecosystems. Once the dissemination routes are known, designing antidissemination drugs and testing their efficacy will become feasible. We discuss in this review how the existing diversity in plasmid genetic structures gives rise to a large diversity in propagation strategies. We would like to propose that, by using an identification methodology based on plasmid mobility types, we can follow the propagation routes of most plasmids in ?-proteobacteria, as well as their cargo genes, in complex ecosystems.

Evolutionarily conserved orthologous families in phages are relatively rare in their prokaryotic hosts

We have identified conserved orthologs in completely sequenced genomes of double-strand DNA phages and arranged them into evolutionary families (phage orthologous groups [POGs]). Using this resource to analyze the collection of known phage genomes, we find that most orthologs are unique in their genomes (having no diverged duplicates [paralogs]), and while many proteins contain multiple domains, the evolutionary recombination of these domains does not appear to be a major factor in evolution of these orthologous families. The number of POGs has been rapidly increasing over the past decade, the percentage of genes in phage genomes that have orthologs in other phages has also been increasing, and the percentage of unknown "ORFans" is decreasing as more proteins find homologs and establish a family. Other properties of phage genomes have remained relatively stable over time, most notably the high fraction of genes that are never or only rarely observed in their cellular hosts. This suggests that despite the renowned ability of phages to transduce cellular genes, these cellular "hitchhiker" genes do not dominate the phage genomic landscape, and a large fraction of the genes in phage genomes maintain an evolutionary trajectory that is distinct from that of the host genes.

A processed noncoding RNA regulates an altruistic bacterial antiviral system

The ≥ 10³⁰ bacteriophages on Earth relentlessly drive adaptive coevolution, forcing the generation of protective mechanisms in their bacterial hosts. One such bacterial phage-resistance system, ToxIN, consists of a protein toxin (ToxN) that is inhibited in vivo by a specific RNA antitoxin (ToxI); however, the mechanisms for this toxicity and inhibition have not been defined. Here we present the crystal structure of the ToxN-ToxI complex from Pectobacterium atrosepticum, determined to 2.75-Å resolution. ToxI is a 36-nucleotide noncoding RNA pseudoknot, and three ToxI monomers bind to three ToxN monomers to generate a trimeric ToxN-ToxI complex. Assembly of this complex is mediated entirely through extensive RNA-protein interactions. Furthermore, a 2'-3' cyclic phosphate at the 3' end of ToxI, and catalytic residues, identify ToxN as an endoRNase that processes ToxI from a repetitive precursor but is regulated by its own catalytic product.

Mycobacteriophages: genes and genomes

Viruses are powerful tools for investigating and manipulating their hosts, but the enormous size and amazing genetic diversity of the bacteriophage population have emerged as something of a surprise. In light of the evident importance of mycobacteria to human health--especially Mycobacterium tuberculosis, which causes tuberculosis--and the difficulties that have plagued their genetic manipulation, mycobacteriophages are especially appealing subjects for discovery, genomic characterization, and manipulation. With more than 70 complete genome sequences available, the mycobacteriophages have provided a wealth of information on the diversity of phages that infect a common bacterial host, revealed the pervasively mosaic nature of phage genome architectures, and identified a huge number of genes of unknown function. Mycobacteriophages have provided key tools for tuberculosis genetics, and new methods for simple construction of mycobacteriophage recombinants will facilitate postgenomic explorations into mycobacteriophage biology.

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.

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.

Characterization of a broad-host-range flagellum-dependent phage that mediates high-efficiency generalized transduction in, and between, Serratia and Pantoea

A phage (PhiOT8) isolated on Serratia sp. ATCC 39006 was shown to be flagellum-dependent, and to mediate generalized transduction with high efficiency (up to 10(-4) transductants per p.f.u.). PhiOT8 was shown to have a broad host range because it also infected a strain of Pantoea agglomerans isolated from the rhizosphere. Transduction of plasmid-borne antibiotic resistance between the two bacterial genera was demonstrated, consistent with purported ecological roles of phages in dissemination of genes between bacterial genera. Serratia sp. ATCC 39006 and P. agglomerans produce a number of interesting secondary metabolites that have potential applications in cancer therapy and biocontrol of fungal infections. PhiOT8 has utility as a powerful functional genomics tool in these bacteria.

Mutagenesis and functional characterization of the RNA and protein components of the toxIN abortive infection and toxin-antitoxin locus of Erwinia

Bacteria are constantly challenged by bacteriophage (phage) infection and have developed multiple adaptive resistance mechanisms. These mechanisms include the abortive infection systems, which promote "altruistic suicide" of an infected cell, protecting the clonal population. A cryptic plasmid of Erwinia carotovora subsp. atroseptica, pECA1039, has been shown to encode an abortive infection system. This highly effective system is active across multiple genera of gram-negative bacteria and against a spectrum of phages. Designated ToxIN, this two-component abortive infection system acts as a toxin-antitoxin module. ToxIN is the first member of a new type III class of protein-RNA toxin-antitoxin modules, of which there are multiple homologues cross-genera. We characterized in more detail the abortive infection phenotype of ToxIN using a suite of Erwinia phages and performed mutagenesis of the ToxI and ToxN components. We determined the minimal ToxI RNA sequence in the native operon that is both necessary and sufficient for abortive infection and to counteract the toxicity of ToxN. Furthermore, site-directed mutagenesis of ToxN revealed key conserved amino acids in this defining member of the new group of toxic proteins. The mechanism of phage activation of the ToxIN system was investigated and was shown to have no effect on the levels of the ToxN protein. Finally, evidence of negative autoregulation of the toxIN operon, a common feature of toxin-antitoxin systems, is presented. This work on the components of the ToxIN system suggests that there is very tight toxin regulation prior to suicide activation by incoming phage.

Exploitation of a new flagellatropic phage of Erwinia for positive selection of bacterial mutants attenuated in plant virulence: towards phage therapy

AIMS

To isolate and characterize novel bacteriophages for the phytopathogen, Erwinia carotovora ssp. atroseptica (Eca), and to isolate phage-resistant mutants attenuated in virulence.

METHODS AND RESULTS

A novel flagellatropic phage was isolated on the potato-rotting bacterial species, Eca, and characterized using electron microscopy and restriction analysis. The phage, named PhiAT1, has an icosahedral head and a long, contractile tail; it belongs to the Myoviridae family. Partial sequencing revealed the presence of genes with homology to those of coliphages T4, T7 and Mu. Phage-resistant transposon mutants of Eca were isolated and studied in vitro for a number of virulence-related phenotypes; only motility was found to be affected. In vivo tuber rotting assays showed that these mutants were attenuated in virulence, presumably because the infection is unable to spread from the initial site of inoculation.

CONCLUSIONS

The Eca flagellum can act as a receptor for PhiAT1 infection, and resistant mutants are enriched for motility and virulence defects.

SIGNIFICANCE AND IMPACT OF THE STUDY

PhiAT1 is the first reported flagellatropic phage found to infect Eca and has enabled further study of the virulence of this economically important phytopathogen.

Genomic islands are dynamic, ancient integrative elements in bacterial evolution

Acquisition of genomic islands plays a central part in bacterial evolution as a mechanism of diversification and adaptation. Genomic islands are non-self-mobilizing integrative and excisive elements that encode diverse functional characteristics but all contain a recombination module comprised of an integrase, associated attachment sites and, in some cases, a recombination directionality factor. Here, we discuss how a group of related genomic islands are evolutionarily ancient elements unrelated to plasmids, phages, integrons and integrative conjugative elements. In addition, we explore the diversity of genomic islands and their insertion sites among Gram-negative bacteria and discuss why they integrate at a limited number of tRNA genes.

The phage abortive infection system, ToxIN, functions as a protein-RNA toxin-antitoxin pair

Various mechanisms exist that enable bacteria to resist bacteriophage infection. Resistance strategies include the abortive infection (Abi) systems, which promote cell death and limit phage replication within a bacterial population. A highly effective 2-gene Abi system from the phytopathogen Erwinia carotovora subspecies atroseptica, designated ToxIN, is described. The ToxIN Abi system also functions as a toxin-antitoxin (TA) pair, with ToxN inhibiting bacterial growth and the tandemly repeated ToxI RNA antitoxin counteracting the toxicity. TA modules are currently divided into 2 classes, protein and RNA antisense. We provide evidence that ToxIN defines an entirely new TA class that functions via a novel protein-RNA mechanism, with analogous systems present in diverse bacteria. Despite the debated role of TA systems, we demonstrate that ToxIN provides viral resistance in a range of bacterial genera against multiple phages. This is the first demonstration of a novel mechanistic class of TA systems and of an Abi system functioning in different bacterial genera, both with implications for the dynamics of phage-bacterial interactions.

Bacteriophage genomics

The past three years have seen an escalation in the number of sequenced bacteriophage genomes with more than 500 now in the NCBI phage database, representing a more than threefold increase since 2005. These span at least 70 different bacterial hosts, with two-thirds of the sequenced genomes of phages representing only eight bacterial hosts. Three key features emerge from the comparative analysis of these genomes. First, they span a very high degree of genetic diversity, suggesting early evolutionary origins. Second, the genome architectures are mosaic, reflecting an unusually high degree of horizontal genetic exchange in their evolution. Third, phage genomes contain a very high proportion of novel genetic sequences of unknown function, and probably represent the largest reservoir of unexplored genes. With an estimated 10(31) bacterial and archael viruses in the biosphere, our view of the virosphere will draw into sharper focus as further bacteriophage genomes are characterized.

Improvisation in evolution of genes and genomes: whose structure is it anyway?

Significant progress has been made in recent years in a variety of seemingly unrelated fields such as sequencing, protein structure prediction, and high-throughput transcriptomics and metabolomics. At the same time, new microscopic models have been developed that made it possible to analyze the evolution of genes and genomes from first principles. The results from these efforts enable, for the first time, a comprehensive insight into the evolution of complex systems and organisms on all scales--from sequences to organisms and populations. Every newly sequenced genome uncovers new genes, families, and folds. Where do these new genes come from? How do gene duplication and subsequent divergence of sequence and structure affect the fitness of the organism? What role does regulation play in the evolution of proteins and folds? Emerging synergism between data and modeling provides first robust answers to these questions.

Unifying classical and molecular taxonomic classification: analysis of the Podoviridae using BLASTP-based tools

We defined phage genera by measuring genome relationships by the numbers of shared homologous/orthologous proteins. We used BLAST-based tools (CoreExtractor.vbs and CoreGenes) to analyze 55 fully sequenced bacteriophage genomes from the NCBI and EBI databases. This approach was first applied to the T7-related phages. Using a cut-off score of 40% homologous proteins, we identified three genera within the T7-related phages, redefined the phi29-related phages, and introduced five novel genera. The T7- and phi29-related phages were given subfamily status and named "Autographivirinae" and "Picovirinae", respectively. Our results confirm and refine the ICTV phage classification, enable elimination of errors in public databases, and provide a straightforward tool for the molecular classification of new phage genomes.