Diversity in lac Operon Regulation among Diverse Escherichia coli Isolates Depends on the Broader Genetic Background but Is Not Explained by Genetic Relatedness.

Transcription of bacterial genes is controlled by the coordinated action of cis- and trans-acting regulators. The activity and mode of action of these regulators can reflect different requirements for gene products in different environments. A well-studied example is the regulatory function that integrates the environmental availability of glucose and lactose to control the Escherichia colilac operon. Most studies of lac operon regulation have focused on a few closely related strains. To determine the range of natural variation in lac regulatory function, we introduced a reporter construct into 23 diverse E. coli strains and measured expression with combinations of inducer concentrations. We found a wide range of regulatory functions. Several functions were similar to the one observed in a reference lab strain, whereas others depended weakly on the presence of cAMP. Some characteristics of the regulatory function were explained by the genetic relatedness of strains, indicating that differences varied on relatively short time scales. The regulatory characteristics explained by genetic relatedness were among those that best predicted the initial growth of strains following transition to a lactose environment, suggesting a role for selection. Finally, we transferred the lac operon, with the lacI regulatory gene, from five natural isolate strains into a reference lab strain. The regulatory function of these hybrid strains revealed the effect of local and global regulatory elements in controlling expression. Together, this work demonstrates that regulatory functions can be varied within a species and that there is variation within a species to best match a function to particular environments.IMPORTANCE The lac operon of Escherichia coli is a classic model for studying gene regulation. This study has uncovered features such as the environmental input logic controlling gene expression, as well as gene expression bistability and hysteresis. Most lac operon studies have focused on a few lab strains, and it is not known how generally those findings apply to the diversity of E. coli strains. We examined the environmental dependence of lac gene regulation in 20 natural isolates of E. coli and found a wide range of regulatory responses. By transferring lac genes from natural isolate strains into a common reference strain, we found that regulation depends on both the lac genes themselves and on the broader genetic background, indicating potential for still-greater regulatory diversity following horizontal gene transfer. Our results reveal that there is substantial natural variation in the regulation of the lac operon and indicate that this variation can be ecologically meaningful.

Lineage space and the propensity of bacterial cells to undergo growth transitions.

The molecular makeup of the offspring of a dividing cell gradually becomes phenotypically decorrelated from the parent cell by noise and regulatory mechanisms that amplify phenotypic heterogeneity. Such regulatory mechanisms form networks that contain thresholds between phenotypes. Populations of cells can be poised near the threshold so that a subset of the population probabilistically undergoes the phenotypic transition. We sought to characterize the diversity of bacterial populations around a growth-modulating threshold via analysis of the effect of non-genetic inheritance, similar to conditions that create antibiotic-tolerant persister cells and other examples of bet hedging. Using simulations and experimental lineage data in Escherichia coli, we present evidence that regulation of growth amplifies the dependence of growth arrest on cellular lineage, causing clusters of related cells undergo growth arrest in certain conditions. Our simulations predict that lineage correlations and the sensitivity of growth to changes in toxin levels coincide in a critical regime. Below the critical regime, the sizes of related growth arrested clusters are distributed exponentially, while in the critical regime clusters sizes are more likely to become large. Furthermore, phenotypic diversity can be nearly as high as possible near the critical regime, but for most parameter values it falls far below the theoretical limit. We conclude that lineage information is indispensable for understanding regulation of cellular growth.

Dynamical predictors of an imminent phenotypic switch in bacteria.

Single cells can stochastically switch across thresholds imposed by regulatory networks. Such thresholds can act as a tipping point, drastically changing global phenotypic states. In ecology and economics, imminent transitions across such tipping points can be predicted using dynamical early warning indicators. A typical example is 'flickering' of a fast variable, predicting a longer-lasting switch from a low to a high state or vice versa. Considering the different timescales between metabolite and protein fluctuations in bacteria, we hypothesized that metabolic early warning indicators predict imminent transitions across a network threshold caused by enzyme saturation. We used stochastic simulations to determine if flickering predicts phenotypic transitions, accounting for a variety of molecular physiological parameters, including enzyme affinity, burstiness of enzyme gene expression, homeostatic feedback, and rates of metabolic precursor influx. In most cases, we found that metabolic flickering rates are robustly peaked near the enzyme saturation threshold. The degree of fluctuation was amplified by product inhibition of the enzyme. We conclude that sensitivity to flickering in fast variables may be a possible natural or synthetic strategy to prepare physiological states for an imminent transition.

