Bistability in a metabolic network underpins the de novo evolution of colony switching in Pseudomonas fluorescens

Evolution of a bistable genetic system in fluctuating and nonfluctuating environments

Epigenetic mechanisms can generate bacterial lineages capable of spontaneously switching between distinct phenotypes. Currently, mathematical models and simulations propose epigenetic switches as a mechanism of adaptation to deal with fluctuating environments. However, bacterial evolution experiments for testing these predictions are lacking. Here, we exploit an epigenetic switch in Salmonella enterica, the opvAB operon, to show clear evidence that OpvAB bistability persists in changing environments but not in stable conditions. Epigenetic control of transcription in the opvAB operon produces OpvABOFF (phage-sensitive) and OpvABON (phage-resistant) cells in a reversible manner and may be interpreted as an example of bet-hedging to preadapt Salmonella populations to the encounter with phages. Our experimental observations and computational simulations illustrate the adaptive value of epigenetic variation as an evolutionary strategy for mutation avoidance in fluctuating environments. In addition, our study provides experimental support to game theory models predicting that phenotypic heterogeneity is advantageous in changing and unpredictable environments.

Evidence for Rho-dependent control of a virulence switch in Acinetobacter baumannii

IMPORTANCE

Acinetobacter baumannii is a significant cause of infections in the healthcare setting. More recently, A. baumannii has been a leading cause of secondary bacterial pneumonia in patients infected with SARS-CoV-2 and the overall frequency of A. baumannii infection increased 78% during the COVID-19 pandemic. A. baumannii can exist in virulent or avirulent subpopulations and this interconversion is mediated by the expression of a family of TetR-type transcriptional regulators. In this study, we demonstrate that Rho is a key regulatory component in the expression of these TetR regulators. Overall, this study is the first to address a role for Rho in A. baumannii and provides additional evidence for the role of Rho in regulating diversity in bacterial subpopulations.

A mathematical model captures the role of adenyl cyclase Cyr1 and guanidine exchange factor Ira2 in creating a growth-to-hyphal bistable switch in Candida albicans

Recent biochemical experiments have indicated that in Candida albicans, a commensal fungal pathogen, the Ras signaling pathway plays a significant role in the yeast-to-hyphal transition; specifically, two enzymes in this pathway, Adenyl Cyclase Cyr1 and GTPase activating protein Ira2, facilitate this transition, in the presence of energy sensor ATP. However, the precise mechanism by which protein interactions between Ira2 and Cyr1 and the energy sensor ATP result in the yeast-to-hyphal transition and create a switch-like process are unknown. We propose a new set of biochemical reaction steps that captures all the essential interactions between Ira2, Cyr1, and ATP in the Ras pathway. With the help of chemical reaction network theory, we demonstrate that this set of biochemical reaction steps results in bistability. Further, bifurcation analysis of the differential equations based on this set of reaction steps supports the existence of a bistable switch, and this switch may act as a checkpoint mechanism for the promotion of growth-to-hyphal transition in C. albicans.

Adaptation to novel spatially-structured environments is driven by the capsule and alters virulence-associated traits

The extracellular capsule is a major virulence factor, but its ubiquity in free-living bacteria with large environmental breadths suggests that it shapes adaptation to novel niches. Yet, how it does so, remains unexplored. Here, we evolve three Klebsiella strains and their capsule mutants in parallel. Their comparison reveals different phenotypic and genotypic evolutionary changes that alter virulence-associated traits. Non-capsulated populations accumulate mutations that reduce exopolysaccharide production and increase biofilm formation and yield, whereas most capsulated populations become hypermucoviscous, a signature of hypervirulence. Hence, adaptation to novel environments primarily occurs by fine-tuning expression of the capsular locus. The same evolutionary conditions selecting for mutations in the capsular gene wzc leading to hypermucoviscosity also result in increased susceptibility to antibiotics by mutations in the ramA regulon. This implies that general adaptive processes outside the host can affect capsule evolution and its role in virulence and infection outcomes may be a by-product of such adaptation.

Phenotypic and genotypic evolution in worrisome Klebsiella spp. is influenced by the capsule. Here the authors show that adaptation outside the host can impact virulence-associated traits, including de novo emergence of hypermucoviscosity.

Evolutionary flexibility in routes to mat formation by Pseudomonas

Many bacteria form mats at the air-liquid interface of static microcosms. These structures typically involve the secretion of exopolysaccharides, the production of which is often controlled by the secondary messenger c-di-GMP. Mechanisms of mat formation have been particularly well characterized in Pseudomonas fluorescens SBW25; stimuli or mutations that increase c-di-GMP production by diguanylate cyclases (WspR, AwsR, and MwsR) result in the secretion of cellulose and mat formation. Here, we characterize and compare mat formation in two close relatives of SBW25: Pseudomonas simiae PICF7 and P. fluorescens A506. We find that PICF7-the strain more closely related to SBW25-can form mats through mutations affecting the activity of the same three diguanylate cyclases as SBW25. However, instead of cellulose, these mutations activate production of the exopolysaccharide Pel. We also provide evidence for at least two further-as yet uncharacterized-routes to mat formation by PICF7. P. fluorescens A506, while retaining the same mutational routes to mat formation as SBW25 and PICF7, preferentially forms mats by a semi-heritable mechanism that culminates in Psl and Pga over-production. Our results demonstrate a high level of evolutionary flexibility in the molecular and structural routes to mat formation, even among close relatives.

Can gene-inactivating mutations lead to evolutionary novelty?

