Sequence tolerance of the phage lambda PRM promoter: implications for evolution of gene regulatory circuitry

Much of the gene regulatory circuitry of phage lambda centers on a complex region called the O(R) region. This approximately 100-bp region is densely packed with regulatory sites, including two promoters and three repressor-binding sites. The dense packing of this region is likely to impose severe constraints on its ability to change during evolution, raising the question of how the specific arrangement of sites and their exact sequences could evolve to their present form. Here we ask whether the sequence of a cis-acting site can be widely varied while retaining its function; if it can, evolution could proceed by a larger number of paths. To help address this question, we developed a lambda cloning vector that allowed us to clone fragments spanning the O(R) region. By using this vector, we carried out intensive mutagenesis of the P(RM) promoter, which drives expression of CI repressor and is activated by CI itself. We made a pool of fragments in which 8 of the 12 positions in the -35 and -10 regions were randomized and cloned this pool into the vector, making a pool of P(RM) variant phage. About 10% of the P(RM) variants were able to lysogenize, suggesting that the lambda regulatory circuitry is compatible with a wide range of P(RM) sequences. Analysis of several of these phages indicated a range of behaviors in prophage induction. Several isolates had induction properties similar to those of the wild type, and their promoters resembled the wild type in their responses to CI. We term this property of different sequences allowing roughly equivalent function "sequence tolerance " and discuss its role in the evolution of gene regulatory circuitry.     

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

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

A comparative study of the immunity region of lambdoid phages including Shiga-toxin-converting phages: molecular basis for cross immunity

Comparison of eight lambdoid phages, including three Shiga-toxin converting phages, has been carried out with respect to the immunity region, especially the recognition helices of their repressor and CRO proteins on the one hand, and operator sequences on the other. Some as yet unassigned components of the regulatory circuits have been inferred by computer search. The cross immunity phenomenon shown by phages VT2-Sa and lambda is explained on the basis of similarity in their sequences. In addition, the similarity of 933W and HK022 in the sequences of their recognition helices of repressor and CRO, on the one hand, and operators, on the other, has led us to predict that they will have identical or similar immunity specificity. This homology has enabled us also to locate the OL (and consequently PL) of phage 933W that has been thought to be non-existent.     

Characterization of the primary immunity region of the Escherichia coli linear plasmid prophage N15.

N15 is the only bacteriophage of Escherichia coli known to lysogenize as a linear plasmid. Clear-plaque mutations lie in at least two regions of the 46-kb genome. We have cloned, sequenced, and characterized the primary immunity region, immB. This region contains a gene, cB, whose product shows homology to lambdoid phage repressors. The cB3 mutation confers thermoinducibility on N15 lysogens, consistent with CB being the primary repressor of N15. Downstream of cB lies the locus of N15 plasmid replication. Upstream of cB lies an operon predicted to encode two products: one homologous to the late repressor of P22 (Cro), the other homologous to the late antiterminator of phi 82 (Q). The Q-like protein is essential for phage development. We show that CB protein regulates the expression of genes that flank the cB gene by binding to DNA at symmetric 16-bp sites. Three sites are clustered upstream of cB and overlap a predicted promoter of the cro and Q-like genes as well as two predicted promoters of cB itself. Two sites downstream of cB overlap a predicted promoter of a plasmid replication gene, repA, consistent with the higher copy number of the mutant, N15cB3. The leader region of repA contains terminators in both orientations and a putative promoter. The organization of these regulatory elements suggests that N15 plasmid replication is controlled not only by CB but also by an antisense RNA and by a balance between termination and antitermination.     

Robustness of a gene regulatory circuit.

Complex interacting systems exhibit system behavior that is often not predictable from the properties of the component parts. We have tested a particular system property, that of robustness. The behavior of a system is termed robust if that behavior is qualitatively normal in the face of substantial changes to the system components. Here we test whether the behavior of the phage lambda gene regulatory circuitry is robust. This circuitry can exist in two alternative patterns of gene expression, and can switch from one regulatory state to the other. These states are stabilized by the action at the O(R) region of two regulatory proteins, CI and Cro, which bind with differential affinities to the O(R)1 and O(R)3 sites, such that each represses the synthesis of the other one. In this work, this pattern of binding was altered by making three mutant phages in which O(R)1 and O(R)3 were identical. These variants had the same qualitative in vivo patterns of gene expression as wild type. We conclude that the behavior of the lambda circuitry is highly robust. Based on these and other results, we propose a two-step pathway, in which robustness plays a key role, for evolution of complex regulatory circuitry.     

Molecular characterization and module composition of P22-related Salmonella phage genomes.

Genomes of newly isolated Salmonella phages were analysed by comparison of their EcoRI restriction patterns and by hybridization. Characteristic hybridization probes from reference phages P22, ES18 and E. coli phage lambda were chosen. Four probes selected from the lysis region examined the dispersal of the lambdoid lysis genes. Other probes characterized were the replication genes and part of the structural genes. The complex immunity region was investigated by means of hybridization as well as biological tests. The results showed the relationship of the isolated phages to the P22 branch of the lambdoid phages and revealed their modular genome organization consisting of different proportions of P22-related sequences. DNA restriction patterns of phages released from Salmonella strains sampled in limited geographical areas were significantly less heterogeneous than those of phages released from the worldwide sampled SARA collection. The use of prophage restriction patterns as a tool for the typing of Salmonellae to support the epidemiologic classification of pathogenic strains is discussed.     

