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.     

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.     

Conservation and diversity in the immunity regions of wild phages with the immunity specificity of phage lambda

The gene regulatory circuitry of phage lambda is among the best-understood circuits. Much of the circuitry centres around the immunity region, which includes genes for two repressors, CI and Cro, and their cis-acting sites. Related phages, termed lambdoid phages, have different immunity regions, but similar regulatory circuitry and genome organization to that of lambda, and show a mosaic organization, arising by recombination between lambdoid phages. We sequenced the immunity regions of several wild phages with the immunity specificity of lambda, both to determine whether natural variation exists in regulation, and to analyse conservation and variability in a region rich in well-studied regulatory elements. CI, Cro and their cis-acting sites are almost identical to those in lambda, implying that regulatory mechanisms controlled by the immunity region are conserved. A segment adjacent to one of the operator regions is also conserved, and may be a novel regulatory element. In most isolates, different alleles of two regulatory proteins (N and CII) flank the immunity region; possibly the lysis-lysogeny decision is more variable among isolates. Extensive mosaicism was observed for several elements flanking the immunity region. Very short sequence elements or microhomologies were also identified. Our findings suggest mechanisms by which fine-scale mosaicism arises.     

Development of genetic circuitry exhibiting toggle switch or oscillatory behavior in Escherichia coli

Analysis of the system design principles of signaling systems requires model systems where all components and regulatory interactions are known. Components of the Lac and Ntr systems were used to construct genetic circuits that display toggle switch or oscillatory behavior. Both devices contain an "activator module" consisting of a modified glnA promoter with lac operators, driving the expression of the activator, NRI. Since NRI activates the glnA promoter, this creates an autoactivated circuit repressible by LacI. The oscillator contains a "repressor module" consisting of the NRI-activated glnK promoter driving LacI expression. This circuitry produced synchronous damped oscillations in turbidostat cultures, with periods much longer than the cell cycle. For the toggle switch, LacI was provided constitutively; the level of active repressor was controlled by using a lacY mutant and varying the concentration of IPTG. This circuitry provided nearly discontinuous expression of activator.     

A novel biotinylated degradable polymer for cell-interactive applications

We previously demonstrated that the lambda system integrated into the host chromosome can overcome the instability encountered in continuous operations of unstable plasmid-based expression vectors. High stability of a cloned gene in a lysogenic state and a high copy number in a lytic state provide cloned-gene stability and overexpression in a two-stage continuous operation. But the expression by the commonly used S- mutant lambda was only twice as high as that of the single copy. To increase the expression in the lambda system, we constructed a Q- mutant lambda vector that can be used in long-term operations such as a two-stage continuous operation. The Q- mutant phage lambda is deficient in the synthesis of proteins involved in cell lysis and lambda DNA packaging, while the S- mutant is deficient in the synthesis of one of two phage proteins required for lysis of the host cell and liberation of the progeny phage. Therefore, it is expected that the replicated Q- lambda DNA containing a cloned gene would not be coated by a phage head and would remain naked for ample expression of the cloned gene and host cells would not lyse easily and consequently would produce larger amounts of cloned-gene products. The beta-galactosidase expression per unit cell by the Q- mutant in a lytic state was about 30 times higher than that in a lysogenic state, while the expression by the commonly used S- mutant in a lytic state was twice as high as that in a lysogenic state. The optimal switching time of the Q- mutant from the lysogenic state to the lytic state for the maximum production of beta-galactosidase was 5.3 h, which corresponds to an early log phase in the batch operation.     

Optimizing genetic circuits by global sensitivity analysis

Artificial genetic circuits are becoming important tools for controlling cellular behavior and studying molecular biosystems. To genetically optimize the properties of complex circuits in a practically feasible fashion, it is necessary to identify the best genes and/or their regulatory components as mutation targets to avoid the mutation experiments being wasted on ineffective regions, but this goal is generally not achievable by current methods. The Random Sampling-High Dimensional Model Representation (RS-HDMR) algorithm is employed in this work as a global sensitivity analysis technique to estimate the sensitivities of the circuit properties with respect to the circuit model parameters, such as rate constants, without knowing the precise parameter values. The sensitivity information can then guide the selection of the optimal mutation targets and thereby reduce the laboratory effort. As a proof of principle, the in vivo effects of 16 pairwise mutations on the properties of a genetic inverter were compared against the RS-HDMR predictions, and the algorithm not only showed good consistency with laboratory results but also revealed useful information, such as different optimal mutation targets for optimizing different circuit properties, not available from previous experiments and modeling.     

Constraint and Contingency in Multifunctional Gene Regulatory Circuits