Evolutionary novelty is difficult to define. It typically involves shifts in organismal or biochemical phenotypes that can be seen as qualitative as well as quantitative changes. In laboratory-based experimental evolution of novel phenotypes and the human domestication of crops, the majority of the mutations that lead to adaptation are loss-of-function mutations that impair or eliminate the function of genes rather than gain-of-function mutations that increase or qualitatively alter the function of proteins. Here, I speculate that easier access to loss-of-function mutations has led them to play a major role in the adaptive radiations that occur when populations have access to many unoccupied ecological niches. I discuss five possible objections to this claim: that genes can only survive if they confer benefits to the organisms that bear them, antagonistic pleiotropy, the importance of pre-existing genetic variation in populations, the danger that adaptation by breaking genes will, over long times, cause organisms to run out of genes, and the recessive nature of most loss-of-function mutations.

Forecasting of phenotypic and genetic outcomes of experimental evolution in Pseudomonas protegens

Experimental evolution with microbes is often highly repeatable under identical conditions, suggesting the possibility to predict short-term evolution. However, it is not clear to what degree evolutionary forecasts can be extended to related species in non-identical environments, which would allow testing of general predictive models and fundamental biological assumptions. To develop an extended model system for evolutionary forecasting, we used previous data and models of the genotype-to-phenotype map from the wrinkly spreader system in Pseudomonas fluorescens SBW25 to make predictions of evolutionary outcomes on different biological levels for Pseudomonas protegens Pf-5. In addition to sequence divergence (78% amino acid and 81% nucleotide identity) for the genes targeted by mutations, these species also differ in the inability of Pf-5 to make cellulose, which is the main structural basis for the adaptive phenotype in SBW25. The experimental conditions were changed compared to the SBW25 system to test if forecasts were extendable to a non-identical environment. Forty-three mutants with increased ability to colonize the air-liquid interface were isolated, and the majority had reduced motility and was partly dependent on the Pel exopolysaccharide as a structural component. Most (38/43) mutations are expected to disrupt negative regulation of the same three diguanylate cyclases as in SBW25, with a smaller number of mutations in promoter regions, including an uncharacterized polysaccharide synthase operon. A mathematical model developed for SBW25 predicted the order of the three main pathways and the genes targeted by mutations, but differences in fitness between mutants and mutational biases also appear to influence outcomes. Mutated regions in proteins could be predicted in most cases (16/22), but parallelism at the nucleotide level was low and mutational hot spot sites were not conserved. This study demonstrates the potential of short-term evolutionary forecasting in experimental populations and provides testable predictions for evolutionary outcomes in other Pseudomonas species.

The role of memory in non-genetic inheritance and its impact on cancer treatment resistance

Intra-tumour heterogeneity is a leading cause of treatment failure and disease progression in cancer. While genetic mutations have long been accepted as a primary mechanism of generating this heterogeneity, the role of phenotypic plasticity is becoming increasingly apparent as a driver of intra-tumour heterogeneity. Consequently, understanding the role of this plasticity in treatment resistance and failure is a key component of improving cancer therapy. We develop a mathematical model of stochastic phenotype switching that tracks the evolution of drug-sensitive and drug-tolerant subpopulations to clarify the role of phenotype switching on population growth rates and tumour persistence. By including cytotoxic therapy in the model, we show that, depending on the strategy of the drug-tolerant subpopulation, stochastic phenotype switching can lead to either transient or permanent drug resistance. We study the role of phenotypic heterogeneity in a drug-resistant, genetically homogeneous population of non-small cell lung cancer cells to derive a rational treatment schedule that drives population extinction and avoids competitive release of the drug-tolerant sub-population. This model-informed therapeutic schedule results in increased treatment efficacy when compared against periodic therapy, and, most importantly, sustained tumour decay without the development of resistance.

Compost spatial heterogeneity promotes evolutionary diversification of a bacterium

Spatial resource heterogeneity is expected to be a key driver for the evolution of diversity. However, direct empirical support for this prediction is limited to studies carried out in simplified laboratory environments. Here, we investigate how altering spatial heterogeneity of potting compost-by the addition of water and mixing-affects the evolutionary diversification of a bacterial species, Pseudomonas fluorescens, that is naturally found in the environment. There was a greater propensity of resource specialists to evolve in the unmanipulated compost, while more generalist phenotypes dominated the compost-water mix. Genomic data were consistent with these phenotypic findings. Competition experiments strongly suggest these results are due to diversifying selection as a result of resource heterogeneity, as opposed to other covariables. Overall, our findings corroborate theoretical and in vitro findings, but in semi-natural, more realistic conditions.

In Silico Characterization and Phylogenetic Distribution of Extracellular Matrix Components in the Model Rhizobacteria Pseudomonas fluorescens F113 and Other Pseudomonads

Biofilms are complex structures that are crucial during host-bacteria interaction and colonization. Bacteria within biofilms are surrounded by an extracellular matrix (ECM) typically composed of proteins, polysaccharides, lipids, and DNA. Pseudomonads contain a variety of ECM components, some of which have been extensively characterized. However, neither the ECM composition of plant-associated pseudomonads nor their phylogenetic distribution within the genus has been so thoroughly studied. In this work, we use in silico methods to describe the ECM composition of Pseudomonas fluorescens F113, a plant growth-promoting rhizobacteria and model for rhizosphere colonization. These components include the polysaccharides alginate, poly-N-acetyl-glucosamine (PNAG) and levan; the adhesins LapA, MapA and PsmE; and the functional amyloids in Pseudomonas. Interestingly, we identified novel components: the Pseudomonas acidic polysaccharide (Pap), whose presence is limited within the genus; and a novel type of Flp/Tad pilus, partially different from the one described in P. aeruginosa. Furthermore, we explored the phylogenetic distribution of the most relevant ECM components in nearly 600 complete Pseudomonas genomes. Our analyses show that Pseudomonas populations contain a diverse set of gene/gene clusters potentially involved in the formation of their ECMs, showing certain commensal versus pathogen lifestyle specialization.

The birth of a bacterial tRNA gene by large-scale, tandem duplication events

