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

The genome sequence of Clostridium botulinum type C neurotoxin-converting phage and the molecular mechanisms of unstable lysogeny

Botulinum neurotoxins (BoNTXs) produced by Clostridium botulinum are among the most poisonous substances known. Of the seven types of BoNTXs, genes for type C1 and D toxins (BoNTX/C1 and D) are carried by bacteriophages. The gene for exoenzyme C3 also resides on these phages. Here, we present the complete genome sequence of c-st, a representative of BoNTX/C1-converting phages. The genome is a linear double-stranded DNA of 185,682 bp with 404-bp terminal direct repeats, the largest known temperate phage genome. We identified 198 potential protein-coding regions, including the genes for production of BoNTX/C1 and exoenzyme C3. Very exceptionally, as a viable bacteriophage, a number of insertion sequences were found on the c-st genome. By analyzing the molecular structure of the c-st genome in lysogens, we also found that it exists as a circular plasmid prophage. These features account for the unstable lysogeny of BoNTX phages, which has historically been called "pseudolysogeny." The PCR scanning analysis of other BoNTX/C1 and D phages based on the c-st sequence further revealed that BoNTX phages comprise a divergent phage family, probably generated by exchanging genomic segments among BoNTX phages and their relatives.

The players in a mutualistic symbiosis: insects, bacteria, viruses, and virulence genes

Aphids maintain mutualistic symbioses involving consortia of coinherited organisms. All possess a primary endosymbiont, Buchnera, which compensates for dietary deficiencies; many also contain secondary symbionts, such as Hamiltonella defensa, which confers defense against natural enemies. Genome sequences of uncultivable secondary symbionts have been refractory to analysis due to the difficulties of isolating adequate DNA samples. By amplifying DNA from hemolymph of infected pea aphids, we obtained a set of genomic sequences of H. defensa and an associated bacteriophage. H. defensa harbors two type III secretion systems, related to those that mediate host cell entry by enteric pathogens. The phage, called APSE-2, is a close relative of the previously sequenced APSE-1 but contains intact homologs of the gene encoding cytolethal distending toxin (cdtB), which interrupts the eukaryotic cell cycle and which is known from a variety of mammalian pathogens. The cdtB homolog is highly expressed, and its genomic position corresponds to that of a homolog of stx (encoding Shiga-toxin) within APSE-1. APSE-2 genomes were consistently abundant in infected pea aphids, and related phages were found in all tested isolates of H. defensa, from numerous insect species. Based on their ubiquity and abundance, these phages appear to be an obligate component of the H. defensa life cycle. We propose that, in these mutualistic symbionts, phage-borne toxin genes provide defense to the aphid host and are a basis for the observed protection against eukaryotic parasites.

Multiple large segment deletion method for Corynebacterium glutamicum

A precise and scarless genome excision method, employing the Cre/loxP system in concert with double-strand break (DSB)-stimulated intramolecular recombination was developed. The DSBs were mediated by the restriction endonuclease, I-SceI. It permitted multiple deletions of independent 14-, 43-, and 10-kb-long genomic regions on the Corynebacterium glutamicum genome. Accuracy of deletion was confirmed by the loss of marker genes, PCR, and sequencing of new genome joints. Eleven, 58, and 4 genes were predicted on the 14-, 43-, and 10-kb deleted regions, respectively. Although the resultant mutant lost a total of 67 kb encoding 73 genes, it still exhibited normal growth under standard laboratory conditions. Such a large segment deletion method in which multiple, successive deletions are possible is useful for genome engineering.

Identification of genes of VSH-1, a prophage-like gene transfer agent of Brachyspira hyodysenteriae

VSH-1 is a mitomycin C-inducible prophage of the anaerobic spirochete Brachyspira hyodysenteriae. Purified VSH-1 virions are noninfectious, contain random 7.5-kb fragments of the bacterial genome, and mediate generalized transduction of B. hyodysenteriae cells. In order to identify and sequence genes of this novel gene transfer agent (GTA), proteins associated either with VSH-1 capsids or with tails were purified by sodium dodecyl sulfate-polyacrylamide gel electrophoresis. The N-terminal amino acid sequences of 11 proteins were determined. Degenerate PCR primers were designed from the amino acid sequences and used to amplify several VSH-1 genes from B. hyodysenteriae strain B204 DNA. A lambda clone library of B. hyodysenteriae B204 DNA was subsequently screened by Southern hybridization methods and used to identify and sequence overlapping DNA inserts containing additional VSH-1 genes. VSH-1 genes spanned 16.3 kb of the B. hyodysenteriae chromosome and were flanked by bacterial genes. VSH-1 identified genes and unidentified, intervening open reading frames were consecutively organized in head (seven genes), tail (seven genes), and lysis (four genes) clusters in the same transcriptional direction. Putative lysis genes encoding endolysin (Lys) and holin proteins were identified from sequence and structural similarities of their translated protein products with GenBank bacteriophage proteins. Recombinant Lys protein hydrolyzed peptidoglycan purified from B. hyodysenteriae cells. The identified VSH-1 genes exceed the DNA capacity of VSH-1 virions and do not encode traditional bacteriophage early functions involved in DNA replication. These genome properties explain the noninfectious nature of VSH-1 virions and further confirm its resemblance to known prophage-like, GTAs of other bacterial species, such as the GTA from Rhodobacter capsulatus. The identification of VSH-1 genes will enable analysis of the regulation of this GTA and should facilitate investigations of VSH-1-like prophages from other Brachyspira species.

Diversity of phage integrases in Enterobacteriaceae: development of markers for environmental analysis of temperate phages

Viruses are the most abundant biological entities in aquatic systems. Temperate bacteriophages have enormous influences on microbial diversity, genetic exchange and bacterial population dynamics. However, development of molecular tools for their detection in the environment has been problematic. The integrase gene is used here as a molecular marker to analyse the diversity of temperate bacteriophages in a population of freshwater bacteria. Interrogation of the GenBank database revealed 32 non-cryptic enteric phage integrase sequences, leading to the development of a suite of 11 degenerate primer sets specific to the extant sequences elucidated. Application of these primer sets to enterobacterial isolates recovered from a freshwater pond and the temperate phages induced from them revealed a number of diverse integrase genes, including novel integrase-like sequences not represented in the databases. This highlights the potential of utilizing the integrase gene family as a marker for phage diversity.

A molecular target for suppression of the evolution of antibiotic resistance: inhibition of the Escherichia coli RecA protein by N(6)-(1-naphthyl)-ADP

We report that N(6)-(1-naphthyl)-ADP inhibits the Escherichia coli RecA protein in vitro. A novel rapid screen identified it as a potent inhibitor of RecA nucleoprotein filament formation, and further characterization established it as an ATP-competitive inhibitor of RecA-catalyzed ATP hydrolysis. This and other inhibitors of RecA activities represent a new approach for understanding the molecular targets and pathways involved in the evolution of antibiotic resistance in bacteria.

