Gene network visualization and quantitative synteny analysis of more than 300 marine T4-like phage scaffolds from the GOS metagenome.

Bacteriophages (phages) are the most abundant biological entities in the biosphere and are the dominant "organisms" in marine environments, exerting an enormous influence on marine microbial populations. Metagenomic projects, such as the Global Ocean Sampling expedition (GOS), have demonstrated the predominance of tailed phages (Caudovirales), particularly T4 superfamily cyanophages (Cyano-T4s), in the marine milieu. Whereas previous metagenomic analyses were limited to gene content information, here we present a comparative analysis of over 300 phage scaffolds assembled from the viral fraction of the GOS data. This assembly permits the examination of synteny (organization) of the genes on the scaffolds and their comparison with the genome sequences from cultured Cyano-T4s. We employ comparative genomics and a novel usage of network visualization software to show that the scaffold phylogenies are similar to those of the traditional marker genes they contain. Importantly, these uncultured metagenomic scaffolds quite closely match the organization of the "core genome" of the known Cyano-T4s. This indicates that the current view of genome architecture in the Cyano-T4s is not seriously biased by being based on a small number of cultured phages, and we can be confident that they accurately reflect the diverse population of such viruses in marine surface waters.

The immense journey of bacteriophage T4--from d'Hérelle to Delbrück and then to Darwin and beyond.

In spite of their importance, the genomics, diversity and evolution of phages and their impact on the biosphere have remained largely unexplored research domains in microbiology. Here, we report on some recent studies with the T4 phage superfamily that shed some new light on these topics.

Plasticity of the gene functions for DNA replication in the T4-like phages.

We have completely sequenced and annotated the genomes of several relatives of the bacteriophage T4, including three coliphages (RB43, RB49 and RB69), three Aeromonas salmonicida phages (44RR2.8t, 25 and 31) and one Aeromonas hydrophila phage (Aeh1). In addition, we have partially sequenced and annotated the T4-like genomes of coliphage RB16 (a close relative of RB43), A. salmonicida phage 65, Acinetobacter johnsonii phage 133 and Vibrio natriegens phage nt-1. Each of these phage genomes exhibited a unique sequence that distinguished it from its relatives, although there were examples of genomes that are very similar to each other. As a group the phages compared here diverge from one another by several criteria, including (a) host range, (b) genome size in the range between approximately 160 kb and approximately 250 kb, (c) content and genetic organization of their T4-like genes for DNA metabolism, (d) mutational drift of the predicted T4-like gene products and their regulatory sites and (e) content of open-reading frames that have no counterparts in T4 or other known organisms (novel ORFs). We have observed a number of DNA rearrangements of the T4 genome type, some exhibiting proximity to putative homing endonuclease genes. Also, we cite and discuss examples of sequence divergence in the predicted sites for protein-protein and protein-nucleic acid interactions of homologues of the T4 DNA replication proteins, with emphasis on the diversity in sequence, molecular form and regulation of the phage-encoded DNA polymerase, gp43. Five of the sequenced phage genomes are predicted to encode split forms of this polymerase. Our studies suggest that the modular construction and plasticity of the T4 genome type and several of its replication proteins may offer resilience to mutation, including DNA rearrangements, and facilitate the adaptation of T4-like phages to different bacterial hosts in nature.

Genome plasticity in the distal tail fiber locus of the T-even bacteriophage: recombination between conserved motifs swaps adhesin specificity.

The adsorption specificity of the T-even phages is determined by the protein sequence near the tip of the long tail fibers. These adhesin sequences are highly variable in both their sequence and specificity for bacterial receptors. The tail fiber adhesin domains are located in different genes in closely related phages of the T-even type. In phage T4, the adhesin sequence is encoded by the C-terminal domain of the large tail fiber gene (gene 37), but in T2, the adhesin is a separate gene product (gene 38) that binds to the tip of T2 tail fibers. Analysis of phage T6 and Ac3 sequences reveals additional variant forms of this locus. The tail fiber host specificity determinants can be exchanged, although the different loci have only limited homology. Chimeric fibers can be created by crossovers either between small homologies within the structural part of the fiber gene or in conserved motifs of the adhesin domain. For example, the T2 adhesin determinants are flanked by G-rich DNA motifs and exchanges involving these sequences can replace the specificity determinants. These features of the distal tail fiber loci genetically link their different forms and can mediate acquisition of diverse host range determinants, including those that allow it to cross species boundaries and infect taxonomically distant hosts.

