Evolutionary relationships among diverse bacteriophages and prophages: all the world's a phage

Elevated abundance of bacteriophage infecting bacteria in soil

Here we report the first direct counts of soil bacteriophage and show that substantial populations of these viruses exist in soil (grand mean = 1.5 x 10(7) g(-1)), at least 350-fold more than the highest numbers estimated from traditional viable plaque counts. Adding pure cultures of a Serratia phage to soil showed that the direct counting methods with electron microscopy developed here underestimated the added phage populations by at least eightfold. So, assuming natural phages were similarly underestimated, virus numbers in soil averaged 1.5 x 10(8) g(-1), which is equivalent to 4% of the total population of bacteria. This high abundance was to some extent confirmed by hybridizing colonies grown on Serratia and Pseudomonas selective media with cocktails of phage infecting these bacteria. This showed that 8.9 and 3.9%, respectively, hybridized with colonies from the two media and confirmed the presence of phage DNA sequences in the cultivable fraction of the natural population. Thus, soil phage, like their aquatic counterparts, are likely to be important in controlling bacterial populations and mediating gene transfer in soil.

Vibrio cholerae phage K139: complete genome sequence and comparative genomics of related phages

In this report, we characterize the complete genome sequence of the temperate phage K139, which morphologically belongs to the Myoviridae phage family (P2 and 186). The prophage genome consists of 33,106 bp, and the overall GC content is 48.9%. Forty-four open reading frames were identified. Homology analysis and motif search were used to assign possible functions for the genes, revealing a close relationship to P2-like phages. By Southern blot screening of a Vibrio cholerae strain collection, two highly K139-related phage sequences were detected in non-O1, non-O139 strains. Combinatorial PCR analysis revealed almost identical genome organizations. One region of variable gene content was identified and sequenced. Additionally, the tail fiber genes were analyzed, leading to the identification of putative host-specific sequence variations. Furthermore, a K139-encoded Dam methyltransferase was characterized.

Bacteriophage HP2 of Haemophilus influenzae

Temperate bacteriophages effect chromosomal evolution of their bacterial hosts, mediating rearrangements and the acquisition of novel genes from other taxa. Although the Haemophilus influenzae genome shows evidence of past phage-mediated lateral transfer, the phages presumed responsible have not been identified. To date, six different H. influenzae phages are known; of these, only the HP1/S2 group, which lyosogenizes exclusively Rd strains (which were originally encapsulated serotype d), is well characterized. Phages in this group are genetically very similar, with a highly conserved set of genes. Because the majority of H. influenzae strains are nonencapsulated (nontypeable), it is important to characterize phages infecting this larger, genetically more diverse group of respiratory pathogens. We have identified and sequenced HP2, a bacteriophage of nontypeable H. influenzae. Although related to the fully sequenced HP1 (and even more so to the partially sequenced S2) and similar in genetic organization, HP2 has a few novel genes and differs in host range; HP2 will not infect or lysogenize Rd strains. Genomic comparisons between HP1/S2 and HP2 suggest recent divergence, with new genes completely replacing old ones at certain loci. Sequence comparisons suggest that H. influenzae phages evolve by recombinational exchange of genes with each other, with cryptic prophages, and with the host chromosome.

The phage N4 virion RNA polymerase catalytic domain is related to single-subunit RNA polymerases

In vitro, bacteriophage N4 virion RNA polymerase (vRNAP) recognizes in vivo sites of transcription initiation on single-stranded templates. N4 vRNAP promoters are comprised of a hairpin structure and conserved sequences. Here, we show that vRNAP consists of a single 3500 amino acid polypeptide, and we define and characterize a transcriptionally active 1106 amino acid domain (mini-vRNAP). Biochemical and genetic characterization of this domain indicates that, despite its peculiar promoter specificity and lack of extensive sequence similarity to other DNA-dependent RNA polymerases, mini-vRNAP is related to the family of T7-like RNA polymerases.

Integration and excision of the Mycobacterium tuberculosis prophage-like element, phiRv1

The genomes of Mycobacterium tuberculosis H37Rv and CDC1551 each contain two prophage-like elements, phiRv1 and phiRv2. The phiRv1 element is not only absent from Mycobacterium bovis BCG but is in different locations within the two sequenced M. tuberculosis genomes; in both cases phiRv1 is inserted into a REP13E12 repeated sequence, which presumably contains the bacterial attachment site, attB, for phiRv1. Although phiRv1 is probably too small to encode infectious phage particles, it may nevertheless have an active integration/excision system and be capable of moving from one chromosomal position to another. We show here that the M. tuberculosis H37Rv phiRv1 element does indeed encode an active site-specific recombination system in which an integrase of the serine recombinase family (Rv1586c) catalyses integration and excision and a small, basic phiRv1-encoded protein (Rv1584c) controls the directionality of re-combination. Integration-proficient plasmid vectors derived from phiRv1 efficiently transform BCG, can utilize four of the seven REP13E12 sites present in BCG as attachment sites, and can occupy more than one site simultaneously.

Imbroglios of viral taxonomy: genetic exchange and failings of phenetic approaches

The practice of classifying organisms into hierarchical groups originated with Aristotle and was codified into nearly immutable biological law by Linnaeus. The heart of taxonomy is the biological species, which forms the foundation for higher levels of classification. Whereas species have long been established among sexual eukaryotes, achieving a meaningful species concept for prokaryotes has been an onerous task and has proven exceedingly difficult for describing viruses and bacteriophages. Moreover, the assembly of viral "species" into higher-order taxonomic groupings has been even more tenuous, since these groupings were based initially on limited numbers of morphological features and more recently on overall genomic similarities. The wealth of nucleotide sequence information that catalyzed a revolution in the taxonomy of free-living organisms necessitates a reevaluation of the concept of viral species, genera, families, and higher levels of classification. Just as microbiologists discarded dubious morphological traits in favor of more accurate molecular yardsticks of evolutionary change, virologists can gain new insight into viral evolution through the rigorous analyses afforded by the molecular phylogenetics of viral genes. For bacteriophages, such dissections of genomic sequences reveal fundamental flaws in the Linnaean paradigm that necessitate a new view of viral evolution, classification, and taxonomy.

The ICESt1 element of Streptococcus thermophilus belongs to a large family of integrative and conjugative elements that exchange modules and change their specificity of integration