Complete genomic sequence of the temperate bacteriophage PhiAT3 isolated from Lactobacillus casei ATCC 393

The complete genomic sequence of a temperate bacteriophage PhiAT3 isolated from Lactobacillus (Lb.) casei ATCC 393 is reported. The phage consists of a linear DNA genome of 39,166 bp, an isometric head of 53 nm in diameter, and a flexible, noncontractile tail of approximately 200 nm in length. The number of potential open reading frames on the phage genome is 53. There are 15 unpaired nucleotides at both 5' ends of the PhiAT3 genome, indicating that the phage uses a cos-site for DNA packaging. The PhiAT3 genome was grouped into five distinct functional clusters: DNA packaging, morphogenesis, lysis, lysogenic/lytic switch, and replication. The amino acid sequences at the NH2-termini of some major proteins were determined. An in vivo integration assay for the PhiAT3 integrase (Int) protein in several lactobacilli was conducted by constructing an integration vector including PhiAT3 int and the attP (int-attP) region. It was found that PhiAT3 integrated at the tRNAArg gene locus of Lactobacillus rhamnosus HN 001, similar to that observed in its native host, Lb. casei ATCC 393.

Comparative genomics and evolution of the tailed-bacteriophages

The number of completely sequenced tailed-bacteriophage genomes that have been published increased to more than 125 last year. The comparison of these genomes has brought their highly mosaic nature into much sharper focus. Furthermore, reports of the complete sequences of about 150 bacterial genomes have shown that the many prophage and parts thereof that reside in these bacterial genomes must comprise a significant fraction of Earth's phage gene pool. These phage and prophage genomes are fertile ground for attempts to deduce the nature of viral evolutionary processes, and such analyses have made it clear that these phage have enjoyed a significant level of horizontal exchange of genetic information throughout their long histories. The strength of these evolutionary deductions rests largely on the extensive knowledge that has accumulated during intensive study into the molecular nature of the life cycles of a few 'model system' phages over the past half century. Recent molecular studies of phages other than these model system phages have made it clear that much remains to be learnt about the variety of lifestyle strategies utilized by the tailed-phage.

Gene order constrains adaptation in bacteriophage T7

The order of genes in the genome is commonly thought to have functional significance for gene regulation and fitness but has not heretofore been tested experimentally. We adapted a bacteriophage T7 variant harboring an ectopically positioned RNA polymerase gene to determine whether it could regain the fitness of the wild type. Two replicate lines maintained the starting gene order and showed only modest recovery of fitness, despite the accumulation of over a dozen mutations. In both lines, a mutation in the early terminator signal is responsible for the majority of the fitness recovery. In a third line, the phage evolved a new gene order, restoring the wild-type position of the RNA polymerase gene but also displacing several other genes to ectopic locations. Due to the recombination, the fitness of this replicate was the highest obtained but it falls short of the wild type adapted to the same growth conditions. The large benefits afforded by the terminator mutation and the recombination are explicable in terms of T7 biology, whereas several mutations with lesser benefits are not easily accounted for. These results support the premise that gene order is important to fitness and that wild-type fitness is not rapidly re-evolved in reorganized genomes.

Xylella and Xanthomonas Mobil'omics

The gamma-proteobacterium Xanthomonadales groups two closely related genera of plant pathogens, Xanthomonas and Xylella. Whole genome sequencing and comparative analyses disclosed a high degree of identity and co-linearity of the chromosome backbone between species and strains. Differences observed are usually clustered into genomic islands, most of which are delimited by genetic mobile elements. Focus is given in this paper to describe which groups of mobile elements are found and what is the relative contribution of these elements to Xanthomonas and Xylella genomes. Insertion sequence (IS) elements have invaded the Xanthomonas genome several times, whereas Xylella is rich in phage-related regions. Also, different plasmids are found inhabiting the bacterial cells studied here. Altogether, these results suggest that the integrative elements such as phages and transposable elements as well as the episomal plasmids are important drivers of the genome evolution of this important group of plant pathogens.

Mobile DNA in obligate intracellular bacteria

The small genomes of obligate intracellular bacteria are often presumed to be impervious to mobile DNA and the fluid genetic processes that drive diversification in free-living bacteria. Categorized by reductive evolution and streamlining, the genomes of some obligate intracellular bacteria manifest striking degrees of stability and gene synteny. However, recent findings from complete genome sequences of obligate intracellular species and their mobile genetic associates favour the abandonment of these wholesale terms for a more complex and tantalizing picture.

Comparative and evolutionary analysis of the bacterial homologous recombination systems

Homologous recombination is a housekeeping process involved in the maintenance of chromosome integrity and generation of genetic variability. Although detailed biochemical studies have described the mechanism of action of its components in model organisms, there is no recent extensive assessment of this knowledge, using comparative genomics and taking advantage of available experimental data on recombination. Using comparative genomics, we assessed the diversity of recombination processes among bacteria, and simulations suggest that we missed very few homologs. The work included the identification of orthologs and the analysis of their evolutionary history and genomic context. Some genes, for proteins such as RecA, the resolvases, and RecR, were found to be nearly ubiquitous, suggesting that the large majority of bacterial genomes are capable of homologous recombination. Yet many genomes show incomplete sets of presynaptic systems, with RecFOR being more frequent than RecBCD/AddAB. There is a significant pattern of co-occurrence between these systems and antirecombinant proteins such as the ones of mismatch repair and SbcB, but no significant association with nonhomologous end joining, which seems rare in bacteria. Surprisingly, a large number of genomes in which homologous recombination has been reported lack many of the enzymes involved in the presynaptic systems. The lack of obvious correlation between the presence of characterized presynaptic genes and experimental data on the frequency of recombination suggests the existence of still-unknown presynaptic mechanisms in bacteria. It also indicates that, at the moment, the assessment of the intrinsic stability or recombination isolation of bacteria in most cases cannot be inferred from the identification of known recombination proteins in the genomes.

Genomes evolve mostly by modifications involving large pieces of genetic material (DNA). Exchanges of chromosome pieces between different organisms as well as intragenomic movements of DNA regions are the result of a process named homologous recombination. The central actor of this process, the RecA protein, is amazingly conserved from bacteria to human. In addition to its role in the generation of genetic variability, homologous recombination is also the guardian of genome integrity, as it acts to repair DNA damage. RecA-catalyzed DNA exchange (synapse) is facilitated by the action of presynaptic enzymes and completed by postsynaptic enzymes (resolvases). In addition, some enzymes counteract RecA. Here, the researchers assess the diversity of recombination proteins among 117 different bacterial species. They find that resolvases are nearly as ubiquitous and as well conserved at the sequence level as RecA. This suggests that the large majority of bacterial genomes are capable of homologous recombination. Presynaptic systems are less ubiquitous, and there is no obvious correlation between their presence and experimental data on the frequency of recombination. However, there is a significant pattern of co-occurrence between these systems and antirecombinant proteins.

Sequence correlations between Cro recognition helices and cognate O(R) consensus half-sites suggest conserved rules of protein-DNA recognition