Organisms differ in the types and numbers of tRNA genes that they carry. While the evolutionary mechanisms behind tRNA gene set evolution have been investigated theoretically and computationally, direct observations of tRNA gene set evolution remain rare. Here, we report the evolution of a tRNA gene set in laboratory populations of the bacterium Pseudomonas fluorescens SBW25. The growth defect caused by deleting the single-copy tRNA gene, serCGA, is rapidly compensated by large-scale (45–290 kb) duplications in the chromosome. Each duplication encompasses a second, compensatory tRNA gene (serTGA) and is associated with a rise in tRNA-Ser(UGA) in the mature tRNA pool. We postulate that tRNA-Ser(CGA) elimination increases the translational demand for tRNA-Ser(UGA), a pressure relieved by increasing serTGA copy number. This work demonstrates that tRNA gene sets can evolve through duplication of existing tRNA genes, a phenomenon that may contribute to the presence of multiple, identical tRNA gene copies within genomes.

Ecological scaffolding and the evolution of individuality

Evolutionary transitions in individuality are central to the emergence of biological complexity. Recent experiments provide glimpses of processes underpinning the transition from single cells to multicellular life and draw attention to the critical role of ecology. Here, we emphasize this ecological dimension and argue that its current absence from theoretical frameworks hampers development of general explanatory solutions. Using mechanistic mathematical models, we show how a minimal ecological structure comprising patchily distributed resources and between-patch dispersal can scaffold Darwinian-like properties on collectives of cells. This scaffolding causes cells to participate directly in the process of evolution by natural selection as if they were members of multicellular collectives, with collectives participating in a death-birth process arising from the interplay between the timing of dispersal events and the rate of resource use by cells. When this timescale is sufficiently long and new collectives are founded by single cells, collectives experience conditions that favour evolution of a reproductive division of labour. Together our simple model makes explicit key events in the major evolutionary transition to multicellularity. It also makes predictions concerning the life history of certain pathogens and serves as an ecological recipe for experimental realization of evolutionary transitions.

Microbial phenotypic heterogeneity in response to a metabolic toxin: Continuous, dynamically shifting distribution of formaldehyde tolerance in Methylobacterium extorquens populations

While microbiologists often make the simplifying assumption that genotype determines phenotype in a given environment, it is becoming increasingly apparent that phenotypic heterogeneity (in which one genotype generates multiple phenotypes simultaneously even in a uniform environment) is common in many microbial populations. The importance of phenotypic heterogeneity has been demonstrated in a number of model systems involving binary phenotypic states (e.g., growth/non-growth); however, less is known about systems involving phenotype distributions that are continuous across an environmental gradient, and how those distributions change when the environment changes. Here, we describe a novel instance of phenotypic diversity in tolerance to a metabolic toxin within wild-type populations of Methylobacterium extorquens, a ubiquitous phyllosphere methylotroph capable of growing on the methanol periodically released from plant leaves. The first intermediate in methanol metabolism is formaldehyde, a potent cellular toxin that is lethal in high concentrations. We have found that at moderate concentrations, formaldehyde tolerance in M. extorquens is heterogeneous, with a cell's minimum tolerance level ranging between 0 mM and 8 mM. Tolerant cells have a distinct gene expression profile from non-tolerant cells. This form of heterogeneity is continuous in terms of threshold (the formaldehyde concentration where growth ceases), yet binary in outcome (at a given formaldehyde concentration, cells either grow normally or die, with no intermediate phenotype), and it is not associated with any detectable genetic mutations. Moreover, tolerance distributions within the population are dynamic, changing over time in response to growth conditions. We characterized this phenomenon using bulk liquid culture experiments, colony growth tracking, flow cytometry, single-cell time-lapse microscopy, transcriptomics, and genome resequencing. Finally, we used mathematical modeling to better understand the processes by which cells change phenotype, and found evidence for both stochastic, bidirectional phenotypic diversification and responsive, directed phenotypic shifts, depending on the growth substrate and the presence of toxin.

Scientists tend to appreciate microbes for their simplicity and predictability: a population of genetically identical cells inhabiting a uniform environment is expected to behave in a uniform way. However, counter-examples to this assumption are frequently being discovered, forcing a re-examination of the relationship between genotype and phenotype. In most such examples, bacterial cells are found to split into two discrete populations, for instance growing and non-growing. Here, we report the discovery of a novel example of microbial phenotypic heterogeneity in which cells are distributed along a gradient of phenotypes, ranging from low to high tolerance of a toxic chemical. Furthermore, we demonstrate that the distribution of phenotypes changes in different growth conditions, and we use mathematical modeling to show that cells may change their phenotype either randomly or in a particular direction in response to the environment. Our work expands our understanding of how a bacterial cell's genome, family history, and environment all contribute to its behavior, with implications for the diverse situations in which we care to understand the growth of any single-celled populations.

Epigenetic switching as a strategy for quick adaptation while attenuating biochemical noise

Epigenetic switches are bistable, molecular systems built from self-reinforcing feedback loops that can spontaneously switch between heritable phenotypes in the absence of DNA mutation. It has been hypothesized that epigenetic switches first evolved as a mechanism of bet-hedging and adaptation, but the evolutionary trajectories and conditions by which an epigenetic switch can outcompete adaptation through genetic mutation remain unknown. Here, we used computer simulations to evolve a mechanistic, biophysical model of a self-activating genetic circuit, which can both adapt genetically through mutation and exhibit epigenetic switching. We evolved these genetic circuits under a fluctuating environment that alternatively selected for low and high protein expression levels. In all tested conditions, the population first evolved by genetic mutation towards a region of genotypes where genetic adaptation can occur faster after each environmental transition. Once in this region, the self-activating genetic circuit can exhibit epigenetic switching, which starts competing with genetic adaptation. We show a trade-off between either minimizing the adaptation time or increasing the robustness of the phenotype to biochemical noise. Epigenetic switching was superior in a fast fluctuating environment because it adapted faster than genetic mutation after an environmental transition, while still attenuating the effect of biochemical noise on the phenotype. Conversely, genetic adaptation was favored in a slowly fluctuating environment because it maximized the phenotypic robustness to biochemical noise during the constant environment between transitions, even if this resulted in slower adaptation. This simple trade-off predicts the conditions and trajectories under which an epigenetic switch evolved to outcompete genetic adaptation, shedding light on possible mechanisms by which bet-hedging strategies might emerge and persist in natural populations.