Genetics of bacteriophage phi 80--a review

The genetic maps of bacteriophage lambda and lambdoid phage phi 80 are compared. The gene organization of phi 80 is very similar to that of lambda, as shown by isolation and characterization of many am, ts and c (clear) mutants of the phage. In general, the essential genes located in the same position on the genetic map of the phages lambda and phi 80 fulfill the same functions. These include the gene clusters coding for the head and tail proteins, genes for DNA synthesis, and the genes controlling lysogeny and late gene expression. The specific regulatory features of phi 80 in relation to the N function of lambda are discussed, but they require further clarification. The two phages differ in immunity specificity, host range, conversion property and temperature sensitivity.     

Structural organization of the ori site of phage P22; comparison with other lambdoid ori sites

The homologous DNA regions of phages P22, lambda and lambdoid coliphages, which code for the amino-terminal portion of genes 18 or O, contain the ori signal. Both the lambdoid and P22 ori regions can be divided into sections, A, B and C. The four direct repeats with internal rotational symmetry of section A in P22 are less regularly organized than in the corresponding signals of the phi 80 and lambda ori sites and show greatest homology to coliphage phi 82. Section B is rich in adenines in the l strand, and section C can be recognized in the P22 ori by the occurrence of overlapping inverted repeats. The latter region is not homologous to the structurally similar section C, 'EcoRI-loop', of the lambdoid coliphages. The results further define the specificity determinants of lambdoid O protein-ori interactions and demonstrate the evolutionary relationship between these functional units.     

Sequence organization of the origins of DNA replication in lambdoid coliphages

We have determined the sequences of the ori region DNA of several phage lambda mutants and hybrids, which shed light on the mechanism of DNA replication in the lambdoid phages. These include the heterologous substitution hybrids lambda rep82:lambda and lambda rep80:lambda, a pseudorevertant of the ori-r93 mutant lambda r93hot5, and the insertion mutant lambda pk35. The ori regions of the three lambdoid phages, lambda, phi 80 and 82, all have repeated sequences, termed iterons, and A . T-rich zones. We note that a similar arrangement of DNA is also found in several other prokaryotic origins of replication. lambda and phi 80 have four iterons, and 82 has five. The origin of lambda r93hot5 is unusual in that contains only three iterons, yet the phage grows normally. Analysis of this mutant indicates that the spacing of iterons is crucial to ori function, whereas their number is not. This argues against the cloverleaf model for lambda ori structure (Hobom et al., 1979). In lambda pk35 the drug resistance element Tn903 is inserted into the "inceptor" (ice) site, proposed to be crucial for lambda replication initiation (Hobom et al., 1979); yet this phage grows normally.     

Sequence analysis of Stx2-converting phage VT2-Sa shows a great divergence in early regulation and replication regions.

In enterohemorrhagic Escherichia coli, Shiga toxin is produced by lysogenic prophages. We have isolated the prophage VT2-Sa that is responsible for production of Shiga toxin type 2 protein, and determined the complete nucleotide sequence of this phage DNA. The entire DNA sequence consisted of 60,942 bp, exhibiting marked similarity to the 933W phage genome. However, several differences were observed in the immunity and replication regions, where cI, cII, cIII, N, cro, O, and P genes were present: Predicted amino acid sequences of N, cI, cro, O and P in the VT2-Sa genome did not show significant similarity to the counterparts of the 933W genome; however its cI showed higher similarity to lambda. Furthermore, O and P closely resembled those of phage HK022. These observations suggest that the various degrees of homology observed in the immunity and replication regions of VT2-Sa could have resulted from frequent recombination events among the lambdoid phages, and that these regions play a key role as a functional unit for phage propagation in competition with other lambdoid phages.     

Complete nucleotide sequence of the prophage VT1-Sakai carrying the Shiga toxin 1 genes of the enterohemorrhagic Escherichia coli O157:H7 strain derived from the Sakai outbreak

Shiga toxins 1 and 2 (Stx1 and Stx2) are encoded by prophages lysogenized in enterohemorrhagic Escherichia coli (EHEC) O157:H7 strains. Lytic growth of the phage particles carrying the stx1 genes (stx1A and stx1B) of the EHEC O157:H7 strain RIMD 0509952, which was derived from the Sakai outbreak in 1996 in Japan, was induced after treatment with mitomycin C, but the plaque formation of the phage was not detected. We have determined the complete nucleotide sequence of the prophage VT1-Sakai. The integration site of the prophage was identified within the yehV gene at 47.7 min on the chromosome. The stx1 genes were downstream of the Q gene in the prophage genome, suggesting that their expression was regulated by the Q protein, the regulator of the late gene expression of the phage, which is similar to that of the stx1 or stx2 genes carried by the lambdoid phages reported previously. The sequences of the N gene and its recognition sites, nutL and nutR, were not homologous to those of the phages carrying the stx genes thus far reported, but they were very similar to those of bacteriophage phi21. The sequences of the repressor proteins, CI and Cro, that regulate expression of the early genes had low similarities with those of the known repressors of other phages, and their operator sequences were different from any sequence reported. These data suggest that multiple genetic recombination among bacteriophages with different immunities took place to generate the prophage VT1-Sakai. Comparison between the sequences of VT1-Sakai and lambda suggests that the ancestor of VT1-Sakai was produced by illegitimate excision, like lambda gal and bio phages.