The O(R) regions from several lambdoid bacteriophages contain the three regulatory sites O(R)1, O(R)2 and O(R)3, to which the Cro and CI proteins can bind. These sites show imperfect dyad symmetry, have similar sequences, and generally lie on the same face of the DNA double helix. We have developed a computational method, which analyzes the O(R) regions of additional phages and predicts the location of these three sites. After tuning the method to predict known O(R) sites accurately, we used it to predict unknown sites, and ultimately compiled a database of 32 known and predicted O(R) binding site sets. We then identified sequences of the recognition helices (RH) for the cognate Cro proteins through manual inspection of multiple sequence alignments. Comparison of Cro RH and consensus O(R) half-site sequences revealed strong one-to-one correlations between two amino acids at each of three RH positions and two bases at each of three half-site positions (H1-->2, H3-->5 and H6-->6). In each of these three cases, one of the two amino acid/base-pairings corresponds to a contact observed in the crystal structure of a lambda Cro/consensus operator complex. The alternate amino acid/base combinations were rationalized using structural models. We suggest that the pairs of amino acid residues act as binary switches that efficiently modulate specificity for different consensus half-site variants during evolution. The observation of structurally reasonable amino acid-to-base correlations suggests that Cro proteins share some common rules of recognition despite their functional and structural diversity.

Complete genome sequence of phiHSIC, a pseudotemperate marine phage of Listonella pelagia

The genome for the marine pseudotemperate member of the Siphoviridae phiHSIC has been sequenced using a combination of linker amplification library construction, restriction digest library construction, and primer walking. phiHSIC enters into a pseudolysogenic relationship with its host, Listonella pelagia, characterized by sigmoidal growth curves producing >10(9) cells/ml and >10(11) phage/ml. The genome (37,966 bp; G+C content, 44%) contained 47 putative open reading frames (ORFs), 17 of which had significant BLASTP hits in GenBank, including a beta subunit of DNA polymerase III, a helicase, a helicase-like subunit of a resolvasome complex, a terminase, a tail tape measure protein, several phage-like structural proteins, and 1 ORF that may assist in host pathogenicity (an ADP ribosyltransferase). The genome was circularly permuted, with no physical ends detected by sequencing or restriction enzyme digestion analysis, and lacked a cos site. This evidence is consistent with a headful packaging mechanism similar to that of Salmonella phage P22 and Shigella phage Sf6. Because none of the phage-like ORFs were closely related to any existing phage sequences in GenBank (i.e., none more than 62% identical and most <25% identical at the amino acid level), phiHSIC is unique among phages that have been sequenced to date. These results further emphasize the need to sequence phages from the marine environment, perhaps the largest reservoir of untapped genetic information.

Large-scale engineering of the Corynebacterium glutamicum genome

The engineering of Corynebacterium glutamicum is important for enhanced production of biochemicals. To construct an improved C. glutamicum genome, we developed a precise genome excision method based on the Cre/loxP recombination system and successfully deleted 11 distinct genomic regions identified by comparative analysis of C. glutamicum genomes. Despite the loss of several predicted open reading frames, the mutant cells exhibited normal growth under standard laboratory conditions. With a total of 250 kb (7.5% of the genome), the 11 genomic regions were loaded with cryptic prophages, transposons, and genes of unknown function which were dispensable for cell growth, indicating recent horizontal acquisitions to the genome. This provides an interesting background for functional genomic studies and can be used in the improvement of cell traits.

Viral metagenomics

Viruses, most of which infect microorganisms, are the most abundant biological entities on the planet. Identifying and measuring the community dynamics of viruses in the environment is complicated because less than one percent of microbial hosts have been cultivated. Also, there is no single gene that is common to all viral genomes, so total uncultured viral diversity cannot be monitored using approaches analogous to ribosomal DNA profiling. Metagenomic analyses of uncultured viral communities circumvent these limitations and can provide insights into the composition and structure of environmental viral communities.

Analysis of the virus population present in equine faeces indicates the presence of hundreds of uncharacterized virus genomes

Virus DNA was isolated from horse faeces and cloned in a sequence-independent fashion. 268 clones were sequenced and 178140 nucleotides of sequence obtained. Statistical analysis suggests the library contains 17560 distinct clones derived from up to 233 different virus genomes. TBLASTX analysis showed that 32% of the clones had significant identity to GenBank entries. Of these 63% were viral; 20% bacterial; 7% archaeal; 6% eukarya; and 5% were related to mobile genetic elements. Fifty-two percent of the virus identities were with Siphoviridae; 26% unclassified phages; 17% Myoviridae; 4% Podoviridae; and one clone (2%) was a vertebrate Orthopoxvirus. Genes coding for predicted virus structural proteins, proteases, glycosidases and nucleic acid-binding proteins were common.

The linear double-stranded DNA of phage Bam35 enters lysogenic host cells, but the late phage functions are suppressed

Bam35, a temperate double-stranded DNA bacteriophage with a 15-kb linear genome, infects gram-positive Bacillus thuringiensis cells. Bam35 morphology and genome organization resemble those of PRD1, a lytic phage infecting gram-negative bacteria. Bam35 and PRD1 have an outer protein coat surrounding a membrane that encloses the viral DNA. We used electrochemical methods to investigate physiological changes of the lysogenic and nonlysogenic hosts during Bam35 DNA entry and host cell lysis. During viral DNA entry, there was an early temporal decrease of membrane voltage associated with K+ efflux that took place when either lysogenic or nonlysogenic hosts were infected. Approximately 40 min postinfection, a second strong K+ efflux was registered that was proposed to be associated with the insertion of holin molecules into the plasma membrane. This phenomenon occurred only when nonlysogenic cells were infected. Lysogenic hosts rarely were observed entering the lytic cycle as demonstrated by thin-section electron microscopy.

Marine phage genomics: what have we learned?

Marine phages are the most abundant and diverse form of life on the planet, and their genomes have been described as the largest untapped reservoir of genomic information. To date, however, the complete genome sequences of only 17 marine phage are known. Nevertheless, these genomes have revealed some interesting features, including the presence of photosynthetic genes in cyanophage and common patterns of genomic organization. Intriguing findings are also being made from studies of the uncultivated marine viral community genome ('metavirome'). The greatest challenge in interpreting the biology of these phages, and for making comparisons with their terrestrial counterparts, is the high proportion of unidentifiable open reading frames (approximately 60%). Future studies are likely to focus on sequencing more marine phage genomes from disparate hosts and diverse environments and on further basic studies of the biology of existing marine phages.

Effect of nutrient addition and environmental factors on prophage induction in natural populations of marine synechococcus species

