@misc{8951,
  abstract     = {Gene expression levels are influenced by multiple coexisting molecular mechanisms. Some of these interactions, such as those of transcription factors and promoters have been studied extensively. However, predicting phenotypes of gene regulatory networks remains a major challenge. Here, we use a well-defined synthetic gene regulatory network to study how network phenotypes depend on local genetic context, i.e. the genetic neighborhood of a transcription factor and its relative position. We show that one gene regulatory network with fixed topology can display not only quantitatively but also qualitatively different phenotypes, depending solely on the local genetic context of its components. Our results demonstrate that changes in local genetic context can place a single transcriptional unit within two separate regulons without the need for complex regulatory sequences. We propose that relative order of individual transcriptional units, with its potential for combinatorial complexity, plays an important role in shaping phenotypes of gene regulatory networks.},
  author       = {Nagy-Staron, Anna A},
  keywords     = {Gene regulatory networks, Gene expression, Escherichia coli, Synthetic Biology},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{Sequences of gene regulatory network permutations for the article "Local genetic context shapes the function of a gene regulatory network"}},
  doi          = {10.15479/AT:ISTA:8951},
  year         = {2020},
}

@misc{7383,
  abstract     = {Organisms cope with change by employing transcriptional regulators. However, when faced with rare environments, the evolution of transcriptional regulators and their promoters may be too slow. We ask whether the intrinsic instability of gene duplication and amplification provides a generic alternative to canonical gene regulation. By real-time monitoring of gene copy number mutations in E. coli, we show that gene duplications and amplifications enable adaptation to fluctuating environments by rapidly generating copy number, and hence expression level, polymorphism. This ‘amplification-mediated gene expression tuning’ occurs on timescales similar to canonical gene regulation and can deal with rapid environmental changes. Mathematical modeling shows that amplifications also tune gene expression in stochastic environments where transcription factor-based schemes are hard to evolve or maintain. The fleeting nature of gene amplifications gives rise to a generic population-level mechanism that relies on genetic heterogeneity to rapidly tune expression of any gene, without leaving any genomic signature.},
  author       = {Grah, Rok},
  keywords     = {Matlab scripts, analysis of microfluidics, mathematical model},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{Matlab scripts for the Paper: Gene Amplification as a Form of Population-Level Gene Expression regulation}},
  doi          = {10.15479/AT:ISTA:7383},
  year         = {2020},
}

@article{7569,
  abstract     = {Genes differ in the frequency at which they are expressed and in the form of regulation used to control their activity. In particular, positive or negative regulation can lead to activation of a gene in response to an external signal. Previous works proposed that the form of regulation of a gene correlates with its frequency of usage: positive regulation when the gene is frequently expressed and negative regulation when infrequently expressed. Such network design means that, in the absence of their regulators, the genes are found in their least required activity state, hence regulatory intervention is often necessary. Due to the multitude of genes and regulators, spurious binding and unbinding events, called “crosstalk”, could occur. To determine how the form of regulation affects the global crosstalk in the network, we used a mathematical model that includes multiple regulators and multiple target genes. We found that crosstalk depends non-monotonically on the availability of regulators. Our analysis showed that excess use of regulation entailed by the formerly suggested network design caused high crosstalk levels in a large part of the parameter space. We therefore considered the opposite ‘idle’ design, where the default unregulated state of genes is their frequently required activity state. We found, that ‘idle’ design minimized the use of regulation and thus minimized crosstalk. In addition, we estimated global crosstalk of S. cerevisiae using transcription factors binding data. We demonstrated that even partial network data could suffice to estimate its global crosstalk, suggesting its applicability to additional organisms. We found that S. cerevisiae estimated crosstalk is lower than that of a random network, suggesting that natural selection reduces crosstalk. In summary, our study highlights a new type of protein production cost which is typically overlooked: that of regulatory interference caused by the presence of excess regulators in the cell. It demonstrates the importance of whole-network descriptions, which could show effects missed by single-gene models.},
  author       = {Grah, Rok and Friedlander, Tamar},
  issn         = {1553-7358},
  journal      = {PLOS Computational Biology},
  number       = {2},
  publisher    = {Public Library of Science},
  title        = {{The relation between crosstalk and gene regulation form revisited}},
  doi          = {10.1371/journal.pcbi.1007642},
  volume       = {16},
  year         = {2020},
}

