Divergence involving global regulatory gene mutations in an Escherichia coli population evolving under phosphate limitation

Serial adaptive laboratory evolution enhances mixed carbon metabolic capacity of Escherichia coli

Microbes have inherent capacities for utilizing various carbon sources, however they often exhibit sub-par fitness due to low metabolic efficiency. To test whether a bacterial strain can optimally utilize multiple carbon sources, Escherichia coli was serially evolved in L-lactate and glycerol. This yielded two end-point strains that evolved first in L-lactate then in glycerol, and vice versa. The end-point strains displayed a universal growth advantage on single and a mixture of adaptive carbon sources, enabled by a concerted action of carbon source-specialists and generalist mutants. The combination of just four variants of glpK, ppsA, ydcI, and rph-pyrE, accounted for more than 80% of end-point strain fitness. In addition, machine learning analysis revealed a coordinated activity of transcriptional regulators imparting condition-specific regulation of gene expression. The effectiveness of the serial adaptive laboratory evolution (ALE) scheme in bioproduction applications was assessed under single and mixed-carbon culture conditions, in which serial ALE strain exhibited superior productivity of acetoin compared to ancestral strains. Together, systems-level analysis elucidated the molecular basis of serial evolution, which hold potential utility in bioproduction applications.

Adaptation of the Mycobacterium tuberculosis transcriptome to biofilm growth

Mycobacterium tuberculosis (M. tb), the causative agent of tuberculosis (TB), is a leading global cause of death from infectious disease. Biofilms are increasingly recognized as a relevant growth form during M. tb infection and may impede treatment by enabling bacterial drug and immune tolerance. M. tb has a complicated regulatory network that has been well-characterized for many relevant disease states, including dormancy and hypoxia. However, despite its importance, our knowledge of the genes and pathways involved in biofilm formation is limited. Here we characterize the biofilm transcriptomes of fully virulent clinical isolates and find that the regulatory systems underlying biofilm growth vary widely between strains and are also distinct from regulatory programs associated with other environmental cues. We used experimental evolution to investigate changes to the transcriptome during adaptation to biofilm growth and found that the application of a uniform selection pressure resulted in loss of strain-to-strain variation in gene expression, resulting in a more uniform biofilm transcriptome. The adaptive trajectories of transcriptomes were shaped by the genetic background of the M. tb population leading to convergence on a sub-lineage specific transcriptome. We identified widespread upregulation of non-coding RNA (ncRNA) as a common feature of the biofilm transcriptome and hypothesize that ncRNA function in genome-wide modulation of gene expression, thereby facilitating rapid regulatory responses to new environments. These results reveal a new facet of the M. tb regulatory system and provide valuable insight into how M. tb adapts to new environments.

RpoS and the bacterial general stress response

SUMMARYThe general stress response (GSR) is a widespread strategy developed by bacteria to adapt and respond to their changing environments. The GSR is induced by one or multiple simultaneous stresses, as well as during entry into stationary phase and leads to a global response that protects cells against multiple stresses. The alternative sigma factor RpoS is the central GSR regulator in E. coli and conserved in most γ-proteobacteria. In E. coli, RpoS is induced under conditions of nutrient deprivation and other stresses, primarily via the activation of RpoS translation and inhibition of RpoS proteolysis. This review includes recent advances in our understanding of how stresses lead to RpoS induction and a summary of the recent studies attempting to define RpoS-dependent genes and pathways.

RpoS role in antibiotic resistance, tolerance and persistence in E. coli natural isolates

BACKGROUND

The intrinsic concentration of RpoS, the second most abundant sigma factor, varies widely across the E. coli species. Bacterial isolates that express high levels of RpoS display high resistance to environmental stresses, such as temperature, pH and osmolarity shifts, but are less nutritional competent, making them less capable of utilising alternative nutrient sources. The role of RpoS in antibiotic resistance and persistence in standard laboratory domesticated strains has been examined in several studies, most demonstrating a positive role for RpoS.

RESULTS

Using disk diffusion assays we examined bacterial resistance to 15 different antibiotics, including β -lactams (penicillins, monobactams, carbapenems and cephalosporins), aminoglycosides, quinolones and anti-folates, in a representative collection of 328 E. coli natural isolates displaying a continuum of different levels of RpoS. There was an overall trend that isolates with higher levels of RpoS were slightly more resistant to these antibiotics. In addition, the effect of RpoS on bacterial tolerance and persistence to 3 different antibiotics - ampicillin, ciprofloxacin and kanamycin was evaluated through time-kill curves. Again, there was a small beneficial effect of RpoS on tolerance and persistence to these antibiotics, but this difference was not statistically significant. Finally, a K-12 strain expressing high levels of RpoS was compared with its isogenic RpoS-null counterpart, and no significant effect of RpoS was found.

CONCLUSION

Based on a representative collection of the species E. coli, RpoS was found to have a very small impact on antibiotic resistance, tolerance, or persistence.

Adaptation of the Mycobacterium tuberculosis transcriptome to biofilm growth

Mycobacterium tuberculosis ( M. tb ), the causative agent of tuberculosis (TB), is a leading global cause of death from infectious disease. Biofilms are increasingly recognized as a relevant growth form during M. tb infection and may impede treatment by enabling bacterial drug and immune tolerance. M. tb has a complicated regulatory network that has been well-characterized for many relevant disease states, including dormancy and hypoxia. However, despite its importance, our knowledge of the genes and pathways involved in biofilm formation is limited. Here we characterize the biofilm transcriptomes of fully virulent clinical isolates and find that the regulatory systems underlying biofilm growth vary widely between strains and are also distinct from regulatory programs associated with other environmental cues. We used experimental evolution to investigate changes to the transcriptome during adaptation to biofilm growth and found that the application of a uniform selection pressure resulted in loss of strain-to-strain variation in gene expression, resulting in a more uniform biofilm transcriptome. The adaptive trajectories of transcriptomes were shaped by the genetic background of the M. tb population leading to convergence on a sub-lineage specific transcriptome. We identified widespread upregulation of non-coding RNA (ncRNA) as a common feature of the biofilm transcriptome and hypothesize that ncRNA function in genome-wide modulation of gene expression, thereby facilitating rapid regulatory responses to new environments. These results reveal a new facet of the M. tb regulatory system and provide valuable insight into how M. tb adapts to new environments.

Importance

Understanding mechanisms of resistance and tolerance in Mycobacterium tuberculosis ( M. tb ) can help us develop new treatments that capitalize on M. tb 's vulnerabilities. Here we used transcriptomics to study both the regulation of biofilm formation in clinical isolates as well as how those regulatory systems adapt to new environments. We find that closely related clinical populations have diverse strategies for growth under biofilm conditions, and that genetic background plays a large role in determining the trajectory of evolution. These results have implications for future treatment strategies that may be informed by our knowledge of the evolutionary constraints on strain(s) from an individual infection. This work provides new information about the mechanisms of biofilm formation in M. tb and outlines a framework for population level approaches for studying bacterial adaptation.

Functional expansion of the natural inorganic phosphorus starvation response system in Escherichia coli

Phosphorus, an indispensable nutrient, plays an essential role in cell composition, metabolism, and signal transduction. When inorganic phosphorus (Pi) is scarce, the Pi starvation response in E. coli is activated to increase phosphorus acquisition and drive the cells into a non-growing state to reduce phosphorus consumption. In the six decades of research history, the initiation, output, and shutdown processes of the Pi starvation response have been extensively studied. Simultaneously, Pi starvation has been used in biosensor development, recombinant protein production, and natural product biosynthesis. In this review, we focus on the output process and the applications of the Pi starvation response that have not been summarized before. Meanwhile, based on the current status of mechanistic studies and applications, we propose practical strategies to develop the natural Pi starvation response into a multifunctional and standardized regulatory system in four aspects, including response threshold, temporal expression, intensity range, and bifunctional regulation, which will contribute to its broader application in more fields such as industrial production, medical analysis, and environmental protection.

A Broad Continuum of E. coli Traits in Nature Associated with the Trade-off Between Self-preservation and Nutritional Competence