A series of experiments were conducted with samples collected in both Tampa Bay and the Gulf of Mexico to assess the impact of nutrient addition on cyanophage induction in natural populations of Synechococcus sp. The samples were virus reduced to decrease the background level of cyanophage and then either left untreated or amended with nitrate, ammonium, urea, or phosphate. Replicate samples were treated with mitomycin C to stimulate cyanophage induction. In five of the nine total experiments performed, cyanophage induction was present in the non-nutrient-amended control samples. Stimulation of cyanophage induction in response to nutrient addition (phosphate) occurred in only one Tampa Bay sample. Nutrient additions caused a decrease in lytic (or control) phage production in three of three offshore stations, in one of three estuarine experiments, and in a lysogenic marine Synechococcus in culture. These results suggest that the process of cyanophage induction as an assay of Synechococcus lysogeny was not inorganically nutrient limited, at least in the samples examined. More importantly, it was observed that the level of cyanophage induction (cyanophage milliliter(-1)) was inversely correlated to Synechococcus and cyanophage abundance. Thus, the intensity of the prophage induction response is defined by ambient population size and cyanophage abundance. This corroborates prior observations that lysogeny in Synechococcus is favored during times of low host abundance.

Effects of diverse environmental conditions on {phi}LC3 prophage stability in Lactococcus lactis

The effects of various growth conditions on spontaneous phiLC3 prophage induction in Lactococcus lactis subsp. cremoris IMN-C1814 was analyzed with a half fraction of a 4(4) factorial experimental design. The four factors included in the study were nutrient availability, acidity, osmolarity, and temperature, each applied at four levels. These environmental factors are related to the fermentation processes in the dairy industry, in which bacteriophage attacks on sensitive starter strains are a constant threat to successful fermentation processes. The frequency of spontaneous phiLC3 induction was determined by quantitative analyses of restored DNA attachment sites (attB) on the bacterial chromosomes in a population of lysogenic cells. Statistical analysis revealed that all four environmental factors tested affected phiLC3 prophage stability and that the environmental factors were involved in interactions (interactions exist when the effect of one factor depends on the level of another factor). The spontaneous phiLC3 induction frequency varied from 0.08 to 1.76%. In general, the induction frequency remained at the same rate or decreased when level 1 to 3 of the four environmental factors was applied. At level 4, which generally gave the least favorable growth conditions, the induction frequency was either unchanged, decreased, or increased, depending on the type of stress. It appeared that the spontaneous induction frequency was independent of the growth behavior of the host. It was the environmental growth conditions that were the decisive factor in induction frequency.

Three Prochlorococcus cyanophage genomes: signature features and ecological interpretations

The oceanic cyanobacteria Prochlorococcus are globally important, ecologically diverse primary producers. It is thought that their viruses (phages) mediate population sizes and affect the evolutionary trajectories of their hosts. Here we present an analysis of genomes from three Prochlorococcus phages: a podovirus and two myoviruses. The morphology, overall genome features, and gene content of these phages suggest that they are quite similar to T7-like (P-SSP7) and T4-like (P-SSM2 and P-SSM4) phages. Using the existing phage taxonomic framework as a guideline, we examined genome sequences to establish “core” genes for each phage group. We found the podovirus contained 15 of 26 core T7-like genes and the two myoviruses contained 43 and 42 of 75 core T4-like genes. In addition to these core genes, each genome contains a significant number of “cyanobacterial” genes, i.e., genes with significant best BLAST hits to genes found in cyanobacteria. Some of these, we speculate, represent “signature” cyanophage genes. For example, all three phage genomes contain photosynthetic genes (psbA, hliP) that are thought to help maintain host photosynthetic activity during infection, as well as an aldolase family gene (talC) that could facilitate alternative routes of carbon metabolism during infection. The podovirus genome also contains an integrase gene (int) and other features that suggest it is capable of integrating into its host. If indeed it is, this would be unprecedented among cultured T7-like phages or marine cyanophages and would have significant evolutionary and ecological implications for phage and host. Further, both myoviruses contain phosphate-inducible genes (phoH and pstS) that are likely to be important for phage and host responses to phosphate stress, a commonly limiting nutrient in marine systems. Thus, these marine cyanophages appear to be variations of two well-known phages—T7 and T4—but contain genes that, if functional, reflect adaptations for infection of photosynthetic hosts in low-nutrient oceanic environments.

An analysis of the genome sequences of three phages capable of infecting marine unicellular cyanobacteria Prochlorococcus reveals they are genetically complex with intriguing adaptations related to their oceanic environment

The generalized transducing Salmonella bacteriophage ES18: complete genome sequence and DNA packaging strategy

The generalized transducing double-stranded DNA bacteriophage ES18 has an icosahedral head and a long noncontractile tail, and it infects both rough and smooth Salmonella enterica strains. We report here the complete 46,900-bp genome nucleotide sequence and provide an analysis of the sequence. Its 79 genes and their organization clearly show that ES18 is a member of the lambda-like (lambdoid) phage group; however, it contains a novel set of genes that program assembly of the virion head. Most of its integration-excision, immunity, Nin region, and lysis genes are nearly identical to those of the short-tailed Salmonella phage P22, while other early genes are nearly identical to Escherichia coli phages lambda and HK97, S. enterica phage ST64T, or a Shigella flexneri prophage. Some of the ES18 late genes are novel, while others are most closely related to phages HK97, lambda, or N15. Thus, the ES18 genome is mosaically related to other lambdoid phages, as is typical for all group members. Analysis of virion DNA showed that it is circularly permuted and about 10% terminally redundant and that initiation of DNA packaging series occurs across an approximately 1-kbp region rather than at a precise location on the genome. This supports a model in which ES18 terminase can move substantial distances along the DNA between recognition and cleavage of DNA destined to be packaged. Bioinformatic analysis of large terminase subunits shows that the different functional classes of phage-encoded terminases can usually be predicted from their amino acid sequence.

Microarray analysis of RpoS-mediated gene expression in Escherichia coli K-12

The alternative sigma factor RpoS controls the expression of many stationary-phase genes in Escherichia coli and other bacteria. Though the RpoS regulon is a large, conserved system that is critical for adaptation to nutrient deprivation and other stresses, it remains incompletely characterized. In this study, we have used oligonucleotide arrays to delineate the transcriptome that is controlled by RpoS during entry into stationary phase of cultures growing in rich medium. The expression of known RpoS-dependent genes was confirmed to be regulated by RpoS, thus validating the use of microarrays for expression analysis. The total number of positively regulated stationary-phase genes was found to be greater than 100. More than 45 new genes were identified as positively controlled by RpoS. Surprisingly, a similar number of genes were found to be negatively regulated by RpoS, and these included almost all genes required for flagellum biosynthesis, genes encoding enzymes of the TCA cycle, and a physically contiguous group of genes located in the Rac prophage region. Negative regulation by RpoS is thus much more extensive than has previously been recognized, and is likely to be an important contributing factor to the competitive growth advantage of rpoS mutants reported in previous studies.

The impact of bacteriophage genomics

The discovery of (bacterio)phages revolutionised microbiology and genetics, while phage research has been integral to answering some of the most fundamental biological questions of the twentieth century. The susceptibility of bacteria to bacteriophage attack can be undesirable in some cases, especially in the dairy industry, but can be desirable in others, for example, the use of bacteriophage therapy to eliminate pathogenic bacteria. The relative ease with which entire bacteriophage genome sequences can now be elucidated has had a profound impact on the study of these bacterial parasites.