@article{7652,
  abstract     = {Organisms cope with change by taking advantage of transcriptional regulators. However, when faced with rare environments, the evolution of transcriptional regulators and their promoters may be too slow. Here, we investigate whether the intrinsic instability of gene duplication and amplification provides a generic alternative to canonical gene regulation. Using real-time monitoring of gene-copy-number mutations in Escherichia coli, we show that gene duplications and amplifications enable adaptation to fluctuating environments by rapidly generating copy-number and, therefore, expression-level polymorphisms. This amplification-mediated gene expression tuning (AMGET) occurs on timescales that are similar to canonical gene regulation and can respond to rapid environmental changes. Mathematical modelling shows that amplifications also tune gene expression in stochastic environments in which transcription-factor-based schemes are hard to evolve or maintain. The fleeting nature of gene amplifications gives rise to a generic population-level mechanism that relies on genetic heterogeneity to rapidly tune the expression of any gene, without leaving any genomic signature.},
  author       = {Tomanek, Isabella and Grah, Rok and Lagator, M. and Andersson, A. M. C. and Bollback, Jonathan P and Tkačik, Gašper and Guet, Calin C},
  issn         = {2397-334X},
  journal      = {Nature Ecology & Evolution},
  number       = {4},
  pages        = {612--625},
  publisher    = {Springer Nature},
  title        = {{Gene amplification as a form of population-level gene expression regulation}},
  doi          = {10.1038/s41559-020-1132-7},
  volume       = {4},
  year         = {2020},
}

@unpublished{7675,
  abstract     = {In prokaryotes, thermodynamic models of gene regulation provide a highly quantitative mapping from promoter sequences to gene expression levels that is compatible with in vivo and in vitro bio-physical measurements. Such concordance has not been achieved for models of enhancer function in eukaryotes. In equilibrium models, it is difficult to reconcile the reported short transcription factor (TF) residence times on the DNA with the high specificity of regulation. In non-equilibrium models, progress is difficult due to an explosion in the number of parameters. Here, we navigate this complexity by looking for minimal non-equilibrium enhancer models that yield desired regulatory phenotypes: low TF residence time, high specificity and tunable cooperativity. We find that a single extra parameter, interpretable as the “linking rate” by which bound TFs interact with Mediator components, enables our models to escape equilibrium bounds and access optimal regulatory phenotypes, while remaining consistent with the reported phenomenology and simple enough to be inferred from upcoming experiments. We further find that high specificity in non-equilibrium models is in a tradeoff with gene expression noise, predicting bursty dynamics — an experimentally-observed hallmark of eukaryotic transcription. By drastically reducing the vast parameter space to a much smaller subspace that optimally realizes biological function prior to inference from data, our normative approach holds promise for mathematical models in systems biology.},
  author       = {Grah, Rok and Zoller, Benjamin and Tkačik, Gašper},
  booktitle    = {bioRxiv},
  publisher    = {Cold Spring Harbor Laboratory},
  title        = {{Normative models of enhancer function}},
  doi          = {10.1101/2020.04.08.029405},
  year         = {2020},
}