A trade-off between reproduction and survival is a characteristic of many organisms. In bacteria, growth is constrained when cellular resources are channelled towards environmental stress protection. At the core of this trade-off in Escherichia coli is RpoS, a sigma factor that diverts transcriptional resources towards general stress resistance. The constancy of RpoS levels in natural isolates is unknown. A uniform RpoS content in E. coli would impart a narrow range of resistance properties to the species, whereas a diverse set of RpoS levels in nature should result in a diverse range of stress susceptibilities. We explore the diversity of trade-off settings and phenotypes by measuring the level of RpoS protein in strains of E. coli cohabiting in a natural environment. Strains from a stream polluted with domestic waste were investigated in monthly samples. Analyses included E. coli phylogroup classification, RpoS protein level, RpoS-dependent stress phenotypes and the sequencing of rpoS mutations. The most striking finding was the continuum of RpoS levels, with a 100-fold range of RpoS amounts consistently found in individuals in the stream. Approximately 1.8% of the sampled strains carried null or non-synonymous mutations in rpoS. The natural isolates also exhibited a broad (>100-fold) range of stress resistance responses. Our results are consistent with the view that a multiplicity of survival-multiplication trade-off settings is a feature of the species E. coli. The phenotypic diversity resulting from the trade-off permits bet-hedging and the adaptation of E. coli strains to a very broad range of environments.

Evolutionary dynamics and structural consequences of de novo beneficial mutations and mutant lineages arising in a constant environment

BACKGROUND

Microbial evolution experiments can be used to study the tempo and dynamics of evolutionary change in asexual populations, founded from single clones and growing into large populations with multiple clonal lineages. High-throughput sequencing can be used to catalog de novo mutations as potential targets of selection, determine in which lineages they arise, and track the fates of those lineages. Here, we describe a long-term experimental evolution study to identify targets of selection and to determine when, where, and how often those targets are hit.

RESULTS

We experimentally evolved replicate Escherichia coli populations that originated from a mutator/nonsense suppressor ancestor under glucose limitation for between 300 and 500 generations. Whole-genome, whole-population sequencing enabled us to catalog 3346 de novo mutations that reached > 1% frequency. We sequenced the genomes of 96 clones from each population when allelic diversity was greatest in order to establish whether mutations were in the same or different lineages and to depict lineage dynamics. Operon-specific mutations that enhance glucose uptake were the first to rise to high frequency, followed by global regulatory mutations. Mutations related to energy conservation, membrane biogenesis, and mitigating the impact of nonsense mutations, both ancestral and derived, arose later. New alleles were confined to relatively few loci, with many instances of identical mutations arising independently in multiple lineages, among and within replicate populations. However, most never exceeded 10% in frequency and were at a lower frequency at the end of the experiment than at their maxima, indicating clonal interference. Many alleles mapped to key structures within the proteins that they mutated, providing insight into their functional consequences.

CONCLUSIONS

Overall, we find that when mutational input is increased by an ancestral defect in DNA repair, the spectrum of high-frequency beneficial mutations in a simple, constant resource-limited environment is narrow, resulting in extreme parallelism where many adaptive mutations arise but few ever go to fixation.

Diversity in E. coli (p)ppGpp Levels and Its Consequences

(p)ppGpp is at the core of global bacterial regulation as it controls growth, the most important aspect of life. It would therefore be expected that at least across a species the intrinsic (basal) levels of (p)ppGpp would be reasonably constant. On the other hand, the historical contingency driven by the selective pressures on bacterial populations vary widely resulting in broad genetic polymorphism. Given that (p)ppGpp controls the expression of many genes including those involved in the bacterial response to environmental challenges, it is not surprising that the intrinsic levels of (p)ppGpp would also vary considerably. In fact, null mutations or less severe genetic polymorphisms in genes associated with (p)ppGpp synthesis and hydrolysis are common. Such variation can be observed in laboratory strains, in natural isolates as well as in evolution experiments. High (p)ppGpp levels result in low growth rate and high tolerance to environmental stresses. Other aspects such as virulence and antimicrobial resistance are also influenced by the intrinsic levels of (p)ppGpp. A case in point is the production of Shiga toxin by certain E. coli strains which is inversely correlated to (p)ppGpp basal level. Conversely, (p)ppGpp concentration is positively correlated to increased tolerance to different antibiotics such as β-lactams, vancomycin, and others. Here we review the variations in intrinsic (p)ppGpp levels and its consequences across the E. coli species.

Adaptive Evolution of Industrial Brewer’s Yeast Strains towards a Snowflake Phenotype

Flocculation or cell aggregation is a well-appreciated characteristic of industrial brewer’s strains, since it allows removal of the cells from the beer in a cost-efficient and environmentally-friendly manner. However, many industrial strains are non-flocculent and genetic interference to increase the flocculation characteristics are not appreciated by the consumers. We applied adaptive laboratory evolution (ALE) to three non-flocculent, industrial Saccharomyces cerevisiae brewer’s strains using small continuous bioreactors (ministats) to obtain an aggregative phenotype, i.e., the “snowflake” phenotype. These aggregates could increase yeast sedimentation considerably. We evaluated the performance of these evolved strains and their produced flavor during lab scale beer fermentations. The small aggregates did not result in a premature sedimentation during the fermentation and did not result in major flavor changes of the produced beer. These results show that ALE could be used to increase the sedimentation behavior of non-flocculent brewer’s strains.

The emergence of adaptive laboratory evolution as an efficient tool for biological discovery and industrial biotechnology

Harnessing the process of natural selection to obtain and understand new microbial phenotypes has become increasingly possible due to advances in culturing techniques, DNA sequencing, bioinformatics, and genetic engineering. Accordingly, Adaptive Laboratory Evolution (ALE) experiments represent a powerful approach both to investigate the evolutionary forces influencing strain phenotypes, performance, and stability, and to acquire production strains that contain beneficial mutations. In this review, we summarize and categorize the applications of ALE to various aspects of microbial physiology pertinent to industrial bioproduction by collecting case studies that highlight the multitude of ways in which evolution can facilitate the strain construction process. Further, we discuss principles that inform experimental design, complementary approaches such as computational modeling that help maximize utility, and the future of ALE as an efficient strain design and build tool driven by growing adoption and improvements in automation.

GREACE-assisted adaptive laboratory evolution in endpoint fermentation broth enhances lysine production by Escherichia coli

BACKGROUND

Late-stage fermentation broth contains high concentrations of target chemicals. Additionally, it contains various cellular metabolites which have leaked from lysed cells, which would exert multifactorial stress to industrial hyperproducers and perturb both cellular metabolism and product formation. Although adaptive laboratory evolution (ALE) has been wildly used to improve stress tolerance of microbial cell factories, single-factor stress condition (i.e. target product or sodium chloride at a high concentration) is currently provided. In order to enhance bacterial stress tolerance to actual industrial production conditions, ALE in late-stage fermentation broth is desired. Genome replication engineering assisted continuous evolution (GREACE) employs mutants of the proofreading element of DNA polymerase complex (DnaQ) to facilitate mutagenesis. Application of GREACE coupled-with selection under stress conditions is expected to accelerate the ALE process.

RESULTS

In this study, GREACE was first modified by expressing a DnaQ mutant KR5-2 using an arabinose inducible promoter on a temperature-sensitive plasmid, which resulted in timed mutagenesis introduction. Using this method, tolerance of a lysine hyperproducer E. coli MU-1 was improved by enriching mutants in a lysine endpoint fermentation broth. Afterwards, the KR5-2 expressing plasmid was cured to stabilize acquired genotypes. By subsequent fermentation evaluation, a mutant RS3 with significantly improved lysine production capacity was selected. The final titer, yield and total amount of lysine produced by RS3 in a 5-L batch fermentation reached 155.0 ± 1.4 g/L, 0.59 ± 0.02 g lysine/g glucose, and 605.6 ± 23.5 g, with improvements of 14.8%, 9.3%, and 16.7%, respectively. Further metabolomics and genomics analyses, coupled with molecular biology studies revealed that mutations SpeB^A302V, AtpB^S165N and SecY^M145V mainly contributed both to improved cell integrity under stress conditions and enhanced metabolic flux into lysine synthesis.

