@article{723,
  abstract     = {Escaping local optima is one of the major obstacles to function optimisation. Using the metaphor of a fitness landscape, local optima correspond to hills separated by fitness valleys that have to be overcome. We define a class of fitness valleys of tunable difficulty by considering their length, representing the Hamming path between the two optima and their depth, the drop in fitness. For this function class we present a runtime comparison between stochastic search algorithms using different search strategies. The (1+1) EA is a simple and well-studied evolutionary algorithm that has to jump across the valley to a point of higher fitness because it does not accept worsening moves (elitism). In contrast, the Metropolis algorithm and the Strong Selection Weak Mutation (SSWM) algorithm, a famous process in population genetics, are both able to cross the fitness valley by accepting worsening moves. We show that the runtime of the (1+1) EA depends critically on the length of the valley while the runtimes of the non-elitist algorithms depend crucially on the depth of the valley. Moreover, we show that both SSWM and Metropolis can also efficiently optimise a rugged function consisting of consecutive valleys.},
  author       = {Oliveto, Pietro and Paixao, Tiago and Pérez Heredia, Jorge and Sudholt, Dirk and Trubenova, Barbora},
  journal      = {Algorithmica},
  number       = {5},
  pages        = {1604 -- 1633},
  publisher    = {Springer},
  title        = {{How to escape local optima in black box optimisation when non elitism outperforms elitism}},
  doi          = {10.1007/s00453-017-0369-2},
  volume       = {80},
  year         = {2018},
}

@article{457,
  abstract     = {Temperate bacteriophages integrate in bacterial genomes as prophages and represent an important source of genetic variation for bacterial evolution, frequently transmitting fitness-augmenting genes such as toxins responsible for virulence of major pathogens. However, only a fraction of bacteriophage infections are lysogenic and lead to prophage acquisition, whereas the majority are lytic and kill the infected bacteria. Unless able to discriminate lytic from lysogenic infections, mechanisms of immunity to bacteriophages are expected to act as a double-edged sword and increase the odds of survival at the cost of depriving bacteria of potentially beneficial prophages. We show that although restriction-modification systems as mechanisms of innate immunity prevent both lytic and lysogenic infections indiscriminately in individual bacteria, they increase the number of prophage-acquiring individuals at the population level. We find that this counterintuitive result is a consequence of phage-host population dynamics, in which restriction-modification systems delay infection onset until bacteria reach densities at which the probability of lysogeny increases. These results underscore the importance of population-level dynamics as a key factor modulating costs and benefits of immunity to temperate bacteriophages},
  author       = {Pleska, Maros and Lang, Moritz and Refardt, Dominik and Levin, Bruce and Guet, Calin C},
  journal      = {Nature Ecology and Evolution},
  number       = {2},
  pages        = {359 -- 366},
  publisher    = {Springer Nature},
  title        = {{Phage-host population dynamics promotes prophage acquisition in bacteria with innate immunity}},
  doi          = {10.1038/s41559-017-0424-z},
  volume       = {2},
  year         = {2018},
}

@article{503,
  abstract     = {Buffers are essential for diluting bacterial cultures for flow cytometry analysis in order to study bacterial physiology and gene expression parameters based on fluorescence signals. Using a variety of constitutively expressed fluorescent proteins in Escherichia coli K-12 strain MG1655, we found strong artifactual changes in fluorescence levels after dilution into the commonly used flow cytometry buffer phosphate-buffered saline (PBS) and two other buffer solutions, Tris-HCl and M9 salts. These changes appeared very rapidly after dilution, and were linked to increased membrane permeability and loss in cell viability. We observed buffer-related effects in several different E. coli strains, K-12, C and W, but not E. coli B, which can be partially explained by differences in lipopolysaccharide (LPS) and outer membrane composition. Supplementing the buffers with divalent cations responsible for outer membrane stability, Mg2+ and Ca2+, preserved fluorescence signals, membrane integrity and viability of E. coli. Thus, stabilizing the bacterial outer membrane is essential for precise and unbiased measurements of fluorescence parameters using flow cytometry.},
  author       = {Tomasek, Kathrin and Bergmiller, Tobias and Guet, Calin C},
  journal      = {Journal of Biotechnology},
  pages        = {40 -- 52},
  publisher    = {Elsevier},
  title        = {{Lack of cations in flow cytometry buffers affect fluorescence signals by reducing membrane stability and viability of Escherichia coli strains}},
  doi          = {10.1016/j.jbiotec.2018.01.008},
  volume       = {268},
  year         = {2018},
}

