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-α. 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.

Stress response genes and their regulators form networks that underlie drug resistance. These networks often have an inherent tradeoff: their expression is costly in the absence of stress, but beneficial in stress. They can quickly emerge in the genomes of infectious microbes and cancer cells, protecting them from treatment. Yet, the evolution of stress resistance networks is not well understood. Here, we use a two‐component synthetic gene circuit integrated into the budding yeast genome to model experimentally the adaptation of a stress response module and its host genome in three different scenarios. In agreement with computational predictions, we find that: (i) intra‐module mutations target and eliminate the module if it confers only cost without any benefit to the cell; (ii) intra‐ and extra‐module mutations jointly activate the module if it is potentially beneficial and confers no cost; and (iii) a few specific mutations repeatedly fine‐tune the module's noisy response if it has excessive costs and/or insufficient benefits. Overall, these findings reveal how the timing and mechanisms of stress response network evolution depend on the environment.

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. Differences in gene regulation underlie many important biological processes and are thought to be important for the adaption of organisms to novel environments. Here we focus on the regulation of a group of well-studied genes, the lac operon, that control the utilization of lactose sugar, and we examine how their regulation changes during the adaptation of populations of Escherichia coli bacteria to environments that differ only in the presence of lactose. We find that lac operon regulation is altered in almost all populations that evolve in the presence of lactose and identify two classes of mutations that explain a large part of this change and that confer significant fitness benefits. Interestingly, our study indicates that other mutations, lying outside of the commonly recognized control region, cause new regulation of the lac operon. Together these findings reinforce the importance of changes in gene regulation during evolution and suggest that the biological basis of these changes can be complex and involve novel interactions between genes.

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