CONCLUSIONS

Our present study indicates that improving a lysine hyperproducer by GREACE-assisted ALE in its stressful living environment is efficient and effective. Accordingly, this is a promising method for improving other valuable chemical hyperproducers.

Quantifying dynamic mechanisms of auto-regulation in Escherichia coli with synthetic promoter in response to varying external phosphate levels

Escherichia coli have developed one of the most efficient regulatory response mechanisms to phosphate starvation. The machinery involves a cascade with a two-component system (TCS) that relays the external signal to the genetic circuit, resulting in a feedback response. Achieving a quantitative understanding of this system has implications in synthetic biology and biotechnology, for example, in applications for wastewater treatment. To this aim, we present a computational model and experimental results with a detailed description of the TCS, consisting of PhoR and PhoB, together with the mechanisms of gene expression. The model is parameterised within the feasible range, and fitted to the dynamic response of our experimental data on PhoB as well as PhoA, the product of this network that is used in alkaline phosphatase production. Deterministic and stochastic simulations with our model predict the regulation dynamics in higher external phosphate concentrations while reproducing the experimental observations. In a cycle of simulations and experimental verification, our model predicts and explores phenotypes with various synthetic promoter designs that can optimise the inorganic phosphate intake in E. coli. Sensitivity analysis demonstrates that the Pho-controlled genes have a significant influence over the phosphate response. Together with experimental findings, our model should thus provide insights for the investigations on engineering new sensors and regulators for living technologies.

Escherichia coli mutation rates and spectra with combinations of environmental limitations

Micro-organisms often face multiple stresses in natural habitats. Individual stresses are well known to influence mutation rates and the spectra of mutational types, but the extent to which multiple stresses affect the genetic variation in populations is unknown. Here we investigate pair-wise combinations of nutritional stresses in Escherichia coli to determine their effect on mutation rates and mutational types. Environmental interactions modified both the rate and spectrum of mutations in double-limited environments, but the effects were not additive or synergistic relative to single stresses. Generally, bacteria in the mixed environments behaved as if one of the two single-stress stimuli was more dominant and the genetic variation seen with every dual limitation was intermediate between known patterns with individual stresses. The composition of mutational types with double stresses was also intermediate between individual stress patterns. At least with mutations, the single stressor results available are reasonable indicators of stress-induced genetic variation in multifaceted natural habitats. With the influence of 11 conditions available on mutational patterns, we can now also see the clustering of mutational types as a function of these environments.

Predicting the evolution of Escherichia coli by a data-driven approach

A tantalizing question in evolutionary biology is whether evolution can be predicted from past experiences. To address this question, we created a coherent compendium of more than 15,000 mutation events for the bacterium Escherichia coli under 178 distinct environmental settings. Compendium analysis provides a comprehensive view of the explored environments, mutation hotspots and mutation co-occurrence. While the mutations shared across all replicates decrease with the number of replicates, our results argue that the pairwise overlapping ratio remains the same, regardless of the number of replicates. An ensemble of predictors trained on the mutation compendium and tested in forward validation over 35 evolution replicates achieves a 49.2 ± 5.8% (mean ± std) precision and 34.5 ± 5.7% recall in predicting mutation targets. This work demonstrates how integrated datasets can be harnessed to create predictive models of evolution at a gene level and elucidate the effect of evolutionary processes in well-defined environments.

How reproducible evolutionary processes are remains an important question in evolutionary biology. Here, the authors compile a compendium of more than 15,000 mutation events for Escherichia coli under 178 distinct environmental settings, and develop an ensemble of predictors to predict evolution at a gene level.

Experimental Design, Population Dynamics, and Diversity in Microbial Experimental Evolution

In experimental evolution, laboratory-controlled conditions select for the adaptation of species, which can be monitored in real time. Despite the current popularity of such experiments, nature's most pervasive biological force was long believed to be observable only on time scales that transcend a researcher's life-span, and studying evolution by natural selection was therefore carried out solely by comparative means. Eventually, microorganisms' propensity for fast evolutionary changes proved us wrong, displaying strong evolutionary adaptations over a limited time, nowadays massively exploited in laboratory evolution experiments. Here, we formulate a guide to experimental evolution with microorganisms, explaining experimental design and discussing evolutionary dynamics and outcomes and how it is used to assess ecoevolutionary theories, improve industrially important traits, and untangle complex phenotypes. Specifically, we give a comprehensive overview of the setups used in experimental evolution. Additionally, we address population dynamics and genetic or phenotypic diversity during evolution experiments and expand upon contributing factors, such as epistasis and the consequences of (a)sexual reproduction. Dynamics and outcomes of evolution are most profoundly affected by the spatiotemporal nature of the selective environment, where changing environments might lead to generalists and structured environments could foster diversity, aided by, for example, clonal interference and negative frequency-dependent selection. We conclude with future perspectives, with an emphasis on possibilities offered by fast-paced technological progress. This work is meant to serve as an introduction to those new to the field of experimental evolution, as a guide to the budding experimentalist, and as a reference work to the seasoned expert.

Different adaptive strategies in E. coli populations evolving under macronutrient limitation and metal ion limitation

BACKGROUND

Adaptive responses to nutrient limitation involve mutations that increase the efficiency of usage or uptake of the limiting nutrient. However, starvation of different nutrients has contrasting effects on physiology, resulting in different evolutionary responses. Most studies performed to understand these evolutionary responses have focused only on macronutrient limitation. Hence our understanding of adaptation under limitation of other forms of nutrients is limited. In this study, we compared the evolutionary response in populations evolving under growth-limiting conditions for a macronutrient and a major cation.

RESULTS

We evolved eight populations of E. coli in nutrient-limited chemostats for 400 generations to identify the genetic basis of the mechanisms involved in efficient usage of two nutrients: nitrogen and magnesium. Our population genomic sequencing work, based on this study and previous work, allowed us to identify targets of selection under these nutrient limiting conditions. Global transcriptional regulators glnGL were targets of selection under nitrogen starvation, while proteins involved in outer-membrane biogenesis (genes from the lpt operon) were targets of selection under magnesium starvation. The protein involved in cell-cycle arrest (yhaV) was a target of selection in both environments. We re-constructed specific mutants to analyze the effect of individual mutations on fitness in nutrient limiting conditions in chemostats and in batch cultures. We further demonstrated that adaptation to nitrogen starvation proceeds via a nutrient specific mechanism, while that to magnesium starvation involves a more general mechanism.

CONCLUSIONS

Our results show two different forms of adaptive strategies under limitation of nutrients that effect cellular physiology in different ways. Adaptation to nitrogen starvation proceeds by upregulation of transcriptional regulator glnG and subsequently of transporter protein amtB, both of which results in increased nitrogen scavenging ability of the cell. On the other hand, adaptation to magnesium starvation proceeds via the restructuring of the cell outer-membrane, allowing magnesium to be redistributed to other biological processes. Also, adaptation to the chemostat environment involves selection for loss of function mutations in genes that under nutrient-limiting conditions interfere with continuous growth.

Genetic Adaptation to Growth Under Laboratory Conditions in Escherichia coli and Salmonella enterica

Experimental evolution under controlled laboratory conditions is becoming increasingly important to address various evolutionary questions, including, for example, the dynamics and mechanisms of genetic adaptation to different growth and stress conditions. In such experiments, mutations typically appear that increase the fitness under the conditions tested (medium adaptation), but that are not necessarily of interest for the specific research question. Here, we have identified mutations that appeared during serial passage of E. coli and S. enterica in four different and commonly used laboratory media and measured the relative competitive fitness and maximum growth rate of 111 genetically re-constituted strains, carrying different single and multiple mutations. Little overlap was found between the mutations that were selected in the two species and the different media, implying that adaptation occurs via different genetic pathways. Furthermore, we show that commonly occurring adaptive mutations can generate undesired genetic variation in a population and reduce the accuracy of competition experiments. However, by introducing media adaptation mutations with large effects into the parental strain that was used for the evolution experiment, the variation (standard deviation) was decreased 10-fold, and it was possible to measure fitness differences between two competitors as small as |s| < 0.001.

