CTX prophages in classical biotype Vibrio cholerae: functional phage genes but dysfunctional phage genomes

CTXphi is a filamentous, lysogenic bacteriophage whose genome encodes cholera toxin, the primary virulence factor produced by Vibrio cholerae. CTX prophages in O1 El Tor and O139 strains of V. cholerae are found within arrays of genetically related elements integrated at a single locus within the V. cholerae large chromosome. The prophages of O1 El Tor and O139 strains generally yield infectious CTXphi. In contrast, O1 classical strains of V. cholerae do not produce CTXphi, although they produce cholera toxin and they contain CTX prophages integrated at two sites. We have identified the second site of CTX prophage integration in O1 classical strains and characterized the classical prophage arrays genetically and functionally. The genes of classical prophages encode functional forms of all of the proteins needed for production of CTXphi. Classical CTX prophages are present either as solitary prophages or as arrays of two truncated, fused prophages. RS1, a genetic element that is closely related to CTXphi and is often interspersed with CTX prophages in El Tor strains, was not detected in classical V. cholerae. Our model for CTXphi production predicts that the CTX prophage arrangements in classical strains will not yield extrachromosomal CTX DNA and thus will not yield virions, and our experimental results confirm this prediction. Thus, failure of O1 classical strains of V. cholerae to produce CTXphi is due to overall deficiencies in the structures of the arrays of classical prophages, rather than to mutations affecting individual CTX prophage genes.     

Identification of the J and K genes in the bacteriophage Mu genome sequence

Bacteriophage Mu was the first transposable phage to be discovered and still serves as the model for a large family of related transposable phages and prophages. The Mu genome sequence is known (NC-000929.1 GI:9633494), but not all of the genes have been assigned to the ORFs in the genome sequence. For this paper, we have sequenced an approximately 3-kb DNA region containing four predicted ORFs, Mup35-Mup38, from lysogens containing amber mutant prophages defective in either the J or the K gene. Amber mutations in prophages with J gene mutations mapped to the Mup36 ORF, and those in the K gene were found in Mup37, identifying the ORFs corresponding to these genes.     

Characterizing spontaneous induction of Stx encoding phages using a selectable reporter system

Shiga toxin (Stx) genes in Stx producing Escherichia coli (STEC) are encoded in prophages of the lambda family, such as H-19B. The subpopulation of STEC lysogens with induced prophages has been postulated to contribute significantly to Stx production and release. To study induced STEC, we developed a selectable in vivo expression technology, SIVET, a reporter system adapted from the RIVET system. The SIVET lysogen has a defective H-19B prophage encoding the TnpR resolvase gene downstream of the phage PR promoter and a cat gene with an inserted tet gene flanked by targets for the TnpR resolvase. Expression of resolvase results in excision of tet, restoring a functional cat gene; induced lysogens survive and are chloramphenicol resistant. Using SIVET we show that: (i) approximately 0.005% of the H-19B lysogens are spontaneously induced per generation during growth in LB. (ii) Variations in cellular physiology (e.g. RecA protein) rather than in levels of expressed repressor explain why members of a lysogen population are spontaneously induced. (iii) A greater fraction of lysogens with stx encoding prophages are induced compared to lysogens with non-Stx encoding prophages, suggesting increased sensitivity to inducing signal(s) has been selected in Stx encoding prophages. (iv) Only a small fraction of the lysogens in a culture spontaneously induce and when the lysogen carries two lambdoid prophages with different repressor/operators, 933W and H-19B, usually both prophages in the same cell are induced.     

Evolution along the parasitism-mutualism continuum determines the genetic repertoire of prophages

Integrated into their bacterial hosts’ genomes, prophage sequences exhibit a wide diversity of length and gene content, from highly degraded cryptic sequences to intact, functional prophages that retain a full complement of lytic-function genes. We apply three approaches—bioinformatics, analytical modelling and computational simulation—to understand the diverse gene content of prophages. In the bioinformatics work, we examine the distributions of over 50,000 annotated prophage genes identified in 1384 prophage sequences, comparing the gene repertoires of intact and incomplete prophages. These data indicate that genes involved in the replication, packaging, and release of phage particles have been preferentially lost in incomplete prophages, while tail fiber, transposase and integrase genes are significantly enriched. Consistent with these results, our mathematical and computational approaches predict that genes involved in phage lytic function are preferentially lost, resulting in shorter prophages that often retain genes that benefit the host. Informed by these models, we offer novel hypotheses for the enrichment of integrase and transposase genes in cryptic prophages. Overall, we demonstrate that functional and cryptic prophages represent a diversity of genetic sequences that evolve along a parasitism-mutualism continuum.     

Prophages and bacterial genomics: what have we learned so far?