Epigenetic switches regulate cell fate decisions during development in multicellular organisms, but their origin predates multicellularity because they are found in viruses, bacteria, and unicellular eukaryotes. It has been suggested that epigenetic switches first evolved as a mechanism of bet-hedging and adaptation to fluctuating environments. To discern the evolutionary pressures that select for epigenetic switches, we used computer simulations to evolve a mechanistic, biophysical model of a self-activating genetic circuit, which can both adapt genetically and exhibit epigenetic switching. Unlike laboratory evolution experiments, this in silico experiment was run many times over a range of evolutionary parameters (population size, selection pressure, mutation step-size, fluctuation frequency) and different model assumptions to uncover statistical regularities in the evolutionary trajectories. Using this computational approach, we could elucidate simple principles that predict the conditions that favor adaptation by epigenetic switching over genetic mutation in a fluctuating environment.

Ribosome Provisioning Activates a Bistable Switch Coupled to Fast Exit from Stationary Phase

Observations of bacteria at the single-cell level have revealed many instances of phenotypic heterogeneity within otherwise clonal populations, but the selective causes, molecular bases, and broader ecological relevance remain poorly understood. In an earlier experiment in which the bacterium Pseudomonas fluorescens SBW25 was propagated under a selective regime that mimicked the host immune response, a genotype evolved that stochastically switched between capsulation states. The genetic cause was a mutation in carB that decreased the pyrimidine pool (and growth rate), lowering the activation threshold of a preexisting but hitherto unrecognized phenotypic switch. Genetic components surrounding bifurcation of UTP flux toward DNA/RNA or UDP-glucose (a precursor of colanic acid forming the capsules) were implicated as key components. Extending these molecular analyses-and based on a combination of genetics, transcriptomics, biochemistry, and mathematical modeling-we show that pyrimidine limitation triggers an increase in ribosome biosynthesis and that switching is caused by competition between ribosomes and CsrA/RsmA proteins for the mRNA transcript of a positively autoregulated activator of colanic acid biosynthesis. We additionally show that in the ancestral bacterium the switch is part of a program that determines stochastic entry into a semiquiescent capsulated state, ensures that such cells are provisioned with excess ribosomes, and enables provisioned cells to exit rapidly from stationary phase under permissive conditions.

Repeated Phenotypic Evolution by Different Genetic Routes in Pseudomonas fluorescens SBW25

Repeated evolution of functionally similar phenotypes is observed throughout the tree of life. The extent to which the underlying genetics are conserved remains an area of considerable interest. Previously, we reported the evolution of colony switching in two independent lineages of Pseudomonas fluorescens SBW25. The phenotypic and genotypic bases of colony switching in the first lineage (Line 1) have been described elsewhere. Here, we deconstruct the evolution of colony switching in the second lineage (Line 6). We show that, as for Line 1, Line 6 colony switching results from an increase in the expression of a colanic acid-like polymer (CAP). At the genetic level, nine mutations occur in Line 6. Only one of these—a nonsynonymous point mutation in the housekeeping sigma factor rpoD—is required for colony switching. In contrast, the genetic basis of colony switching in Line 1 is a mutation in the metabolic gene carB. A molecular model has recently been proposed whereby the carB mutation increases capsulation by redressing the intracellular balance of positive (ribosomes) and negative (RsmAE/CsrA) regulators of a positive feedback loop in capsule expression. We show that Line 6 colony switching is consistent with this model; the rpoD mutation generates an increase in ribosomal gene expression, and ultimately an increase in CAP expression.

Shortsighted Evolution Constrains the Efficacy of Long-Term Bet Hedging

To survive unpredictable environmental change, many organisms adopt bet-hedging strategies that are initially costly but provide a long-term fitness benefit. The temporal extent of these deferred fitness benefits determines whether bet-hedging organisms can survive long enough to realize them. In this article, we examine a model of microbial bet hedging in which there are two paths to extinction: unpredictable environmental change and demographic stochasticity. In temporally correlated environments, these drivers of extinction select for different switching strategies. Rapid phenotype switching ensures survival in the face of unpredictable environmental change, while slower-switching organisms become extinct. However, when both switching strategies are present in the same population, then demographic stochasticity-enforced by a limited population size-leads to extinction of the faster-switching organism. As a result, we find a novel form of evolutionary suicide whereby selection in a fluctuating environment can favor bet-hedging strategies that ultimately increase the risk of extinction. Population structures with multiple subpopulations and dispersal can reduce the risk of extinction from unpredictable environmental change and shift the balance so as to facilitate the evolution of slower-switching organisms.

Genetic exchanges are more frequent in bacteria encoding capsules