@phdthesis{7680,
  abstract     = {Proteins and their complex dynamic interactions regulate cellular mechanisms from sensing and transducing extracellular signals, to mediating genetic responses, and sustaining or changing cell morphology. To manipulate these protein-protein interactions (PPIs) that govern the behavior and fate of cells, synthetically constructed, genetically encoded tools provide the means to precisely target proteins of interest (POIs), and control their subcellular localization and activity in vitro and in vivo. Ideal synthetic tools react to an orthogonal cue, i.e. a trigger that does not activate any other endogenous process, thereby allowing manipulation of the POI alone.
In optogenetics, naturally occurring photosensory domain from plants, algae and bacteria are re-purposed and genetically fused to POIs. Illumination with light of a specific wavelength triggers a conformational change that can mediate PPIs, such as dimerization or oligomerization. By using light as a trigger, these tools can be activated with high spatial and temporal precision, on subcellular and millisecond scales. Chemogenetic tools consist of protein domains that recognize and bind small molecules. By genetic fusion to POIs, these domains can mediate PPIs upon addition of their specific ligands, which are often synthetically designed to provide highly specific interactions and exhibit good bioavailability.
Most optogenetic tools to mediate PPIs are based on well-studied photoreceptors responding to red, blue or near-UV light, leaving a striking gap in the green band of the visible light spectrum. Among both optogenetic and chemogenetic tools, there is an abundance of methods to induce PPIs, but tools to disrupt them require UV illumination, rely on covalent linkage and subsequent enzymatic cleavage or initially result in protein clustering of unknown stoichiometry.
This work describes how the recently structurally and photochemically characterized green-light responsive cobalamin-binding domains (CBDs) from bacterial transcription factors were re-purposed to function as a green-light responsive optogenetic tool. In contrast to previously engineered optogenetic tools, CBDs do not induce PPI, but rather confer a PPI already upon expression, which can be rapidly disrupted by illumination. This was employed to mimic inhibition of constitutive activity of a growth factor receptor, and successfully implement for cell signalling in mammalian cells and in vivo to rescue development in zebrafish. This work further describes the development and application of a chemically induced de-dimerizer (CDD) based on a recently identified and structurally described bacterial oxyreductase. CDD forms a dimer upon expression in absence of its cofactor, the flavin derivative F420. Safety and of domain expression and ligand exposure are demonstrated in vitro and in vivo in zebrafish. The system is further applied to inhibit cell signalling output from a chimeric receptor upon F420 treatment.
CBDs and CDD expand the repertoire of synthetic tools by providing novel mechanisms of mediating PPIs, and by recognizing previously not utilized cues. In the future, they can readily be combined with existing synthetic tools to functionally manipulate PPIs in vitro and in vivo.},
  author       = {Kainrath, Stephanie},
  issn         = {2663-337X},
  pages        = {98},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{Synthetic tools for optogenetic and chemogenetic inhibition of cellular signals}},
  doi          = {10.15479/AT:ISTA:7680},
  year         = {2020},
}

@article{196,
  abstract     = {The abelian sandpile serves as a model to study self-organized criticality, a phenomenon occurring in biological, physical and social processes. The identity of the abelian group is a fractal composed of self-similar patches, and its limit is subject of extensive collaborative research. Here, we analyze the evolution of the sandpile identity under harmonic fields of different orders. We show that this evolution corresponds to periodic cycles through the abelian group characterized by the smooth transformation and apparent conservation of the patches constituting the identity. The dynamics induced by second and third order harmonics resemble smooth stretchings, respectively translations, of the identity, while the ones induced by fourth order harmonics resemble magnifications and rotations. Starting with order three, the dynamics pass through extended regions of seemingly random configurations which spontaneously reassemble into accentuated patterns. We show that the space of harmonic functions projects to the extended analogue of the sandpile group, thus providing a set of universal coordinates identifying configurations between different domains. Since the original sandpile group is a subgroup of the extended one, this directly implies that it admits a natural renormalization. Furthermore, we show that the harmonic fields can be induced by simple Markov processes, and that the corresponding stochastic dynamics show remarkable robustness over hundreds of periods. Finally, we encode information into seemingly random configurations, and decode this information with an algorithm requiring minimal prior knowledge. Our results suggest that harmonic fields might split the sandpile group into sub-sets showing different critical coefficients, and that it might be possible to extend the fractal structure of the identity beyond the boundaries of its domain. },
  author       = {Lang, Moritz and Shkolnikov, Mikhail},
  issn         = {1091-6490},
  journal      = {Proceedings of the National Academy of Sciences},
  number       = {8},
  pages        = {2821--2830},
  publisher    = {National Academy of Sciences},
  title        = {{Harmonic dynamics of the Abelian sandpile}},
  doi          = {10.1073/pnas.1812015116},
  volume       = {116},
  year         = {2019},
}