The 34,734-bp element ICESt1 from Streptococcus thermophilus CNRZ368 is site-specifically integrated into the 3(') end of the gene fda. ICESt1 encodes integrative functions and putative transfer functions. Six proteins of the putative conjugative system of ICESt1 are related to those encoded by the conjugative transposon Tn916 from Enterococcus faecalis. A comparison of these proteins with those encoded by the complete or partial genome sequences of various low G+C bacteria including Bacillus subtilis, Clostridium difficile, E. faecalis, Listeria monocytogenes, Staphylococcus aureus, and Streptococcus mutans revealed the presence of numerous putative site-specific integrative conjugative elements and/or conjugative transposons within these genomes. Sequence comparisons revealed that these elements possess a modular structure and that exchanges of unrelated or distantly related modules and genes have occurred between these elements, and also plasmids and prophages. These exchanges have probably led to modifications in the site specificity of integration of these elements. Therefore, a distinction between low specificity integrative conjugative elements (i.e., conjugative transposons) and site-specific integrative conjugative elements does not appear to be relevant. We propose to call all the conjugative elements that excise by site-specific recombination and integrate by recombination between a specific site of a circular intermediate and another site, "Integrative and Conjugative Elements" (ICEs), irrespective of the integration specificity.

Genomic analysis of Clostridium perfringens bacteriophage phi3626, which integrates into guaA and possibly affects sporulation

Two temperate viruses, phi3626 and phi8533, have been isolated from lysogenic Clostridium perfringens strains. Phage phi3626 was chosen for detailed analysis and was inspected by electron microscopy, protein profiling, and host range determination. For the first time, the nucleotide sequence of a bacteriophage infecting Clostridium species was determined. The virus belongs to the Siphoviridae family of the tailed phages, the order Caudovirales. Its genome consists of a linear double-stranded DNA molecule of 33,507 nucleotides, with invariable 3'-protruding cohesive ends of nine residues. Fifty open reading frames were identified, which are organized in three major life cycle-specific gene clusters. The genes required for lytic development show an opposite orientation and arrangement compared to the lysogeny control region. A function could be assigned to 19 gene products, based upon bioinformatic analyses, N-terminal amino acid sequencing, or experimental evidence. These include DNA-packaging proteins, structural components, a dual lysis system, a putative lysogeny switch, and proteins that are involved in replication, recombination, and modification of phage DNA. The presence of genes encoding a putative sigma factor related to sporulation-dependent sigma factors and a putative sporulation-dependent transcription regulator suggests a possible interaction of phi3626 with onset of sporulation in C. perfringens. We found that the phi3626 attachment site attP lies in a noncoding region immediately downstream of int. Integration of the viral genome occurs into the bacterial attachment site attB, which is located within the 3' end of a guaA homologue. This essential housekeeping gene is functionally independent of the integration status, due to reconstitution of its terminal codons by phage sequence.

One of two copies of the gene for the activatable shiga toxin type 2d in Escherichia coli O91:H21 strain B2F1 is associated with an inducible bacteriophage

Shiga toxin (Stx) types 1 and 2 are encoded within intact or defective temperate bacteriophages in Stx-producing Escherichia coli (STEC), and expression of these toxins is linked to bacteriophage induction. Among Stx2 variants, only stx(2e) from one human STEC isolate has been reported to be carried within a toxin-converting phage. In this study, we examined the O91:H21 STEC isolate B2F1, which carries two functional alleles for the potent activatable Stx2 variant toxin, Stx2d, for the presence of Stx2d-converting bacteriophages. We first constructed mutants of B2F1 that produced one or the other Stx2d toxin and found that the mutant that produced only Stx2d1 made less toxin than the Stx2d2-producing mutant. Consistent with that result, the Stx2d1-producing mutant was attenuated in a streptomycin-treated mouse model of STEC infection. When the mutants were treated with mitomycin C to promote bacteriophage induction, Vero cell cytotoxicity was elevated only in extracts of the Stx2d1-producing mutant. Additionally, when mice were treated with ciprofloxacin, an antibiotic that induces the O157:H7 Stx2-converting phage, the animals were more susceptible to the Stx2d1-producing mutant. Moreover, an stx(2d1)-containing lysogen was isolated from plaques on strain DH5alpha that had been exposed to lysates of the mutant that produced Stx2d1 only, and supernatants from that lysogen transformed with a plasmid encoding RecA were cytotoxic when the lysogen was induced with mitomycin C. Finally, electron-microscopic examination of extracts from the Stx2d1-producing mutant showed hexagonal particles that resemble the prototypic Stx2-converting phage 933W. Together these observations provide strong evidence that expression of Stx2d1 is bacteriophage associated. We conclude that despite the sequence similarity of the stx(2d1)- and stx(2d2)-flanking regions in B2F1, Stx2d1 expression is repressed within the context of its toxin-converting phage while Stx2d2 expression is independent of phage induction.

The Phage Proteomic Tree: a genome-based taxonomy for phage

There are approximately 10(31) phage in the biosphere, making them the most abundant biological entities on the planet. Despite their great numbers and ubiquitous presence, very little is known about phage biodiversity, biogeography, or phylogeny. Information is limited, in part, because the current ICTV taxonomical system is based on culturing phage and measuring physical parameters of the free virion. No sequence-based taxonomic systems have previously been established for phage. We present here the "Phage Proteomic Tree," which is based on the overall similarity of 105 completely sequenced phage genomes. The Phage Proteomic Tree places phage relative to both their near neighbors and all other phage included in the analysis. This method groups phage into taxa that predicts several aspects of phage biology and highlights genetic markers that can be used for monitoring phage biodiversity. We propose that the Phage Proteomic Tree be used as the basis of a genome-based taxonomical system for phage.

Activation of prophage eib genes for immunoglobulin-binding proteins by genes from the IbrAB genetic island of Escherichia coli ECOR-9

Four distinct Escherichia coli immunoglobulin-binding (eib) genes, each of which encodes a surface-exposed protein that binds immunoglobulins in a nonimmune manner, are carried by separate prophages in E. coli reference (ECOR) strain ECOR-9. Each eib gene was transferred to test E. coli strains, both in the form of multicopy recombinant plasmids and as lysogenized prophage. The derived lysogens express little or no Eib protein, in sharp contrast to the parental lysogen, suggesting that ECOR-9 has an expression-enhancing activity that the derived lysogens lack. Supporting this hypothesis, we cloned from ECOR-9 overlapping genes, ibrA and ibrB (designation is derived from "immunoglobulin-binding regulator"), which together activated eib expression in the derived lysogens. The proteins encoded by ibrA and ibrB are very similar to uncharacterized proteins encoded by genes of Salmonella enterica serovar Typhi and E. coli O157:H7 (in a prophage-like element of the Sakai strain and in two O islands of strain EDL933). The genomic segment containing ibrA and ibrB has been designated the IbrAB island. It contains regions of homology to the Shiga toxin-converting prophage, Stx2, as well as genes homologous to phage antirepressor genes. The left boundary between the IbrAB island and the chromosomal framework is located near min 35.8 of the E. coli K-12 genome. Homology to IbrAB was found in certain other ECOR strains, including the other five eib-positive strains and most strains of the phylogenetic group B2. Sequencing of a 1.1-kb portion of ibrAB revealed that the other eib-positive strains diverge by </=0.1% from ECOR-9, whereas eib-negative ECOR-47 diverges by 16%.

Burkholderia thailandensis E125 harbors a temperate bacteriophage specific for Burkholderia mallei

Burkholderia thailandensis is a nonpathogenic gram-negative bacillus that is closely related to Burkholderia mallei and Burkholderia pseudomallei. We found that B. thailandensis E125 spontaneously produced a bacteriophage, termed phiE125, which formed turbid plaques in top agar containing B. mallei ATCC 23344. We examined the host range of phiE125 and found that it formed plaques on B. mallei but not on any other bacterial species tested, including B. thailandensis and B. pseudomallei. Examination of the bacteriophage by transmission electron microscopy revealed an isometric head and a long noncontractile tail. B. mallei NCTC 120 and B. mallei DB110795 were resistant to infection with phiE125 and did not produce lipopolysaccharide (LPS) O antigen due to IS407A insertions in wbiE and wbiG, respectively. wbiE was provided in trans on a broad-host-range plasmid to B. mallei NCTC 120, and it restored LPS O-antigen production and susceptibility to phiE125. The 53,373-bp phiE125 genome contained 70 genes, an IS3 family insertion sequence (ISBt3), and an attachment site (attP) encompassing the 3' end of a proline tRNA (UGG) gene. While the overall genetic organization of the phiE125 genome was similar to lambda-like bacteriophages and prophages, it also possessed a novel cluster of putative replication and lysogeny genes. The phiE125 genome encoded an adenine and a cytosine methyltransferase, and purified bacteriophage DNA contained both N6-methyladenine and N4-methylcytosine. The results presented here demonstrate that phiE125 is a new member of the lambda supergroup of Siphoviridae that may be useful as a diagnostic tool for B. mallei.

A novel bacteriophage-encoded RNA polymerase binding protein inhibits transcription initiation and abolishes transcription termination by host RNA polymerase

Xp10 is a lytic bacteriophage of Xanthomonas oryzae, a Gram-negative bacterium that causes rice blight. We purified an Xp10 protein, p7, that binds to and inhibits X. oryzae RNA polymerase (RNAP). P7 is a novel 73 amino acid-long protein; it does not bind to and hence does not affect transcription by Escherichia coli RNAP. Analysis of E. coli/X. oryzae RNAP hybrids locates the p7 binding site to the largest X. oryzae RNAP subunit, beta'. Binding of p7 to X. oryzae RNAP holoenzyme prevents large conformational change that places the sigma subunit region 4 into the correct position for interaction with the -35 promoter element. As a result, open promoter complex formation on the -10/-35 class promoters is inhibited. Inhibition of promoter complex formation on the extended -10 class promoters is less efficient. The p7 protein also abolishes factor-independent transcription termination by X. oryzae RNAP by preventing the release of nascent RNA at terminators. Further physiological and mechanistic studies of this novel transcription factor should provide additional insights into its biological role and the processes of promoter recognition and transcription termination.

Complete nucleotide sequence and likely recombinatorial origin of bacteriophage T3

We report the complete genome sequence (38,208 bp) of bacteriophage T3 and provide a bioinformatic comparative analysis with other completely sequenced members of the T7 group of phages. This comparison suggests that T3 has evolved from a recombinant between a T7-like coliphage and a yersiniophage. To assess this, recombination between T7 and the Yersinia enterocolitica serotype O:3 phage phiYeO3-12 was accomplished in vivo; coliphage progeny from this cross were selected that had many biological properties of T3. This represents the first experimentally observed recombination between lytic phages whose normal hosts are different bacterial genera.

Evolution of viral structure

Viruses vastly outnumber their host cells and must present a huge selective pressure. It is also becoming evident that only a small percent of the eukaryotic genome codes for molecules involved in cellular structures and functions, and that much of the remainder may have a viral origin. Viruses clearly play a central role in the biosphere, but how is this viral world organized? Classification was originally based on virus morphology and the particular host infected, but now there is an increasing trend to rely on sequence information. The type of genome (e.g., RNA or DNA, single- or double-stranded) provides fundamental classification criteria, while sequence comparisons can provide fine mapping for closely related viruses. However, it is currently very difficult to identify long-range evolutionary relationships. We present here a different approach, based on the idea that each virus has an innate "self." When the structures and functions characteristic of this "self" are identified, then they uncover relationships beyond those accessible from sequence information alone. The new approach is illustrated by sketching some possible viral lineages. We propose that urviruses were present before the division of cellular life into its current domains, and that the viral world has lineages that can be traced back to the root of the universal tree of life.

Horizontal gene transfer in microbial genome evolution

Horizontal gene transfer is the collective name for processes that permit the exchange of DNA among organisms of different species. Only recently has it been recognized as a significant contribution to inter-organismal gene exchange. Traditionally, it was thought that microorganisms evolved clonally, passing genes from mother to daughter cells with little or no exchange of DNA among diverse species. Studies of microbial genomes, however, have shown that genomes contain genes that are closely related to a number of different prokaryotes, sometimes to phylogenetically very distantly related ones. (Doolittle et al., 1990, J. Mol. Evol. 31, 383-388; Karlin et al., 1997, J. Bacteriol. 179, 3899-3913; Karlin et al., 1998, Annu. Rev. Genet. 32, 185-225; Lawrence and Ochman, 1998, Proc. Natl. Acad. Sci. USA 95, 9413-9417; Rivera et al., 1998, Proc. Natl. Acad. Sci. USA 95, 6239-6244; Campbell, 2000, Theor. Popul. Biol. 57 71-77; Doolittle, 2000, Sci. Am. 282, 90-95; Ochman and Jones, 2000, Embo. J. 19, 6637-6643; Boucher et al. 2001, Curr. Opin., Microbiol. 4, 285-289; Wang et al., 2001, Mol. Biol. Evol. 18, 792-800). Whereas prokaryotic and eukaryotic evolution was once reconstructed from a single 16S ribosomal RNA (rRNA) gene, the analysis of complete genomes is beginning to yield a different picture of microbial evolution, one that is wrought with the lateral movement of genes across vast phylogenetic distances. (Lane et al., 1988, Methods Enzymol. 167, 138-144; Lake and Rivera, 1996, Proc. Natl. Acad. Sci. USA 91, 2880-2881; Lake et al., 1999, Science 283, 2027-2028).

Bacteriophages: evolution of the majority

The dsDNA-tailed bacteriophages are probably the largest evolving group in the Biosphere and they are arguably very ancient. Comparative examination of genomes indicates that the hallmark of phage evolution is horizontal exchange of sequences. This is accomplished, first, by rampant non-homologous recombination between different genomes and, second, by reassortment of the variant sequences so created through homologous recombination. The comparative analysis suggests mechanisms by which new genes can be added to phage genomes and by which genes with novel functions may be assembled from parts. Horizontal exchange of sequences occurs most frequently among closely related phages, but it also extends across the entire global population at lower frequency. Bacteriophages also have probable ancestral connections with viruses of eukaryotes and archaea.

Snapshot of the genome of the pseudo-T-even bacteriophage RB49

RB49 is a virulent bacteriophage that infects Escherichia coli. Its virion morphology is indistinguishable from the well-known T-even phage T4, but DNA hybridization indicated that it was phylogenetically distant from T4 and thus it was classified as a pseudo-T-even phage. To further characterize RB49, we randomly sequenced small fragments corresponding to about 20% of the approximately 170-kb genome. Most of these nucleotide sequences lacked sufficient homology to T4 to be detected in an NCBI BlastN analysis. However, when translated, about 70% of them encoded proteins with homology to T4 proteins. Among these sequences were the numerous components of the virion and the phage DNA replication apparatus. Mapping the RB49 genes revealed that many of them had the same relative order found in the T4 genome. The complete nucleotide sequence was determined for the two regions of RB49 genome that contain most of the genes involved in DNA replication. This sequencing revealed that RB49 has homologues of all the essential T4 replication genes, but, as expected, their sequences diverged considerably from their T4 homologues. Many of the nonessential T4 genes are absent from RB49 and have been replaced by unknown sequences. The intergenic sequences of RB49 are less conserved than the coding sequences, and in at least some cases, RB49 has evolved alternative regulatory strategies. For example, an analysis of transcription in RB49 revealed a simpler pattern of regulation than in T4, with only two, rather than three, classes of temporally controlled promoters. These results indicate that RB49 and T4 have diverged substantially from their last common ancestor. The different T4-type phages appear to contain a set of common genes that can be exploited differently, by means of plasticity in the regulatory sequences and the precise choice of a large group of facultative genes.

Genomic sequence and evolution of marine cyanophage P60: a new insight on lytic and lysogenic phages

The genome of cyanophage P60, a lytic virus which infects marine Synechococcus WH7803, was completely sequenced. The P60 genome contained 47,872 bp with 80 potential open reading frames that were mostly similar to the genes found in lytic phages like T7, phi-YeO3-12, and SIO1. The DNA replication system, consisting of primase-helicase and DNA polymerase, appeared to be more conserved in podoviruses than in siphoviruses and myoviruses, suggesting that DNA replication genes could be the critical elements for lytic phages. Strikingly high sequence similarities in the regions coding for nucleotide metabolism were found between cyanophage P60 and marine unicellular cyanobacteria.

Genomic subtraction to identify and characterize sequences of Shiga toxin-producing Escherichia coli O91:H21

To identify Shiga toxin-producing Escherichia coli genes associated with severe human disease, a genomic subtraction technique was used with hemolytic-uremic syndrome-associated O91:H21 strain CH014 and O6:H10 bovine strains. The method was adapted to the Shiga toxin-producing E. coli genome: three rounds of subtraction were used to isolate DNA fragments specific to strain CH014. The fragments were characterized by genetic support analysis, sequencing, and hybridization to the genome of a collection of Shiga toxin-producing E. coli strains. A total of 42 fragments were found, 19 of which correspond to previously identified unique DNA sequences in the enterohemorrhagic E. coli EDL933 reference strain, including 7 fragments corresponding to prophage sequences and others encoding candidate virulence factors, such a SepA homolog protein and a fimbrial usher protein. In addition, the subtraction procedure yielded plasmid-related sequences from Shigella flexneri and enteropathogenic and Shiga toxin-producing E. coli virulence plasmids. We found that lateral gene transfer is extensive in strain CH014, and we discuss the role of genomic mobile elements, especially bacteriophages, in the evolution and possible transfer of virulence determinants.

Comparative analysis of the genomes of the temperate bacteriophages phi 11, phi 12 and phi 13 of Staphylococcus aureus 8325

The genomes of the three temperate bacteriophages contained in the chromosome of Staphylococcus aureus 8325 have been extracted from the sequence database and analyzed. phi 11, phi 12 and phi 13 are members of the same lytic group but different serogroups and consequently co-habitate the same host cell. Their genomes are approximately 42 kb to 45 kb and contain about 90 ORFs of at least 50 codons. Of these, about 50 have similarities to known genes or to genes of other staphylococcal phages. Each of the phages clusters within a homology group that share large regions of sequence identity while intergroup homology is comparatively low. The arrangement of genes on the chromosomes of the three phages is similar and consistent with current modular theory of phage gene organization. The replicated genomes appear to be packaged by different mechanisms. Phage phi 11 and phi 12 have been found to contain sequences consistent with pac-site phages while phi 13 has sequences consistent with cos-site phages. The attBsite for phi 11 is located in an intergenic region of the S. aureus chromosome while phi 12 and phi 13 integrate into specific genes. The phi 12 att-site is within an unknown gene, but the phi 13 att-site is within the beta-toxin gene. In contrast to the other two phages, phi 13 also introduces the staphylokinase gene (sak) and a second gene related to expression of fib.

Alterations of the portal protein, gpB, of bacteriophage lambda suppress mutations in cosQ, the site required for termination of DNA packaging

The cosQ site of bacteriophage lambda is required for DNA packaging termination. Previous studies have shown that cosQ mutations can be suppressed in three ways: by a local suppressor within cosQ, an increase in the length of the lambda chromosome, and missense mutations affecting the prohead's portal protein, gpB. In the present work, revertants of a set of lethal cosQ mutants were screened for suppressors. Seven new cosQ suppressors affected gene B, which encodes the portal protein of the prohead. All seven were allele-nonspecific suppressors of cosQ mutations. Experiments with several phages having two cosQ suppressors showed that the suppression effects were additive. Furthermore, these double suppressors had minimal effects on the growth of cosQ(+) phages. These trans-acting suppressors affecting the portal protein are proposed to allow the mutant cosQ site to be more efficiently recognized, due to the slowing of the rate of translocation.

Evolution of the herpesviruses

The complete genome sequences of 26 herpesvirus species, some represented by more than one strain, are currently deposited in the public databases. Their combined length totals over four million base pairs. Several additional genomes are in the pipeline, and a host of partial sequence information is also available. Consequently, researchers have a very detailed picture of the genetic content of herpesviruses and their relationships to each other. This review of the herpesvirus phylogenetic tree explains what is currently known about the evolution of this important virus family by proceeding from the twigs, along the branches, boughs and trunk and towards the root. The account focuses naturally on places where the grasp is secure, but also ventures where the bark is slippery and looks out on regions where footholds have not yet been established.

C-terminal domains of Listeria monocytogenes bacteriophage murein hydrolases determine specific recognition and high-affinity binding to bacterial cell wall carbohydrates

Listeria monocytogenes phage endolysins Ply118 and Ply500 share a unique enzymatic activity and specifically hydrolyse Listeria cells at the completion of virus multiplication in order to release progeny phage. With the aim of determining the molecular basis for the lytic specificity of these enzymes, we have elucidated their domain structure and examined the function of their unrelated and unique C-terminal cell wall binding domains (CBDs). Analysis of deletion mutants showed that both domains are needed for lytic activity. Fusions of CBDs with green fluorescent protein (GFP) demonstrated that the C-terminal 140 amino acids of Ply500 and the C-terminal 182 residues of Ply118 are necessary and sufficient to direct the murein hydrolases to the bacterial cell wall. CBD500 was able to target GFP to the surface of Listeria cells belonging to serovar groups 4, 5 and 6, resulting in an even staining of the entire cell surface. In contrast, the CBD118 hybrid bound to a ligand predominantly present at septal regions and cell poles, but only on cells of serovars 1/2, 3 and 7. Non-covalent binding to surface carbohydrate ligands occurred in a rapid, saturation-dependent manner. We measured 4 x 104 and 8 x 104 binding sites for CBD118 and CBD500 respectively. Surface plasmon resonance analysis revealed unexpected high molecular affinity constants for the CBD-ligand interactions, corresponding to nanomolar affinities. In conclusion, we show that the CBDs are responsible for targeting the phage endolysins to their substrates and function to confer recognition specificity on the proteins. As the CBD sequences contain no repeats and lack all known sequence motifs for anchoring of proteins to the bacterial cell, we conclude that they use unique structural motifs for specific association with the surface of Gram-positive bacteria.

Sequence analysis of the lactococcal bacteriophage bIL170: insights into structural proteins and HNH endonucleases in dairy phages

The complete 31754 bp genome of bIL170, a virulent bacteriophage of Lactococcus lactis belonging to the 936 group, was analysed. Sixty-four ORFs were predicted and the function of 16 of them was assigned by significant homology to proteins in databases. Three putative homing endonucleases of the HNH family were found in the early region. An HNH endonuclease with zinc-binding motif was identified in the late cluster, potentially being part of the same functional module as terminase. Three putative structural proteins were analysed in detail and show interesting features among dairy phages. Notably, gpl12 (putative fibre) and gpl20 (putative baseplate protein) of bIL170 are related by at least one of their domains to a number of multi-domain proteins encoded by lactococcal or streptococcal phages. A 110- to 150-aa-long hypervariable domain flanked by two conserved motifs of about 20 aa was identified. The analysis presented here supports the participation of some of these proteins in host-range determination and suggests that specific adsorption to the host may involve a complex multi-component system. Divergences in the genome of phages of the 936 group, that may have important biological properties, were noted. Insertions/deletions of units of one or two ORFs were the main source of divergence in the early clusters of the two entirely sequenced phages, bIL170 and sk1. An exchange of fragments probably affected the regions containing the putative origin of replication. It led to the absence in bIL170 of the direct repeats recognized in sk1 and to the presence of different ORFs in the ori region. Shuffling of protein domains affected the endolysin (putative cell-wall binding part), as well as gpl12 and gpl20.

HF2: a double-stranded DNA tailed haloarchaeal virus with a mosaic genome

HF2 is a haloarchaeal virus infecting two Halorubrum species (Family Halobacteriaceae). It is lytic, has a head-and-tail morphology and belongs to the Myoviridae (contractile tails). The linear double-stranded DNA genome was sequenced and found to be 77 670 bp in length, with a mol% G+C of 55.8. A total of 121 likely open reading frames (ORFs) were identified, of which 37 overlapped at start and stop codons. The predicted proteins were usually acidic (average pI of 4.8), and less than about 12% of them had homologues in the sequence databases. Four complete tRNA-like sequences (tRNA-Arg, -Asx, -Pro and -Tyr) and an incomplete tRNA-Thr were detected. A transcription map showed that most of the genome was transcribed and that the synthesis of transcripts occurred in a highly organized and reproducible pattern over a 5 h infection cycle. Transcripts often spanned multiple ORFs, suggesting that viral genes were organized into operons. The predicted ORF and observed transcript directions matched well and showed that transcription is mainly directed inwards from the genome termini, meeting at about 45-48 kb, and this was also a turning point in a cumulative GC-skew plot. The low point in cumulative GC-skew, near the left end, was a region rich in short repeats and lacking ORFs, which is likely to be an origin of replication. The HF2 genome is a mosaic of components from widely different sources, demonstrating clearly that viruses of haloarchaea, like their bacteriophage counterparts, are vectors for the exchange and transmission of genetic material between wide taxonomic distances, even across domains.

Complete genomic sequence of SfV, a serotype-converting temperate bacteriophage of Shigella flexneri

Bacteriophage SfV is a temperate serotype-converting phage of Shigella flexneri. SfV encodes the factors involved in type V O-antigen modification, and the serotype conversion and integration-excision modules of the phage have been isolated and characterized. We now report on the complete sequence of the SfV genome (37,074 bp). A total of 53 open reading frames were predicted from the nucleotide sequence, and analysis of the corresponding proteins was used to construct a functional map. The general organization of the genes in the SfV genome is similar to that of bacteriophage lambda, and numerous features of the sequence are described. The superinfection immunity system of SfV includes a lambda-like repression system and a P4-like transcription termination mechanism. Sequence analysis also suggests that SfV encodes multiple DNA methylases, and experiments confirmed that orf-41 encodes a Dam methylase. Studies conducted to determine if the phage-encoded methylase confers host DNA methylation showed that the two S. flexneri strains analyzed encode their own Dam methylase. Restriction mapping and sequence analysis revealed that the phage genome has cos sites at the termini. The tail assembly and structural genes of SfV show homology to those of phage Mu and Mu-like prophages in the genome of Escherichia coli O157:H7 and Haemophilus influenzae. Significant homology (30% of the genome in total) between sections of the early, regulatory, and structural regions of the SfV genome and the e14 and KpLE1 prophages in the E. coli K-12 genome were noted, suggesting that these three phages have common evolutionary origins.

Filamentous phage active on the gram-positive bacterium Propionibacterium freudenreichii

We present the first description of a single-stranded DNA filamentous phage able to replicate in a gram-positive bacterium. Phage B5 infects Propionibacterium freudenreichii and has a genome consisting of 5,806 bases coding for 10 putative open reading frames. The organization of the genome is very similar to the organization of the genomes of filamentous phages active on gram-negative bacteria. The putative coat protein exhibits homology with the coat proteins of phages PH75 and Pf3 active on Thermus thermophilus and Pseudomonas aeruginosa, respectively. B5 is, therefore, evolutionarily related to the filamentous phages active on gram-negative bacteria.

The nucleotide sequence of Shiga toxin (Stx) 2e-encoding phage phiP27 is not related to other Stx phage genomes, but the modular genetic structure is conserved

In this study we determined the complete nucleotide sequence of Shiga toxin 2e-encoding bacteriophage phi P27, isolated from the Shiga toxin-producing Escherichia coli patient isolate 2771/97. phi P27 is integrated as a prophage in the chromosomal yecE gene. This integration generates identity segments of attL and attR sites with lengths of 11 nucleotides. The integrated prophage genome has a size of 42,575 bp. We identified 58 open reading frames (ORFs), each with a length of >150 nucleotides. The deduced proteins of 44 ORFs showed significant homologies to other proteins present in sequence databases, whereas 14 putative proteins did not. For 29 proteins, we could deduce a putative function. Most of these are related to the basic phage propagation cycle. The phi P27 genome represents a mosaic composed of genetic elements which are obviously derived from related and unrelated phages. We identified five short linker sequences of 22 to 151 bp in the phi P27 sequence which have also been detected in a couple of other lambdoid phages. These linkers are located between functional modules in the phage genome and are thought to play a role in genetic recombination. Although the overall DNA sequence of phi P27 is not highly related to other known phages, the data obtained demonstrate a typical lambdoid genome structure.

Bacteriophage Mu genome sequence: analysis and comparison with Mu-like prophages in Haemophilus, Neisseria and Deinococcus

We report the complete 36,717 bp genome sequence of bacteriophage Mu and provide an analysis of the sequence, both with regard to the new genes and other genetic features revealed by the sequence itself and by a comparison to eight complete or nearly complete Mu-like prophage genomes found in the genomes of a diverse group of bacteria. The comparative studies confirm that members of the Mu-related family of phage genomes are genetically mosaic with respect to each other, as seen in other groups of phages such as the phage lambda-related group of phages of enteric hosts and the phage L5-related group of mycobacteriophages. Mu also possesses segments of similarity, typically gene-sized, to genomes of otherwise non-Mu-like phages. The comparisons show that some well-known features of the Mu genome, including the invertible segment encoding tail fiber sequences, are not present in most members of the Mu genome sequence family examined here, suggesting that their presence may be relatively volatile over evolutionary time. The head and tail-encoding structural genes of Mu have only very weak similarity to the corresponding genes of other well-studied phage types. However, these weak similarities, and in some cases biochemical data, can be used to establish tentative functional assignments for 12 of the head and tail genes. These assignments are strongly supported by the fact that the order of gene functions assigned in this way conforms to the strongly conserved order of head and tail genes established in a wide variety of other phages. We show that the Mu head assembly scaffolding protein is encoded by a gene nested in-frame within the C-terminal half of another gene that encodes the putative head maturation protease. This is reminiscent of the arrangement established for phage lambda.

The genome of bacteriophage phiKZ of Pseudomonas aeruginosa

Bacteriophage phiKZ is a giant virus that efficiently infects Pseudomonas aeruginosa strains pathogenic to human and, therefore, it is attractive for phage therapy. We present here the complete phiKZ genome sequence and a preliminary analysis of its genome structure. The 280,334 bp genome is a linear, circularly permutated and terminally redundant, A+T-rich double-stranded DNA molecule. The phiKZ DNA has no detectable sequence homology to other viruses and microorganisms, and it does not contain NotI, PstI, SacI, SmaI, XhoI, and XmaIII endonuclease restriction sites. The genome has 306 open reading frames (ORFs) varying in size from 50 to 2237 amino acid residues. According to the orientation of transcription, ORFs are apparently organized into clusters and most have a clockwise direction. The phiKZ genome also encodes six tRNAs specific for Met (AUG), Asn (AAC), Asp (GAC), Leu (TTA), Thr (ACA), and Pro (CCA). A putative promoter sequence containing a TATATTAC block was identified. Most potential stem-loop transcription terminators contain the tetranucleotide UUCG loops. Some genes may be assigned as phage-encoded RNA polymerase subunits. Only 59 phiKZ gene products exhibit similarity to proteins of known function from a diversity of organisms. Most of these conserved gene products, such as dihydrofolate reductase, ribonucleoside diphosphate reductase, thymidylate synthase, thymidylate kinase, and deoxycytidine triphosphate deaminase are involved in nucleotide metabolism. However, no virus-encoded DNA polymerase, DNA replication-associated proteins, or single-stranded DNA-binding protein were found based on amino acid homology, and they may therefore be strongly divergent from known homologous proteins. Fifteen phiKZ gene products show homology to proteins of pathogenic organisms, including Mycobacterium tuberculosis, Haemophilus influenzae, Listeria sp., Rickettsia prowazakeri, and Vibrio cholerae that must be considered before using this phage as a therapeutic agent. The phiKZ coat contains at least 40 polypeptides, and several proteins are cleaved during virus assembly in a way similar to phage T4. Eleven phiKZ-encoded polypeptides are related to proteins of other bacteriphages that infect a variety of hosts. Among these are four gene products that contain a putative intron-encoded endonuclease harboring the H-N-H motif common to many double-stranded DNA phages. These observations provide evidence that phages infecting diverse hosts have had access to a common genetic pool. However, limited homology on the DNA and protein levels indicates that bacteriophage phiKZ represents an evolutionary distinctive branch of the Myoviridae family.

Evolution of lambdoid replication modules

Comparison of the putative iteron-binding proteins of lambdoid phages allows us to propose that in the case of lambdoid replication modules, the units on which natural selection acts do not coincide with the open reading frames. Rather, the first replication gene is split into two segments, and its 3' part (corresponding to the C-terminal domain of the iteron-binding protein) forms one unit with the second gene. We also propose from the phylogenetic analysis of phage-encoded homologs of E. coli DnaB and DnaC, that the recombination with the host sequences is not frequent. Accessory ATP-ases for helicase loading (E. coli DnaC homologs) may not be universal replication proteins. Our analysis may suggest that the bacterial helicase loaders might be of phage origin. The comparison of DnaC homologs of enterobacteria and enterobacterial phages supports the experimental data on residues important in interaction with DnaB. We propose that construction of plasmids carrying the replication origins of lambdoid prophages could be useful not only in further research on DNA replication but also on the role of these prophages in shuttling genes for bacterial virulence. The phage replication sequences could be also useful for identification of clinical enterobacterial isolates.

Temperate phages in Salmonella enterica serovar Typhimurium: implications for epidemiology

Salmonella enterica serovar Typhimurium is the most common Salmonella serovar isolated from humans in Australia. The most common definitive phage types (DT) include 9, 64 and 135. Induction of lysogenic phages from DT 64 with mitomycin C followed by cesium chloride gradient purification, resulted in separation of two populations of phage particles. DNA extracted from these particles that was digested with SmaI showed two distinct patterns of banding. Transmission electron microscopy showed that both phage particles belong to the podovirus family of the C1 morphotype. One of the phages, ST64T is capable of mediating both generalized transduction and bacteriophage type conversion. Crude phage lysate induced from S. Typhimurium DT 64 was capable of phage type conversion. S. Typhimurium DT 9 was converted to DT 64 and DT 135 was converted to DT 16. S. Typhimurium DT 41 was also converted to DT 29. Amplified-fragment length polymorphism revealed differences between the original isolates and the convertants. Phage type conversion raises the question of the stability of the bacterial phage types in natural settings and the possibility of its occurrence during an outbreak scenario.

Integration sites for genetic elements in prokaryotic tRNA and tmRNA genes: sublocation preference of integrase subfamilies

Most classical integrases of prokaryotic genetic elements specify integration into tRNA or tmRNA genes. Sequences shared between element and host integration sites suggest that crossover can occur at any of three sublocations within a tRNA gene, two with flanking symmetry (anticodon-loop and T-loop tDNA) and the third at the asymmetric 3' end of the gene. Integrase phylogeny matches this classification: integrase subfamilies use exclusively either the symmetric sublocations or the asymmetric sublocation, although tRNA genes of several different aminoacylation identities may be used within any subfamily. These two familial sublocation preferences imply two modes by which new integration site usage evolves. The tmRNA gene has been adopted as an integration site in both modes, and its distinctive structure imposes some constraints on proposed evolutionary mechanisms.

Microviridae, a family divided: isolation, characterization, and genome sequence of phiMH2K, a bacteriophage of the obligate intracellular parasitic bacterium Bdellovibrio bacteriovorus

A novel single-stranded DNA phage, phiMH2K, of Bdellovibrio bacteriovorus was isolated, characterized, and sequenced. This phage is a member of the Microviridae, a family typified by bacteriophage phiX174. Although B. bacteriovorus and Escherichia coli are both classified as proteobacteria, phiMH2K is only distantly related to phiX174. Instead, phiMH2K exhibits an extremely close relationship to the Microviridae of Chlamydia in both genome organization and encoded proteins. Unlike the double-stranded DNA bacteriophages, for which a wide spectrum of diversity has been observed, the single-stranded icosahedral bacteriophages appear to fall into two distinct subfamilies. These observations suggest that the mechanisms driving single-stranded DNA bacteriophage evolution are inherently different from those driving the evolution of the double-stranded bacteriophages.

Phage genomics: small is beautiful

The Age of Genomics dawned only gradually for bacteriophages. It was 1977 when the genome of phage phi X174 was published and 1983 when the "large" genome of phage lambda hit the streets. More recently, the pace has quickened, so that we now have over 100 complete phage genomes and can expect thousands in a very few years. These sequences have been marvelously informative for the biology of the individual phages, but with the advent of high volume sequencing technology, the real excitement for phage biology is that it is now possible to analyze the sequences together and thereby address--for the first time at whole genome resolution--a set of fundamental biological questions related to populations: What is the structure of the global phage population? What are its dynamics? How do phages evolve? This is Comparative Genomics with a capital "C".

Cryptic plasmids of Mycobacterium avium: Tn552 to the rescue

Plasmids have been described in almost all bacterial species analysed and have proven to be essential genetic tools. In many bacteria these extrachromosomal DNAs are cryptic with no known markers or function, which makes their characterization and genetic exploitation extremely difficult. Here we describe a system that will allow the rescue of any circular DNA (plasmid or phage) using an in vitro transposition system to deliver both a selectable marker (kanamycin) and an Escherichia coli plasmid origin of replication. In this study, we demonstrate the rescue of four cryptic plasmids from the opportunistic pathogen Mycobacterium avium. To evaluate the host range of the rescued plasmids, we have examined their ability to be propagated in Mycobacterium smegmatis and Mycobacterium bovis BCG, and their compatibility with other mycobacterial plasmids. In addition, we use a library of transposon insertions to sequence one plasmid, pVT2, and to begin a genetic analysis of plasmid genes. Using this approach, we identified a putative conjugative relaxase, suggesting this myco-bacterial plasmid is transferable, and three genes required for plasmid establishment and replication.

Structure and function of the Escherichia coli RecE protein, a member of the RecB nuclease domain family

The RecB subunit of the Escherichia coli RecBCD enzyme has both helicase and nuclease activities. The helicase function was localized to an N-terminal domain, whereas the nuclease activity was found in a C-terminal domain. Recent analysis has uncovered a group of proteins that have weak amino acid sequence similarity to the RecB nuclease domain and that are proposed to constitute a family of related proteins (Aravind, L., Walker, D. R., and Koonin, E. V. (1999) Nucleic Acids Res. 27, 1223-1242). One is the E. coli RecE protein (exonuclease VIII), an ATP-independent exonuclease that degrades the 5'-terminated strand of double-stranded DNA. We have made mutations in several residues of RecE that align with the critical residues of RecB, and we find that the mutations reduce or abolish the nuclease activity of RecE but do not affect the enzyme binding to linear double-stranded DNA. Proteolysis experiments with subtilisin show that a stable 34-kilodalton C-terminal domain that contains these critical residues has nuclease activity, whereas no stable proteolytic fragments accumulate from the N-terminal portion of RecE. These results show that RecE has a nuclease domain and active site that are similar to RecB, despite the very weak sequence similarity between the two proteins. These similarities support the hypothesis that the nuclease domains of the two proteins are evolutionarily related.

Ancestral divergence, genome diversification, and phylogeographic variation in subpopulations of sorbitol-negative, beta-glucuronidase-negative enterohemorrhagic Escherichia coli O157

The O157:H7 lineage of enterohemorrhagic Escherichia coli is a geographically disseminated complex of highly related genotypes that share common ancestry. The common clone that is found worldwide carries several markers of events in its evolution, including markers for acquisition of virulence genes and loss of physiological characteristics, such as sorbitol fermentation ability and beta-glucuronidase production. Populations of variants that are distinct with respect to motility and the sorbitol and beta-glucuronidase markers appear to have diverged at several points along the inferred evolutionary pathway. In addition to these variants, distinct subpopulations of the contemporary non-sorbitol-fermenting, beta-glucuronidase-negative O157:H7 clone were recently detected among bovine and human clinical isolates in the United States by using high-resolution genome comparison. In order to determine if these recently described subpopulations were derived from a regional or ancestral divergence event, we used octamer-based genome scanning, marker sorting, and DNA sequence analysis to examine their phylogenetic relationship to populations of non-sorbitol-fermenting, beta-glucuronidase negative O157:H7 and O157:H- strains from Australia. The inferred phylogeny is consistent with the hypothesis that subpopulations on each continent resulted from geographic spread of an ancestral divergence event and subsequent expansion of distinct subpopulations. Marker sorting and DNA sequence analyses identified sets of monophyletic markers consistent with the pattern of divergence and demonstrated that phylogeographic variation occurred through emergence of regional subclones and concentration of regional polymorphisms among distinct subpopulations. DNA sequence analysis of representative polyphyletic markers showed that genome diversity accrued through random drift and bacteriophage-mediated events.

Detection of homologous recombination among bacteriophage P2 relatives

Sequencing of five late genes from 18 isolates of P2-like bacteriophages showed that these are at least 96% identical to the genes of phage P2. A maximum-parsimony phylogenetic analysis of these genes showed excess homoplasy of a magnitude three to six times higher than that expected. Examination of the distribution of the number of homoplasies at parsimoniously informative sites and incompatibility matrices of such sites revealed a pattern typical for extensive recombination. It has been shown that phage P2 probably incorporated some functionally complete genes or gene modules by recombination with other phages or with different hosts, but homologous recombination within genes has previously not been shown. In this paper we demonstrate that homologous recombination between P2-like bacteriophages occurs randomly at multiple breakpoints in five late genes. The rate of recombination is high but, since some phages were sampled decades apart and in different parts of the world, this has to be viewed on an evolutionary time scale. The applicability of different methods used for detection of recombination breakpoints and estimation of rates of recombination in bacteriophages is discussed.

Shiga-toxin-converting bacteriophages

Shiga toxins (Stx) comprise a family of potent cytotoxins that are involved in severe human disease. Stx are mainly produced by Escherichia coli isolated from human and nonhuman sources, and by Shigella dysenteriae type 1. The genes encoding Stx are thought to be generally encoded in the genome of lambdoid prophages (Stx-converting bacteriophages; Stx phages). They share a unique position in the late region of the phage genome downstream of the late promoter PR'. This location suggests that expression of stx is controlled by a Q-like antiterminator. Therefore, induction of Stx-converting prophages appears to trigger increased production of Stx. Following induction, Stx phages can be transduced in vivo and in vitro into other bacteria. Stx phages play an important role in the expression of Stx and in lateral gene transfer and are therefore a contribution to the emergence of new Stx-producing E. coli (STEC) variants.

Diversification of Escherichia coli genomes: are bacteriophages the major contributors?

Determination of the genome sequence of enterohemorrhagic Escherichia coli O157 Sakai and genomic comparison with the laboratory strain K-12 has revealed that the two strains share a highly conserved 4.1-Mb sequence and that each also contains a large amount of strain-specific sequence. The analysis also revealed the presence of a surprisingly large number of prophages in O157, most of which are lambda-like phages that resemble each other. Based on these results, we discuss how the E. coli strains have diverged from a common ancestral strain, and how bacteriophages contributed to this process. We also describe possible mechanisms by which O157 acquired many closely related phages, and raise the possibility that such bacteria might function as 'phage factories', releasing a variety of chimeric or mosaic phages into the environment.

Comparative genomics reveals close genetic relationships between phages from dairy bacteria and pathogenic Streptococci: evolutionary implications for prophage-host interactions

The genome of the highly pathogenic M1 serotype Streptococcus pyogenes isolate SF370 contains eight prophage elements. Only prophage SF370.1 could be induced by mitomycin C treatment. Prophage SF370.3 showed a 33.5-kb-long genome that closely resembled the genome organization of the cos-site temperate Siphovirus r1t infecting the dairy bacterium Lactococcus lactis. The two-phage genomes shared between 60 and 70% nucleotide sequence identity over the DNA packaging, head and tail genes. Analysis of the SF370.3 genome revealed mutations in the replisome organizer gene that may prevent the induction of the prophage. The mutated phage replication gene was closely related to a virulence marker identified in recently emerged M3 serotype S. pyogenes strains in Japan. This observation suggests that prophage genes confer selective advantage to the lysogenic host. SF370.3 encodes a hyaluronidase and a DNase that may facilitate the spreading of S. pyogenes through tissue planes of its human host. Prophage SF370.2 showed a 43-kb-long genome that closely resembled the genome organization of pac-site temperate Siphoviridae infecting the dairy bacteria S. thermophilus and L. lactis. Over part of the structural genes, the similarity between SF370.2 and S. thermophilus phage O1205 extended to the nucleotide sequence level. SF370.2 showed two probable inactivating mutations: one in the replisome organizer gene and another in the gene encoding the portal protein. Prophage SF370.2 also encodes a hyaluronidase and in addition two very likely virulence factors: prophage-encoded toxins acting as superantigens that may contribute to the immune deregulation observed during invasive streptococcal infections. The superantigens are encoded between the phage lysin and the right attachment site of the prophage genome. The genes were nearly sequence identical with a DNA segment in S. equi, suggesting horizontal gene transfer. The trend for prophage genome inactivation was even more evident for the remaining five prophage sequences that showed massive losses of prophage DNA. In these prophage remnants only 13-0.3 kb of putative prophage DNA was detected. We discuss the genomics data from S. pyogenes strain SF370 within the framework of Darwinian coevolution of prophages and lysogenic bacteria and suggest elements of genetic cooperation and elements of an arms race in this host-parasite relationship.

A conserved genetic module that encodes the major virion components in both the coliphage T4 and the marine cyanophage S-PM2

Sequence analysis of a 10-kb region of the genome of the marine cyanomyovirus S-PM2 reveals a homology to coliphage T4 that extends as a contiguous block from gene (g)18 to g23. The order of the S-PM2 genes in this region is similar to that of T4, but there are insertions and deletions of small ORFs of unknown function. In T4, g18 codes for the tail sheath, g19, the tail tube, g20, the head portal protein, g21, the prohead core protein, g22, a scaffolding protein, and g23, the major capsid protein. Thus, the entire module that determines the structural components of the phage head and contractile tail is conserved between T4 and this cyanophage. The significant differences in the morphology of these phages must reflect the considerable divergence of the amino acid sequence of their homologous virion proteins, which uniformly exceeds 50%. We suggest that their enormous diversity in the sea could be a result of genetic shuffling between disparate phages mediated by such commonly shared modules. These conserved sequences could facilitate genetic exchange by providing partially homologous substrates for recombination between otherwise divergent phage genomes. Such a mechanism would thus expand the pool of phage genes accessible by recombination to all those phages that share common modules.

Behavior of restriction-modification systems as selfish mobile elements and their impact on genome evolution

Restriction-modification (RM) systems are composed of genes that encode a restriction enzyme and a modification methylase. RM systems sometimes behave as discrete units of life, like viruses and transposons. RM complexes attack invading DNA that has not been properly modified and thus may serve as a tool of defense for bacterial cells. However, any threat to their maintenance, such as a challenge by a competing genetic element (an incompatible plasmid or an allelic homologous stretch of DNA, for example) can lead to cell death through restriction breakage in the genome. This post-segregational or post-disturbance cell killing may provide the RM complexes (and any DNA linked with them) with a competitive advantage. There is evidence that they have undergone extensive horizontal transfer between genomes, as inferred from their sequence homology, codon usage bias and GC content difference. They are often linked with mobile genetic elements such as plasmids, viruses, transposons and integrons. The comparison of closely related bacterial genomes also suggests that, at times, RM genes themselves behave as mobile elements and cause genome rearrangements. Indeed some bacterial genomes that survived post-disturbance attack by an RM gene complex in the laboratory have experienced genome rearrangements. The avoidance of some restriction sites by bacterial genomes may result from selection by past restriction attacks. Both bacteriophages and bacteria also appear to use homologous recombination to cope with the selfish behavior of RM systems. RM systems compete with each other in several ways. One is competition for recognition sequences in post-segregational killing. Another is super-infection exclusion, that is, the killing of the cell carrying an RM system when it is infected with another RM system of the same regulatory specificity but of a different sequence specificity. The capacity of RM systems to act as selfish, mobile genetic elements may underlie the structure and function of RM enzymes.

Two recombination-dependent DNA replication pathways of bacteriophage T4, and their roles in mutagenesis and horizontal gene transfer

Two major pathways of recombination-dependent DNA replication, "join-copy" and "join-cut-copy," can be distinguished in phage T4: join-copy requires only early and middle genes, but two late proteins, endonuclease VII and terminase, are uniquely important in the join-cut-copy pathway. In wild-type T4, timing of these pathways is integrated with the developmental program and related to transcription and packaging of DNA. In primase mutants, which are defective in origin-dependent lagging-strand DNA synthesis, the late pathway can bypass the lack of primers for lagging-strand DNA synthesis. The exquisitely regulated synthesis of endo VII, and of two proteins from its gene, explains the delay of recombination-dependent DNA replication in primase (as well as topoisomerase) mutants, and the temperature-dependence of the delay. Other proteins (e.g., the single-stranded DNA binding protein and the products of genes 46 and 47) are important in all recombination pathways, but they interact differently with other proteins in different pathways. These homologous recombination pathways contribute to evolution because they facilitate acquisition of any foreign DNA with limited sequence homology during horizontal gene transfer, without requiring transposition or site-specific recombination functions. Partial heteroduplex repair can generate what appears to be multiple mutations from a single recombinational intermediate. The resulting sequence divergence generates barriers to formation of viable recombinants. The multiple sequence changes can also lead to erroneous estimates in phylogenetic analyses.

Role of genomic typing in taxonomy, evolutionary genetics, and microbial epidemiology

Currently, genetic typing of microorganisms is widely used in several major fields of microbiological research. Taxonomy, research aimed at elucidation of evolutionary dynamics or phylogenetic relationships, population genetics of microorganisms, and microbial epidemiology all rely on genetic typing data for discrimination between genotypes. Apart from being an essential component of these fundamental sciences, microbial typing clearly affects several areas of applied microbiological research. The epidemiological investigation of outbreaks of infectious diseases and the measurement of genetic diversity in relation to relevant biological properties such as pathogenicity, drug resistance, and biodegradation capacities are obvious examples. The diversity among nucleic acid molecules provides the basic information for all fields described above. However, researchers in various disciplines tend to use different vocabularies, a wide variety of different experimental methods to monitor genetic variation, and sometimes widely differing modes of data processing and interpretation. The aim of the present review is to summarize the technological and fundamental concepts used in microbial taxonomy, evolutionary genetics, and epidemiology. Information on the nomenclature used in the different fields of research is provided, descriptions of the diverse genetic typing procedures are presented, and examples of both conceptual and technological research developments for Escherichia coli are included. Recommendations for unification of the different fields through standardization of laboratory techniques are made.

Detecting anomalous gene clusters and pathogenicity islands in diverse bacterial genomes

A gene in a genome is defined as putative alien (pA) if its codon usage difference from the average gene exceeds a high threshold and codon usage differences from ribosomal protein genes, chaperone genes and protein-synthesis-processing factors are also high. pA gene clusters in bacterial genomes are relevant for detecting genomic islands (GIs), including pathogenicity islands (PAIs). Four other analyses appropriate to this task are G+C genome variation (the standard method); genomic signature divergences (dinucleotide bias); extremes of codon bias; and anomalies of amino acid usage. For example, the cagA domain of Helicobacter pylori is highly deviant in its genome signature and codon bias from the rest of the genome. Using these methods we can detect two potential PAIs in the Neisseria meningitidis genome, which contain hemagglutinin and/or hemolysin-related genes. Additionally, G+C variation and genome signature differences of the Mycobacterium tuberculosis genome indicate two pA gene clusters.

Evolutionary Role of Restriction/Modification Systems as Revealed by Comparative Genome Analysis

Type II restriction modification systems (RMSs) have been regarded either as defense tools or as molecular parasites of bacteria. We extensively analyzed their evolutionary role from the study of their impact in the complete genomes of 26 bacteria and 35 phages in terms of palindrome avoidance. This analysis reveals that palindrome avoidance is not universally spread among bacterial species and that it does not correlate with taxonomic proximity. Palindrome avoidance is also not universal among bacteriophage, even when their hosts code for RMSs, and depends strongly on the genetic material of the phage. Interestingly, palindrome avoidance is intimately correlated with the infective behavior of the phage. We observe that the degree of palindrome and restriction site avoidance is significantly and consistently less important in phages than in their bacterial hosts. This result brings to the fore a larger selective load for palindrome and restriction site avoidance on the bacterial hosts than on their infecting phages. It is then consistent with a view where type II RMSs are considered as parasites possibly at the verge of mutualism. As a consequence, RMSs constitute a nontrivial third player in the host–parasite relationship between bacteria and phages.

Control of directionality in integrase-mediated recombination: examination of recombination directionality factors (RDFs) including Xis and Cox proteins

Similarity between the DNA substrates and products of integrase-mediated site-specific recombination reactions results in a single recombinase enzyme being able to catalyze both the integration and excision reactions. The control of directionality in these reactions is achieved through a class of small accessory factors that favor one reaction while interfering with the other. These proteins, which we will refer to collectively as recombination directionality factors (RDFs), play architectural roles in reactions catalyzed by their cognate recombinases and have been identified in conjunction with both tyrosine and serine integrases. Previously identified RDFs are typically small, basic and have diverse amino acid sequences. A subset of RDFs, the cox genes, also function as transcriptional regulators. We present here a compilation of all the known RDF proteins as well as those identified through database mining that we predict to be involved in conferring recombination directionality. Analysis of this group of proteins shows that they can be grouped into distinct sub-groups based on their sequence similarities and that they are likely to have arisen from several independent evolutionary lineages. This compilation will prove useful in recognizing new proteins that confer directionality upon site-specific recombination reactions encoded by plasmids, transposons, phages and prophages.