Capsules allow bacteria to colonize novel environments, to withstand numerous stresses, and to resist antibiotics. Yet, even though genetic exchanges with other cells should be adaptive under such circumstances, it has been suggested that capsules lower the rates of homologous recombination and horizontal gene transfer. We analysed over one hundred pan-genomes and thousands of bacterial genomes for the evidence of an association between genetic exchanges (or lack thereof) and the presence of a capsule system. We found that bacteria encoding capsules have larger pan-genomes, higher rates of horizontal gene transfer, and higher rates of homologous recombination in their core genomes. Accordingly, genomes encoding capsules have more plasmids, conjugative elements, transposases, prophages, and integrons. Furthermore, capsular loci are frequent in plasmids, and can be found in prophages. These results are valid for Bacteria, independently of their ability to be naturally transformable. Since we have shown previously that capsules are commonly present in nosocomial pathogens, we analysed their co-occurrence with antibiotic resistance genes. Genomes encoding capsules have more antibiotic resistance genes, especially those encoding efflux pumps, and they constitute the majority of the most worrisome nosocomial bacteria. We conclude that bacteria with capsule systems are more genetically diverse and have fast-evolving gene repertoires, which may further contribute to their success in colonizing novel niches such as humans under antibiotic therapy.

Previous works showed that bacteria encoding capsules are better colonizers and are dominant in most environments suggesting a positive role for capsules in the genetic diversification of bacteria. Yet, it has been repeatedly suggested, based almost exclusively studies in few model species, that such bacteria are less diverse and engage in fewer genetic exchanges. Here, we reverse the current paradigm and show that bacteria encoding capsules have larger and more diverse gene repertoires, which change faster by horizontal gene transfer and recombination. Our study alters the traditional view of the capsule as a barrier to gene flow and raises novel questions about the role of capsules in bacterial adaptation.

Experimental Evolution as a High-Throughput Screen for Genetic Adaptations

Experimental evolution is a method in which populations of organisms, often microbes, are founded by one or more ancestors of known genotype and then propagated under controlled conditions to study the evolutionary process. These evolving populations are influenced by all population genetic forces, including selection, mutation, drift, and recombination, and the relative contributions of these forces may be seen as mysterious. Here, I describe why the outcomes of experimental evolution should be viewed with greater certainty because the force of selection typically dominates. Importantly, any mutant rising rapidly to high frequency in large populations must have acquired adaptive traits in the selective environment. Sequencing the genomes of these mutants can identify genes or pathways that contribute to an adaptation. I review the logic and simple mathematics why this evolve-and-resequence approach is a powerful way to find the mutations or mutation combinations that best increase fitness in any new environment.

Interplay Between Capsule Expression and Uracil Metabolism in Streptococcus pneumoniae D39

Pyrimidine nucleotides play an important role in the biosynthesis of activated nucleotide sugars (NDP-sugars). NDP-sugars are the precursors of structural polysaccharides in bacteria, including capsule, which is a major virulence factor of the human pathogen S. pneumoniae. In this work, we identified a spontaneous non-reversible mutant of strain D39 that displayed a non-producing capsule phenotype. Whole-genome sequencing analysis of this mutant revealed several non-synonymous single base modifications, including in genes of the de novo synthesis of pyrimidines and in the -10 box of capsule operon promoter (Pcps). By directed mutagenesis we showed that the point mutation in Pcps was solely responsible for the drastic decrease in capsule expression. We also demonstrated that D39 subjected to uracil deprivation shows increased biomass and decreased Pcps activity and capsule amounts. Importantly, Pcps expression is further decreased by mutating the first gene of the de novo synthesis of pyrimidines, carA. In contrast, the absence of uracil from the culture medium showed no effect on the spontaneous mutant strain. Co-cultivation of the wild-type and the mutant strain indicated a competitive advantage of the spontaneous mutant (non-producing capsule) in medium devoid of uracil. We propose a model in that uracil may act as a signal for the production of different capsule amounts in S. pneumoniae.

The evolutionary advantage of heritable phenotypic heterogeneity

Phenotypic plasticity is an evolutionary driving force in diverse biological processes, including the adaptive immune system, the development of neoplasms, and the persistence of pathogens despite drug pressure. It is essential, therefore, to understand the evolutionary advantage of an allele that confers on cells the ability to express a range of phenotypes. Here, we study the fate of a new mutation that allows the expression of multiple phenotypic states, introduced into a finite population of individuals that can express only a single phenotype. We show that the advantage of such a mutation depends on the degree of phenotypic heritability between generations, called phenotypic memory. We analyze the fixation probability of the phenotypically plastic allele as a function of phenotypic memory, the variance of expressible phenotypes, the rate of environmental changes, and the population size. We find that the fate of a phenotypically plastic allele depends fundamentally on the environmental regime. In constant environments, plastic alleles are advantageous and their fixation probability increases with the degree of phenotypic memory. In periodically fluctuating environments, by contrast, there is an optimum phenotypic memory that maximizes the probability of the plastic allele's fixation. This same optimum memory also maximizes geometric mean fitness, in steady state. We interpret these results in the context of previous studies in an infinite-population framework. We also discuss the implications of our results for the design of therapies that can overcome persistence and, indirectly, drug resistance.

Nonselective Bottlenecks Control the Divergence and Diversification of Phase-Variable Bacterial Populations

Phase variation occurs in many pathogenic and commensal bacteria and is a major generator of genetic variability. A putative advantage of phase variation is to counter reductions in variability imposed by nonselective bottlenecks during transmission. Genomes of Campylobacter jejuni, a widespread food-borne pathogen, contain multiple phase-variable loci whose rapid, stochastic variation is generated by hypermutable simple sequence repeat tracts. These loci can occupy a vast number of combinatorial expression states (phasotypes) enabling populations to rapidly access phenotypic diversity. The imposition of nonselective bottlenecks can perturb the relative frequencies of phasotypes, changing both within-population diversity and divergence from the initial population. Using both in vitro testing of C. jejuni populations and a simple stochastic simulation of phasotype change, we observed that single-cell bottlenecks produce output populations of low diversity but with bimodal patterns of either high or low divergence. Conversely, large bottlenecks allow divergence only by accumulation of diversity, while interpolation between these extremes is observed in intermediary bottlenecks. These patterns are sensitive to the genetic diversity of initial populations but stable over a range of mutation rates and number of loci. The qualitative similarities of experimental and in silico modeling indicate that the observed patterns are robust and applicable to other systems where localized hypermutation is a defining feature. We conclude that while phase variation will maintain bacterial population diversity in the face of intermediate bottlenecks, narrow transmission-associated bottlenecks could produce host-to-host variation in bacterial phenotypes and hence stochastic variation in colonization and disease outcomes.