@article{6717,
  abstract     = {With the recent publication by Silpe and Bassler (2019), considering phage detection of a bacterial quorum-sensing (QS) autoinducer, we now have as many as five examples of phage-associated intercellular communication (Table 1). Each potentially involves ecological inferences by phages as to concentrations of surrounding phage-infected or uninfected bacteria. While the utility of phage detection of bacterial QS molecules may at first glance appear to be straightforward, we suggest in this commentary that the underlying ecological explanation is unlikely to be simple.},
  author       = {Igler, Claudia and Abedon, Stephen T.},
  journal      = {Frontiers in Microbiology},
  publisher    = {Frontiers},
  title        = {{Commentary: A host-produced quorum-sensing autoinducer controls a phage lysis-lysogeny decision}},
  doi          = {10.3389/fmicb.2019.01171},
  volume       = {10},
  year         = {2019},
}

@article{6784,
  abstract     = {Mathematical models have been used successfully at diverse scales of biological organization, ranging from ecology and population dynamics to stochastic reaction events occurring between individual molecules in single cells. Generally, many biological processes unfold across multiple scales, with mutations being the best studied example of how stochasticity at the molecular scale can influence outcomes at the population scale. In many other contexts, however, an analogous link between micro- and macro-scale remains elusive, primarily due to the challenges involved in setting up and analyzing multi-scale models. Here, we employ such a model to investigate how stochasticity propagates from individual biochemical reaction events in the bacterial innate immune system to the ecology of bacteria and bacterial viruses. We show analytically how the dynamics of bacterial populations are shaped by the activities of immunity-conferring enzymes in single cells and how the ecological consequences imply optimal bacterial defense strategies against viruses. Our results suggest that bacterial populations in the presence of viruses can either optimize their initial growth rate or their population size, with the first strategy favoring simple immunity featuring a single restriction modification system and the second strategy favoring complex bacterial innate immunity featuring several simultaneously active restriction modification systems.},
  author       = {Ruess, Jakob and Pleska, Maros and Guet, Calin C and Tkačik, Gašper},
  issn         = {1553-7358},
  journal      = {PLoS Computational Biology},
  number       = {7},
  publisher    = {Public Library of Science},
  title        = {{Molecular noise of innate immunity shapes bacteria-phage ecologies}},
  doi          = {10.1371/journal.pcbi.1007168},
  volume       = {15},
  year         = {2019},
}

@misc{7016,
  abstract     = {Organisms cope with change by employing transcriptional regulators. However, when faced with rare environments, the evolution of transcriptional regulators and their promoters may be too slow. We ask whether the intrinsic instability of gene duplication and amplification provides a generic alternative to canonical gene regulation. By real-time monitoring of gene copy number mutations in E. coli, we show that gene duplications and amplifications enable adaptation to fluctuating environments by rapidly generating copy number, and hence expression level, polymorphism. This ‘amplification-mediated gene expression tuning’ occurs on timescales similar to canonical gene regulation and can deal with rapid environmental changes. Mathematical modeling shows that amplifications also tune gene expression in stochastic environments where transcription factor-based schemes are hard to evolve or maintain. The fleeting nature of gene amplifications gives rise to a generic population-level mechanism that relies on genetic heterogeneity to rapidly tune expression of any gene, without leaving any genomic signature.},
  author       = {Tomanek, Isabella},
  keywords     = {Escherichia coli, gene amplification, galactose, DOG, experimental evolution, Illumina sequence data, FACS data, microfluidics data},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{Data for the paper "Gene amplification as a form of population-level gene expression regulation"}},
  doi          = {10.15479/AT:ISTA:7016},
  year         = {2019},
}