Complete genomic sequence of bacteriophage B3, a Mu-like phage of Pseudomonas aeruginosa

Bacteriophage B3 is a transposable phage of Pseudomonas aeruginosa. In this report, we present the complete DNA sequence and annotation of the B3 genome. DNA sequence analysis revealed that the B3 genome is 38,439 bp long with a G+C content of 63.3%. The genome contains 59 proposed open reading frames (ORFs) organized into at least three operons. Of these ORFs, the predicted proteins from 41 ORFs (68%) display significant similarity to other phage or bacterial proteins. Many of the predicted B3 proteins are homologous to those encoded by the early genes and head genes of Mu and Mu-like prophages found in sequenced bacterial genomes. Only two of the predicted B3 tail proteins are homologous to other well-characterized phage tail proteins; however, several Mu-like prophages and transposable phage D3112 encode approximately 10 highly similar proteins in their predicted tail gene regions. Comparison of the B3 genomic organization with that of Mu revealed evidence of multiple genetic rearrangements, the most notable being the inversion of the proposed B3 immunity/early gene region, the loss of Mu-like tail genes, and an extreme leftward shift of the B3 DNA modification gene cluster. These differences illustrate and support the widely held view that tailed phages are genetic mosaics arising by the exchange of functional modules within a diverse genetic pool.

Movement of viruses between biomes

Viruses are abundant in all known ecosystems. In the present study, we tested the possibility that viruses from one biome can successfully propagate in another. Viral concentrates were prepared from different near-shore marine sites, lake water, marine sediments, and soil. The concentrates were added to microcosms containing dissolved organic matter as a food source (after filtration to allow 100-kDa particles to pass through) and a 3% (vol/vol) microbial inoculum from a marine water sample (after filtration through a 0.45-microm-pore-size filter). Virus-like particle abundances were then monitored using direct counting. Viral populations from lake water, marine sediments, and soil were able to replicate when they were incubated with the marine microbes, showing that viruses can move between different ecosystems and propagate. These results imply that viruses can laterally transfer DNA between microbes in different biomes.

Complete genome sequence of the industrial bacterium Bacillus licheniformis and comparisons with closely related Bacillus species

The complete sequence of the Bacillus licheniformis ATCC 14580 genome was determined, revealing 4,208 predicted protein-coding genes, 7 rRNA operons and 72 tRNA genes.

Background

Bacillus licheniformis is a Gram-positive, spore-forming soil bacterium that is used in the biotechnology industry to manufacture enzymes, antibiotics, biochemicals and consumer products. This species is closely related to the well studied model organism Bacillus subtilis, and produces an assortment of extracellular enzymes that may contribute to nutrient cycling in nature.

Results

We determined the complete nucleotide sequence of the B. licheniformis ATCC 14580 genome which comprises a circular chromosome of 4,222,336 base-pairs (bp) containing 4,208 predicted protein-coding genes with an average size of 873 bp, seven rRNA operons, and 72 tRNA genes. The B. licheniformis chromosome contains large regions that are colinear with the genomes of B. subtilis and Bacillus halodurans, and approximately 80% of the predicted B. licheniformis coding sequences have B. subtilis orthologs.

Conclusions

Despite the unmistakable organizational similarities between the B. licheniformis and B. subtilis genomes, there are notable differences in the numbers and locations of prophages, transposable elements and a number of extracellular enzymes and secondary metabolic pathway operons that distinguish these species. Differences include a region of more than 80 kilobases (kb) that comprises a cluster of polyketide synthase genes and a second operon of 38 kb encoding plipastatin synthase enzymes that are absent in the B. licheniformis genome. The availability of a completed genome sequence for B. licheniformis should facilitate the design and construction of improved industrial strains and allow for comparative genomics and evolutionary studies within this group of Bacillaceae.

Phages and the evolution of bacterial pathogens: from genomic rearrangements to lysogenic conversion

Comparative genomics demonstrated that the chromosomes from bacteria and their viruses (bacteriophages) are coevolving. This process is most evident for bacterial pathogens where the majority contain prophages or phage remnants integrated into the bacterial DNA. Many prophages from bacterial pathogens encode virulence factors. Two situations can be distinguished: Vibrio cholerae, Shiga toxin-producing Escherichia coli, Corynebacterium diphtheriae, and Clostridium botulinum depend on a specific prophage-encoded toxin for causing a specific disease, whereas Staphylococcus aureus, Streptococcus pyogenes, and Salmonella enterica serovar Typhimurium harbor a multitude of prophages and each phage-encoded virulence or fitness factor makes an incremental contribution to the fitness of the lysogen. These prophages behave like "swarms" of related prophages. Prophage diversification seems to be fueled by the frequent transfer of phage material by recombination with superinfecting phages, resident prophages, or occasional acquisition of other mobile DNA elements or bacterial chromosomal genes. Prophages also contribute to the diversification of the bacterial genome architecture. In many cases, they actually represent a large fraction of the strain-specific DNA sequences. In addition, they can serve as anchoring points for genome inversions. The current review presents the available genomics and biological data on prophages from bacterial pathogens in an evolutionary framework.

Cloning Serratia entomophila antifeeding genes--a putative defective prophage active against the grass grub Costelytra zealandica

Serratia entomophila and Serratia proteamaculans (Enterobacteriaceae) cause amber disease in the grass grub Costelytra zealandica (Coleoptera: Scarabaeidae), an important pasture pest in New Zealand. Larval disease symptoms include cessation of feeding, clearance of the gut, amber coloration, and eventual death. A 155-kb plasmid, pADAP, carries the genes sepA, sepB, and sepC, which are essential for production of amber disease symptoms. Transposon insertions in any of the sep genes in pADAP abolish gut clearance but not cessation of feeding, indicating the presence of an antifeeding gene(s) elsewhere on pADAP. Based on deletion analysis of pADAP and subsequent sequence data, a 47-kb clone was constructed, which when placed in either an Escherichia coli or a Serratia background exerted strong antifeeding activity and often led to rapid death of the infected grass grub larvae. Sequence data show that the antifeeding component is part of a large gene cluster that may form a defective prophage and that six potential members of this prophage are present in Photorhabdus luminescens subsp. laumondii TTO1, a species which also has sep gene homologues.

Transfer of photosynthesis genes to and from Prochlorococcus viruses