Experimental evolution of diverse Escherichia coli metabolic mutants identifies genetic loci for convergent adaptation of growth rate

Cell growth is determined by substrate availability and the cell’s metabolic capacity to assimilate substrates into building blocks. Metabolic genes that determine growth rate may interact synergistically or antagonistically, and can accelerate or slow growth, depending on genetic background and environmental conditions. We evolved a diverse set of Escherichia coli single-gene deletion mutants with a spectrum of growth rates and identified mutations that generally increase growth rate. Despite the metabolic differences between parent strains, mutations that enhanced growth largely mapped to core transcription machinery, including the β and β’ subunits of RNA polymerase (RNAP) and the transcription elongation factor, NusA. The structural segments of RNAP that determine enhanced growth have been previously implicated in antibiotic resistance and in the control of transcription elongation and pausing. We further developed a computational framework to characterize how the transcriptional changes that occur upon acquisition of these mutations affect growth rate across strains. Our experimental and computational results provide evidence for cases in which RNAP mutations shift the competitive balance between active transcription and gene silencing. This study demonstrates that mutations in specific regions of RNAP are a convergent adaptive solution that can enhance the growth rate of cells from distinct metabolic states.

The loss of a metabolic function caused by gene deletion can be compensated, in certain cases, by the concurrent mutation of a second gene. Whether such gene pairs share a local chemical or regulatory relationship or interact via a non-local mechanism has implications for the co-evolution of genetic changes, development of alternatives to gene therapy, and the design of combination antimicrobial therapies that select against resistance. Yet, we lack a comprehensive knowledge of adaptive responses to metabolic mutations, and our understanding of the mechanisms underlying genetic rescue remains limited. We present results of a laboratory evolution approach that has the potential to address both challenges, showing that mutations in specific regions of RNA polymerase enhance growth rates of distinct mutant strains of Escherichia coli with a spectrum of growth defects. Several of these adaptive mutations are deleterious when engineered directly into the original wild-type strain under alternative cultivation conditions, and thus have epistatic rescue properties when paired with the corresponding primary metabolic gene deletions. Our combination of adaptive evolution, directed genetic engineering, and mathematical analysis of transcription and growth rate distinguishes between rescue interactions that are specific or non-specific to a particular deletion. Our study further supports a model for RNA polymerase as a locus of convergent adaptive evolution from different sub-optimal metabolic starting points.

Conditional Function of Autoaggregative Protein Cah and Common cah Mutations in Shiga Toxin-Producing Escherichia coli

Cah is a calcium-binding autotransporter protein involved in autoaggregation and biofilm formation. Although cah is widespread in Shiga toxin-producing Escherichia coli (STEC), we detected mutations in cah at a frequency of 31.3% in this pathogen. In STEC O157:H7 supershedder strain SS17, a large deletion results in a smaller coding sequence, encoding a protein lacking the C-terminal 71 amino acids compared with Cah in STEC O157:H7 strain EDL933. We examined the function of Cah in biofilm formation and host colonization to better understand the selective pressures for cah mutations. EDL933-Cah played a conditional role in biofilm formation in vitro: it enhanced E. coli DH5α biofilm formation on glass surfaces under agitated culture conditions that prevented autoaggregation but inhibited biofilm formation under hydrostatic conditions that facilitated autoaggregation. This function appeared to be strain dependent since Cah-mediated biofilm formation was diminished when an EDL933 cah gene was expressed in SS17. Deletion of cah in EDL933 enhanced bacterial attachment to spinach leaves and altered the adherence pattern of EDL933 to bovine recto-anal junction squamous epithelial (RSE) cells. In contrast, in trans expression of EDL933 cah in SS17 increased its attachment to leaf surfaces, and in DH5α, it enhanced its adherence to RSE cells. Hence, the ecological function of Cah appears to be modulated by environmental conditions and other bacterial strain-specific properties. Considering the prevalence of cah in STEC and its role in attachment and biofilm formation, cah mutations might be selected in ecological niches in which inactivation of Cah would result in an increased fitness in STEC during colonization of plants or animal hosts.IMPORTANCE Shiga toxin-producing Escherichia coli (STEC) harbors genes encoding diverse adhesins, and many of these are known to play an important role in bacterial attachment and host colonization. We demonstrated here that the autotransporter protein Cah confers on E. coli DH5α cells a strong autoaggregative phenotype that is inversely correlated with its ability to form biofilms and plays a strain-specific role in plant and animal colonization by STEC. Although cah is widespread in the STEC population, we detected a mutation rate of 31.3% in cah, which is similar to that reported for rpoS and fimH The formation of cell aggregates due to increased bacterium-to-bacterium interactions may be disadvantageous to bacterial populations under conditions that favor a planktonic state in STEC. Therefore, a loss-of-function mutation in cah is likely a selective trait in STEC when autoaggregative properties become detrimental to bacterial cells and may contribute to the adaptability of STEC to fluctuating environments.

Modulation of Global Transcriptional Regulatory Networks as a Strategy for Increasing Kanamycin Resistance of the Translational Elongation Factor-G Mutants in Escherichia coli

Evolve and resequence experiments have provided us a tool to understand bacterial adaptation to antibiotics. In our previous work, we used short-term evolution to isolate mutants resistant to the ribosome targeting antibiotic kanamycin, and reported that Escherichia coli develops low cost resistance to kanamycin via different point mutations in the translation Elongation Factor-G (EF-G). Furthermore, we had shown that the resistance of EF-G mutants could be increased by second site mutations in the genes rpoD/cpxA/topA/cyaA Mutations in three of these genes had been discovered in earlier screens for aminoglycoside resistance. In this work, we expand our understanding of these second site mutations, the goal being to understand how these mutations affect the activities of the mutated gene products to confer resistance. We show that the mutation in cpxA most likely results in an active Cpx stress response. Further evolution of an EF-G mutant in a higher concentration of kanamycin than what was used in our previous experiments identified the cpxA locus as a primary target for a significant increase in resistance. The mutation in cyaA results in a loss of catalytic activity and probably results in resistance via altered CRP function. Despite a reduction in cAMP levels, the CyaA^N600Y mutant has a transcriptome indicative of increased CRP activity, pointing to an unknown role for CyaA and / or cAMP in gene expression. From the transcriptomes of double and single mutants, we describe the epistasis between the mutation in EF-G and these second site mutations. We show that the large scale transcriptomic changes in the topoisomerase I (FusA^A608E-TopA^S180L) mutant likely result from increased negative supercoiling in the cell. Finally, genes with known roles in aminoglycoside resistance were present among the misregulated genes in the mutants.

The fitness costs and benefits of antibiotic resistance in drug-free microenvironments encountered in the human body

The relationship between bacterial drug resistance and growth fitness is a contentious topic, but some antibiotic resistance mutations clearly have a fitness cost in the laboratory. Whether these costs translate into deleterious effects in natural habitats is less certain however. Previously, fitness effects of resistance mutations were mostly characterized in nutrient-rich, fast-growth conditions, which bacteria rarely encounter in natural habitats. Carbon, phosphate, iron or oxygen limitations are conditions met by bacterial pathogens in various compartments of the human body. Here, we measured the fitness of four different rpoB mutations commonly found in rifampicin-resistant bacterial isolates. The fitness properties and the emergence of these and other alleles were studied in Escherichia coli populations growing under nutrient excess and in four different nutrient-limited states. Consistent with previous findings, all four mutations exhibited deleterious fitness effects under nutrient-rich conditions. In stark contrast, we found positive or neutral fitness effects under nutrient-limited conditions. Two particular rpoB alleles had a remarkable fitness increase under phosphate limitation and these alleles arose to high frequencies specifically under phosphate limitation. These findings suggest that it is not meaningful to draw general conclusions on fitness costs without considering bacterial microenvironments in humans and other animals.