The genome of the pseudo T-even bacteriophages, a diverse group that resembles T4.

Polymerase chain reaction analysis of a large collection of bacteriophages with T-even morphology revealed four phages that are distantly related to all the others. The genomes of these pseudo T-even phages hybridized under stringent conditions to only a limited portion of the T4 genome that encodes virus head, head-to-tail joining and contractile tail genes. Except for this region, no extensive hybridization was detected between most pairs of the different pseudo T-even genomes. Sequencing of this conserved region of the pseudo T-even phage RB49 revealed substantial nucleotide sequence divergence from T4 (approximately 30% to 40%), and random genomic sequencing of this phage indicated that more than a third of its sequences had no detectable homology to T4. Among those sequences related to the T-even genes were virion structural components including the constituents of the phage base plate. Only a few sequences had homology to T4 early functions; these included ribonucleotide diphosphatase reductase, DNA ligase and the large subunit of DNA topoisomerase. The genomes of the pseudo T-even phage were digested by restriction enzymes that are unable to digest the T-even DNAs which contain glucosylated hydroxymethyl-cytosine residues. This suggests that only limited nucleotide modifications must be present in the pseudo T-even genomes. Conservation of much of the morphogenetic region of these diverse phage genomes may reflect particularly strong sequence constraints on these gene products. However, other explanations are considered, including the possibility that the various morphogenetic segments were acquired by the pseudo T-even genomes by modular evolution. These results support the notion that phage evolution may proceed within a network of both closely and distantly related genomes.

Bacteriophage T4 host range is expanded by duplications of a small domain of the tail fiber adhesin.

The adsorption specificity of T4 is determined by the tip of the gene 37 tail fibers which bind to receptors on the bacterial surface. T4 infects only Escherichia coli and closely related Shigella species, but rare host range mutants can be isolated that infect Yersinia pseudotuberculosis I, an evolutionally distant bacterium. Some of these mutations result in amino acid residue substitutions in the C-terminal portion of gene 37, but others involve unequal exchanges between a series of sequence motifs (His boxes) in the same region. The duplication or mutational alteration of this segment apparently suffices for phage adsorption to a Yersinia receptor. It is suggested that recombination between the His box sequences can generate diversity in phage host range by shuffling receptor recognition domains.

Sense and antisense transcription of bacteriophage T4 gene 32. Processing and stability of the mRNAs.

Analysis of bacteriophage T4 gene 32 transcription has revealed a multiplicity of mRNAs. In plasmids, gene 32 is expressed primarily from a strong promoter that is shut off after phage infection. In a wild-type infection, gene 32 is initially transcribed from prereplicative polycistronic and monocistronic promoters; subsequently, a monocistronic late mRNA predominates. This transcript, as well as a post-transcriptionally processed product of the earlier mRNA, can be stable. The eventual degradation of the stable mRNAs is temporally regulated by the phage. Finally, the transcription termination region of gene 32 can function as an antisense promoter both in vitro and in vivo.

Non-nearest neighbor effects on the thermodynamics of unfolding of a model mRNA pseudoknot.