Gene regulatory circuits drive the development, physiology, and behavior of organisms from bacteria to humans. The phenotypes or functions of such circuits are embodied in the gene expression patterns they form. Regulatory circuits are typically multifunctional, forming distinct gene expression patterns in different embryonic stages, tissues, or physiological states. Any one circuit with a single function can be realized by many different regulatory genotypes. Multifunctionality presumably constrains this number, but we do not know to what extent. We here exhaustively characterize a genotype space harboring millions of model regulatory circuits and all their possible functions. As a circuit''s number of functions increases, the number of genotypes with a given number of functions decreases exponentially but can remain very large for a modest number of functions. However, the sets of circuits that can form any one set of functions becomes increasingly fragmented. As a result, historical contingency becomes widespread in circuits with many functions. Whether a circuit can acquire an additional function in the course of its evolution becomes increasingly dependent on the function it already has. Circuits with many functions also become increasingly brittle and sensitive to mutation. These observations are generic properties of a broad class of circuits and independent of any one circuit genotype or phenotype.     

The neural circuits and synaptic mechanisms underlying motor initiation in C. elegans

C. elegans is widely used to dissect how neural circuits and genes generate behavior. During locomotion, worms initiate backward movement to change locomotion direction spontaneously or in response to sensory cues; however, the underlying neural circuits are not well defined. We applied a multidisciplinary approach to map neural circuits in freely behaving worms by integrating functional imaging, optogenetic interrogation, genetic manipulation, laser ablation, and electrophysiology. We found that a disinhibitory circuit and a stimulatory circuit together promote initiation of backward movement and that circuitry dynamics is differentially regulated by sensory cues. Both circuits require glutamatergic transmission but depend on distinct glutamate receptors. This dual mode of motor initiation control is found in mammals, suggesting that distantly related organisms with anatomically distinct nervous systems may adopt similar strategies for motor control. Additionally, our studies illustrate how a multidisciplinary approach facilitates dissection of circuit and synaptic mechanisms underlying behavior in a genetic model organism.     

Stimulation of the lysogenization of Pseudomonas aeruginosa cells by phage transposon B39 induced by plasmid Rms163

The influence of Rms163 plasmid on lysogenization of Pseudomonas aeruginosa cells by B39 phage was studied. Plasmid Rms163 was shown to increase the frequency of lysogenization of PAO1 cells 7-8 times. C-mutants of B39 phage were isolated. According to complementation test, c-mutants were distributed into two groups--cI and cII/III. The product of cI is essential for establishment and maintenance of lysogenic state, cII/cIII product being only necessary for establishment of lysogenization. The mutants with special characteristics were isolated: B39cx1 phage carries a mutation which seems to be located on a regulatory site essential for establishment of lysogenic state. The region of the B39 genome responsible for interaction with Rms163 plasmid was mapped. Possible mechanisms of Rms163 plasmid interference with transposable B39 phage are discussed.     

Neural model of the genetic network

Many cell control processes consist of networks of interacting elements that affect the state of each other over time. Such an arrangement resembles the principles of artificial neural networks, in which the state of a particular node depends on the combination of the states of other neurons. The lambda bacteriophage lysis/lysogeny decision circuit can be represented by such a network. It is used here as a model for testing the validity of a neural approach to the analysis of genetic networks. The model considers multigenic regulation including positive and negative feedback. It is used to simulate the dynamics of the lambda phage regulatory system; the results are compared with experimental observation. The comparison proves that the neural network model describes behavior of the system in full agreement with experiments; moreover, it predicts its function in experimentally inaccessible situations and explains the experimental observations. The application of the principles of neural networks to the cell control system leads to conclusions about the stability and redundancy of genetic networks and the cell functionality. Reverse engineering of the biochemical pathways from proteomics and DNA micro array data using the suggested neural network model is discussed.     

Analysis of noise in quorum sensing

Noise may play a pivotal role in gene circuit functionality, as demonstrated for the genetic switch in the bacterial phage lambda. Like the lambda switch, bacterial quorum sensing (QS) systems operate within a population and contain a bistable switching element, making it likely that noise plays a functional role in QS circuit operation. Therefore, a detailed analysis of the noise behavior of QS systems is needed. We have developed a set of tools generally applicable to the analysis of gene circuits, with an emphasis on investigations in the frequency domain (FD), that we apply here to the QS system in the marine bacterium Vibrio fischeri. We demonstrate that a tight coupling between exact stochastic simulation and FD analysis provides insights into the structure/function relationships in the QS circuit. Furthermore, we argue that a noise analysis is incomplete without consideration of the power spectral densities (PSDs) of the important molecular output signals. As an example we consider reversible reactions in the QS circuit, and show through analysis and exact stochastic simulation that these circuits make significant and dynamic modifications to the noise spectra. In particular, we demonstrate a "whitening" effect, which occurs as the noise is processed through these reversible reactions.     