Comparative genomics gives us a new window into phage-host interactions and their evolutionary implications. Here we report the presence of genes central to oxygenic photosynthesis in the genomes of three phages from two viral families (Myoviridae and Podoviridae) that infect the marine cyanobacterium Prochlorococcus. The genes that encode the photosystem II core reaction center protein D1 (psbA), and a high-light-inducible protein (HLIP) (hli) are present in all three genomes. Both myoviruses contain additional hli gene types, and one of them encodes the second photosystem II core reaction center protein D2 (psbD), whereas the other encodes the photosynthetic electron transport proteins plastocyanin (petE) and ferredoxin (petF). These uninterrupted, full-length genes are conserved in their amino acid sequence, suggesting that they encode functional proteins that may help maintain photosynthetic activity during infection. Phylogenetic analyses show that phage D1, D2, and HLIP proteins cluster with those from Prochlorococcus, indicating that they are of cyanobacterial origin. Their distribution among several Prochlorococcus clades further suggests that the genes encoding these proteins were transferred from host to phage multiple times. Phage HLIPs cluster with multicopy types found exclusively in Prochlorocococus, suggesting that phage may be mediating the expansion of the hli gene family by transferring these genes back to their hosts after a period of evolution in the phage. These gene transfers are likely to play a role in the fitness landscape of hosts and phages in the surface oceans.

The chromosome of Shigella flexneri bacteriophage Sf6: complete nucleotide sequence, genetic mosaicism, and DNA packaging

Shigella flexneri temperate bacteriophage Sf6 is of interest in part because its prophage expresses the oac gene that alters the antigenic properties of the surface O-antigen polysaccharide of its host bacterium. We have determined the complete sequence of its 39,044 bp genome. The sequence shows that Sf6 is a member of the canonical lambdoid phage group, and like other phages of this type has a highly mosaic genome. It has chromosomal regions that encode proteins >80% identical with at least 15 different previously characterized lambdoid phages and prophages, but 43% of the genome, including the virion assembly genes, is homologous to the genome of one phage, HK620. An analysis of the nucleotide differences between Sf6 and HK620 indicates that even these similar regions are highly mosaic. This mosaicism suggests ways in which the virion structural proteins might interact with each other. The Sf6 early operons are arranged like a typical lambdoid phage, with "boundary sequences" often found between functional modules in the "metabolic" genome domain. By virtue of high degree of similarity in the encoding genes and their DNA target sites, we predict that the integrase, early transcription anti-terminator, CI and Cro repressors, and CII protein of Sf6 have DNA binding specificities very similar to the homologous proteins encoded by phages HK620, lambda, 434 and P22, respectively. The late operon contains two tRNA genes. The Sf6 terminase genes are unusual. Analysis of in vivo initiation of the DNA packaging series showed that the Sf6 apparatus that recognizes DNA for packaging appears to cleave DNA for initiation of packaging series at many sites within a large region of about 1800 bp that includes a possible pac site. This is unlike previously characterized phage packaging mechanisms.

Burkholderia cenocepacia phage BcepMu and a family of Mu-like phages encoding potential pathogenesis factors

We have isolated BcepMu, a Mu-like bacteriophage whose host range includes human pathogenic Burkholderia cenocepacia (formally B. cepacia genomovar III) isolates, and determined its complete 36748 bp genomic sequence. Like enteric bacteriophage Mu, the BcepMu genomic DNA is flanked by variable host sequences, a result of transposon-mediated replication. The BcepMu genome encodes 53 proteins, including capsid assembly components related to those of Mu, and tail sheath and tube proteins related to those of bacteriophage P2. Seventeen of the BcepMu genes were demonstrated to encode homotypic interacting domains by using a cI fusion system. Most BcepMu genes have close homologs to prophage elements present in the two published Salmonella typhi genomes, and in the database sequences of Photorhabdus luminescens, and Chromobacterium violaceum. These prophage elements, designated SalMu, PhotoMu and ChromoMu, respectively, are collinear with BcepMu through nearly their entire lengths and show only limited mosaicism, despite the divergent characters of their hosts. The BcepMu family of Mu-like phages has a number of notable differences from Mu. Most significantly, the critical left end region of BcepMu is inverted with respect to Mu, and the BcepMu family of transposases is clearly of a distinct lineage with different molecular requirements at the transposon ends. Interestingly, a survey of 33 B.cepacia complex strains indicated that the BcepMu prophage is widespread in human pathogenic B.cenocepacia ET12 lineage isolates, but not in isolates from the PHDC or Midwest lineages. Identified members of the BcepMu family all contain a gene possibly involved in bacterial pathogenicity, a homolog of the type-two-secretion component exeA, but only BcepMu also carries a lipopolysaccharide modification acyltransferase which may also contribute a pathogenicity factor.

Mu-like prophage strong gyrase site sequences: analysis of properties required for promoting efficient mu DNA replication

The bacteriophage Mu genome contains a centrally located strong gyrase site (SGS) that is required for efficient prophage replication. To aid in studying the unusual properties of the SGS, we sought other gyrase sites that might be able to substitute for the SGS in Mu replication. Five candidate sites were obtained by PCR from Mu-like prophage sequences present in Escherichia coli O157:H7 Sakai, Haemophilus influenzae Rd, Salmonella enterica serovar Typhi CT18, and two strains of Neisseria meningitidis. Each of the sites was used to replace the natural Mu SGS to form recombinant prophages, and the effects on Mu replication and host lysis were determined. The site from the E. coli prophage supported markedly enhanced replication and host lysis over that observed with a Mu derivative lacking the SGS, those from the N. meningitidis prophages allowed a small enhancement, and the sites from the Haemophilus and Salmonella prophages gave none. Each of the candidate sites was cleaved specifically by E. coli DNA gyrase both in vitro and in vivo. Supercoiling assays performed in vitro, with the five sites or the Mu SGS individually cloned into a pUC19 reporter plasmid, showed that the Mu SGS and the E. coli or N. meningitidis sequences allowed an enhancement of processive, gyrase-dependent supercoiling, whereas the H. influenzae or Salmonella serovar Typhi sequences did not. While consistent with a requirement for enhanced processivity of supercoiling for a site to function in Mu replication, these data suggest that other factors are also important. The relevance of these observations to an understanding of the function of the SGS is discussed.

Microarray analysis of transposition targets in Escherichia coli: the impact of transcription

Transposable elements have influenced the genetic and physical composition of all modern organisms. Defining how different transposons select target sites is critical for understanding the biochemical mechanism of this type of recombination and the impact of mobile genes on chromosome structure and function. Phage Mu replicates in Gram-negative bacteria using an extremely efficient transposition reaction. Replicated copies are excised from the chromosome and packaged into virus particles. Each viral genome plus several hundred base pairs of host DNA covalently attached to the prophage right end is packed into a virion. To study Mu transposition preferences, we used DNA microarray technology to measure the abundance of >4,000 Escherichia coli genes in purified Mu phage DNA. Insertion hot- and cold-spot genes were found throughout the genome, reflecting >1,000-fold variation in utilization frequency. A moderate preference was observed for genes near the origin compared to terminus of replication. Large biases were found at hot and cold spots, which often include several consecutive genes. Efficient transcription of genes had a strong negative influence on transposition. Our results indicate that local chromosome structure is more important than DNA sequence in determining Mu target-site selection.

Diversity and population structure of a near-shore marine-sediment viral community