@inproceedings{7147,
  abstract     = {The expression of a gene is characterised by its transcription factors and the function processing them. If the transcription factors are not affected by gene products, the regulating function is often represented as a combinational logic circuit, where the outputs (product) are determined by current input values (transcription factors) only, and are hence independent on their relative arrival times. However, the simultaneous arrival of transcription factors (TFs) in genetic circuits is a strong assumption, given that the processes of transcription and translation of a gene into a protein introduce intrinsic time delays and that there is no global synchronisation among the arrival times of different molecular species at molecular targets.

In this paper, we construct an experimentally implementable genetic circuit with two inputs and a single output, such that, in presence of small delays in input arrival, the circuit exhibits qualitatively distinct observable phenotypes. In particular, these phenotypes are long lived transients: they all converge to a single value, but so slowly, that they seem stable for an extended time period, longer than typical experiment duration. We used rule-based language to prototype our circuit, and we implemented a search for finding the parameter combinations raising the phenotypes of interest.

The behaviour of our prototype circuit has wide implications. First, it suggests that GRNs can exploit event timing to create phenotypes. Second, it opens the possibility that GRNs are using event timing to react to stimuli and memorise events, without explicit feedback in regulation. From the modelling perspective, our prototype circuit demonstrates the critical importance of analysing the transient dynamics at the promoter binding sites of the DNA, before applying rapid equilibrium assumptions.},
  author       = {Guet, Calin C and Henzinger, Thomas A and Igler, Claudia and Petrov, Tatjana and Sezgin, Ali},
  booktitle    = {17th International Conference on Computational Methods in Systems Biology},
  isbn         = {9783030313036},
  issn         = {1611-3349},
  location     = {Trieste, Italy},
  pages        = {155--187},
  publisher    = {Springer Nature},
  title        = {{Transient memory in gene regulation}},
  doi          = {10.1007/978-3-030-31304-3_9},
  volume       = {11773},
  year         = {2019},
}

@phdthesis{6371,
  abstract     = {Decades of studies have revealed the mechanisms of gene regulation in molecular detail. We make use of such well-described regulatory systems to explore how the molecular mechanisms of protein-protein and protein-DNA interactions shape the dynamics and evolution of gene regulation. 

i) We uncover how the biophysics of protein-DNA binding determines the potential of regulatory networks to evolve and adapt, which can be captured using a simple mathematical model. 
ii) The evolution of regulatory connections can lead to a significant amount of crosstalk between binding proteins. We explore the effect of crosstalk on gene expression from a target promoter, which seems to be modulated through binding competition at non-specific DNA sites. 
iii) We investigate how the very same biophysical characteristics as in i) can generate significant fitness costs for cells through global crosstalk, meaning non-specific DNA binding across the genomic background. 
iv) Binding competition between proteins at a target promoter is a prevailing regulatory feature due to the prevalence of co-regulation at bacterial promoters. However, the dynamics of these systems are not always straightforward to determine even if the molecular mechanisms of regulation are known. A detailed model of the biophysical interactions reveals that interference between the regulatory proteins can constitute a new, generic form of system memory that records the history of the input signals at the promoter. 

We demonstrate how the biophysics of protein-DNA binding can be harnessed to investigate the principles that shape and ultimately limit cellular gene regulation. These results provide a basis for studies of higher-level functionality, which arises from the underlying regulation.   
},
  author       = {Igler, Claudia},
  issn         = {2663-337X},
  keywords     = {gene regulation, biophysics, transcription factor binding, bacteria},
  pages        = {152},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{On the nature of gene regulatory design - The biophysics of transcription factor binding shapes gene regulation}},
  doi          = {10.15479/AT:ISTA:6371},
  year         = {2019},
}

@article{6465,
  abstract     = {Tight control over protein degradation is a fundamental requirement for cells to respond rapidly to various stimuli and adapt to a fluctuating environment. Here we develop a versatile, easy-to-handle library of destabilizing tags (degrons) for the precise regulation of protein expression profiles in mammalian cells by modulating target protein half-lives in a predictable manner. Using the well-established tetracycline gene-regulation system as a model, we show that the dynamics of protein expression can be tuned by fusing appropriate degron tags to gene regulators. Next, we apply this degron library to tune a synthetic pulse-generating circuit in mammalian cells. With this toolbox we establish a set of pulse generators with tailored pulse lengths and magnitudes of protein expression. This methodology will prove useful in the functional roles of essential proteins, fine-tuning of gene-expression systems, and enabling a higher complexity in the design of synthetic biological systems in mammalian cells.},
  author       = {Chassin, Hélène and Müller, Marius and Tigges, Marcel and Scheller, Leo and Lang, Moritz and Fussenegger, Martin},
  issn         = {20411723},
  journal      = {Nature Communications},
  number       = {1},
  publisher    = {Springer Nature},
  title        = {{A modular degron library for synthetic circuits in mammalian cells}},
  doi          = {10.1038/s41467-019-09974-5},
  volume       = {10},
  year         = {2019},
}