Evolutionary implications of Liebig's law of the minimum: Selection under low concentrations of two nonsubstitutable nutrients

Interactions between different axes of an organism's niche determine the evolutionary trajectory of a population. An extreme case of these interactions is predicted from ecological theory in Liebig's law of the minimum. This law states that in environments where multiple nutrients are in relatively low concentrations, only one nutrient will affect the growth of the organism. This implies that the evolutionary response of the population would be dictated by the most growth-limiting nutrient. Alternatively, it is possible that an initial adaptation to the most limiting nutrient results in other nutrients present in low concentration affecting the evolutionary dynamics of the population. To test these hypotheses, we conducted twelve evolution experiments in chemostats using Escherichia coli populations: four under nitrogen limitation, four under magnesium limitation, and four in which both nitrogen and magnesium are in low concentrations. In the last environment, only magnesium seems to limit growth (Low Nitrogen Magnesium Limited environment, LNML). We observe a decrease in nitrogen concentration in the LNML environment over the course of our evolution experiment indicating that nitrogen might become limiting in these environments. Genetic reconstruction results show that clones adapted to magnesium limitation have genes involved in nitrogen starvation, that is, glnG (nitrogen starvation transcriptional regulator) and amtB (transport protein) to be upregulated only in the LNML environment as compared to magnesium-limiting environments. Together, our results highlights that in low-nutrient environments, adaptation to the growth-limiting nutrient results in other nutrients at low concentrations to play a role in the evolutionary dynamics of the population.

An evolutionary optimization of a rhodopsin-based phototrophic metabolism in Escherichia coli

Background

The expression of the Gloeobacter rhodopsin (GR) in a chemotrophic Escherichia coli enables the light-driven phototrophic energy generation. Adaptive laboratory evolution has been used for acquiring desired phenotype of microbial cells and for the elucidation of basic mechanism of molecular evolution. To develop an optimized strain for the artificially acquired phototrophic metabolism, an ancestral E. coli expressing GR was adaptively evolved in a chemostat reactor with constant illumination and limited glucose conditions. This study was emphasized at an unexpected genomic mutation contributed to the improvement of microbial performance.

Results

During the chemostat culture, increase of cell size was observed, which were distinguished from that of the typical rod-shaped ancestral cells. A descendant ET5 strain was randomly isolated from the chemostat culture at 88-days. The phototrophic growth and the light-induced proton pumping of the ET5 strain were twofold and eightfold greater, respectively, than those of the ancestral E. coli strain. Single point mutation of C1082A at dgcQ gene (encoding diguanylate cyclase, also known as the yedQ gene) in the chromosome of ET5 strain was identified from whole genome sequencing analysis. An ancestral E. coli complemented with the same dgcQ mutation from the ET5 was repeated the subsequently enhancements of light-driven phototrophic growth and proton pumping. Intracellular c-di-GMP, the product of the diguanylate cyclase (dgcQ), of the descendant ET5 strain was suddenly increased while that of the ancestral strain was negligible.

Conclusions

Newly acquired phototrophic metabolism of E. coli was further improved via adaptive laboratory evolution by the rise of a point mutation on a transmembrane cell signaling protein followed by increase of signal molecule that eventually led an increase proton pumping and phototrophic growth.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-017-0725-6) contains supplementary material, which is available to authorized users.

Inactivation of Transcriptional Regulators during Within-Household Evolution of Escherichia coli

We analyzed the within-household evolution of two household-associated Escherichia coli strains from pandemic clonal group ST131-H30, using isolates recovered from five individuals within two families, each of which had a distinct strain. Family 1's strain was represented by a urine isolate from the index patient (older sister) with recurrent cystitis and a blood isolate from her younger sister with fatal urosepsis. Family 2's strain was represented by a urine isolate from the index patient (father) with pyelonephritis and renal abscesses, blood and kidney drainage isolates from the daughter with emphysematous pyelonephritis, and urine and fecal isolates from the mother with cystitis. Collectively, the several variants of each family's strain had accumulated a total of 8 (family 1) and 39 (family 2) point mutations; no two isolates were identical. Of the 47 total mutations, 36 resulted in amino acid changes or truncation of coded proteins. Fourteen such mutations (39%) targeted genes encoding transcriptional regulators, and 9 (25%) involved DNA-binding transcription factors (TFs), which significantly exceeded the relative contribution of TF genes to the isolates' genomes (∼6%). At least one-half of the transcriptional regulator mutations were inactivating, based on phenotypic and/or transcriptional analysis. In particular, inactivating mutations in the global regulator LrhA (repressor of type 1 fimbriae and flagella) occurred in the blood isolates from both households and increased the virulence of E. coli strains in a murine sepsis model. The results indicate that E. coli undergoes adaptive evolution between and/or within hosts, generating subpopulations with distinctive phenotypes and virulence potential.IMPORTANCE The clonal evolution of bacterial strains associated with interhost transmission is poorly understood. We characterized the genome sequences of clonal descendants of two Escherichia coli strains, recovered at different time points from multiple individuals within two households who had different types of urinary tract infection. We found evidence that the E. coli strains underwent extensive mutational diversification between and within these individuals, driven disproportionately by inactivation of transcriptional regulators. In urosepsis isolates, the mutations observed in the global regulator LrhA increased bacterial virulence in a murine sepsis model. Our findings help in understanding the adaptive dynamics and strategies of E. coli during short-term natural evolution.

A shifting mutational landscape in 6 nutritional states: Stress-induced mutagenesis as a series of distinct stress input-mutation output relationships

Environmental stresses increase genetic variation in bacteria, plants, and human cancer cells. The linkage between various environments and mutational outcomes has not been systematically investigated, however. Here, we established the influence of nutritional stresses commonly found in the biosphere (carbon, phosphate, nitrogen, oxygen, or iron limitation) on both the rate and spectrum of mutations in Escherichia coli. We found that each limitation was associated with a remarkably distinct mutational profile. Overall mutation rates were not always elevated, and nitrogen, iron, and oxygen limitation resulted in major spectral changes but no net increase in rate. Our results thus suggest that stress-induced mutagenesis is a diverse series of stress input-mutation output linkages that is distinct in every condition. Environment-specific spectra resulted in the differential emergence of traits needing particular mutations in these settings. Mutations requiring transpositions were highest under iron and oxygen limitation, whereas base-pair substitutions and indels were highest under phosphate limitation. The unexpected diversity of input-output effects explains some important phenomena in the mutational biases of evolving genomes. The prevalence of bacterial insertion sequence transpositions in the mammalian gut or in anaerobically stored cultures is due to environmentally determined mutation availability. Likewise, the much-discussed genomic bias towards transition base substitutions in evolving genomes can now be explained as an environment-specific output. Altogether, our conclusion is that environments influence genetic variation as well as selection.

CHO Cells adapted to inorganic phosphate limitation show higher growth and higher pyruvate carboxylase flux in phosphate replete conditions

Inorganic phosphate (P_i ) is an essential ion involved in diverse cellular processes including metabolism. Changes in cellular metabolism upon long term adaptation to P_i limitation have been reported in E. coli. Given the essential role of P_i , adaptation to P_i limitation may also result in metabolic changes in animal cells. In this study, we have adapted CHO cells producing recombinant IgG to limiting P_i conditions for 75 days. Not surprisingly, adapted cells showed better survival under P_i limitation. Here, we report the finding that such cells also showed better growth characteristics compared to control in batch culture replete with P_i (higher peak density and integral viable cell density), accompanied by a lower specific oxygen uptake rate and cytochrome oxidase activity towards the end of exponential phase. Surprisingly, the adapted cells grew to a lower peak density under glucose limitation. This suggests long term P_i limitation may lead to selection for an altered metabolism with higher dependence on glucose availability for biomass assimilation compared to control. Steady state U-¹³ C glucose labeling experiments suggest that adapted cells have a higher pyruvate carboxylase flux. Consistent with this observation, supplementation with aspartate abolished the peak density difference whereas supplementation with serine did not abolish the difference. This supports the hypothesis that cell growth in the adapted culture might be higher due to a higher pyruvate carboxylase flux. Decreased fitness under carbon limitation and mutations in the sucABCD operon has been previously reported in E. coli upon long term adaptation to P_i limitation, suggestive of a similarity in cellular response among such diverse species. © 2017 American Institute of Chemical Engineers Biotechnol. Prog., 33:749-758, 2017.