Cellular Growth Arrest and Persistence from Enzyme Saturation.

Metabolic efficiency depends on the balance between supply and demand of metabolites, which is sensitive to environmental and physiological fluctuations, or noise, causing shortages or surpluses in the metabolic pipeline. How cells can reliably optimize biomass production in the presence of metabolic fluctuations is a fundamental question that has not been fully answered. Here we use mathematical models to predict that enzyme saturation creates distinct regimes of cellular growth, including a phase of growth arrest resulting from toxicity of the metabolic process. Noise can drive entry of single cells into growth arrest while a fast-growing majority sustains the population. We confirmed these predictions by measuring the growth dynamics of Escherichia coli utilizing lactose as a sole carbon source. The predicted heterogeneous growth emerged at high lactose concentrations, and was associated with cell death and production of antibiotic-tolerant persister cells. These results suggest how metabolic networks may balance costs and benefits, with important implications for drug tolerance.

Survival of Phenotypic Information during Cellular Growth Transitions.

Phenotypic memory can predispose cells to physiological outcomes, contribute to heterogeneity in cellular populations, and allow computation of environmental features, such as nutrient gradients. In bacteria and synthetic circuits in general, memory can often be set by protein concentrations: because of the relative stability of proteins, the degradation rate is often dominated by the growth rate, and inheritance is a significant factor. Cells can then be primed to respond to events that recur with frequencies faster than the time to eliminate proteins. Protein memory can be extended if cells reach extremely low growth rates or no growth. Here we characterize the necessary time scales for different quantities of protein memory, measured as relative entropy (Kullback-Leibler divergence), for a variety of cellular growth arrest transition dynamics. We identify a critical manifold in relative protein degradation/growth arrest time scales where information is, in principle, preserved indefinitely because proteins are trapped at a concentration determined by the competing time scales as long as nongrowth-mediated protein degradation is negligible. We next asked what characteristics of growth arrest dynamics and initial protein distributions best preserve or eliminate information about previous environments. We identified that sharp growth arrest transitions with skewed initial protein distributions optimize flexibility, with information preservation and minimal cost of residual protein. As a result, a nearly memoryless regime, corresponding to a form of bet-hedging, may be an optimal strategy for storage of information by protein concentrations in growth-arrested cells.

Amino acid positions subject to multiple coevolutionary constraints can be robustly identified by their eigenvector network centrality scores.

As proteins evolve, amino acid positions key to protein structure or function are subject to mutational constraints. These positions can be detected by analyzing sequence families for amino acid conservation or for coevolution between pairs of positions. Coevolutionary scores are usually rank-ordered and thresholded to reveal the top pairwise scores, but they also can be treated as weighted networks. Here, we used network analyses to bypass a major complication of coevolution studies: For a given sequence alignment, alternative algorithms usually identify different, top pairwise scores. We reconciled results from five commonly-used, mathematically divergent algorithms (ELSC, McBASC, OMES, SCA, and ZNMI), using the LacI/GalR and 1,6-bisphosphate aldolase protein families as models. Calculations used unthresholded coevolution scores from which column-specific properties such as sequence entropy and random noise were subtracted; "central" positions were identified by calculating various network centrality scores. When compared among algorithms, network centrality methods, particularly eigenvector centrality, showed markedly better agreement than comparisons of the top pairwise scores. Positions with large centrality scores occurred at key structural locations and/or were functionally sensitive to mutations. Further, the top central positions often differed from those with top pairwise coevolution scores: instead of a few strong scores, central positions often had multiple, moderate scores. We conclude that eigenvector centrality calculations reveal a robust evolutionary pattern of constraints-detectable by divergent algorithms--that occur at key protein locations. Finally, we discuss the fact that multiple patterns coexist in evolutionary data that, together, give rise to emergent protein functions.

Interplay of gene expression noise and ultrasensitive dynamics affects bacterial operon organization.

