@misc{9847, abstract = {information on culture conditions, phage mutagenesis, verification and lysate preparation; Raw data}, author = {Pleska, Maros and Guet, Calin C}, publisher = {The Royal Society}, title = {{Supplementary materials and methods; Full data set from effects of mutations in phage restriction sites during escape from restriction–modification}}, doi = {10.6084/m9.figshare.5633917.v1}, year = {2017}, } @misc{9845, abstract = {Estimates of 13 C-arabinose and 2 H-glucose uptake from the fractions of heavy isotopes measured in single cells}, author = {Nikolic, Nela and Schreiber, Frank and Dal Co, Alma and Kiviet, Daniel and Bergmiller, Tobias and Littmann, Sten and Kuypers, Marcel and Ackermann, Martin}, publisher = {Public Library of Science}, title = {{Mathematical model}}, doi = {10.1371/journal.pgen.1007122.s017}, year = {2017}, } @misc{9849, abstract = {This text provides additional information about the model, a derivation of the analytic results in Eq (4), and details about simulations of an additional parameter set.}, author = {Lukacisinova, Marta and Novak, Sebastian and Paixao, Tiago}, publisher = {Public Library of Science}, title = {{Modelling and simulation details}}, doi = {10.1371/journal.pcbi.1005609.s001}, year = {2017}, } @misc{9850, abstract = {In this text, we discuss how a cost of resistance and the possibility of lethal mutations impact our model.}, author = {Lukacisinova, Marta and Novak, Sebastian and Paixao, Tiago}, publisher = {Public Library of Science}, title = {{Extensions of the model}}, doi = {10.1371/journal.pcbi.1005609.s002}, year = {2017}, } @misc{9846, author = {Nikolic, Nela and Schreiber, Frank and Dal Co, Alma and Kiviet, Daniel and Bergmiller, Tobias and Littmann, Sten and Kuypers, Marcel and Ackermann, Martin}, publisher = {Public Library of Science}, title = {{Supplementary methods}}, doi = {10.1371/journal.pgen.1007122.s016}, year = {2017}, } @misc{9851, abstract = {Based on the intuitive derivation of the dynamics of SIM allele frequency pM in the main text, we present a heuristic prediction for the long-term SIM allele frequencies with χ > 1 stresses and compare it to numerical simulations.}, author = {Lukacisinova, Marta and Novak, Sebastian and Paixao, Tiago}, publisher = {Public Library of Science}, title = {{Heuristic prediction for multiple stresses}}, doi = {10.1371/journal.pcbi.1005609.s003}, year = {2017}, } @misc{9852, abstract = {We show how different combination strategies affect the fraction of individuals that are multi-resistant.}, author = {Lukacisinova, Marta and Novak, Sebastian and Paixao, Tiago}, publisher = {Public Library of Science}, title = {{Resistance frequencies for different combination strategies}}, doi = {10.1371/journal.pcbi.1005609.s004}, year = {2017}, } @misc{9844, author = {Nikolic, Nela and Schreiber, Frank and Dal Co, Alma and Kiviet, Daniel and Bergmiller, Tobias and Littmann, Sten and Kuypers, Marcel and Ackermann, Martin}, publisher = {Public Library of Science}, title = {{Source data for figures and tables}}, doi = {10.1371/journal.pgen.1007122.s018}, year = {2017}, } @article{561, abstract = {Restriction–modification systems are widespread genetic elements that protect bacteria from bacteriophage infections by recognizing and cleaving heterologous DNA at short, well-defined sequences called restriction sites. Bioinformatic evidence shows that restriction sites are significantly underrepresented in bacteriophage genomes, presumably because bacteriophages with fewer restriction sites are more likely to escape cleavage by restriction–modification systems. However, how mutations in restriction sites affect the likelihood of bacteriophage escape is unknown. Using the bacteriophage l and the restriction–modification system EcoRI, we show that while mutation effects at different restriction sites are unequal, they are independent. As a result, the probability of bacteriophage escape increases with each mutated restriction site. Our results experimentally support the role of restriction site avoidance as a response to selection imposed by restriction–modification systems and offer an insight into the events underlying the process of bacteriophage escape.}, author = {Pleska, Maros and Guet, Calin C}, issn = {1744-9561}, journal = {Biology Letters}, number = {12}, publisher = {The Royal Society}, title = {{Effects of mutations in phage restriction sites during escape from restriction–modification}}, doi = {10.1098/rsbl.2017.0646}, volume = {13}, year = {2017}, } @phdthesis{202, abstract = {Restriction-modification (RM) represents the simplest and possibly the most widespread mechanism of self/non-self discrimination in nature. In order to provide bacteria with immunity against bacteriophages and other parasitic genetic elements, RM systems rely on a balance between two enzymes: the restriction enzyme, which cleaves non-self DNA at specific restriction sites, and the modification enzyme, which tags the host’s DNA as self and thus protects it from cleavage. In this thesis, I use population and single-cell level experiments in combination with mathematical modeling to study different aspects of the interplay between RM systems, bacteria and bacteriophages. First, I analyze how mutations in phage restriction sites affect the probability of phage escape – an inherently stochastic process, during which phages accidently get modified instead of restricted. Next, I use single-cell experiments to show that RM systems can, with a low probability, attack the genome of their bacterial host and that this primitive form of autoimmunity leads to a tradeoff between the evolutionary cost and benefit of RM systems. Finally, I investigate the nature of interactions between bacteria, RM systems and temperate bacteriophages to find that, as a consequence of phage escape and its impact on population dynamics, RM systems can promote acquisition of symbiotic bacteriophages, rather than limit it. The results presented here uncover new fundamental biological properties of RM systems and highlight their importance in the ecology and evolution of bacteria, bacteriophages and their interactions.}, author = {Pleska, Maros}, issn = {2663-337X}, pages = {126}, publisher = {Institute of Science and Technology Austria}, title = {{Biology of restriction-modification systems at the single-cell and population level}}, doi = {10.15479/AT:ISTA:th_916}, year = {2017}, } @article{1351, abstract = {The behaviour of gene regulatory networks (GRNs) is typically analysed using simulation-based statistical testing-like methods. In this paper, we demonstrate that we can replace this approach by a formal verification-like method that gives higher assurance and scalability. We focus on Wagner’s weighted GRN model with varying weights, which is used in evolutionary biology. In the model, weight parameters represent the gene interaction strength that may change due to genetic mutations. For a property of interest, we synthesise the constraints over the parameter space that represent the set of GRNs satisfying the property. We experimentally show that our parameter synthesis procedure computes the mutational robustness of GRNs—an important problem of interest in evolutionary biology—more efficiently than the classical simulation method. We specify the property in linear temporal logic. We employ symbolic bounded model checking and SMT solving to compute the space of GRNs that satisfy the property, which amounts to synthesizing a set of linear constraints on the weights.}, author = {Giacobbe, Mirco and Guet, Calin C and Gupta, Ashutosh and Henzinger, Thomas A and Paixao, Tiago and Petrov, Tatjana}, issn = {00015903}, journal = {Acta Informatica}, number = {8}, pages = {765 -- 787}, publisher = {Springer}, title = {{Model checking the evolution of gene regulatory networks}}, doi = {10.1007/s00236-016-0278-x}, volume = {54}, year = {2017}, } @article{1336, abstract = {Evolutionary algorithms (EAs) form a popular optimisation paradigm inspired by natural evolution. In recent years the field of evolutionary computation has developed a rigorous analytical theory to analyse the runtimes of EAs on many illustrative problems. Here we apply this theory to a simple model of natural evolution. In the Strong Selection Weak Mutation (SSWM) evolutionary regime the time between occurrences of new mutations is much longer than the time it takes for a mutated genotype to take over the population. In this situation, the population only contains copies of one genotype and evolution can be modelled as a stochastic process evolving one genotype by means of mutation and selection between the resident and the mutated genotype. The probability of accepting the mutated genotype then depends on the change in fitness. We study this process, SSWM, from an algorithmic perspective, quantifying its expected optimisation time for various parameters and investigating differences to a similar evolutionary algorithm, the well-known (1+1) EA. We show that SSWM can have a moderate advantage over the (1+1) EA at crossing fitness valleys and study an example where SSWM outperforms the (1+1) EA by taking advantage of information on the fitness gradient.}, author = {Paixao, Tiago and Pérez Heredia, Jorge and Sudholt, Dirk and Trubenova, Barbora}, issn = {01784617}, journal = {Algorithmica}, number = {2}, pages = {681 -- 713}, publisher = {Springer}, title = {{Towards a runtime comparison of natural and artificial evolution}}, doi = {10.1007/s00453-016-0212-1}, volume = {78}, year = {2017}, } @article{1084, abstract = {BceRS and PsdRS are paralogous two-component systems in Bacillus subtilis controlling the response to antimicrobial peptides. In the presence of extracellular bacitracin and nisin, respectively, the two response regulators (RRs) bind their target promoters, PbceA or PpsdA, resulting in a strong up-regulation of target gene expression and ultimately antibiotic resistance. Despite high sequence similarity between the RRs BceR and PsdR and their known binding sites, no cross-regulation has been observed between them. We therefore investigated the specificity determinants of PbceA and PpsdA that ensure the insulation of these two paralogous pathways at the RR–promoter interface. In vivo and in vitro analyses demonstrate that the regulatory regions within these two promoters contain three important elements: in addition to the known (main) binding site, we identified a linker region and a secondary binding site that are crucial for functionality. Initial binding to the high-affinity, low-specificity main binding site is a prerequisite for the subsequent highly specific binding of a second RR dimer to the low-affinity secondary binding site. In addition to this hierarchical cooperative binding, discrimination requires a competition of the two RRs for their respective binding site mediated by only slight differences in binding affinities.}, author = {Fang, Chong and Nagy-Staron, Anna A and Grafe, Martin and Heermann, Ralf and Jung, Kirsten and Gebhard, Susanne and Mascher, Thorsten}, issn = { 0950382X}, journal = {Molecular Microbiology}, number = {1}, pages = {16 -- 31}, publisher = {Wiley-Blackwell}, title = {{Insulation and wiring specificity of BceR like response regulators and their target promoters in Bacillus subtilis}}, doi = {10.1111/mmi.13597}, volume = {104}, year = {2017}, } @article{954, abstract = {Understanding the relation between genotype and phenotype remains a major challenge. The difficulty of predicting individual mutation effects, and particularly the interactions between them, has prevented the development of a comprehensive theory that links genotypic changes to their phenotypic effects. We show that a general thermodynamic framework for gene regulation, based on a biophysical understanding of protein-DNA binding, accurately predicts the sign of epistasis in a canonical cis-regulatory element consisting of overlapping RNA polymerase and repressor binding sites. Sign and magnitude of individual mutation effects are sufficient to predict the sign of epistasis and its environmental dependence. Thus, the thermodynamic model offers the correct null prediction for epistasis between mutations across DNA-binding sites. Our results indicate that a predictive theory for the effects of cis-regulatory mutations is possible from first principles, as long as the essential molecular mechanisms and the constraints these impose on a biological system are accounted for.}, author = {Lagator, Mato and Paixao, Tiago and Barton, Nicholas H and Bollback, Jonathan P and Guet, Calin C}, issn = {2050084X}, journal = {eLife}, publisher = {eLife Sciences Publications}, title = {{On the mechanistic nature of epistasis in a canonical cis-regulatory element}}, doi = {10.7554/eLife.25192}, volume = {6}, year = {2017}, } @article{1007, abstract = {A nonlinear system possesses an invariance with respect to a set of transformations if its output dynamics remain invariant when transforming the input, and adjusting the initial condition accordingly. Most research has focused on invariances with respect to time-independent pointwise transformations like translational-invariance (u(t) -> u(t) + p, p in R) or scale-invariance (u(t) -> pu(t), p in R>0). In this article, we introduce the concept of s0-invariances with respect to continuous input transformations exponentially growing/decaying over time. We show that s0-invariant systems not only encompass linear time-invariant (LTI) systems with transfer functions having an irreducible zero at s0 in R, but also that the input/output relationship of nonlinear s0-invariant systems possesses properties well known from their linear counterparts. Furthermore, we extend the concept of s0-invariances to second- and higher-order s0-invariances, corresponding to invariances with respect to transformations of the time-derivatives of the input, and encompassing LTI systems with zeros of multiplicity two or higher. Finally, we show that nth-order 0-invariant systems realize – under mild conditions – nth-order nonlinear differential operators: when excited by an input of a characteristic functional form, the system’s output converges to a constant value only depending on the nth (nonlinear) derivative of the input.}, author = {Lang, Moritz and Sontag, Eduardo}, issn = {0005-1098}, journal = {Automatica}, pages = {46 -- 55}, publisher = {International Federation of Automatic Control}, title = {{Zeros of nonlinear systems with input invariances}}, doi = {10.1016/j.automatica.2017.03.030}, volume = {81C}, year = {2017}, } @misc{5564, abstract = {Compressed Fastq files with whole-genome sequencing data of IS-wt strain D and clones from four evolved populations (A11, C08, C10, D08). Information on this data collection is available in the Methods Section of the primary publication.}, author = {Steinrück, Magdalena and Guet, Calin C}, publisher = {Institute of Science and Technology Austria}, title = {{Fastq files for "Complex chromosomal neighborhood effects determine the adaptive potential of a gene under selection"}}, doi = {10.15479/AT:ISTA:65}, year = {2017}, } @misc{5560, abstract = {This repository contains the data collected for the manuscript "Biased partitioning of the multi-drug efflux pump AcrAB-TolC underlies long-lived phenotypic heterogeneity". The data is compressed into a single archive. Within the archive, different folders correspond to figures of the main text and the SI of the related publication. Data is saved as plain text, with each folder containing a separate readme file describing the format. Typically, the data is from fluorescence microscopy measurements of single cells growing in a microfluidic "mother machine" device, and consists of relevant values (primarily arbitrary unit or normalized fluorescence measurements, and division times / growth rates) after raw microscopy images have been processed, segmented, and their features extracted, as described in the methods section of the related publication.}, author = {Bergmiller, Tobias and Andersson, Anna M and Tomasek, Kathrin and Balleza, Enrique and Kiviet, Daniel and Hauschild, Robert and Tkacik, Gasper and Guet, Calin C}, keywords = {single cell microscopy, mother machine microfluidic device, AcrAB-TolC pump, multi-drug efflux, Escherichia coli}, publisher = {Institute of Science and Technology Austria}, title = {{Biased partitioning of the multi-drug efflux pump AcrAB-TolC underlies long-lived phenotypic heterogeneity}}, doi = {10.15479/AT:ISTA:53}, year = {2017}, } @article{665, abstract = {The molecular mechanisms underlying phenotypic variation in isogenic bacterial populations remain poorly understood.We report that AcrAB-TolC, the main multidrug efflux pump of Escherichia coli, exhibits a strong partitioning bias for old cell poles by a segregation mechanism that is mediated by ternary AcrAB-TolC complex formation. Mother cells inheriting old poles are phenotypically distinct and display increased drug efflux activity relative to daughters. Consequently, we find systematic and long-lived growth differences between mother and daughter cells in the presence of subinhibitory drug concentrations. A simple model for biased partitioning predicts a population structure of long-lived and highly heterogeneous phenotypes. This straightforward mechanism of generating sustained growth rate differences at subinhibitory antibiotic concentrations has implications for understanding the emergence of multidrug resistance in bacteria.}, author = {Bergmiller, Tobias and Andersson, Anna M and Tomasek, Kathrin and Balleza, Enrique and Kiviet, Daniel and Hauschild, Robert and Tkacik, Gasper and Guet, Calin C}, issn = {00368075}, journal = {Science}, number = {6335}, pages = {311 -- 315}, publisher = {American Association for the Advancement of Science}, title = {{Biased partitioning of the multidrug efflux pump AcrAB TolC underlies long lived phenotypic heterogeneity}}, doi = {10.1126/science.aaf4762}, volume = {356}, year = {2017}, } @article{1028, abstract = {Optogenetics and photopharmacology provide spatiotemporally precise control over protein interactions and protein function in cells and animals. Optogenetic methods that are sensitive to green light and can be used to break protein complexes are not broadly available but would enable multichromatic experiments with previously inaccessible biological targets. Herein, we repurposed cobalamin (vitamin B12) binding domains of bacterial CarH transcription factors for green-light-induced receptor dissociation. In cultured cells, we observed oligomerization-induced cell signaling for the fibroblast growth factor receptor 1 fused to cobalamin-binding domains in the dark that was rapidly eliminated upon illumination. In zebrafish embryos expressing fusion receptors, green light endowed control over aberrant fibroblast growth factor signaling during development. Green-light-induced domain dissociation and light-inactivated receptors will critically expand the optogenetic toolbox for control of biological processes.}, author = {Kainrath, Stephanie and Stadler, Manuela and Gschaider-Reichhart, Eva and Distel, Martin and Janovjak, Harald L}, issn = {14337851}, journal = {Angewandte Chemie - International Edition}, number = {16}, pages = {4608--4611}, publisher = {Wiley-Blackwell}, title = {{Green-light-induced inactivation of receptor signaling using cobalamin-binding domains}}, doi = {10.1002/anie.201611998}, volume = {56}, year = {2017}, } @article{704, abstract = {How the organization of genes on a chromosome shapes adaptation is essential for understanding evolutionary paths. Here, we investigate how adaptation to rapidly increasing levels of antibiotic depends on the chromosomal neighborhood of a drug-resistance gene inserted at different positions of the Escherichia coli chromosome. Using a dual-fluorescence reporter that allows us to distinguish gene amplifications from other up-mutations, we track in real-time adaptive changes in expression of the drug-resistance gene. We find that the relative contribution of several mutation types differs systematically between loci due to properties of neighboring genes: essentiality, expression, orientation, termination, and presence of duplicates. These properties determine rate and fitness effects of gene amplification, deletions, and mutations compromising transcriptional termination. Thus, the adaptive potential of a gene under selection is a system-property with a complex genetic basis that is specific for each chromosomal locus, and it can be inferred from detailed functional and genomic data.}, author = {Steinrück, Magdalena and Guet, Calin C}, issn = {2050084X}, journal = {eLife}, publisher = {eLife Sciences Publications}, title = {{Complex chromosomal neighborhood effects determine the adaptive potential of a gene under selection}}, doi = {10.7554/eLife.25100}, volume = {6}, year = {2017}, }