Transmission and within-host spread of pathogenic organisms are associated with selective and nonselective bottlenecks that significantly reduced population diversity. In several bacterial pathogens, hypermutable mechanisms have evolved that mediate high-frequency reversible switching of specific phenotypes, such as surface structures, and hence counteract bottleneck-associated reductions in population diversity. Here, we investigated how combinations of hypermutable simple sequence repeats interact with nonselective bottlenecks by using a stochastic computer model and experimental data for Campylobacter jejuni, a food-borne pathogen. We find that bottleneck size qualitatively alters the output populations, with large bottlenecks maintaining population diversity while small bottlenecks produce dramatic shifts in the prevalence of particular variants. We conclude that narrow bottlenecks are capable of producing host-to-host variation in repeat-controlled bacterial phenotypes, leading to a potential for stochastic person-to-person variations in disease outcome for C. jejuni and other organisms with similar hypermutable mechanisms.

When sensing is gambling: An experimental system reveals how plasticity can generate tunable bet-hedging strategies

Genotypes can persist in unpredictable environments by "hedging their bets" and producing diverse phenotypes. Theoretical studies have shown that the phenotypic variability needed for a bet-hedging strategy can be generated by factors either inside or outside an organism. However, sensing the environment and bet hedging are frequently treated as distinct evolutionary strategies. Furthermore, nearly all empirical studies of the molecular underpinnings of bet-hedging strategies to date have focused on internal sources of variability. We took a synthetic approach and constructed an experimental system where a phenotypic trade-off is mediated by actively sensing a cue present in the environment. We show that active sensing can generate a diversified bet-hedging strategy. Mutations affecting the norm of reaction to the cue alter the diversification strategy, indicating that bet hedging by active sensing is evolvable. Our results indicate that a broader class of biological systems should be considered as potential examples of bet-hedging strategies, and that research into the structure of environmental variability is needed to distinguish bet-hedging strategies from adaptive plasticity.

The Shift of an Intestinal "Microbiome" to a "Pathobiome" Governs the Course and Outcome of Sepsis Following Surgical Injury

Sepsis following surgical injury remains a growing and worrisome problem following both emergent and elective surgery. Although early resuscitation efforts and prompt antibiotic therapy have improved outcomes in the first 24 to 48  h, late onset sepsis is now the most common cause of death in modern intensive care units. This time shift may be, in part, a result of prolonged exposure of the host to the stressors of critical illness which, over time, erode the health promoting intestinal microbiota and allow for virulent pathogens to predominate. Colonizing pathogens can then subvert the immune system and contribute to the deterioration of the host response. Here, we posit that novel approaches integrating the molecular, ecological, and evolutionary dynamics of the evolving gut microbiome/pathobiome during critical illness are needed to understand and prevent the late onset sepsis that develops following prolonged critical illness.

Evolutionary convergence in experimental Pseudomonas populations

Model microbial systems provide opportunity to understand the genetic bases of ecological traits, their evolution, regulation and fitness contributions. Experimental populations of Pseudomonas fluorescens rapidly diverge in spatially structured microcosms producing a range of surface-colonising forms. Despite divergent molecular routes, wrinkly spreader (WS) niche specialist types overproduce a cellulosic polymer allowing mat formation at the air-liquid interface and access to oxygen. Given the range of ways by which cells can form mats, such phenotypic parallelism is unexpected. We deleted the cellulose-encoding genes from the ancestral genotype and asked whether this mutant could converge on an alternate phenotypic solution. Two new traits were discovered. The first involved an exopolysaccharide encoded by pgaABCD that functions as cell-cell glue similar to cellulose. The second involved an activator of an amidase (nlpD) that when defective causes cell chaining. Both types form mats, but were less fit in competition with cellulose-based WS types. Surprisingly, diguanylate cyclases linked to cellulose overexpression underpinned evolution of poly-beta-1,6-N-acetyl-d-glucosamine (PGA)-based mats. This prompted genetic analyses of the relationships between the diguanylate cyclases WspR, AwsR and MwsR, and both cellulose and PGA. Our results suggest that c-di-GMP regulatory networks may have been shaped by evolution to accommodate loss and gain of exopolysaccharide modules facilitating adaptation to new environments.

Stochasticity in the Genotype-Phenotype Map: Implications for the Robustness and Persistence of Bet-Hedging

Nongenetic variation in phenotypes, or bet-hedging, has been observed as a driver of drug resistance in both bacterial infections and cancers. Here, we study how bet-hedging emerges in genotype-phenotype (GP) mapping through a simple interaction model: a molecular switch. We use simple chemical reaction networks to implement stochastic switches that map gene products to phenotypes, and investigate the impact of structurally distinct mappings on the evolution of phenotypic heterogeneity. Bet-hedging naturally emerges within this model, and is robust to evolutionary loss through mutations to both the expression of individual genes, and to the network itself. This robustness explains an apparent paradox of bet-hedging-why does it persist in environments where natural selection necessarily acts to remove it? The structure of the underlying molecular mechanism, itself subject to selection, can slow the evolutionary loss of bet-hedging to ensure a survival mechanism against environmental catastrophes even when they are rare. Critically, these properties, taken together, have profound implications for the use of treatment-holidays to combat bet-hedging-driven resistant disease, as the efficacy of breaks from treatment will ultimately be determined by the structure of the GP mapping.