Characterization of acetic acid-detoxifying Escherichia coli evolved under phosphate starvation conditions

BACKGROUND

During prolonged incubation of Escherichia coli K-12 in batch culture under aerobic, phosphate (Pi) starvation conditions, excess glucose is converted into acetic acid, which may trigger cell death. Following serial cultures, we isolated five evolved strains in two populations that survived prolonged incubation.

METHODS

We sequenced the genomes of the ancestral and evolved strains, and determined the effects of the genetic changes, tested alone and in combination, on characteristic phenotypes in pure and in mixed cultures.

RESULTS

Evolved strains used two main strategies: (1) the constitutive expression of the Trk- and Kdp-dependent K(+) transport systems, and (2) the inactivation of the ArcA global regulator. Both processes helped to maintain a residual activity of the tricarboxylic acid cycle, which decreased the production of acetic acid and eventually allowed its re-consumption. Evolved strains acquired a few additional genetic changes besides the trkH, kdpD and arcA mutations, which might increase the scavenging of organophosphates (phnE (+), lapB, and rseP) and the resistance to oxidative (rsxC) and acetic acid stresses (e14(-)/icd (+)).

CONCLUSIONS

Evolved strains rapidly acquired mutations (phnE (+) lapB rpoS trkH and phnE (+) rseP kdpD) that were globally beneficial to growth on glucose and organophosphates, but detrimental to long-term viability. The spread of these mutant strains might give the ancestral strain time to accumulate up to five genetic changes (phnE (+) arcA rsxC crfC e14(-)/icd (+)), which allowed growth on glucose and organophosphates, and provided a long-term survival. The latter strain, which expressed several mechanisms of protection against endogenous and exogenous stresses, might provide a platform for producing toxic recombinant proteins and chemicals during prolonged incubation under aerobic, Pi starvation conditions.

Chromobacterium violaceum adaptation to low-phosphate conditions

Chromobacterium violaceum is a free-living bacterium that inhabits low-nutrient environments such as the Amazon basin. Bacteria respond to phosphate (Pi) shortage by expressing a range of genes involved in Pi uptake and assimilation, known as the PHO regulon. Several PHO regulon genes have been annotated in the genome of C. violaceum. Here we show that C. violaceum is extremely well adapted to low-Pi conditions. Remarkably, this bacterium is able to grow in media containing only traces of Pi. The PHO regulon genes are induced upon Pi depletion, but the bacteria continued to grow under these conditions. Unlike other Proteobacteria hitherto analyzed, neither PstS nor PhoU play a role in the repression of the PHO regulon under Pi excess.

An Environmental Shiga Toxin-Producing Escherichia coli O145 Clonal Population Exhibits High-Level Phenotypic Variation That Includes Virulence Traits

Shiga toxin-producing Escherichia coli (STEC) serotype O145 is one of the major non-O157 serotypes associated with severe human disease. Here we examined the genetic diversity, population structure, virulence potential, and antimicrobial resistance profiles of environmental O145 strains recovered from a major produce production region in California. Multilocus sequence typing analyses revealed that sequence type 78 (ST-78), a common ST in clinical strains, was the predominant genotype among the environmental strains. Similarly, all California environmental strains belonged to H28, a common H serotype in clinical strains. Although most environmental strains carried an intact fliC gene, only one strain retained swimming motility. Diverse stx subtypes were identified, including stx1a, stx2a, stx2c, and stx2e. Although no correlation was detected between the stx genotype and Stx1 production, high Stx2 production was detected mainly in strains carrying stx2a only and was correlated positively with the cytotoxicity of Shiga toxin. All environmental strains were capable of producing enterohemolysin, whereas only 10 strains were positive for anaerobic hemolytic activity. Multidrug resistance appeared to be common, as nearly half of the tested O145 strains displayed resistance to at least two different classes of antibiotics. The core virulence determinants of enterohemorrhagic E. coli were conserved in the environmental STEC O145 strains; however, there was large variation in the expression of virulence traits among the strains that were highly related genotypically, implying a trend of clonal divergence. Several cattle isolates exhibited key virulence traits comparable to those of the STEC O145 outbreak strains, emphasizing the emergence of hypervirulent strains in agricultural environments.

Bacterial transcriptome reorganization in thermal adaptive evolution

Background

Evolution optimizes a living system at both the genome and transcriptome levels. Few studies have investigated transcriptome evolution, whereas many studies have explored genome evolution in experimentally evolved cells. However, a comprehensive understanding of evolutionary mechanisms requires knowledge of how evolution shapes gene expression. Here, we analyzed Escherichia coli strains acquired during long-term thermal adaptive evolution.

Results

Evolved and ancestor Escherichia coli cells were exponentially grown under normal and high temperatures for subsequent transcriptome analysis. We found that both the ancestor and evolved cells had comparable magnitudes of transcriptional change in response to heat shock, although the evolutionary progression of their expression patterns during exponential growth was different at either normal or high temperatures. We also identified inverse transcriptional changes that were mediated by differences in growth temperatures and genotypes, as well as negative epistasis between genotype—and heat shock-induced transcriptional changes. Principal component analysis revealed that transcriptome evolution neither approached the responsive state at the high temperature nor returned to the steady state at the regular temperature. We propose that the molecular mechanisms of thermal adaptive evolution involve the optimization of steady-state transcriptomes at high temperatures without disturbing the heat shock response.

Conclusions

Our results suggest that transcriptome evolution works to maintain steady-state gene expression during constrained differentiation at various evolutionary stages, while also maintaining responsiveness to environmental stimuli and transcriptome homeostasis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1999-x) contains supplementary material, which is available to authorized users.

Evolutionary Transition from Pathogenicity to Commensalism: Global Regulator Mutations Mediate Fitness Gains through Virulence Attenuation

Symbiotic interactions are indispensable for metazoan function, but their origin and evolution remain elusive. We use a controlled evolution experiment to demonstrate the emergence of novel commensal interactions between Pseudomonas aeruginosa, an initially pathogenic bacterium, and a metazoan host, Caenorhabditis elegans. We show that commensalism evolves through loss of virulence, because it provides bacteria with a double fitness advantage: Increased within-host fitness and a larger host population to infect. Commensalism arises irrespective of host immune status, as the adaptive path in immunocompromised C. elegans knockouts does not differ from that in wild type. Dissection of temporal dynamics of genomic adaptation for 125 bacterial populations reveals highly parallel evolution of incipient commensalism across independent biological replicates. Adaptation is mainly achieved through frame shift mutations in the global regulator lasR and nonsynonymous point mutations in the polymerase gene rpoB that arise early in evolution. Genetic knockouts of lasR not only corroborate its role in virulence attenuation but also show that further mutations are necessary for the fully commensal phenotype. The evolutionary transition from pathogenicity to commensalism as we observe here is facilitated by mutations in global regulators such as lasR, because few genetic changes cause pleiotropic effects across the genome with large phenotypic effects. Finally, we found that nucleotide diversity increased more quickly in bacteria adapting to immunocompromised hosts than in those adapting to immunocompetent hosts. Nevertheless, the outcome of evolution was comparable across host types. Commensalism can thus evolve independently of host immune state solely as a side-effect of bacterial adaptation to novel hosts.

Determination of single nucleotide variants in Escherichia coli DH5α by using short-read sequencing