Viruses, most of which are phage, are extremely abundant in marine sediments, yet almost nothing is known about their identity or diversity. We present the metagenomic analysis of an uncultured near-shore marine-sediment viral community. Three-quarters of the sequences in the sample were not related to anything previously reported. Among the sequences that could be identified, the majority belonged to double-stranded DNA phage. Temperate phage were more common than lytic phage, suggesting that lysogeny may be an important lifestyle for sediment viruses. Comparisons between the sediment sample and previously sequenced seawater viral communities showed that certain phage phylogenetic groups were abundant in all marine viral communities, while other phage groups were under-represented or absent. This 'marineness' suggests that marine phage are derived from a common set of ancestors. Several independent mathematical models, based on the distribution of overlapping shotgun sequence fragments from the library, were used to show that the diversity of the viral community was extremely high, with at least 10(4) viral genotypes per kilogram of sediment and a Shannon index greater than 9 nats. Based on these observations we propose that marine-sediment viral communities are one of the largest unexplored reservoirs of sequence space on the planet.

The impact of prophages on bacterial chromosomes

Prophages were automatically localized in sequenced bacterial genomes by a simple semantic script leading to the identification of 190 prophages in 115 investigated genomes. The distribution of prophages with respect to presence or absence in a given bacterial species, the location and orientation of the prophages on the replichore was not homogeneous. In bacterial pathogens, prophages are particularly prominent. They frequently encoded virulence genes and were major contributors to the genetic individuality of the strains. However, some commensal and free-living bacteria also showed prominent prophage contributions to the bacterial genomes. Lysogens containing multiple sequence-related prophages can experience rearrangements of the bacterial genome across prophages, leading to prophages with new gene constellations. Transfer RNA genes are the preferred chromosomal integration sites, and a number of prophages also carry tRNA genes. Prophage integration into protein coding sequences can lead to either gene disruption or new proteins. The phage repressor, immunity and lysogenic conversion genes are frequently transcribed from the prophage. The expression of the latter is sometimes integrated into control circuits linking prophages, the lysogenic bacterium and its animal host. Prophages are apparently as easily acquired as they are lost from the bacterial chromosome. Fixation of prophage genes seems to be restricted to those with functions that have been co-opted by the bacterial host.

Genomic diversity of Burkholderia pseudomallei clinical isolates: subtractive hybridization reveals a Burkholderia mallei-specific prophage in B. pseudomallei 1026b

Burkholderia pseudomallei is the etiologic agent of the disease melioidosis and is a category B biological threat agent. The genomic sequence of B. pseudomallei K96243 was recently determined, but little is known about the overall genetic diversity of this species. Suppression subtractive hybridization was employed to assess the genetic variability between two distinct clinical isolates of B. pseudomallei, 1026b and K96243. Numerous mobile genetic elements, including a temperate bacteriophage designated phi1026b, were identified among the 1026b-specific suppression subtractive hybridization products. Bacteriophage phi1026b was spontaneously produced by 1026b, and it had a restricted host range, infecting only Burkholderia mallei. It possessed a noncontractile tail, an isometric head, and a linear 54,865-bp genome. The mosaic nature of the phi1026b genome was revealed by comparison with bacteriophage phiE125, a B. mallei-specific bacteriophage produced by Burkholderia thailandensis. The phi1026b genes for DNA packaging, tail morphogenesis, host lysis, integration, and DNA replication were nearly identical to the corresponding genes in phiE125. On the other hand, phi1026b genes involved in head morphogenesis were similar to head morphogenesis genes encoded by Pseudomonas putida and Pseudomonas aeruginosa bacteriophages. Consistent with this observation, immunogold electron microscopy demonstrated that polyclonal antiserum against phiE125 reacted with the tail of phi1026b but not with the head. The results presented here suggest that B. pseudomallei strains are genetically heterogeneous and that bacteriophages are major contributors to the genomic diversity of this species. The bacteriophage characterized in this study may be a useful diagnostic tool for differentiating B. pseudomallei and B. mallei, two closely related biological threat agents.

Genome of staphylococcal phage K: a new lineage of Myoviridae infecting gram-positive bacteria with a low G+C content

Phage K is a polyvalent phage of the Myoviridae family which is active against a wide range of staphylococci. Phage genome sequencing revealed a linear DNA genome of 127,395 bp, which carries 118 putative open reading frames. The genome is organized in a modular form, encoding modules for lysis, structural proteins, DNA replication, and transcription. Interestingly, the structural module shows high homology to the structural module from Listeria phage A511, suggesting intergenus horizontal transfer. In addition, phage K exhibits the potential to encode proteins necessary for its own replisome, including DNA ligase, primase, helicase, polymerase, RNase H, and DNA binding proteins. Phage K has a complete absence of GATC sites, making it insensitive to restriction enzymes which cleave this sequence. Three introns (lys-I1, pol-I2, and pol-I3) encoding putative endonucleases were located in the genome. Two of these (pol-I2 and pol-I3) were found to interrupt the DNA polymerase gene, while the other (lys-I1) interrupts the lysin gene. Two of the introns encode putative proteins with homology to HNH endonucleases, whereas the other encodes a 270-amino-acid protein which contains two zinc fingers (CX(2)CX(22)CX(2)C and CX(2)CX(23)CX(2)C). The availability of the genome of this highly virulent phage, which is active against infective staphylococci, should provide new insights into the biology and evolution of large broad-spectrum polyvalent phages.

Annotation of the pRhico plasmid of Azospirillum brasilense reveals its role in determining the outer surface composition

The plant growth-promoting soil bacterium Azospirillum brasilense enhances growth of economically important crops, such as wheat, corn and rice. In order to improve plant growth, a close bacterial association with the plant roots is needed. Genes encoded on a 90-MDa plasmid, denoted pRhico plasmid, present in A. brasilense Sp7, play an important role in plant root interaction. Sequencing, annotation and in silico analysis of this 90-MDa plasmid revealed the presence of a large collection of genes encoding enzymes involved in surface polysaccharide biosynthesis. Analysis of the 90-MDa plasmid genome provided evidence for its essential role in the viability of the bacterial cell.

The prophages of Lactobacillus johnsonii NCC 533: comparative genomics and transcription analysis

Two non-inducible, but apparently complete prophages were identified in the genome of the sequenced Lactobacillus johnsonii strain NCC 533. The 38- and 40-kb-long prophages Lj928 and Lj965 represent distinct lineages of Sfi11-like pac-site Siphoviridae unrelated at the DNA sequence level. The deduced structural proteins from Lj928 demonstrated aa sequence identity with Lactococcus lactis phage TP901-1, while Lj965 shared sequence links with Streptococcus thermophilus phage O1205. With the exception of tRNA genes, inserted between DNA replication and DNA packaging genes, the transcription of the prophage was restricted to the genome segments near both attachment sites. Transcribed genes unrelated to phage functions were inserted between the phage repressor and integrase genes; one group of genes shared sequence relatedness with a mobile DNA element in Staphylococcus aureus. A short, but highly transcribed region was located between the phage lysin and right attachment site; it lacked a protein-encoding function in one prophage.