The upstream autoregulatory mRNA leader sequence of gene 32 of 17 T-even and related bacteriophages folds into a simple tertiary structural motif, a hairpin-type RNA pseudoknot. In phage T4, the pseudoknot is contained within 28 contiguous nucleotides which adopt a pseudocontinuous helical structure derived from two coaxially stacked helical stems of four (stem 1) and seven (stem 2) base-pairs connected by two inequivalent single-stranded loops of five and one nucleotide(s). These two loops cross the minor and major grooves of stems 1 and 2, respectively. In this study, the equilibrium unfolding pathway of a 35-nucleotide RNA fragment corresponding to the wild-type and sequence variants of the T4 gene 32 mRNA has been determined through analysis of dual-wave-length, equilibrium thermal melting profiles via application of a van't Hoff model based on multiple sequential, two-state transitions. The melting profile of the wild-type RNA is well-described by two sequential melting transitions over a wide range of magnesium concentration. Compensatory base-pair substitutions incorporated into helical stems 1 and 2 were used to assign the first low enthalpy, moderate tm melting transition to the denaturation of the short three to four base-pair stem 1, followed by unfolding of the larger seven base-pair stem 2. We find that loop 1 substitution mutants (A10 to G10, C10, U10 or GA10) strikingly uncouple the melting of stems 1 and 2, with the U10 substitution and the GA10 loop expansion more destabilizing than the G10 and C10 substitutions. A significant increase in the extent of cleavage by RNase T1 following the conserved G26 (the 3' nucleotide in loop 2) in the U10, G10, and GA10 mutants suggests that an altered helix-helix junction region in this mutant may be responsible, at least in part, for this uncoupling. In addition to a modest destabilization of stem 2, the major effect of deletion or nucleotide substitution in the 3' single-stranded tail is a destabilization of stem 1, a non-nearest neighbor tertiary structural effect, which may well be transmitted through an altered loop 1-core helix interaction. In contrast, truncation of the 5' tail has no effect on the stability of the molecule.

Polynucleotide ligase in bacteriophage T4D recombination.

Following infection of E. coli B with ligase-deficient rII bacteriophage T4D recombination between linked markers is increased 4.2 fold and heterozygote frequency increased 2.3 fold. In such infection recombination occurs at a rapid rate for an extended period. This is in contrast to the time course of recombination observed in wild-type, lysis-inhibited, or lysis-defective (gene t defective) infection. In all of these cases recombination under standard cross conditions occurs early in the vegetative cycle. The increased recombination in ligase-deficient rII infection is reduced in a bacterial strain which produces greater than normal levels of host ligase. These results indicate that ligase has a crucial role not only in the replication of DNA but also in recombination. The level of ligase may determine whether DNA replication occurs with or without concomitant recombination.

A modified pBR322 vector with improved properties for the cloning, recovery, and sequencing of blunt-ended DNA fragments.

The construction of a plasmid vector which facilitates the cloning and recovery of blunt-ended DNA fragments is described. This plasmid, called pHP34, differs from pBR322 by a 10-bp insertion which introduces a unique SmaI site immediately flanked by two EcoRI sites. Blunt-ended DNA fragments cloned in the SmaI site can be recovered by digestion with EcoRI. Small cloned fragments can be chemically sequenced using a strategy which does not require their purification. The use of a plasmid related to pHP34 for in vitro mutagenesis by the insertion of a DNA linker fragment conferring an antibiotic resistance is also discussed.

In vitro insertional mutagenesis with a selectable DNA fragment.

A new method for in vitro insertional mutagenesis of genes cloned in Escherichia coli is presented. This simple procedure combines the advantages of in vitro DNA linker mutagenesis with those of in vivo transposition mutagenesis. It makes use of the omega fragment, a 2.0-kb DNA segment consisting of an antibiotic resistance gene (the Smr/Spcr gene of the R100.1 plasmid) flanked by short inverted repeats carrying transcription and translation termination signals and synthetic polylinkers. The omega fragment is inserted into a linearized plasmid by in vitro ligation, and the recombinant DNA molecules are selected by their resistance to streptomycin and spectinomycin. The omega fragment terminates RNA and protein synthesis prematurely, thus allowing the definition and mapping of both transcription and translation units. Because of the symmetrical structure of omega, the same effect is obtained with insertions in either orientation. The antibiotic resistance gene can be subsequently excised from the mutated molecules, leaving behind its flanking restriction site(s).

Omega mutagenesis in gram-negative bacteria: a selectable interposon which is strongly polar in a wide range of bacterial species.

We have used the 2.0-kb DNA fragment omega [Prentki and Krisch, Gene 29 (1984) 303-313] to mutagenize in vitro a broad-host-range plasmid carrying the entire meta-cleavage pathway of the Pseudomonas putida TOL plasmid pWW0. The mutant plasmids were subsequently introduced by conjugal mobilization into a variety of Gram-negative bacteria. The omega fragment carries a selectable marker (aadA+; SpcR/SmR), which is expressed in all species tested, as well as flanking transcription and translation termination signals and synthetic polylinkers. Expression of the plasmid-borne catechol 2,3-dioxygenase (C23O) gene, situated downstream from the site of omega insertion, was substantially reduced in all strains tested. The transcription terminators originally cloned from bacteriophage T4 gene 32, are apparently functional in a wide range of hosts. Insertional mutagenesis with the omega 'interposon' can thus be used in a wide variety of species, with the advantages of a positive selection for the presence of the fragment, the termination of RNA and protein synthesis beyond the site of insertion, and genetic stability of the resulting mutation.