Multiple roles for a novel RND-type efflux system in Acinetobacter baumannii AB5075

Colony opacity phase variation in Acinetobacter baumannii strain AB5075 is regulated by a reversible high‐frequency switch. Transposon mutagenesis was used to generate mutations that decreased the opaque to translucent switch and a gene encoding a predicted periplasmic membrane fusion component of a resistance–nodulation–cell division (RND)‐type efflux system was isolated. This gene was designated arpA and immediately downstream was a gene designated arpB that encodes a predicted membrane transporter of RND‐type systems. A nonpolar, in‐frame deletion in arpA resulted in a 70‐fold decrease in the opaque to translucent switch. An arpB::Tc mutant exhibited a 769‐fold decrease in the opaque to translucent switch. However, the translucent to opaque switch was largely unchanged in both the arpA and arpB mutants. The arpA and arpB mutants also exhibited increased surface motility in the opaque form and the arpB mutant exhibited increased susceptibility to aminoglycosides. The arpA and arpB mutants were both attenuated in a Galleria mellonella model of virulence. A divergently transcribed TetR‐type regulator ArpR was capable of repressing the arpAB operon when this TetR regulator was overexpressed. The arpR gene was also involved in regulating the opaque to translucent switch as an in‐frame arpR mutation decreased this switch by 1,916‐fold.

Effect of Shear Stress on Pseudomonas aeruginosa Isolated from the Cystic Fibrosis Lung

Chronic colonization of the lungs by Pseudomonas aeruginosa is one of the major causes of morbidity and mortality in cystic fibrosis (CF) patients. To gain insights into the characteristic biofilm phenotype of P. aeruginosa in the CF lungs, mimicking the CF lung environment is critical. We previously showed that growth of the non-CF-adapted P. aeruginosa PAO1 strain in a rotating wall vessel, a device that simulates the low fluid shear (LS) conditions present in the CF lung, leads to the formation of in-suspension, self-aggregating biofilms. In the present study, we determined the phenotypic and transcriptomic changes associated with the growth of a highly adapted, transmissible P. aeruginosa CF strain in artificial sputum medium under LS conditions. Robust self-aggregating biofilms were observed only under LS conditions. Growth under LS conditions resulted in the upregulation of genes involved in stress response, alginate biosynthesis, denitrification, glycine betaine biosynthesis, glycerol metabolism, and cell shape maintenance, while genes involved in phenazine biosynthesis, type VI secretion, and multidrug efflux were downregulated. In addition, a number of small RNAs appeared to be involved in the response to shear stress. Finally, quorum sensing was found to be slightly but significantly affected by shear stress, resulting in higher production of autoinducer molecules during growth under high fluid shear (HS) conditions. In summary, our study revealed a way to modulate the behavior of a highly adapted P. aeruginosa CF strain by means of introducing shear stress, driving it from a biofilm lifestyle to a more planktonic lifestyle.

Biofilm formation by Pseudomonas aeruginosa is one of the hallmarks of chronic cystic fibrosis (CF) lung infections. The biofilm matrix protects this bacterium from antibiotics as well as from the immune system. Hence, the prevention or reversion of biofilm formation is believed to have a great impact on treatment of chronic P. aeruginosa CF lung infections. In the present study, we showed that it is possible to modulate the behavior of a highly adapted transmissible P. aeruginosa CF isolate at both the transcriptomic and phenotypic levels by introducing shear stress in a CF-like environment, driving it from a biofilm to a planktonic lifestyle. Consequently, the results obtained in this study are of great importance with regard to therapeutic applications that introduce shear stress in the lungs of CF patients.

Phenotypic Evolution With and Beyond Genome Evolution

DNA does not make phenotypes on its own. In this volume entitled "Genes and Phenotypic Evolution," the present review draws the attention on the process of phenotype construction-including development of multicellular organisms-and the multiple interactions and feedbacks between DNA, organism, and environment at various levels and timescales in the evolutionary process. First, during the construction of an individual's phenotype, DNA is recruited as a template for building blocks within the cellular context and may in addition be involved in dynamical feedback loops that depend on the environmental and organismal context. Second, in the production of phenotypic variation among individuals, stochastic, environmental, genetic, and parental sources of variation act jointly. While in controlled laboratory settings, various genetic and environmental factors can be tested one at a time or in various combinations, they cannot be separated in natural populations because the environment is not controlled and the genotype can rarely be replicated. Third, along generations, genotype and environment each have specific properties concerning the origin of their variation, the hereditary transmission of this variation, and the evolutionary feedbacks. Natural selection acts as a feedback from phenotype and environment to genotype. This review integrates recent results and concrete examples that illustrate these three points. Although some themes are shared with recent calls and claims to a new conceptual framework in evolutionary biology, the viewpoint presented here only means to add flesh to the standard evolutionary synthesis.

Lineage Tracking for Probing Heritable Phenotypes at Single-Cell Resolution

Determining the phenotype and genotype of single cells is central to understand microbial evolution. DNA sequencing technologies allow the detection of mutants at high resolution, but similar approaches for phenotypic analyses are still lacking. We show that a drop-based millifluidic system enables the detection of heritable phenotypic changes in evolving bacterial populations. At time intervals, cells were sampled and individually compartmentalized in 100 nL drops. Growth through 15 generations was monitored using a fluorescent protein reporter. Amplification of heritable changes–via growth–over multiple generations yields phenotypically distinct clusters reflecting variation relevant for evolution. To demonstrate the utility of this approach, we follow the evolution of Escherichia coli populations during 30 days of starvation. Phenotypic diversity was observed to rapidly increase upon starvation with the emergence of heritable phenotypes. Mutations corresponding to each phenotypic class were identified by DNA sequencing. This scalable lineage-tracking technology opens the door to large-scale phenotyping methods with special utility for microbiology and microbial population biology.