Crystal Structure of the Mor Protein of Bacteriophage Mu, a Member of the Mor/C Family of Transcription Activators

Transcription from the middle promoter, Pm, of bacteriophage Mu requires the phage-encoded activator protein Mor and bacterial RNA polymerase. Mor is a sequence-specific DNA-binding protein that mediates transcription activation through its interactions with the C-terminal domains of the alpha and sigma subunits of bacterial RNA polymerase. Here we present the first structure for a member of the Mor/C family of transcription activators, the crystal structure of Mor to 2.2-A resolution. Each monomer of the Mor dimer is composed of two domains, the N-terminal dimerization domain and C-terminal DNA-binding domain, which are connected by a linker containing a beta strand. The N-terminal dimerization domain has an unusual mode of dimerization; helices alpha1 and alpha2 of both monomers are intertwined to form a four-helix bundle, generating a hydrophobic core that is further stabilized by antiparallel interactions between the two beta strands. Mutational analysis of key leucine residues in helix alpha1 demonstrated a role for this hydrophobic core in protein solubility and function. The C-terminal domain has a classical helix-turn-helix DNA-binding motif that is located at opposite ends of the elongated dimer. Since the distance between the two helix-turn-helix motifs is too great to allow binding to two adjacent major grooves of the 16-bp Mor-binding site, we propose that conformational changes in the protein and DNA will be required for Mor to interact with the DNA. The highly conserved glycines flanking the beta strand may act as pivot points, facilitating the conformational changes of Mor, and the DNA may be bent.

In vivo recombineering of bacteriophage lambda by PCR fragments and single-strand oligonucleotides

We demonstrate that the bacteriophage lambda Red functions efficiently recombine linear DNA or single-strand oligonucleotides (ss-oligos) into bacteriophage lambda to create specific changes in the viral genome. Point mutations, deletions, and gene replacements have been created. While recombineering with oligonucleotides, we encountered other mutations accompanying the desired point mutational change. DNA sequence analysis suggests that these unwanted mutations are mainly frameshift deletions introduced during oligonucleotide synthesis.

Secondary Structure Switching in Cro Protein Evolution

We report the solution structure of the Cro protein from bacteriophage P22. Comparisons of its sequence and structure to those of lambda Cro strongly suggest an alpha-to-beta secondary structure switching event during Cro evolution. The folds of P22 Cro and lambda Cro share a three alpha helix fragment comprising the N-terminal half of the domain. However, P22 Cro's C terminus folds as two helices, while lambda Cro's folds as a beta hairpin. The all-alpha fold found for P22 Cro appears to be ancestral, since it also occurs in cI proteins, which are anciently duplicated paralogues of Cro. PSI-BLAST and transitive homology analyses strongly suggest that the sequences of P22 Cro and lambda Cro are globally homologous despite encoding different folds. The alpha+beta fold of lambda Cro therefore likely evolved from its all-alpha ancestor by homologous secondary structure switching, rather than by nonhomologous replacement of both sequence and structure.

The pKO2 linear plasmid prophage of Klebsiella oxytoca

Temperate bacteriophages with plasmid prophages are uncommon in nature, and of these only phages N15 and PY54 are known to have a linear plasmid prophage with closed hairpin telomeres. We report here the complete nucleotide sequence of the 51,601-bp Klebsiella oxytoca linear plasmid pKO2, and we demonstrate experimentally that it is also a prophage. We call this bacteriophage phiKO2. An analysis of the 64 predicted phiKO2 genes indicate that it is a fairly close relative of phage N15; they share a mosaic relationship that is typical of different members of double-stranded DNA tailed-phage groups. Although the head, tail shaft, and lysis genes are not recognizably homologous between these phages, other genes such as the plasmid partitioning, replicase, prophage repressor, and protelomerase genes (and their putative targets) are so similar that we predict that they must have nearly identical DNA binding specificities. The phiKO2 virion is unusual in that its phage lambda-like tails have an exceptionally long (3,433 amino acids) central tip tail fiber protein. The phiKO2 genome also carries putative homologues of bacterial dinI and umuD genes, both of which are involved in the host SOS response. We show that these divergently transcribed genes are regulated by LexA protein binding to a single target site that overlaps both promoters.

Analysis of the lambdoid prophage element e14 in the E. coli K-12 genome

BACKGROUND

Many sequenced bacterial genomes harbor phage-like elements or cryptic prophages. These elements have been implicated in pathogenesis, serotype conversion and phage immunity. The e14 element is a defective lambdoid prophage element present at 25 min in the E. coli K-12 genome. This prophage encodes important functional genes such as lit (T4 exclusion), mcrA (modified cytosine restriction activity) and pin (recombinase).

RESULTS

Bioinformatic analysis of the e14 prophage sequence shows the modular nature of the e14 element which shares a large part of its sequence with the Shigella flexneri phage SfV. Based on this similarity, the regulatory region including the repressor and Cro proteins and their binding sites were identified. The protein product of b1149 was found to be a fusion of a replication protein and a terminase. The genes b1143, b1151 and b1152 were identified as putative pseudogenes. A number of duplications of the stfE tail fibre gene of the e14 are seen in plasmid p15B. A protein based comparative approach using the COG database as a starting point helped detect lambdoid prophage like elements in a representative set of completely sequenced genomes.

CONCLUSIONS

The e14 element was characterized for the function of its encoded genes, the regulatory regions, replication origin and homology with other phage and bacterial sequences. Comparative analysis at nucleotide and protein levels suggest that a number of important phage related functions are missing in the e14 genome including parts of the early left operon, early right operon and late operon. The loss of these genes is the result of at least three major deletions that have occurred on e14 since its integration. A comparative protein level approach using the COG database can be effectively used to detect defective lambdoid prophage like elements in bacterial genomes.

The Bacillus thuringiensis linear double-stranded DNA phage Bam35, which is highly similar to the Bacillus cereus linear plasmid pBClin15, has a prophage state

Bam35, a 15-kbp double-stranded DNA phage, infects Bacillus thuringiensis. Recently, sequencing of the related Bacillus cereus revealed a 15.1-kbp linear plasmid, pBClin15. We show that pBClin15 closely resembles Bam35 and demonstrate conversion of Bam35 to a prophage. This state is common, as several B. thuringiensis strains release Bam35-related viruses.

Bacteriophage genomics

Comparative genomic studies of bacteriophages, especially the tailed phages, together with environmental studies, give a dramatic new picture of the size, genetic structure and dynamics of this population. Sequence comparisons reveal some of the detailed mechanisms by which these viruses evolve and influence the evolution of their bacterial and archaeal hosts. We see rampant horizontal exchange of sequences among genomes, mediated by both homologous and nonhomologous recombination. High frequency exchange among phages occupying similar ecological niches leads to a high degree of mosaic diversity in local populations. Horizontal exchange also takes place at lower frequency across the entire span of phage sequence space.