@misc{5569,
  abstract     = {Nela Nikolic, Tobias Bergmiller, Alexandra Vandervelde, Tanino G. Albanese, Lendert Gelens, and Isabella Moll (2018)
“Autoregulation of mazEF expression underlies growth heterogeneity in bacterial populations” Nucleic Acids Research, doi: 10.15479/AT:ISTA:74;
microscopy experiments by Tobias Bergmiller; image and data analysis by Nela Nikolic.},
  author       = {Bergmiller, Tobias and Nikolic, Nela},
  keywords     = {microscopy, microfluidics},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{Time-lapse microscopy data}},
  doi          = {10.15479/AT:ISTA:74},
  year         = {2018},
}

@misc{5585,
  abstract     = {Mean repression values and standard error of the mean are given for all operator mutant libraries.},
  author       = {Igler, Claudia and Lagator, Mato and Tkacik, Gasper and Bollback, Jonathan P and Guet, Calin C},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{Data for the paper Evolutionary potential of transcription factors for gene regulatory rewiring}},
  doi          = {10.15479/AT:ISTA:108},
  year         = {2018},
}

@article{5984,
  abstract     = {G-protein-coupled receptors (GPCRs) form the largest receptor family, relay environmental stimuli to changes in cell behavior and represent prime drug targets. Many GPCRs are classified as orphan receptors because of the limited knowledge on their ligands and coupling to cellular signaling machineries. Here, we engineer a library of 63 chimeric receptors that contain the signaling domains of human orphan and understudied GPCRs functionally linked to the light-sensing domain of rhodopsin. Upon stimulation with visible light, we identify activation of canonical cell signaling pathways, including cAMP-, Ca2+-, MAPK/ERK-, and Rho-dependent pathways, downstream of the engineered receptors. For the human pseudogene GPR33, we resurrect a signaling function that supports its hypothesized role as a pathogen entry site. These results demonstrate that substituting unknown chemical activators with a light switch can reveal information about protein function and provide an optically controlled protein library for exploring the physiology and therapeutic potential of understudied GPCRs.},
  author       = {Morri, Maurizio and Sanchez-Romero, Inmaculada and Tichy, Alexandra-Madelaine and Kainrath, Stephanie and Gerrard, Elliot J. and Hirschfeld, Priscila and Schwarz, Jan and Janovjak, Harald L},
  issn         = {2041-1723},
  journal      = {Nature Communications},
  number       = {1},
  publisher    = {Springer Nature},
  title        = {{Optical functionalization of human class A orphan G-protein-coupled receptors}},
  doi          = {10.1038/s41467-018-04342-1},
  volume       = {9},
  year         = {2018},
}

@article{161,
  abstract     = {Which properties of metabolic networks can be derived solely from stoichiometry? Predictive results have been obtained by flux balance analysis (FBA), by postulating that cells set metabolic fluxes to maximize growth rate. Here we consider a generalization of FBA to single-cell level using maximum entropy modeling, which we extend and test experimentally. Specifically, we define for Escherichia coli metabolism a flux distribution that yields the experimental growth rate: the model, containing FBA as a limit, provides a better match to measured fluxes and it makes a wide range of predictions: on flux variability, regulation, and correlations; on the relative importance of stoichiometry vs. optimization; on scaling relations for growth rate distributions. We validate the latter here with single-cell data at different sub-inhibitory antibiotic concentrations. The model quantifies growth optimization as emerging from the interplay of competitive dynamics in the population and regulation of metabolism at the level of single cells.},
  author       = {De Martino, Daniele and Mc, Andersson Anna and Bergmiller, Tobias and Guet, Calin C and Tkacik, Gasper},
  journal      = {Nature Communications},
  number       = {1},
  publisher    = {Springer Nature},
  title        = {{Statistical mechanics for metabolic networks during steady state growth}},
  doi          = {10.1038/s41467-018-05417-9},
  volume       = {9},
  year         = {2018},
}