Transcription factor levels enable metabolic diversification of single cells of environmental bacteria

Transcriptional noise is a necessary consequence of the molecular events that drive gene expression in prokaryotes. However, some environmental microorganisms that inhabit polluted sites, for example, the m-xylene degrading soil bacterium Pseudomonas putida mt-2 seem to have co-opted evolutionarily such a noise for deploying a metabolic diversification strategy that allows a cautious exploration of new chemical landscapes. We have examined this phenomenon under the light of deterministic and stochastic models for activation of the main promoter of the master m-xylene responsive promoter of the system (Pu) by its cognate transcriptional factor (XylR). These analyses consider the role of co-factors for Pu activation and determinants of xylR mRNA translation. The model traces the onset and eventual disappearance of the bimodal distribution of Pu activity along time to the growth-phase dependent abundance of XylR itself, that is, very low in exponentially growing cells and high in stationary. This tenet was validated by examining the behaviour of a Pu-GFP fusion in a P. putida strain in which xylR expression was engineered under the control of an IPTG-inducible system. This work shows how a relatively simple regulatory scenario (for example, growth-phase dependent expression of a limiting transcription factor) originates a regime of phenotypic diversity likely to be advantageous in competitive environmental settings.

The Identification of Genes Important in Pseudomonas syringae pv. phaseolicola Plant Colonisation Using In Vitro Screening of Transposon Libraries

The bacterial plant pathogen Pseudomonas syringae pv. phaseolicola (Pph) colonises the surface of common bean plants before moving into the interior of plant tissue, via wounds and stomata. In the intercellular spaces the pathogen proliferates in the apoplastic fluid and forms microcolonies (biofilms) around plant cells. If the pathogen can suppress the plant’s natural resistance response, it will cause halo blight disease. The process of resistance suppression is fairly well understood, but the mechanisms used by the pathogen in colonisation are less clear. We hypothesised that we could apply in vitro genetic screens to look for changes in motility, colony formation, and adhesion, which are proxies for infection, microcolony formation and cell adhesion. We made transposon (Tn) mutant libraries of Pph strains 1448A and 1302A and found 106/1920 mutants exhibited alterations in colony morphology, motility and biofilm formation. Identification of the insertion point of the Tn identified within the genome highlighted, as expected, a number of altered motility mutants bearing mutations in genes encoding various parts of the flagellum. Genes involved in nutrient biosynthesis, membrane associated proteins, and a number of conserved hypothetical protein (CHP) genes were also identified. A mutation of one CHP gene caused a positive increase in in planta bacterial growth. This rapid and inexpensive screening method allows the discovery of genes important for in vitro traits that can be correlated to roles in the plant interaction.

Individuality, phenotypic differentiation, dormancy and 'persistence' in culturable bacterial systems: commonalities shared by environmental, laboratory, and clinical microbiology

For bacteria, replication mainly involves growth by binary fission. However, in a very great many natural environments there are examples of phenotypically dormant, non-growing cells that do not replicate immediately and that are phenotypically 'nonculturable' on media that normally admit their growth. They thereby evade detection by conventional culture-based methods. Such dormant cells may also be observed in laboratory cultures and in clinical microbiology. They are usually more tolerant to stresses such as antibiotics, and in clinical microbiology they are typically referred to as 'persisters'. Bacterial cultures necessarily share a great deal of relatedness, and inclusive fitness theory implies that there are conceptual evolutionary advantages in trading a variation in growth rate against its mean, equivalent to hedging one's bets. There is much evidence that bacteria exploit this strategy widely. We here bring together data that show the commonality of these phenomena across environmental, laboratory and clinical microbiology. Considerable evidence, using methods similar to those common in environmental microbiology, now suggests that many supposedly non-communicable, chronic and inflammatory diseases are exacerbated (if not indeed largely caused) by the presence of dormant or persistent bacteria (the ability of whose components to cause inflammation is well known). This dormancy (and resuscitation therefrom) often reflects the extent of the availability of free iron. Together, these phenomena can provide a ready explanation for the continuing inflammation common to such chronic diseases and its correlation with iron dysregulation. This implies that measures designed to assess and to inhibit or remove such organisms (or their access to iron) might be of much therapeutic benefit.

Individuality, phenotypic differentiation, dormancy and 'persistence' in culturable bacterial systems: commonalities shared by environmental, laboratory, and clinical microbiology

For bacteria, replication mainly involves growth by binary fission. However, in a very great many natural environments there are examples of phenotypically dormant, non-growing cells that do not replicate immediately and that are phenotypically 'nonculturable' on media that normally admit their growth. They thereby evade detection by conventional culture-based methods. Such dormant cells may also be observed in laboratory cultures and in clinical microbiology. They are usually more tolerant to stresses such as antibiotics, and in clinical microbiology they are typically referred to as 'persisters'. Bacterial cultures necessarily share a great deal of relatedness, and inclusive fitness theory implies that there are conceptual evolutionary advantages in trading a variation in growth rate against its mean, equivalent to hedging one's bets. There is much evidence that bacteria exploit this strategy widely. We here bring together data that show the commonality of these phenomena across environmental, laboratory and clinical microbiology. Considerable evidence, using methods similar to those common in environmental microbiology, now suggests that many supposedly non-communicable, chronic and inflammatory diseases are exacerbated (if not indeed largely caused) by the presence of dormant or persistent bacteria (the ability of whose components to cause inflammation is well known). This dormancy (and resuscitation therefrom) often reflects the extent of the availability of free iron. Together, these phenomena can provide a ready explanation for the continuing inflammation common to such chronic diseases and its correlation with iron dysregulation. This implies that measures designed to assess and to inhibit or remove such organisms (or their access to iron) might be of much therapeutic benefit.