Bacterial genome nucleotide sequences are being completed at a rapid and increasing rate. Integrated virus genomes (prophages) are common in such genomes. Fifty-one of the 82 such genomes published to date carry prophages, and these contain 230 recognizable putative prophages. Prophages can constitute as much as 10-20% of a bacterium's genome and are major contributors to differences between individuals within species. Many of these prophages appear to be defective and are in a state of mutational decay. Prophages, including defective ones, can contribute important biological properties to their bacterial hosts. Therefore, if we are to comprehend bacterial genomes fully, it is essential that we are able to recognize accurately and understand their prophages from nucleotide sequence analysis. Analysis of the evolution of prophages can shed light on the evolution of both bacteriophages and their hosts. Comparison of the Rac prophages in the sequenced genomes of three Escherichia coli strains and the Pnm prophages in two Neisseria meningitidis strains suggests that some prophages can lie in residence for very long times, perhaps millions of years, and that recombination events have occurred between related prophages that reside at different locations in a bacterium's genome. In addition, many genes in defective prophages remain functional, so a significant portion of the temperate bacteriophage gene pool resides in prophages.     

Expression of Pseudomonas aeruginosa phage transposons in Pseudomonas putida PpG1 cells. II. Zygotic induction--a necessary condition for the formation of defective lysogens

The transfer of hybrid plasmid RP4::PT (where PT is the genome of a transposable phage specific for Pseudomonas aeruginosa) into recipient cells of P. putida strain PpG1 occurs with the same frequency as into P. aeruginosa, the homologous host for PT. Approximately 1/3 of all PpG1 exconjugants carrying RP4 markers lost the capability to produce viable PT phage. In contrast, in a cross with homologous recipient P. aeruginosa all exconjugant clones contained nondefective prophages in the hybrid plasmids. Zygotic induction is an obligatory condition for detection of PpG1 exconjugants with defective phages. The defective prophages in RP4::PT hybrid plasmids have deletions of different size; the other carry mutations indistinguishable from point mutations in an essential phage gene. Some of deletions also cover plasmid genes. At least some of the defective prophages, including deleted ones, have arisen in the recipient cells of P. putida after transfer of the hybrid plasmid.     

Rapid membrane protein topology prediction

State-of-the-art methods for topology of α-helical membrane proteins are based on the use of time-consuming multiple sequence alignments obtained from PSI-BLAST or other sources. Here, we examine if it is possible to use the consensus of topology prediction methods that are based on single sequences to obtain a similar accuracy as the more accurate multiple sequence-based methods. Here, we show that TOPCONS-single performs better than any of the other topology prediction methods tested here, but ~6% worse than the best method that is utilizing multiple sequence alignments. AVAILABILITY AND IMPLEMENTATION: TOPCONS-single is available as a web server from http://single.topcons.net/ and is also included for local installation from the web site. In addition, consensus-based topology predictions for the entire international protein index (IPI) is available from the web server and will be updated at regular intervals.     

Vienna RNA secondary structure server

The Vienna RNA secondary structure server provides a web interface to the most frequently used functions of the Vienna RNA software package for the analysis of RNA secondary structures. It currently offers prediction of secondary structure from a single sequence, prediction of the consensus secondary structure for a set of aligned sequences and the design of sequences that will fold into a predefined structure. All three services can be accessed via the Vienna RNA web server at http://rna.tbi.univie.ac.at/.     

Prophages integrating into prophages: A mechanism to accumulate type III secretion effector genes and duplicate Shiga toxin-encoding prophages in Escherichia coli

Bacteriophages (or phages) play major roles in the evolution of bacterial pathogens via horizontal gene transfer. Multiple phages are often integrated in a host chromosome as prophages, not only carrying various novel virulence-related genetic determinants into host bacteria but also providing various possibilities for prophage-prophage interactions in bacterial cells. In particular, Escherichia coli strains such as Shiga toxin (Stx)-producing E. coli (STEC) and enteropathogenic E. coli (EPEC) strains have acquired more than 10 prophages (up to 21 prophages), many of which encode type III secretion system (T3SS) effector gene clusters. In these strains, some prophages are present at a single locus in tandem, which is usually interpreted as the integration of phages that use the same attachment (att) sequence. Here, we present phages integrating into T3SS effector gene cluster-associated loci in prophages, which are widely distributed in STEC and EPEC. Some of the phages integrated into prophages are Stx-encoding phages (Stx phages) and have induced the duplication of Stx phages in a single cell. The identified attB sequences in prophage genomes are apparently derived from host chromosomes. In addition, two or three different attB sequences are present in some prophages, which results in the generation of prophage clusters in various complex configurations. These phages integrating into prophages represent a medically and biologically important type of inter-phage interaction that promotes the accumulation of T3SS effector genes in STEC and EPEC, the duplication of Stx phages in STEC, and the conversion of EPEC to STEC and that may be distributed in other types of E. coli strains as well as other prophage-rich bacterial species.     

FrameD: A flexible program for quality check and gene prediction in prokaryotic genomes and noisy matured eukaryotic sequences