Escherichia coli DH5α is a common laboratory strain that provides an important platform for routine use in cloning and synthetic biology applications. Many synthetic circuits have been constructed and successfully expressed in E. coli DH5α; however, its genome sequence has not been determined yet. Here, we determined E. coli DH5α genome sequence and identified genetic mutations that affect its phenotypic functions by using short-read sequencing. The sequencing results clearly described the genotypes of E. coli DH5α, which aid in further studies using the strain. Additionally, we observed 105 single nucleotide variants (SNVs), 83% of which were detected in protein-coding regions compared to the parental strain E. coli DH1. Interestingly, 23% of the protein-coding regions have mutations in their amino acid residues, whose biological functions were categorized into two-component systems, peptidoglycan biosynthesis and lipopolysaccharide biosynthesis. These results underscore the advantages of E. coli DH5α, which tolerates the components of transformation buffer and expresses foreign plasmids efficiently. Moreover, these SNVs were also observed in the commercially available strain. These data provide the genetic information of E. coli DH5α for its future application in metabolic engineering and synthetic biology.

Improved feature-based prediction of SNPs in human cytochrome P450 enzymes

Single nucleotide polymorphisms (SNPs) make up the most common form of mutations in human cytochrome P450 enzymes family, and have the potential to bring with different drug responses or specific diseases in individual patients. Here, based on machine learning technology, we aim to explore an effective set of sequence-based features for improving prediction of SNPs by using support vector machine algorithms. The features are derived from the target residues and flanking protein sequences, such as amino acid types, sequences composition, physicochemical properties, position-specific scoring matrix, phylogenetic entropy and the number of possible codons of target residues. In order to deal with the imbalance data with a majority of non-SNPs and a minority of SNPs, a preprocessing strategy based on fuzzy set theory was applied to the datasets. Our final model achieves the performance of 93.8% in sensitivity, 88.8% in specificity, 91.3% in accuracy and 0.971 of AUC value, which is significantly higher than the previous DNA sequence-based or protein sequence-based methods. Furthermore, our study also suggested the roles of individual features for prediction of SNPs. The most important features consist of the amino acid type, the number of available codons, position-specific scoring matrix and phylogenetic entropy. The improved model will be a promising tool for SNP predictions, and assist in the research of genome mutation and personalized prescriptions.

Simple phenotypic sweeps hide complex genetic changes in populations

Changes in allele frequencies and the fixation of beneficial mutations are central to evolution. The precise relationship between mutational and phenotypic sweeps is poorly described however, especially when multiple alleles are involved. Here, we investigate these relationships in a bacterial population over 60 days in a glucose-limited chemostat in a large population. High coverage metagenomic analysis revealed a disconnection between smooth phenotypic sweeps and the complexity of genetic changes in the population. Phenotypic adaptation was due to convergent evolution and involved soft sweeps by 7–26 highly represented alleles of several genes in different combinations. Allele combinations spread from undetectably low baselines, indicating that minor subpopulations provide the basis of most innovations. A hard sweep was also observed, involving a single combination of rpoS, mglD, malE, sdhC, and malT mutations sweeping to greater than 95% of the population. Other mutant genes persisted but at lower abundance, including hfq, consistent with its demonstrated frequency-dependent fitness under glucose limitation. Other persistent, newly identified low-frequency mutations were in the aceF, galF, ribD and asm genes, in noncoding regulatory regions, three large indels and a tandem duplication; these were less affected by fluctuations involving more dominant mutations indicating separate evolutionary paths. Our results indicate a dynamic subpopulation structure with a minimum of 42 detectable mutations maintained over 60 days. We also conclude that the massive population-level mutation supply in combination with clonal interference leads to the soft sweeps observed, but not to the exclusion of an occasional hard sweep.

The functional basis of adaptive evolution in chemostats

Two of the central problems in biology are determining the molecular basis of adaptive evolution and understanding how cells regulate their growth. The chemostat is a device for culturing cells that provides great utility in tackling both of these problems: it enables precise control of the selective pressure under which organisms evolve and it facilitates experimental control of cell growth rate. The aim of this review is to synthesize results from studies of the functional basis of adaptive evolution in long-term chemostat selections using Escherichia coli and Saccharomyces cerevisiae. We describe the principle of the chemostat, provide a summary of studies of experimental evolution in chemostats, and use these studies to assess our current understanding of selection in the chemostat. Functional studies of adaptive evolution in chemostats provide a unique means of interrogating the genetic networks that control cell growth, which complements functional genomic approaches and quantitative trait loci (QTL) mapping in natural populations. An integrated approach to the study of adaptive evolution that accounts for both molecular function and evolutionary processes is critical to advancing our understanding of evolution. By renewing efforts to integrate these two research programs, experimental evolution in chemostats is ideally suited to extending the functional synthesis to the study of genetic networks.

The enduring utility of continuous culturing in experimental evolution

Studying evolution in the laboratory provides a means of understanding the processes, dynamics and outcomes of adaptive evolution in precisely controlled and readily replicated conditions. The advantages of experimental evolution are maximized when the selection is well defined, which enables linking genotype, phenotype and fitness. One means of maintaining a defined selection is continuous culturing: chemostats enable the study of adaptive evolution in constant nutrient-limited environments, whereas cells in turbidostats evolve in constant nutrient abundance. Although the experimental effort required for continuous culturing is considerable relative to the experimental simplicity of serial batch culture, the opposite is true of the environments they produce: continuous culturing results in simplified and invariant conditions whereas serially diluted batch cultures are complex and dynamic. The comparative simplicity of the selective environment that is unique to continuous culturing provides an ideal experimental system for addressing key questions in adaptive evolution.

The rpoS gene is predominantly inactivated during laboratory storage and undergoes source-sink evolution in Escherichia coli species

The rpoS gene codes for an alternative RNA polymerase sigma factor, which acts as a general regulator of the stress response. Inactivating alleles of rpoS in collections of natural Escherichia coli isolates have been observed at very variable frequencies, from less than 1% to more than 70% of strains. rpoS is easily inactivated in nutrient-deprived environments such as stab storage, which makes it difficult to determine the true frequency of rpoS inactivation in nature. We studied the evolutionary history of rpoS and compared it to the phylogenetic history of bacteria in two collections of 82 human commensal and extraintestinal E. coli strains. These strains were representative of the phylogenetic diversity of the species and differed only by their storage conditions. In both collections, the phylogenetic histories of rpoS and of the strains were congruent, indicating that horizontal gene transfer had not occurred at the rpoS locus, and rpoS was under strong purifying selection, with a ratio of the nonsynonymous mutation rate (Ka) to the synonymous substitution rate (Ks) substantially smaller than 1. Stab storage was associated with a high frequency of inactivating alleles, whereas almost no amino acid sequence variation was observed in RpoS in the collection studied directly after isolation of the strains from the host. Furthermore, the accumulation of variations in rpoS was typical of source-sink dynamics. In conclusion, rpoS is rarely inactivated in natural E. coli isolates within their mammalian hosts, probably because such strains rapidly become evolutionary dead ends. Our data should encourage bacteriologists to freeze isolates immediately and to avoid the use of stab storage.

Multihost experimental evolution of the pathogen Ralstonia solanacearum unveils genes involved in adaptation to plants

Ralstonia solanacearum, the causal agent of a lethal bacterial wilt plant disease, infects an unusually wide range of hosts. These hosts can further be split into plants where R. solanacearum is known to cause disease (original hosts) and those where this bacterium can grow asymptomatically (distant hosts). Moreover, this pathogen is able to adapt to many plants as supported by field observations reporting emergence of strains with enlarged pathogenic properties. To investigate the genetic bases of host adaptation, we conducted evolution experiments by serial passages of a single clone of the pathogen on three original and two distant hosts over 300 bacterial generations and then analyzed the whole-genome of nine evolved clones. Phenotypic analysis of the evolved clones showed that the pathogen can increase its fitness on both original and distant hosts although the magnitude of fitness increase was greater on distant hosts. Only few genomic modifications were detected in evolved clones compared with the ancestor but parallel evolutionary changes in two genes were observed in independent evolved populations. Independent mutations in the regulatory gene efpR were selected for in three populations evolved on beans, a distant host. Reverse genetic approaches confirmed that these mutations were associated with fitness gain on bean plants. This work provides a first step toward understanding the within-host evolutionary dynamics of R. solanacearum during infection and identifying bacterial genes subjected to in planta selection. The discovery of EfpR as a determinant conditioning host adaptation of the pathogen illustrates how experimental evolution coupled with whole-genome sequencing is a potent tool to identify novel molecular players involved in central life-history traits.