@misc{9810,
  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},
  publisher    = {Public Library of Science},
  title        = {{Numerical data used in figures}},
  doi          = {10.1371/journal.pbio.2005971.s008},
  year         = {2018},
}

@article{438,
  abstract     = {The MazF toxin sequence-specifically cleaves single-stranded RNA upon various stressful conditions, and it is activated as a part of the mazEF toxin–antitoxin module in Escherichia coli. Although autoregulation of mazEF expression through the MazE antitoxin-dependent transcriptional repression has been biochemically characterized, less is known about post-transcriptional autoregulation, as well as how both of these autoregulatory features affect growth of single cells during conditions that promote MazF production. Here, we demonstrate post-transcriptional autoregulation of mazF expression dynamics by MazF cleaving its own transcript. Single-cell analyses of bacterial populations during ectopic MazF production indicated that two-level autoregulation of mazEF expression influences cell-to-cell growth rate heterogeneity. The increase in growth rate heterogeneity is governed by the MazE antitoxin, and tuned by the MazF-dependent mazF mRNA cleavage. Also, both autoregulatory features grant rapid exit from the stress caused by mazF overexpression. Time-lapse microscopy revealed that MazF-mediated cleavage of mazF mRNA leads to increased temporal variability in length of individual cells during ectopic mazF overexpression, as explained by a stochastic model indicating that mazEF mRNA cleavage underlies temporal fluctuations in MazF levels during stress.},
  author       = {Nikolic, Nela and Bergmiller, Tobias and Vandervelde, Alexandra and Albanese, Tanino and Gelens, Lendert and Moll, Isabella},
  journal      = {Nucleic Acids Research},
  number       = {6},
  pages        = {2918--2931},
  publisher    = {Oxford University Press},
  title        = {{Autoregulation of mazEF expression underlies growth heterogeneity in bacterial populations}},
  doi          = {10.1093/nar/gky079},
  volume       = {46},
  year         = {2018},
}

@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},
}

@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{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{696,
  abstract     = {Mutator strains are expected to evolve when the availability and effect of beneficial mutations are high enough to counteract the disadvantage from deleterious mutations that will inevitably accumulate. As the population becomes more adapted to its environment, both availability and effect of beneficial mutations necessarily decrease and mutation rates are predicted to decrease. It has been shown that certain molecular mechanisms can lead to increased mutation rates when the organism finds itself in a stressful environment. While this may be a correlated response to other functions, it could also be an adaptive mechanism, raising mutation rates only when it is most advantageous. Here, we use a mathematical model to investigate the plausibility of the adaptive hypothesis. We show that such a mechanism can be mantained if the population is subjected to diverse stresses. By simulating various antibiotic treatment schemes, we find that combination treatments can reduce the effectiveness of second-order selection on stress-induced mutagenesis. We discuss the implications of our results to strategies of antibiotic therapy.},
  author       = {Lukacisinova, Marta and Novak, Sebastian and Paixao, Tiago},
  issn         = {1553734X},
  journal      = {PLoS Computational Biology},
  number       = {7},
  publisher    = {Public Library of Science},
  title        = {{Stress induced mutagenesis: Stress diversity facilitates the persistence of mutator genes}},
  doi          = {10.1371/journal.pcbi.1005609},
  volume       = {13},
  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},
}