@article{138,
  abstract     = {Autoregulation is the direct modulation of gene expression by the product of the corresponding gene. Autoregulation of bacterial gene expression has been mostly studied at the transcriptional level, when a protein acts as the cognate transcriptional repressor. A recent study investigating dynamics of the bacterial toxin–antitoxin MazEF system has shown how autoregulation at both the transcriptional and post-transcriptional levels affects the heterogeneity of Escherichia coli populations. Toxin–antitoxin systems hold a crucial but still elusive part in bacterial response to stress. This perspective highlights how these modules can also serve as a great model system for investigating basic concepts in gene regulation. However, as the genomic background and environmental conditions substantially influence toxin activation, it is important to study (auto)regulation of toxin–antitoxin systems in well-defined setups as well as in conditions that resemble the environmental niche.},
  author       = {Nikolic, Nela},
  journal      = {Current Genetics},
  number       = {1},
  pages        = {133--138},
  publisher    = {Springer},
  title        = {{Autoregulation of bacterial gene expression: lessons from the MazEF toxin–antitoxin system}},
  doi          = {10.1007/s00294-018-0879-8},
  volume       = {65},
  year         = {2019},
}

@misc{9786,
  author       = {Ruess, Jakob and Pleska, Maros and Guet, Calin C and Tkačik, Gašper},
  publisher    = {Public Library of Science},
  title        = {{Supporting text and results}},
  doi          = {10.1371/journal.pcbi.1007168.s001},
  year         = {2019},
}

@phdthesis{26,
  abstract     = {Expression of genes is a fundamental molecular phenotype that is subject to evolution by different types of mutations. Both the rate and the effect of mutations may depend on the DNA sequence context of a particular gene or a particular promoter sequence. In this thesis I investigate the nature of this dependence using simple genetic systems in Escherichia coli. With these systems I explore the evolution of constitutive gene expression from random starting sequences at different loci on the chromosome and at different locations in sequence space. First, I dissect chromosomal neighborhood effects that underlie locus-dependent differences in the potential of a gene under selection to become more highly expressed. Next, I find that the effects of point mutations in promoter sequences are dependent on sequence context, and that an existing energy matrix model performs poorly in predicting relative expression of unrelated sequences. Finally, I show that a substantial fraction of random sequences contain functional promoters and I present an extended thermodynamic model that predicts promoter strength in full sequence space. Taken together, these results provide new insights and guides on how to integrate information on sequence context to improve our qualitative and quantitative understanding of bacterial gene expression, with implications for rapid evolution of drug resistance, de novo evolution of genes, and horizontal gene transfer.},
  author       = {Steinrück, Magdalena},
  issn         = {2663-337X},
  pages        = {109},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{The influence of sequence context on the evolution of bacterial gene expression}},
  doi          = {10.15479/AT:ISTA:th1059},
  year         = {2018},
}

@article{305,
  abstract     = {The hanging-drop network (HDN) is a technology platform based on a completely open microfluidic network at the bottom of an inverted, surface-patterned substrate. The platform is predominantly used for the formation, culturing, and interaction of self-assembled spherical microtissues (spheroids) under precisely controlled flow conditions. Here, we describe design, fabrication, and operation of microfluidic hanging-drop networks.},
  author       = {Misun, Patrick and Birchler, Axel and Lang, Moritz and Hierlemann, Andreas and Frey, Olivier},
  journal      = {Methods in Molecular Biology},
  pages        = {183 -- 202},
  publisher    = {Springer},
  title        = {{Fabrication and operation of microfluidic hanging drop networks}},
  doi          = {10.1007/978-1-4939-7792-5_15},
  volume       = {1771},
  year         = {2018},
}