Identification and characterization of mutations in Escherichia coli that selectively influence the growth of hybrid lambda bacteriophages carrying the immunity region of bacteriophage P22.

Mutations in two Escherichia coli genes, sipA and sipB, result in a specific inhibition of the growth of certain hybrid lambdoid bacteriophages, lambda immP22, that have the early regulatory regions and adjacent genes from bacteriophage P22. The sipB391 mutation maps near minute 56 and exerts the strongest inhibitory effect on the growth of the hybrid phages. The sipA1 mutation maps near minute 72 and plays an auxiliary role: enhancing the action of sipB391. Such a role is not limited to sipA1, since there is a similar enhancement by the nusA1 and nusE71 mutations. The Sip-imposed restriction on the growth of lambda immP22 phages is not observed if the phage carries a mutation in the c1 gene. Perhaps this reflects the fact that the c1 product regulates phage DNA replication and is a major determinant in the decision governing whether the phage takes the lytic or lysogenic pathway. Consistent with this idea is the observation that lambda immP22 DNA replication is severely inhibited in bacteria carrying the sipB391 mutation. It is suggested that sip mutations exaggerate the normal role of c1 in limiting lytic growth. This causes a failure in the expression of sufficient amounts of some or all of the lytic gene products required for phage growth.     

The P1 ParA protein and its ATPase activity play a direct role in the segregation of plasmid copies to daughter cells

The P1 ParA protein is an ATPase that recognizes the parA promoter region where it acts to autoregulate the P1 parA–parB operon. The ParB protein is essential for plasmid partition and recognizes the cis -acting partition site parS . The regulatory role of ParA is also essential because a controlled level of ParB protein is critical for partition. However, we show that this regulatory activity is not the only role for ParA in partition. Efficient partition can be achieved without autoregulation as long as Par protein levels are kept within a range of low values. The properties of ParA mutants in these conditions showed that ParA is essential for some critical step in the partition process that is independent of par operon regulation. The putative nucleotide-binding site for the ParA ATPase was identified and disrupted by mutation. The resulting mutant was substantially defective for autoregulation and completely inactive for partition in a system in which the need for autoregulation is abolished. Thus, the ParA nucleotide-binding site appears to be necessary both for the repressor activity of ParA and for some essential step in the partition process itself. We propose that the nucleotide-bound form of the enzyme adopts a configuration that favours binding to the operator, but that the ATPase activity of ParA is required for some energetic step in partition of the plasmid copies to daughter cells.     

Biological pattern generation: the cellular and computational logic of networks in motion

In 1900, Ramón y Cajal advanced the neuron doctrine, defining the neuron as the fundamental signaling unit of the nervous system. Over a century later, neurobiologists address the circuit doctrine: the logic of the core units of neuronal circuitry that control animal behavior. These are circuits that can be called into action for perceptual, conceptual, and motor tasks, and we now need to understand whether there are coherent and overriding principles that govern the design and function of these modules. The discovery of central motor programs has provided crucial insight into the logic of one prototypic set of neural circuits: those that generate motor patterns. In this review, I discuss the mode of operation of these pattern generator networks and consider the neural mechanisms through which they are selected and activated. In addition, I will outline the utility of computational models in analysis of the dynamic actions of these motor networks.     

Plasmid formation from bacteriophage lambda in Escherichia coli

Summary Among the survivors of Escherichia coli derivatives infected with phage c1 or vir that are unable to establish ordinal lysogeny, clones arise which perpetuate the nondefective phage genome. When the bacteria bears a mutation(s) that makes the cell tolerant to the phage multiplication, such clones appear readily.The bacteria- complex was studied genetically and chemically, and it was concluded that the intact phage genomes, about two to four circular copies per bacterial chromosome, are perpetuated in bacterial cytoplasm as plasmids or in lysogenic state in cytoplasm.Several lines of evidence suggests that the phage genome in the lysogenic state in cytoplasm is under a different regulatory system from that in the normal prophage state on chromosome.     

Induction and repression of the Drosophila Sgs-3 glue gene are mediated by distinct sequences in the proximal promoter.

The normal developmental expression of the Drosophila salivary gland secretion protein gene Sgs-3 requires the interaction of a distal and proximal regulatory element. A deletion/replacement analysis of the proximal promoter in stably transformed lines shows that induction of an Sgs-3/Adh fusion gene is normal if sequences from +10 to -50 are replaced by those of the hsp70 gene. Sequences between -98 and -50 are necessary for this expression but there is internal redundancy within this region as two distinct upstream sequences of 18 and 22 bp respectively are sufficient for stage- and tissue-specific expression, albeit at reduced levels. A point mutation at -53 eliminates the ecdysone-mediated repression of the Sgs-3 promoter at pupariation. We report mosaicisms of expression within the salivary gland for a number of stably transformed lines.