@article{735,
  abstract     = {Cell-cell contact formation constitutes an essential step in evolution, leading to the differentiation of specialized cell types. However, remarkably little is known about whether and how the interplay between contact formation and fate specification affects development. Here, we identify a positive feedback loop between cell-cell contact duration, morphogen signaling, and mesendoderm cell-fate specification during zebrafish gastrulation. We show that long-lasting cell-cell contacts enhance the competence of prechordal plate (ppl) progenitor cells to respond to Nodal signaling, required for ppl cell-fate specification. We further show that Nodal signaling promotes ppl cell-cell contact duration, generating a positive feedback loop between ppl cell-cell contact duration and cell-fate specification. Finally, by combining mathematical modeling and experimentation, we show that this feedback determines whether anterior axial mesendoderm cells become ppl or, instead, turn into endoderm. Thus, the interdependent activities of cell-cell signaling and contact formation control fate diversification within the developing embryo.},
  author       = {Barone, Vanessa and Lang, Moritz and Krens, Gabriel and Pradhan, Saurabh and Shamipour, Shayan and Sako, Keisuke and Sikora, Mateusz K and Guet, Calin C and Heisenberg, Carl-Philipp J},
  issn         = {15345807},
  journal      = {Developmental Cell},
  number       = {2},
  pages        = {198 -- 211},
  publisher    = {Cell Press},
  title        = {{An effective feedback loop between cell-cell contact duration and morphogen signaling determines cell fate}},
  doi          = {10.1016/j.devcel.2017.09.014},
  volume       = {43},
  year         = {2017},
}

@article{538,
  abstract     = {Optogenetik und Photopharmakologie ermöglichen präzise räumliche und zeitliche Kontrolle von Proteinwechselwirkung und -funktion in Zellen und Tieren. Optogenetische Methoden, die auf grünes Licht ansprechen und zum Trennen von Proteinkomplexen geeignet sind, sind nichtweitläufig verfügbar, würden jedoch mehrfarbige Experimente zur Beantwortung von biologischen Fragestellungen ermöglichen. Hier demonstrieren wir die Verwendung von Cobalamin(Vitamin B12)-bindenden Domänen von bakteriellen CarH-Transkriptionsfaktoren zur Grünlicht-induzierten Dissoziation von Rezeptoren. Fusioniert mit dem Fibroblasten-W achstumsfaktor-Rezeptor 1 führten diese im Dunkeln in kultivierten Zellen zu Signalaktivität durch Oligomerisierung, welche durch Beleuchten umgehend aufgehoben wurde. In Zebrafischembryonen, die einen derartigen Rezeptor exprimieren, ermöglichte grünes Licht die Kontrolle über abnormale Signalaktivität während der Embryonalentwicklung. },
  author       = {Kainrath, Stephanie and Stadler, Manuela and Gschaider-Reichhart, Eva and Distel, Martin and Janovjak, Harald L},
  journal      = {Angewandte Chemie},
  number       = {16},
  pages        = {4679 -- 4682},
  publisher    = {Wiley},
  title        = {{Grünlicht-induzierte Rezeptorinaktivierung durch Cobalamin-bindende Domänen}},
  doi          = {10.1002/ange.201611998},
  volume       = {129},
  year         = {2017},
}

@article{541,
  abstract     = {While we have good understanding of bacterial metabolism at the population level, we know little about the metabolic behavior of individual cells: do single cells in clonal populations sometimes specialize on different metabolic pathways? Such metabolic specialization could be driven by stochastic gene expression and could provide individual cells with growth benefits of specialization. We measured the degree of phenotypic specialization in two parallel metabolic pathways, the assimilation of glucose and arabinose. We grew Escherichia coli in chemostats, and used isotope-labeled sugars in combination with nanometer-scale secondary ion mass spectrometry and mathematical modeling to quantify sugar assimilation at the single-cell level. We found large variation in metabolic activities between single cells, both in absolute assimilation and in the degree to which individual cells specialize in the assimilation of different sugars. Analysis of transcriptional reporters indicated that this variation was at least partially based on cell-to-cell variation in gene expression. Metabolic differences between cells in clonal populations could potentially reduce metabolic incompatibilities between different pathways, and increase the rate at which parallel reactions can be performed.},
  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},
  issn         = {15537390},
  journal      = {PLoS Genetics},
  number       = {12},
  publisher    = {Public Library of Science},
  title        = {{Cell-to-cell variation and specialization in sugar metabolism in clonal bacterial populations}},
  doi          = {10.1371/journal.pgen.1007122},
  volume       = {13},
  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},
}

@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},
}

@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},
}