Bacterial chromosomes are organized into polycistronic cotranscribed operons, but the evolutionary pressures maintaining them are unclear. We hypothesized that operons alter gene expression noise characteristics, resulting in selection for or against maintaining operons depending on network architecture. Mathematical models for 6 functional classes of network modules showed that three classes exhibited decreased noise and 3 exhibited increased noise with same-operon cotranscription of interacting proteins. Noise reduction was often associated with a decreased chance of reaching an ultrasensitive threshold. Stochastic simulations of the lac operon demonstrated that the predicted effects of transcriptional coupling hold for a complex network module. We employed bioinformatic analysis to find overrepresentation of noise-minimizing operon organization compared with randomized controls. Among constitutively expressed physically interacting protein pairs, higher coupling frequencies appeared at lower expression levels, where noise effects are expected to be dominant. Our results thereby suggest an important role for gene expression noise, in many cases interacting with an ultrasensitive switch, in maintaining or selecting for operons in bacterial chromosomes.

Adaptive evolution of the lactose utilization network in experimentally evolved populations of Escherichia coli.

Adaptation to novel environments is often associated with changes in gene regulation. Nevertheless, few studies have been able both to identify the genetic basis of changes in regulation and to demonstrate why these changes are beneficial. To this end, we have focused on understanding both how and why the lactose utilization network has evolved in replicate populations of Escherichia coli. We found that lac operon regulation became strikingly variable, including changes in the mode of environmental response (bimodal, graded, and constitutive), sensitivity to inducer concentration, and maximum expression level. In addition, some classes of regulatory change were enriched in specific selective environments. Sequencing of evolved clones, combined with reconstruction of individual mutations in the ancestral background, identified mutations within the lac operon that recapitulate many of the evolved regulatory changes. These mutations conferred fitness benefits in environments containing lactose, indicating that the regulatory changes are adaptive. The same mutations conferred different fitness effects when present in an evolved clone, indicating that interactions between the lac operon and other evolved mutations also contribute to fitness. Similarly, changes in lac regulation not explained by lac operon mutations also point to important interactions with other evolved mutations. Together these results underline how dynamic regulatory interactions can be, in this case evolving through mutations both within and external to the canonical lactose utilization network.

Non-transcriptional regulatory processes shape transcriptional network dynamics.

Information about the extra- or intracellular environment is often captured as biochemical signals that propagate through regulatory networks. These signals eventually drive phenotypic changes, typically by altering gene expression programmes in the cell. Reconstruction of transcriptional regulatory networks has given a compelling picture of bacterial physiology, but transcriptional network maps alone often fail to describe phenotypes. Cellular response dynamics are ultimately determined by interactions between transcriptional and non-transcriptional networks, with dramatic implications for physiology and evolution. Here, we provide an overview of non-transcriptional interactions that can affect the performance of natural and synthetic bacterial regulatory networks.

Bistable responses in bacterial genetic networks: designs and dynamical consequences.

A key property of living cells is their ability to react to stimuli with specific biochemical responses. These responses can be understood through the dynamics of underlying biochemical and genetic networks. Evolutionary design principles have been well studied in networks that display graded responses, with a continuous relationship between input signal and system output. Alternatively, biochemical networks can exhibit bistable responses so that over a range of signals the network possesses two stable steady states. In this review, we discuss several conceptual examples illustrating network designs that can result in a bistable response of the biochemical network. Next, we examine manifestations of these designs in bacterial master-regulatory genetic circuits. In particular, we discuss mechanisms and dynamic consequences of bistability in three circuits: two-component systems, sigma-factor networks, and a multistep phosphorelay. Analyzing these examples allows us to expand our knowledge of evolutionary design principles networks with bistable responses.

Adaptable functionality of transcriptional feedback in bacterial two-component systems.

A widespread mechanism of bacterial signaling occurs through two-component systems, comprised of a sensor histidine kinase (SHK) and a transcriptional response regulator (RR). The SHK activates RR by phosphorylation. The most common two-component system structure involves expression from a single operon, the transcription of which is activated by its own phosphorylated RR. The role of this feedback is poorly understood, but it has been associated with an overshooting kinetic response and with fast recovery of previous interrupted signaling events in different systems. Mathematical models show that overshoot is only attainable with negative feedback that also improves response time. Our models also predict that fast recovery of previous interrupted signaling depends on high accumulation of SHK and RR, which is more likely in a positive feedback regime. We use Monte Carlo sampling of the parameter space to explore the range of attainable model behaviors. The model predicts that the effective feedback sign can change from negative to positive depending on the signal level. Variations in two-component system architectures and parameters may therefore have evolved to optimize responses in different bacterial lifestyles. We propose a conceptual model where low signal conditions result in a responsive system with effectively negative feedback while high signal conditions with positive feedback favor persistence of system output.