A bacteriophage T4 expression cassette that functions efficiently in a wide range of gram-negative bacteria.

We have constructed a derivative of the broad-host-range vector RSF1010. This plasmid, p alpha omega, contains an expression cassette derived from bacteriophage T4 gene 32, into which we have inserted the coding sequence for the xylE enzyme (C2,3O) of the TOL plasmid pWWO. The composite plasmid, p alpha xylE omega, was transferred by conjugal mobilisation into a variety of Gram-negative bacteria (Agrobacter, Paracoccus, Erwinia, Pseudomonas, Rhizobium and Xanthomonas). High levels of C2,3O activity were found in almost all of the extracts. Polyacrylamide gel electrophoresis of these extracts revealed a prominent protein band at 35 kDa whose identity as the C2,3O gene product was confirmed by immunoblotting. We have mapped the 5' ends of the gene 32/xylE hybrid transcripts. In all of the Gram-negative bacteria, the proximal P2 promoter is the most efficient promoter in the cassette. In most of the strains a weaker and more distal promoter activity (Pl) was also detected. In both uninfected and phage-infected Escherichia coli cells, the transcript produced from this promoter is processed at a specific site upstream from the gene 32 start codon. The same processing occurred in all the bacterial species examined. The decay of the hybrid xylE transcript has been analyzed in E. coli and Erwinia, and in both strains this mRNA was among the most stable.

A plasmid expression vector that permits stabilization of both mRNAs and proteins encoded by the cloned genes.

Two new expression vectors have been constructed to take advantage of several useful properties of bacteriophage T4-infected Escherichia coli. These plasmids, pRDB8 and pRDB9, contain the promoter region and start codon of T4 gene 32, a contiguous multiple cloning site (MCS), and translation and transcription termination signals. DNA fragments inserted into the MCS are transcribed and translated at a high level in both uninfected and phage T4-infected cells. Furthermore, the extreme stability of the hybrid mRNA after infection permits the specific biosynthetic labeling of the protein encoded by the cloned gene. In addition, the cloned gene product is stabilized, since the host-mediated degradation of foreign proteins is inhibited by phage infection. The properties of this expression system were demonstrated with the constant region of a rabbit immunoglobulin lambda light chain (C lambda) gene. Although proteolytic degradation of the C lambda fusion protein was rapid in uninfected cells, degradation was blocked in phage-infected cells and the protein accumulated in greater amounts.

Omegon-Km: a transposable element designed for in vivo insertional mutagenesis and cloning of genes in gram-negative bacteria.

To combine the features of the omega interposons with the advantages of in vivo transposition mutagenesis, we have constructed an artificial transposon, called Omegon-Km. The Omegon-Km transposon is carried on the plasmid pJFF350 which can be conjugally mobilized into a broad range of Gram-negative bacteria. Omegon-Km is flanked, in inverted orientation, by synthetic 28-bp repeats derived from the ends of IS1. In addition, each end of Omegon-Km has the very efficient transcription and translation terminators of the omega interposon. Internally, Omegon-Km carries the selectable kanamycin (Km)-neomycin resistance gene (alph A) which is expressed well in many Gram-negative bacteria. The IS1 transposition functions are located on the donor plasmid but external to Omegon-Km. Thus, insertions of Omegon-Km are very stable because they lack the capacity for further transposition. Omegon-Km mutagenesis is performed by conjugal transfer of pJFF350 from Escherichia coli into any Gram-negative recipient strain in which this plasmid is unable to replicate. Those cells which have had a transposition event are selected by their resistance to Km. Very high frequencies of Omegon-Km transposition were observed in Pseudomonas putida. Preliminary experiments with other Gram-negative soil and water bacteria (Rhizobium leguminosarum, Paracoccus denitrificans) yielded mutants at reasonable levels. The presence of an E. coli-specific origin of replication (ori) within Omegon-Km allows the rapid and easy cloning, in E. coli, of the nucleotide sequences flanking the site of the transposition event.