@article{82,
  abstract     = {In experimental cultures, when bacteria are mixed with lytic (virulent) bacteriophage, bacterial cells resistant to the phage commonly emerge and become the dominant population of bacteria. Following the ascent of resistant mutants, the densities of bacteria in these simple communities become limited by resources rather than the phage. Despite the evolution of resistant hosts, upon which the phage cannot replicate, the lytic phage population is most commonly maintained in an apparently stable state with the resistant bacteria. Several mechanisms have been put forward to account for this result. Here we report the results of population dynamic/evolution experiments with a virulent mutant of phage Lambda, λVIR, and Escherichia coli in serial transfer cultures. We show that, following the ascent of λVIR-resistant bacteria, λVIRis maintained in the majority of cases in maltose-limited minimal media and in all cases in nutrient-rich broth. Using mathematical models and experiments, we show that the dominant mechanism responsible for maintenance of λVIRin these resource-limited populations dominated by resistant E. coli is a high rate of either phenotypic or genetic transition from resistance to susceptibility—a hitherto undemonstrated mechanism we term &quot;leaky resistance.&quot; We discuss the implications of leaky resistance to our understanding of the conditions for the maintenance of phage in populations of bacteria—their “existence conditions.”.},
  author       = {Chaudhry, Waqas and Pleska, Maros and Shah, Nilang and Weiss, Howard and Mccall, Ingrid and Meyer, Justin and Gupta, Animesh and Guet, Calin C and Levin, Bruce},
  journal      = {PLoS Biology},
  number       = {8},
  publisher    = {Public Library of Science},
  title        = {{Leaky resistance and the conditions for the existence of lytic bacteriophage}},
  doi          = {10.1371/journal.pbio.2005971},
  volume       = {16},
  year         = {2018},
}

@article{19,
  abstract     = {Bacteria regulate genes to survive antibiotic stress, but regulation can be far from perfect. When regulation is not optimal, mutations that change gene expression can contribute to antibiotic resistance. It is not systematically understood to what extent natural gene regulation is or is not optimal for distinct antibiotics, and how changes in expression of specific genes quantitatively affect antibiotic resistance. Here we discover a simple quantitative relation between fitness, gene expression, and antibiotic potency, which rationalizes our observation that a multitude of genes and even innate antibiotic defense mechanisms have expression that is critically nonoptimal under antibiotic treatment. First, we developed a pooled-strain drug-diffusion assay and screened Escherichia coli overexpression and knockout libraries, finding that resistance to a range of 31 antibiotics could result from changing expression of a large and functionally diverse set of genes, in a primarily but not exclusively drug-specific manner. Second, by synthetically controlling the expression of single-drug and multidrug resistance genes, we observed that their fitness-expression functions changed dramatically under antibiotic treatment in accordance with a log-sensitivity relation. Thus, because many genes are nonoptimally expressed under antibiotic treatment, many regulatory mutations can contribute to resistance by altering expression and by activating latent defenses.},
  author       = {Palmer, Adam and Chait, Remy P and Kishony, Roy},
  issn         = {0737-4038},
  journal      = {Molecular Biology and Evolution},
  number       = {11},
  pages        = {2669 -- 2684},
  publisher    = {Oxford University Press},
  title        = {{Nonoptimal gene expression creates latent potential for antibiotic resistance}},
  doi          = {10.1093/molbev/msy163},
  volume       = {35},
  year         = {2018},
}

@article{67,
  abstract     = {Gene regulatory networks evolve through rewiring of individual components—that is, through changes in regulatory connections. However, the mechanistic basis of regulatory rewiring is poorly understood. Using a canonical gene regulatory system, we quantify the properties of transcription factors that determine the evolutionary potential for rewiring of regulatory connections: robustness, tunability and evolvability. In vivo repression measurements of two repressors at mutated operator sites reveal their contrasting evolutionary potential: while robustness and evolvability were positively correlated, both were in trade-off with tunability. Epistatic interactions between adjacent operators alleviated this trade-off. A thermodynamic model explains how the differences in robustness, tunability and evolvability arise from biophysical characteristics of repressor–DNA binding. The model also uncovers that the energy matrix, which describes how mutations affect repressor–DNA binding, encodes crucial information about the evolutionary potential of a repressor. The biophysical determinants of evolutionary potential for regulatory rewiring constitute a mechanistic framework for understanding network evolution.},
  author       = {Igler, Claudia and Lagator, Mato and Tkacik, Gasper and Bollback, Jonathan P and Guet, Calin C},
  journal      = {Nature Ecology and Evolution},
  number       = {10},
  pages        = {1633 -- 1643},
  publisher    = {Nature Publishing Group},
  title        = {{Evolutionary potential of transcription factors for gene regulatory rewiring}},
  doi          = {10.1038/s41559-018-0651-y},
  volume       = {2},
  year         = {2018},
}