Synergy between individual TNF-dependent functions determines granuloma performance for controlling Mycobacterium tuberculosis infection.

Mycobacterium tuberculosis is one of the world's most deadly human pathogens; an integrated understanding of how it successfully survives in its host is crucial to developing new treatment strategies. One notable characteristic of infection with M. tuberculosis is the formation of granulomas, aggregates of immune cells whose structure and function may reflect success or failure of the host to contain infection. One central regulator of host responses to infection, including granuloma formation, is the pleiotropic cytokine TNF-alpha. Experimental work has characterized roles for TNF in macrophage activation; regulation of apoptosis; chemokine and cytokine production; and regulation of cellular recruitment via transendothelial migration. Separating the effects of these functions is presently difficult or impossible in vivo. To this end, we applied a computational model to understand specific roles of TNF in control of tuberculosis in a single granuloma. In the model, cells are represented as discrete entities on a spatial grid responding to environmental stimuli by following programmed rules determined from published experimental studies. Simulated granulomas emerge as a result of these rules. After confirming the importance of TNF in this model, we assessed the effects of individual TNF functions. The model predicts that multiple TNF activities contribute to control of infection within the granuloma, with macrophage activation as a key effector mechanism for controlling bacterial growth. Results suggest that bacterial numbers are a strong contributing factor to granuloma structure with TNF. Finally, TNF-dependent apoptosis may reduce inflammation at the cost of impairing mycobacterial clearance.

The timing of TNF and IFN-gamma signaling affects macrophage activation strategies during Mycobacterium tuberculosis infection.

During most infections, the population of immune cells known as macrophages are key to taking up and killing bacteria as an integral part of the immune response. However, during infection with Mycobacterium tuberculosis (Mtb), host macrophages serve as the preferred environment for mycobacterial growth. Further, killing of Mtb by macrophages is impaired unless they become activated. Activation is induced by stimulation from bacterial antigens and inflammatory cytokines derived from helper T cells. The key macrophage-activating cytokines in Mtb infection are tumor necrosis factor-alpha (TNF) and interferon (IFN)-gamma. Due to differences in cellular sources and secretion pathways for TNF and IFN-gamma, the possibility of heterogeneous cytokine distributions exists, suggesting that the timing of macrophage activation from these signals may affect activation kinetics and thus impact the outcome of Mtb infection. Here we use a mathematical model to show that negative feedback from production of nitric oxide (the key mediator of mycobacterial killing) that typically optimizes macrophage responses to activating stimuli may reduce effective killing of Mtb. Statistical sensitivity analysis predicts that if TNF and IFN-gamma signals precede infection, the level of negative feedback may have a strong effect on how effectively macrophages kill Mtb. However, this effect is relaxed when IFN-gamma or TNF+IFN-gamma signals are received coincident with infection. Under these conditions, the model suggests that negative feedback induces fast responses and an initial overshoot of nitric oxide production for given doses of TNF and IFN-gamma, favoring killing of Mtb. Together, our results suggest that direct entry of macrophages into a granuloma site (and not distal to it) from lung vascular sources represents a preferred host strategy for mycobacterial control. We examine implications of these results in establishment of latent Mtb infection.

Requirement for multiple activation signals by anti-inflammatory feedback in macrophages.

Pathogen killing is one of the primary roles of macrophages, utilizing potent effectors such as nitric oxide (NO) and involving other cellular machinery including iron regulatory apparatus. Macrophages become strongly activated upon receipt of appropriate signaling with cytokines and pathogen-derived endotoxins. However, they must resist activation in the absence of decisive signaling due to the energetic demands of activation coupled with the toxic nature of effector molecules to surrounding tissues. We have developed a mathematical model of the modular biochemical network of macrophages involved with activation, pathogen killing and iron regulation. This model requires synergistic interaction of multiple activation signals to overcome the quiescent state. To achieve a trade-off between macrophage quiescence and activation, strong activation signals are modulated via negative regulation by NO. In this way a single activation signal is insufficient for complete activation. In addition, our results suggest that iron regulation is usually controlled by activation signals. However, under conditions of partial macrophage activation, exogenous iron levels play a key role in regulating NO production. This model will be useful for evaluating macrophage control of intracellular pathogens in addition to the biochemical mechanisms examined here.