Direct PCR sequencing of the ndd gene of bacteriophage T4: identification of a product involved in bacterial nucleoid disruption.

The rapid disruption of the Escherichia coli nucleoid after T4 infection requires the activity of the phage-encoded ndd gene. We have genetically identified the sequence encoding ndd. Determination of the sequence of a 2.5-kb segment including ndd closed the last significant gap in the sequence of the T4 genome. This analysis was performed on PCR-amplified fragments that were purified by gel-exclusion chromatography and then submitted to linear amplification cycle sequencing. This technology permitted sequence comparison of two ndd mutants (ndd44 and ndd98) with the wild-type gene. The analysis of ndd from six bacteriophages of the T-even family indicated that the protein encoded by this nonessential gene is surprisingly conserved.

Bacterial poly(A) polymerase: an enzyme that modulates RNA stability.

We have constructed a strain that overexpresses E coli poly(A) polymerase (PAP I). The recombinant protein was soluble, and a partially purified extract had high levels of poly(A) polymerising activity. An antiserum raised against the overexpressed PAP I has permitted two types of analysis: the identification of other E coli proteins that may interact with PAP I, and the search for PAP I-like proteins in other bacteria. Immunoprecipitation experiments suggest that PAP I is associated with a 48-kDa protein. This protein remains to be identified. Western blotting using the antiserum against E coli PAP I revealed related proteins in a variety of Gram-negative bacteria and in B subtilis. A comparison of the E coli protein with putative poly(A) polymerases recently identified in H influenza and B subtilis showed highly conserved sequences in the amino terminal and central portions of the proteins that may be important for enzyme activity.

The episodic evolution of fibritin: traces of ancient global environmental alterations may remain in the genomes of T4-like phages.

The evolutionary adaptation of bacteriophages to their environment is achieved by alterations of their genomes involving a combination of both point mutations and lateral gene transfer. A phylogenetic analysis of a large set of collar fiber protein (fibritin) loci from diverse T4-like phages indicates that nearly all the modular swapping involving the C-terminal domain of this gene occurred in the distant past and has since ceased. In phage T4, this fibritin domain encodes the sequence that mediates both the attachment of the long tail fibers to the virion and also controls, in an environmentally sensitive way, the phage's ability to infect its host bacteria. Subsequent to its distant period of modular exchange, the evolution of fibritin has proceeded primarily by the slow vertical divergence mechanism. We suggest that ancient and sudden changes in the environment forced the T4-like phages to alter fibritin's mode of action or function. The genome's response to such episodes of rapid environmental change could presumably only be achieved quickly enough by employing the modular evolution mechanism. A phylogenetic analysis of the fibritin locus reveals the possible traces of such events within the T4 superfamily's genomes.

The gp38 adhesins of the T4 superfamily: a complex modular determinant of the phage's host specificity.

The tail fiber adhesins are the primary determinants of host range in the T4-type bacteriophages. Among the indispensable virion components, the sequences of the long tail fiber genes and their associated adhesins are among the most variable. The predominant form of the adhesin in the T4-type phages is not even the version of the gene encoded by T4, the archetype of the superfamily, but rather a small unrelated protein (gp38) encoded by closely related phages such as T2 and T6. This gp38 adhesin has a modular design: its N-terminal attachment domain binds at the tip of the tail fiber, whereas the C-terminal specificity domain determines its host receptor affinity. This specificity domain has a series of four hypervariable segments (HVSs) that are separated by a set of highly conserved glycine-rich motifs (GRMs) that apparently form the domain's conserved structural core. The role of gp38's various components was examined by a comparative analysis of a large series of gp38 adhesins from T-even superfamily phages with differing host specificities. A deletion analysis revealed that the individual HVSs and GRMs are essential to the T6 adhesin's function and suggests that these different components all act in synergy to mediate adsorption. The evolutionary advantages of the modular design of the adhesin involving both conserved structural elements and multiple independent and easily interchanged specificity determinants are discussed.