We describe FrameD, a program that predicts coding regions in prokaryotic and matured eukaryotic sequences. Initially targeted at gene prediction in bacterial GC rich genomes, the gene model used in FrameD also allows to predict genes in the presence of frameshifts and partially undetermined sequences which makes it also very suitable for gene prediction and frameshift correction in unfinished sequences such as EST and EST cluster sequences. Like recent eukaryotic gene prediction programs, FrameD also includes the ability to take into account protein similarity information both in its prediction and its graphical output. Its performances are evaluated on different bacterial genomes. The web site (http://genopole.toulouse.inra.fr/bioinfo/FrameD/FD) allows direct prediction, sequence correction and translation and the ability to learn new models for new organisms.     

AGenDA: gene prediction by comparative sequence analysis

Comparative sequence analysis is a powerful approach to identify functional elements in genomic sequences. Herein, we describe AGenDA (Alignment-based GENe Detection Algorithm), a novel method for gene prediction that is based on long-range alignment of syntenic regions in eukaryotic genome sequences. Local sequence homologies identified by the DIALIGN program are searched for conserved splice signals to define potential protein-coding exons; these candidate exons are then used to assemble complete gene structures. The performance of our method was tested on a set of 105 human-mouse sequence pairs. These test runs showed that sensitivity and specificity of AGenDA are comparable with the best gene- prediction program that is currently available. However, since our method is based on a completely different type of input information, it can detect genes that are not detectable by standard methods and vice versa. Thus, our approach seems to be a useful addition to existing gene-prediction programs. Availability: DIALIGN is available through the Bielefeld Bioinformatics Server (BiBiServ) at http://bibiserv.techfak.uni-bielefeld.de/dialign/ The gene-prediction program AGenDA described in this paper will be available through the BiBiServ or MIPS web server at http://mips.gsf.de.     

Prophage Tracer: precisely tracing prophages in prokaryotic genomes using overlapping split-read alignment

The life cycle of temperate phages includes a lysogenic cycle stage when the phage integrates into the host genome and becomes a prophage. However, the identification of prophages that are highly divergent from known phages remains challenging. In this study, by taking advantage of the lysis-lysogeny switch of temperate phages, we designed Prophage Tracer, a tool for recognizing active prophages in prokaryotic genomes using short-read sequencing data, independent of phage gene similarity searching. Prophage Tracer uses the criterion of overlapping split-read alignment to recognize discriminative reads that contain bacterial (attB) and phage (attP) att sites representing prophage excision signals. Performance testing showed that Prophage Tracer could predict known prophages with precise boundaries, as well as novel prophages. Two novel prophages, dsDNA and ssDNA, encoding highly divergent major capsid proteins, were identified in coral-associated bacteria. Prophage Tracer is a reliable data mining tool for the identification of novel temperate phages and mobile genetic elements. The code for the Prophage Tracer is publicly available at https://github.com/WangLab-SCSIO/Prophage_Tracer.     

Genome engineering reveals large dispensable regions in Bacillus subtilis

Bacterial genomes contain 250 to 500 essential genes, as suggested by single gene disruptions and theoretical considerations. If this view is correct, the remaining nonessential genes of an organism, such as Bacillus subtilis, have been acquired during evolution in its perpetually changing ecological niches. Notably, approximately 47% of the approximately 4,100 genes of B. subtilis belong to paralogous gene families in which several members have overlapping functions. Thus, essential gene functions will outnumber essential genes. To answer the question to what extent the most recently acquired DNA contributes to the life of B. subtilis under standard laboratory growth conditions, we initiated a "reconstruction" of the B. subtilis genome by removing prophages and AT-rich islands. Stepwise deletion of two prophages (SPbeta, PBSX), three prophage-like regions, and the largest operon of B. subtilis (pks) resulted in a genome reduction of 7.7% and elimination of 332 genes. The resulting strain was phenotypically characterized by metabolic flux analysis, proteomics, and specific assays for protein secretion, competence development, sporulation, and cell motility. We show that genome engineering is a feasible strategy for functional analysis of large gene clusters, and that removal of dispensable genomic regions may pave the way toward an optimized Bacillus cell factory.     

Distribution and function of prophage phiRv1 and phiRv2 among Mycobacterium tuberculosis complex

Mycobacterium tuberculosis complex (MTBC) is notorious for causing diseases, such as tuberculosis. Tuberculosis caused by M. tuberculosis remains a global public health concern. Two prophages, phiRv1 and phiRv2, can be found among most MTBC genomes. However, no precise functions have been assigned for the two prophages. In this paper, to find out the function of these two prophages, the distribution and function of phiRv1 and phiRv2 in MTBC genomes were analyzed from multiple omics data. We found that complex insertion, deletion, and reorganization appeared on the locus of two prophages in MTBC genomes; some genes of the two prophages can be translated and are functional from proteomic data; the expression of other prophage genes, such as Rv1577c, Rv2650c, Rv2652c, Rv2659c, and Rv2658c, can vary with environmental stresses and might enhance the fitness of MTBC. These data will facilitate our in-depth understanding of their function.