Characterization of mutations in the PAS domain of the EvgS sensor kinase selected by laboratory evolution for acid resistance in Escherichia coli

Laboratory-based evolution and whole-genome sequencing can link genotype and phenotype. We used evolution of acid resistance in exponential phase Escherichia coli to study resistance to a lethal stress. Iterative selection at pH 2.5 generated five populations that were resistant to low pH in early exponential phase. Genome sequencing revealed multiple mutations, but the only gene mutated in all strains was evgS, part of a two-component system that has already been implicated in acid resistance. All these mutations were in the cytoplasmic PAS domain of EvgS, and were shown to be solely responsible for the resistant phenotype, causing strong upregulation at neutral pH of genes normally induced by low pH. Resistance to pH 2.5 in these strains did not require the transporter GadC, or the sigma factor RpoS. We found that EvgS-dependent constitutive acid resistance to pH 2.5 was retained in the absence of the regulators GadE or YdeO, but was lost if the oxidoreductase YdeP was also absent. A deletion in the periplasmic domain of EvgS abolished the response to low pH, but not the activity of the constitutive mutants. On the basis of these results we propose a model for how EvgS may become activated by low pH.

Natural rpoS mutations contribute to population heterogeneity in Escherichia coli O157:H7 strains linked to the 2006 US spinach-associated outbreak

We previously reported significantly different acid resistance between curli variants derived from the same Escherichia coli O157:H7 strain, although the curli fimbriae were not associated with this phenotypic divergence. Here we investigated the underlying molecular mechanism by examining the genes encoding the common transcriptional regulators of curli biogenesis and acid resistance. rpoS null mutations were detected in all curli-expressing variants of the 2006 spinach-associated outbreak strains, whereas a wild-type rpoS was present in all curli-deficient variants. Consequently curli-expressing variants were much more sensitive to various stress challenges than curli-deficient variants. This loss of general stress fitness appeared solely to be the result of rpoS mutation since the stress resistances could be restored in curli-expressing variants by a functional rpoS. Comparative transcriptomic analyses between the curli variants revealed a large number of differentially expressed genes, characterized by the enhanced expression of metabolic genes in curli-expressing variants, but a marked decrease in transcription of genes related to stress resistances. Unlike the curli-expressing variants of the 1993 US hamburger-associated outbreak strains (Applied Environmental Microbiology 78: 7706-7719), all curli-expressing variants of the 2006 spinach-associated outbreak strains carry a functional rcsB gene, suggesting an alternative mechanism governing intra-strain phenotypic divergence in E. coli O157:H7.

The effect of the rpoSam allele on gene expression and stress resistance in Escherichia coli

The RNA polymerase associated with RpoS transcribes many genes related to stationary phase and stress survival in Escherichia coli. The DNA sequence of rpoS exhibits a high degree of polymorphism. A C to T transition at position 99 of the rpoS ORF, which results in a premature amber stop codon often found in E. coli strains. The rpoSam mutant expresses a truncated and partially functional RpoS protein. Here, we present new evidence regarding rpoS polymorphism in common laboratory E. coli strains. One out of the six tested strains carries the rpoSam allele, but expressed a full-length RpoS protein owing to the presence of an amber supressor mutation. The rpoSam allele was transferred to a non-suppressor background and tested for RpoS level, stress resistance and for the expression of RpoS and sigma70-dependent genes. Overall, the rpoSam strain displayed an intermediate phenotype regarding stress resistance and the expression of σ(S)-dependent genes when compared to the wild-type rpoS(+) strain and to the rpoS null mutant. Surprisingly, overexpression of rpoSam had a differential effect on the expression of the σ(70)-dependent genes phoA and lacZ that, respectively, encode the enzymes alkaline phosphatase and β-galactosidase. The former was enhanced while the latter was inhibited by high levels of RpoSam.

Evolutionary history and genetic parallelism affect correlated responses to evolution

We investigated the relationship between genomic and phenotypic evolution among replicate populations of Escherichia coli evolved for 1000 generations in four different environments. By resequencing evolved genomes, we identified parallel changes in genes encoding transcription regulators within and between environments. Depending on both the environment and the altered gene, genetic parallelism at the gene level involved mutations that affected identical codons, protein domains or were widely distributed across the gene. Evolved clones were characterized by parallel phenotypic changes in their respective evolution environments but also in the three alternative environments. Phenotypic parallelism was high for clones that evolved in the same environment, even in the absence of genetic parallelism. By contrast, clones that evolved in different environments revealed a higher parallelism in correlated responses when they shared mutated genes. Altogether, this work shows that after an environmental change or the colonization of a new habitat, similar ecological performance might be expected from individuals that share mutated genes or that experienced similar past selective pressures.

Molecular specificity, convergence and constraint shape adaptive evolution in nutrient-poor environments

One of the central goals of evolutionary biology is to explain and predict the molecular basis of adaptive evolution. We studied the evolution of genetic networks in Saccharomyces cerevisiae (budding yeast) populations propagated for more than 200 generations in different nitrogen-limiting conditions. We find that rapid adaptive evolution in nitrogen-poor environments is dominated by the de novo generation and selection of copy number variants (CNVs), a large fraction of which contain genes encoding specific nitrogen transporters including PUT4, DUR3 and DAL4. The large fitness increases associated with these alleles limits the genetic heterogeneity of adapting populations even in environments with multiple nitrogen sources. Complete identification of acquired point mutations, in individual lineages and entire populations, identified heterogeneity at the level of genetic loci but common themes at the level of functional modules, including genes controlling phosphatidylinositol-3-phosphate metabolism and vacuole biogenesis. Adaptive strategies shared with other nutrient-limited environments point to selection of genetic variation in the TORC1 and Ras/PKA signaling pathways as a general mechanism underlying improved growth in nutrient-limited environments. Within a single population we observed the repeated independent selection of a multi-locus genotype, comprised of the functionally related genes GAT1, MEP2 and LST4. By studying the fitness of individual alleles, and their combination, as well as the evolutionary history of the evolving population, we find that the order in which these mutations are acquired is constrained by epistasis. The identification of repeatedly selected variation at functionally related loci that interact epistatically suggests that gene network polymorphisms (GNPs) may be a frequent outcome of adaptive evolution. Our results provide insight into the mechanistic basis by which cells adapt to nutrient-limited environments and suggest that knowledge of the selective environment and the regulatory mechanisms important for growth and survival in that environment greatly increase the predictability of adaptive evolution.

We studied adaptive evolution in different nitrogen-limited environments using long-term selection of asexually reproducing Saccharomyces cerevisiae populations in chemostats. Using next generation sequencing and DNA microarrays, we identified all acquired genetic variation associated with increased fitness, in both individual lineages and entire populations. We find that amplification alleles that include nutrient transporter genes specific to the molecular form of the nitrogen present in the environment are a common mechanism underlying increased fitness. In addition, we identified a general strategy for adaptation to nitrogen-limited environments that entails remodeling of phospholipid biogenesis required for producing important cellular components including vacuoles and autophagosomes. More general strategies for adaptation to nutrient-limited environments point to a role for re-wiring of signaling pathways that coordinate cell growth with nutrient availability. We reconstructed the evolutionary dynamics of a population evolving in ammonium-limited conditions and find that a multi-locus genotype is repeatedly selected within the population and constrained by epistasis. We propose that this genotype constitutes a “gene network polymorphism (GNP),” which may be a common outcome of adaptive evolution. Our study suggests that when the selective pressure is understood the molecular basis of adaptive evolution in large microbial populations may be predicted with reasonable precision.

A metabolic trade-off between phosphate and glucose utilization in Escherichia coli

Using full genome sequencing and metabolomics we show that adaptation to chronic nutrient starvation reduces metabolic flexibility inEscherichia coli.