Adaptive mgl-regulatory mutations and genetic diversity evolving in glucose-limited Escherichia coli populations

Glucose consumption rate-dependent transcriptome profiling of Escherichia coli provides insight on performance as microbial factories

Background

The modification of glucose import capacity is an engineering strategy that has been shown to improve the characteristics of Escherichia coli as a microbial factory. A reduction in glucose import capacity can have a positive effect on production strain performance, however, this is not always the case. In this study, E. coli W3110 and a group of four isogenic derivative strains, harboring single or multiple deletions of genes encoding phosphoenolpyruvate:sugar phosphotransferase system (PTS)-dependent transporters as well as non-PTS transporters were characterized by determining their transcriptomic response to reduced glucose import capacity.

Results

These strains were grown in bioreactors with M9 mineral salts medium containing 20 g/L of glucose, where they displayed specific growth rates ranging from 0.67 to 0.27 h⁻¹, and specific glucose consumption rates (qs) ranging from 1.78 to 0.37 g/g h. RNA-seq analysis revealed a transcriptional response consistent with carbon source limitation among all the mutant strains, involving functions related to transport and metabolism of alternate carbon sources and characterized by a decrease in genes encoding glycolytic enzymes and an increase in gluconeogenic functions. A total of 107 and 185 genes displayed positive and negative correlations with qs, respectively. Functions displaying positive correlation included energy generation, amino acid biosynthesis, and sugar import.

Conclusion

Changes in gene expression of E. coli strains with impaired glucose import capacity could be correlated with qs values and this allowed an inference of the physiological state of each mutant. In strains with lower qs values, a gene expression pattern is consistent with energy limitation and entry into the stationary phase. This physiological state could explain why these strains display a lower capacity to produce recombinant protein, even when they show very low rates of acetate production. The comparison of the transcriptomes of the engineered strains employed as microbial factories is an effective approach for identifying favorable phenotypes with the potential to improve the synthesis of biotechnological products.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-022-01909-y.

A mutational hotspot that determines highly repeatable evolution can be built and broken by silent genetic changes

Mutational hotspots can determine evolutionary outcomes and make evolution repeatable. Hotspots are products of multiple evolutionary forces including mutation rate heterogeneity, but this variable is often hard to identify. In this work, we reveal that a near-deterministic genetic hotspot can be built and broken by a handful of silent mutations. We observe this when studying homologous immotile variants of the bacteria Pseudomonas fluorescens, AR2 and Pf0-2x. AR2 resurrects motility through highly repeatable de novo mutation of the same nucleotide in >95% lines in minimal media (ntrB A289C). Pf0-2x, however, evolves via a number of mutations meaning the two strains diverge significantly during adaptation. We determine that this evolutionary disparity is owed to just 6 synonymous variations within the ntrB locus, which we demonstrate by swapping the sites and observing that we are able to both break (>95% to 0%) and build (0% to 80%) a deterministic mutational hotspot. Our work reveals a key role for silent genetic variation in determining adaptive outcomes.

Genetic interaction network has a very limited impact on the evolutionary trajectories in continuous culture-grown populations of yeast

Background

The impact of genetic interaction networks on evolution is a fundamental issue. Previous studies have demonstrated that the topology of the network is determined by the properties of the cellular machinery. Functionally related genes frequently interact with one another, and they establish modules, e.g., modules of protein complexes and biochemical pathways. In this study, we experimentally tested the hypothesis that compensatory evolutionary modifications, such as mutations and transcriptional changes, occur frequently in genes from perturbed modules of interacting genes.

Results

Using Saccharomyces cerevisiae haploid deletion mutants as a model, we investigated two modules lacking COG7 or NUP133, which are evolutionarily conserved genes with many interactions. We performed laboratory evolution experiments with these strains in two genetic backgrounds (with or without additional deletion of MSH2), subjecting them to continuous culture in a non-limiting minimal medium. Next, the evolved yeast populations were characterized through whole-genome sequencing and transcriptome analyses. No obvious compensatory changes resulting from inactivation of genes already included in modules were identified. The supposedly compensatory inactivation of genes in the evolved strains was only rarely observed to be in accordance with the established fitness effect of the genetic interaction network. In fact, a substantial majority of the gene inactivations were predicted to be neutral in the experimental conditions used to determine the interaction network. Similarly, transcriptome changes during continuous culture mostly signified adaptation to growth conditions rather than compensation of the absence of the COG7, NUP133 or MSH2 genes. However, we noticed that for genes whose inactivation was deleterious an upregulation of transcription was more common than downregulation.

Conclusions

Our findings demonstrate that the genetic interactions and the modular structure of the network described by others have very limited effects on the evolutionary trajectory following gene deletion of module elements in our experimental conditions and has no significant impact on short-term compensatory evolution. However, we observed likely compensatory evolution in functionally related (albeit non-interacting) genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01830-9.

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.

Evolution of an Escherichia coli PTS- strain: a study of reproducibility and dynamics of an adaptive evolutive process

Adaptive laboratory evolution (ALE) has been used to study and solve pressing questions about evolution, especially for the study of the development of mutations that confer increased fitness during evolutionary processes. In this contribution, we investigated how the evolutionary process conducted with the PTS^- mutant of Escherichia coli PB11 in three parallel batch cultures allowed the restoration of rapid growth with glucose as the carbon source. The significant findings showed that genomic sequence analysis of a set of newly evolved mutants isolated from ALE experiments 2-3 developed some essential mutations, which efficiently improved the fast-growing phenotypes throughout different fitness landscapes. Regulator galR was the target of several mutations such as SNPs, partial and total deletions, and insertion of an IS1 element and thus indicated the relevance of a null mutation of this gene in the adaptation of the evolving population of PB11 during the parallel ALE experiments. These mutations resulted in the selection of MglB and GalP as the primary glucose transporters by the evolving population, but further selection of at least a second adaptive mutation was also necessary. We found that mutations in the yfeO, rppH, and rng genes improved the fitness advantage of evolving PTS^- mutants and resulted in amplification of leaky activity in Glk for glucose phosphorylation and upregulation of glycolytic and other growth-related genes. Notably, we determined that these mutations appeared and were fixed in the evolving populations between 48 and 72 h of cultivation, which resulted in the selection of fast-growing mutants during one ALE experiments in batch cultures of 80 h duration.Key points• ALE experiments selected evolved mutants through different fitness landscapes in which galR was the target of different mutations: SNPs, deletions, and insertion of IS.• Key mutations in evolving mutants appeared and fixed at 48-72 h of cultivation.• ALE experiments led to increased understanding of the genetics of cellular adaptation to carbon source limitation.

Chasing bacterial chassis for metabolic engineering: a perspective review from classical to non-traditional microorganisms

The last few years have witnessed an unprecedented increase in the number of novel bacterial species that hold potential to be used for metabolic engineering. Historically, however, only a handful of bacteria have attained the acceptance and widespread use that are needed to fulfil the needs of industrial bioproduction - and only for the synthesis of very few, structurally simple compounds. One of the reasons for this unfortunate circumstance has been the dearth of tools for targeted genome engineering of bacterial chassis, and, nowadays, synthetic biology is significantly helping to bridge such knowledge gap. Against this background, in this review, we discuss the state of the art in the rational design and construction of robust bacterial chassis for metabolic engineering, presenting key examples of bacterial species that have secured a place in industrial bioproduction. The emergence of novel bacterial chassis is also considered at the light of the unique properties of their physiology and metabolism, and the practical applications in which they are expected to outperform other microbial platforms. Emerging opportunities, essential strategies to enable successful development of industrial phenotypes, and major challenges in the field of bacterial chassis development are also discussed, outlining the solutions that contemporary synthetic biology-guided metabolic engineering offers to tackle these issues.

Contingency and determinism in evolution: Replaying life’s tape

Historical processes display some degree of “contingency,” meaning their outcomes are sensitive to seemingly inconsequential events that can fundamentally change the future. Contingency is what makes historical outcomes unpredictable. Unlike many other natural phenomena, evolution is a historical process. Evolutionary change is often driven by the deterministic force of natural selection, but natural selection works upon variation that arises unpredictably through time by random mutation, and even beneficial mutations can be lost by chance through genetic drift. Moreover, evolution has taken place within a planetary environment with a particular history of its own. This tension between determinism and contingency makes evolutionary biology a kind of hybrid between science and history. While philosophers of science examine the nuances of contingency, biologists have performed many empirical studies of evolutionary repeatability and contingency. Here, we review the experimental and comparative evidence from these studies. Replicate populations in evolutionary “replay” experiments often show parallel changes, especially in overall performance, although idiosyncratic outcomes show that the particulars of a lineage’s history can affect which of several evolutionary paths is taken. Comparative biologists have found many notable examples of convergent adaptation to similar conditions, but quantification of how frequently such convergence occurs is difficult. On balance, the evidence indicates that evolution tends to be surprisingly repeatable among closely related lineages, but disparate outcomes become more likely as the footprint of history grows deeper. Ongoing research on the structure of adaptive landscapes is providing additional insight into the interplay of fate and chance in the evolutionary process.

Effect of Global Regulators RpoS and Cyclic-AMP/CRP on the Catabolome and Transcriptome of Escherichia coli K12 during Carbon- and Energy-Limited Growth

For heterotrophic microbes, limited availability of carbon and energy sources is one of the major nutritional factors restricting the rate of growth in most ecosystems. Physiological adaptation to this hunger state requires metabolic versatility which usually involves expression of a wide range of different catabolic pathways and of high-affinity carbon transporters; together, this allows for simultaneous utilization of mixtures of carbonaceous compounds at low concentrations. In Escherichia coli the stationary phase sigma factor RpoS and the signal molecule cAMP are the major players in the regulation of transcription under such conditions; however, their interaction is still not fully understood. Therefore, during growth of E. coli in carbon-limited chemostat culture at different dilution rates, the transcriptomes, expression of periplasmic proteins and catabolomes of strains lacking one of these global regulators, either rpoS or adenylate cyclase (cya), were compared to those of the wild-type strain. The inability to synthesize cAMP exerted a strong negative influence on the expression of alternative carbon source uptake and degradation systems. In contrast, absence of RpoS increased the transcription of genes belonging to high-affinity uptake systems and central metabolism, presumably due to reduced competition of σ^D with σ^S. Phenotypical analysis confirmed this observation: The ability to respire alternative carbon substrates and to express periplasmic high-affinity binding proteins was eliminated in cya and crp mutants, while these properties were not affected in the rpoS mutant. As expected, transcription of numerous stress defence genes was negatively affected by the rpoS knock-out mutation. Interestingly, several genes of the RpoS stress response regulon were also down-regulated in the cAMP-negative strain indicating a coordinated global regulation. The results demonstrate that cAMP is crucial for catabolic flexibility during slow, carbon-limited growth, whereas RpoS is primarily involved in the regulation of stress response systems necessary for the survival of this bacterium under hunger conditions.

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.

Physiological impact of transposable elements encoding DDE transposases in the environmental adaptation of Streptococcus agalactiae

We have referenced and described Streptococcus agalactiae transposable elements encoding DDE transposases. These elements belonged to nine families of insertion sequences (ISs) and to a family of conjugative transposons (TnGBSs). An overview of the physiological impact of the insertion of all these elements is provided. DDE-transposable elements affect S. agalactiae in a number of aspects of its capability to adapt to various environments and modulate the expression of several virulence genes, the scpB–lmB genomic region and the genes involved in capsule expression and haemolysin transport being the targets of several different mobile elements. The referenced mobile elements modify S. agalactiae behaviour by transferring new gene(s) to its genome, by modifying the expression of neighbouring genes at the integration site or by promoting genomic rearrangements. Transposition of some of these elements occurs in vivo, suggesting that by dynamically regulating some adaptation and/or virulence genes, they improve the ability of S. agalactiae to reach different niches within its host and ensure the ‘success’ of the infectious process.

Ex uno plures: clonal reinforcement drives evolution of a simple microbial community

A major goal of genetics is to define the relationship between phenotype and genotype, while a major goal of ecology is to identify the rules that govern community assembly. Achieving these goals by analyzing natural systems can be difficult, as selective pressures create dynamic fitness landscapes that vary in both space and time. Laboratory experimental evolution offers the benefit of controlling variables that shape fitness landscapes, helping to achieve both goals. We previously showed that a clonal population of E. coli experimentally evolved under continuous glucose limitation gives rise to a genetically diverse community consisting of one clone, CV103, that best scavenges but incompletely utilizes the limiting resource, and others, CV101 and CV116, that consume its overflow metabolites. Because this community can be disassembled and reassembled, and involves cooperative interactions that are stable over time, its genetic diversity is sustained by clonal reinforcement rather than by clonal interference. To understand the genetic factors that produce this outcome, and to illuminate the community's underlying physiology, we sequenced the genomes of ancestral and evolved clones. We identified ancestral mutations in intermediary metabolism that may have predisposed the evolution of metabolic interdependence. Phylogenetic reconstruction indicates that the lineages that gave rise to this community diverged early, as CV103 shares only one Single Nucleotide Polymorphism with the other evolved clones. Underlying CV103's phenotype we identified a set of mutations that likely enhance glucose scavenging and maintain redox balance, but may do so at the expense of carbon excreted in overflow metabolites. Because these overflow metabolites serve as growth substrates that are differentially accessible to the other community members, and because the scavenging lineage shares only one SNP with these other clones, we conclude that this lineage likely served as an “engine” generating diversity by creating new metabolic niches, but not the occupants themselves.

The variability of natural systems makes it difficult to deduce how organisms' genotypes manifest as phenotypes, and how communities of interacting organisms arise. Using laboratory experimental evolution we can control this variation. We previously showed that a population of E. coli that originated from a single clone and was cultured in the presence of a single limiting resource, evolves into a stable, three-membered community, wherein one clone excretes metabolites that the others utilize as carbon sources. To discern the genetic factors at work in producing this outcome and to illuminate the community's physiology, we sequenced the genomes of the ancestral and evolved clones. We identified in the ancestor mutations that may have predisposed evolution of cross-feeding. We found that the lineages which gave rise to the community diverged early on, and that the numerically dominant lineage that best scavenges limiting glucose does so as a result of adaptive mutations that enhance glucose uptake but favor fermentative over respiratory pathways, resulting in overflow metabolites. Because this clone produces secondary resources that sustain other community members, and because it shares with them only one mutation, we conclude that it is an “engine” generating diversity by creating new niches, but not the occupants themselves.

Starvation induces phenotypic diversification and convergent evolution in Vibrio vulnificus

Starvation is a common stress experienced by bacteria living in natural environments and the ability to adapt to and survive intense stress is of paramount importance for any bacterial population. A series of starvation experiments were conducted using V. vulnificus 93U204 in phosphate-buffered saline and seawater. The starved population entered the death phase during the first week and approximately 1% of cells survived. After that the population entered a long-term stationary phase, and could survive for years. Starvation-induced diversification (SID) of phenotypes was observed in starved populations and phenotypic variants (PVs) appeared in less than 8 days. The cell density, rather than the population size, had a major effect on the extent of SID. SID was also observed in strain YJ016, where it evolved at a faster pace. PVs appeared to emerge in a fixed order: PV with reduced motility, PV with reduced proteolytic activity, and PV with reduced hemolytic activity. All of the tested PVs had growth advantages in the stationary phase phenotypes and increased fitness compared with 93U204 cells in co-culture competition experiments, which indicates that they had adapted to starvation. We also found that SID occurred in natural seawater with a salinity of 1%–3%, so this mechanism may facilitate bacterial adaptation in natural environments.

A case of adaptation through a mutation in a tandem duplication during experimental evolution in Escherichia coli

Background

DNA duplications constitute important precursors for genome variation. Here we analyzed an unequal duplication harboring a beneficial mutation that may provide alternative evolutionary outcomes.

Results

We characterized this evolutionary event during experimental evolution for only 100 generations of an Escherichia coli strain under glucose limitation within chemostats. By combining Insertion Sequence based Restriction Length Polymorphism experiments, pulsed field gel electrophoresis and two independent genome re-sequencing experiments, we identified an evolved lineage carrying a 180 kb duplication of the 46’ region of the E. coli chromosome. This evolved duplication revealed a heterozygous state, with one copy harboring a 2668 bp deletion that included part of the ogrK gene and both the yegR and yegS genes. By genetically manipulating ancestral and evolved strains, we showed that the single yegS inactivation was sufficient to confer a frequency dependent fitness increase under the chemostat selective conditions in both the ancestor and evolved genetic contexts, implying that the duplication itself was not a direct fitness contributor. Nonetheless, the heterozygous duplicated state was relatively stable in the conditions prevailing during evolution in chemostats, in striking contrast to non selective conditions in which the duplication resolved at high frequency into either its ancestral or deleted copy.

Conclusions

Our results suggest that the duplication state may constitute a second order selection process providing higher evolutionary potential. Moreover, its heterozygous nature may provide differential evolutionary opportunities in alternating environments. Our results also highlighted how careful analyses of whole genome data are needed to identify such complex rearrangements.

Epistatic interactions determine the mutational pathways and coexistence of lineages in clonal Escherichia coli populations

Understanding how diversity emerges in a single niche is not fully understood. Rugged fitness landscapes and epistasis between beneficial mutations could explain coexistence among emerging lineages. To provide an experimental test of this notion, we investigated epistasis among four pleiotropic mutations in rpoS, mglD, malT, and hfq present in two coexisting lineages that repeatedly fixed in experimental populations of Escherichia coli. The mutations were transferred into the ancestral background individually or in combination of double or triple alleles. The combined competitive fitness of two or three beneficial mutations from the same lineage was consistently lower than the sum of the competitive fitness of single mutants--a clear indication of negative epistasis within lineages. We also found sign epistasis (i.e., the combined fitness of two beneficial mutations lower than the ancestor), not only from two different lineages (i.e., hfq and rpoS) but also from the same lineage (i.e., mglD and malT). The sign epistasis between loci of different lineages indeed indicated a rugged fitness landscape, providing an epistatic explanation for the coexistence of distinct rpoS and hfq lineages in evolving populations. The negative and sign epistasis between beneficial mutations within the same lineage can further explain the order of mutation acquisition.

Can you sequence ecology? Metagenomics of adaptive diversification

Few areas of science have benefited more from the expansion in sequencing capability than the study of microbial communities. Can sequence data, besides providing hypotheses of the functions the members possess, detect the evolutionary and ecological processes that are occurring? For example, can we determine if a species is adapting to one niche, or if it is diversifying into multiple specialists that inhabit distinct niches? Fortunately, adaptation of populations in the laboratory can serve as a model to test our ability to make such inferences about evolution and ecology from sequencing. Even adaptation to a single niche can give rise to complex temporal dynamics due to the transient presence of multiple competing lineages. If there are multiple niches, this complexity is augmented by segmentation of the population into multiple specialists that can each continue to evolve within their own niche. For a known example of parallel diversification that occurred in the laboratory, sequencing data gave surprisingly few obvious, unambiguous signs of the ecological complexity present. Whereas experimental systems are open to direct experimentation to test hypotheses of selection or ecological interaction, the difficulty in "seeing ecology" from sequencing for even such a simple system suggests translation to communities like the human microbiome will be quite challenging. This will require both improved empirical methods to enhance the depth and time resolution for the relevant polymorphisms and novel statistical approaches to rigorously examine time-series data for signs of various evolutionary and ecological phenomena within and between species.

Adaptation and heterogeneity of Escherichia coli MC1000 growing in complex environments

In a study aiming to assess bacterial evolution in complex growth media, we evaluated the long-term adaptive response of Escherichia coli MC1000 in Luria-Bertani (LB) medium. Seven parallel populations were founded and followed over 150 days in sequential batch cultures under three different oxygen conditions (defined environments), and 19 evolved forms were isolated. The emergence of forms with enhanced fitness was evident in competition experiments of all evolved forms versus the ancestral strain. The evolved forms were then subjected to phenotypic and genomic analyses relative to the ancestor. Profound changes were found in their phenotypes as well as whole-genome sequences. Interestingly, considerable heterogeneity was found at the intrapopulational level. However, consistently occurring parallel adaptive responses were found across all populations. The evolved forms all contained a mutation in galR, a repressor of the galactose operon. Concomitantly, the new forms revealed enhanced growth on galactose as well as galactose-containing disaccharides. This response was likely driven by the LB medium.

Synchronous waves of failed soft sweeps in the laboratory: remarkably rampant clonal interference of alleles at a single locus

It has increasingly been recognized that adapting populations of microbes contain not one, but many lineages continually arising and competing at once. This process, termed "clonal interference," alters the rate and dynamics of adaptation and biases winning mutations toward those with the largest selective effect. Here we uncovered a dramatic example of clonal interference between multiple similar mutations occurring at the same locus within replicate populations of Methylobacterium extorquens AM1. Because these mutational events involved the transposition of an insertion sequence into a narrow window of a single gene, they were both readily detectable at low frequencies and could be distinguished due to differences in insertion sites. This allowed us to detect up to 17 beneficial alleles of this type coexisting in a single population. Despite conferring a large selective benefit, the majority of these alleles rose and then fell in frequency due to other lineages emerging that were more fit. By comparing allele-frequency dynamics to the trajectories of fitness gains by these populations, we estimated the fitness values of the genotypes that contained these mutations. Collectively across all populations, these alleles arose upon backgrounds with a wide range of fitness values. Within any single population, however, multiple alleles tended to rise and fall synchronously during a single wave of multiple genotypes with nearly identical fitness values. These results suggest that alleles of large benefit arose repeatedly in failed "soft sweeps" during narrow windows of adaptation due to the combined effects of epistasis and clonal interference.

Mutation accumulation and fitness in mutator subpopulations of Escherichia coli

Bacterial populations in clinical and laboratory settings contain a significant proportion of mutants with elevated mutation rates (mutators). Mutators have a particular advantage when multiple beneficial mutations are needed for fitness, as in antibiotic resistance. Nevertheless, high mutation rates potentially lead to increasing numbers of deleterious mutations and subsequently to the decreased fitness of mutators. To test how fitness changed with mutation accumulation, genome sequencing and fitness assays of nine Escherichia coli mutY mutators were undertaken in an evolving chemostat population at three time points. Unexpectedly, the fitness in members of the mutator subpopulation became constant despite a growing number of mutations over time. To test if the accumulated mutations affected fitness, we replaced each of the known beneficial mutations with wild-type alleles in a mutator isolate. We found that the other 25 accumulated mutations were not deleterious. Our results suggest that isolates with deleterious mutations are eliminated by competition in a continuous culture, leaving mutators with mostly neutral mutations. Interestingly, the mutator-non-mutator balance in the population reversed after the fitness plateau of mutators was reached, suggesting that the mutator-non-mutator ratio in populations has more to do with competition between members of the population than the accumulation of deleterious mutations.

The basis of antagonistic pleiotropy in hfq mutations that have opposite effects on fitness at slow and fast growth rates

Mutations beneficial in one environment may cause costs in different environments, resulting in antagonistic pleiotropy. Here, we describe a novel form of antagonistic pleiotropy that operates even within the same environment, where benefits and deleterious effects exhibit themselves at different growth rates. The fitness of hfq mutations in Escherichia coli affecting the RNA chaperone involved in small-RNA regulation is remarkably sensitive to growth rate. E. coli populations evolving in chemostats under nutrient limitation acquired beneficial mutations in hfq during slow growth (0.1 h(-1)) but not in populations growing sixfold faster. Four identified hfq alleles from parallel populations were beneficial at 0.1 h(-1) and deleterious at 0.6 h(-1). The hfq mutations were beneficial, deleterious or neutral at an intermediate growth rate (0.5 h(-1)) and one changed from beneficial to deleterious within a 36 min difference in doubling time. The benefit of hfq mutations was due to the greater transport of limiting nutrient, which diminished at higher growth rates. The deleterious effects of hfq mutations at 0.6 h(-1) were less clear, with decreased viability a contributing factor. The results demonstrate distinct pleiotropy characteristics in the alleles of the same gene, probably because the altered residues in Hfq affected the regulation of expression of different genes in distinct ways. In addition, these results point to a source of variation in experimental measurement of the selective advantage of a mutation; estimates of fitness need to consider variation in growth rate impacting on the magnitude of the benefit of mutations and on their fitness distributions.

Glucose transport in Escherichia coli mutant strains with defects in sugar transport systems

In Escherichia coli, several systems are known to transport glucose into the cytoplasm. The main glucose uptake system under batch conditions is the glucose phosphoenolpyruvate:carbohydrate phosphotransferase system (glucose PTS), but the mannose PTS and the galactose and maltose transporters also can translocate glucose. Mutant strains which lack the enzyme IIBC (EIIBC) protein of the glucose PTS have been investigated previously because their lower rate of acetate formation offers advantages in industrial applications. Nevertheless, a systematic study to analyze the impact of the different glucose uptake systems has not been undertaken. Specifically, how the bacteria cope with the deletion of the major glucose uptake system and which alternative transporters react to compensate for this deficit have not been studied in detail. Therefore, a series of mutant strains were analyzed in aerobic and anaerobic batch cultures, as well as glucose-limited continuous cultivations. Deletion of EIIBC disturbs glucose transport severely in batch cultures; cyclic AMP (cAMP)-cAMP receptor protein (CRP) levels rise, and induction of the mgl operon occurs. Nevertheless, Mgl activity is not essential for growth of these mutants, since deletion of this transporter did not affect the growth rate; the activities of the remaining transporters seem to be sufficient. Under conditions of glucose limitation, mgl is upregulated 23-fold compared to levels for growth under glucose excess. Despite the strong induction of mgl upon glucose limitation, deletion of this transport system did not lead to further changes. Although the galactose transporters are often regarded as important for glucose uptake at micromolar concentrations, the glucose as well as mannose PTS might be sufficient for growth at this relatively low dilution rate.

Evolutionary insight from whole-genome sequencing of experimentally evolved microbes

Experimental evolution (EE) combined with whole-genome sequencing (WGS) has become a compelling approach to study the fundamental mechanisms and processes that drive evolution. Most EE-WGS studies published to date have used microbes, owing to their ease of propagation and manipulation in the laboratory and relatively small genome sizes. These experiments are particularly suited to answer long-standing questions such as: How many mutations underlie adaptive evolution, and how are they distributed across the genome and through time? Are there general rules or principles governing which genes contribute to adaptation, and are certain kinds of genes more likely to be targets than others? How common is epistasis among adaptive mutations, and what does this reveal about the variety of genetic routes to adaptation? How common is parallel evolution, where the same mutations evolve repeatedly and independently in response to similar selective pressures? Here, we summarize the significant findings of this body of work, identify important emerging trends and propose promising directions for future research. We also outline an example of a computational pipeline for use in EE-WGS studies, based on freely available bioinformatics tools.

Optimization of gene expression through divergent mutational paths

Adaptation under similar selective pressure often leads to comparable phenotypes. A longstanding question is whether such phenotypic repeatability entails similar (parallelism) or different genotypic changes (convergence). To better understand this, we characterized mutations that optimized expression of a plasmid-borne metabolic pathway during laboratory evolution of a bacterium. Expressing these pathway genes was essential for growth but came with substantial costs. Starting from overexpression, replicate populations founded by this bacterium all evolved to reduce expression. Despite this phenotypic repetitiveness, the underlying mutational spectrum was highly diverse. Analysis of these plasmid mutations identified three distinct means to modulate gene expression: (1) reducing the gene copy number, (2) lowering transcript stability, and (3) integration of the pathway-bearing plasmid into the host genome. Our study revealed diverse molecular changes beneath convergence to a simple phenotype. This complex genotype-phenotype mapping presents a challenge to inferring genetic evolution based solely on phenotypic changes.

Microbial laboratory evolution in the era of genome-scale science

Advances in DNA sequencing, high-throughput technologies, and genetic manipulation systems have enabled empirical studies of the molecular and genomic bases of adaptive evolution. This review discusses key insights learned from direct observation of the evolution process.

Laboratory evolution studies provide fundamental biological insight through direct observation of the evolution process. They not only enable testing of evolutionary theory and principles, but also have applications to metabolic engineering and human health. Genome-scale tools are revolutionizing studies of laboratory evolution by providing complete determination of the genetic basis of adaptation and the changes in the organism's gene expression state. Here, we review studies centered on four central themes of laboratory evolution studies: (1) the genetic basis of adaptation; (2) the importance of mutations to genes that encode regulatory hubs; (3) the view of adaptive evolution as an optimization process; and (4) the dynamics with which laboratory populations evolve.

Stock culture heterogeneity rather than new mutational variation complicates short-term cell physiology studies of Escherichia coli K-12 MG1655 in continuous culture

Nutrient-limited continuous cultures in chemostats have been used to study microbial cell physiology for over 60 years. Genome instability and genetic heterogeneity are possible uncontrolled factors in continuous cultivation experiments. We investigated these issues by using high-throughput (HT) DNA sequencing to characterize samples from different phases of a glucose-limited accelerostat (A-stat) experiment with Escherichia coli K-12 MG1655 and a duration regularly used in cell physiology studies (20 generations of continuous cultivation). Seven consensus mutations from the reference sequence and five subpopulations characterized by different mutations were detected in the HT-sequenced samples. This genetic heterogeneity was confirmed to result from the stock culture by Sanger sequencing. All the subpopulations in which allele frequencies increased (betA, cspG/cspH, glyA) during the experiment were also present at the end of replicate A-stats, indicating that no new subpopulations emerged during our experiments. The fact that ~31 % of the cells in our initial cultures obtained directly from a culture stock centre were mutants raises concerns that even if cultivations are started from single colonies, there is a significant chance of picking a mutant clone with an altered phenotype. Our results show that current HT DNA sequencing technology allows accurate subpopulation analysis and demonstrates that a glucose-limited E. coli K-12 MG1655 A-stat experiment with a duration of tens of generations is suitable for studying cell physiology and collecting quantitative data for metabolic modelling without interference from new mutations.

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

Many of the important changes in evolution are regulatory in nature. Sequenced bacterial genomes point to flexibility in regulatory circuits but we do not know how regulation is remodeled in evolving bacteria. Here, we study the regulatory changes that emerge in populations evolving under controlled conditions during experimental evolution of Escherichia coli in a phosphate-limited chemostat culture. Genomes were sequenced from five clones with different combinations of phenotypic properties that coexisted in a population after 37 days. Each of the distinct isolates contained a different mutation in 1 of 3 highly pleiotropic regulatory genes (hfq, spoT, or rpoS). The mutations resulted in dissimilar proteomic changes, consistent with the documented effects of hfq, spoT, and rpoS mutations. The different mutations do share a common benefit, however, in that the mutations each redirect cellular resources away from stress responses that are redundant in a constant selection environment. The hfq mutation lowers several individual stress responses as well the small RNA–dependent activation of rpoS translation and hence general stress resistance. The spoT mutation reduces ppGpp levels, decreasing the stringent response as well as rpoS expression. The mutations in and upstream of rpoS resulted in partial or complete loss of general stress resistance. Our observations suggest that the degeneracy at the core of bacterial stress regulation provides alternative solutions to a common evolutionary challenge. These results can explain phenotypic divergence in a constant environment and also how evolutionary jumps and adaptive radiations involve altered gene regulation.

Beneficial mutations and the dynamics of adaptation in asexual populations

We discuss the dynamics of adaptive evolution in asexual (clonal) populations. The classical 'periodic selection' model of clonal evolution assumed that beneficial mutations are very rare and therefore substitute unfettered into populations as occasional, isolated events. Newer models allow for the possibility that beneficial mutations are sufficiently common to coexist and compete for fixation within populations. Experimental evolution studies in microbes provide empirical support for stochastic models in which both selection and mutation are strong effects and clones compete for fixation; however, the relative importance of competition among clones bearing mutations of different selective effects versus competition among clones bearing multiple mutations remains unresolved. We provide some new theoretical results, moreover, suggesting that population dynamics consistent with the periodic selection model can arise even in a deterministic model that can accommodate a very high beneficial mutation rate.

Genomic identification of a novel mutation in hfq that provides multiple benefits in evolving glucose-limited populations of Escherichia coli

Beneficial mutations in diversifying glucose-limited Escherichia coli populations are mostly unidentified. The genome of an evolved isolate with multiple differences from that of the ancestor was fully assembled. Remarkably, a single mutation in hfq was responsible for the multiple benefits under glucose limitation through changes in at least five regulation targets.

Fluctuating selection: the perpetual renewal of adaptation in variable environments

Darwin insisted that evolutionary change occurs very slowly over long periods of time, and this gradualist view was accepted by his supporters and incorporated into the infinitesimal model of quantitative genetics developed by R. A. Fisher and others. It dominated the first century of evolutionary biology, but has been challenged in more recent years both by field surveys demonstrating strong selection in natural populations and by quantitative trait loci and genomic studies, indicating that adaptation is often attributable to mutations in a few genes. The prevalence of strong selection seems inconsistent, however, with the high heritability often observed in natural populations, and with the claim that the amount of morphological change in contemporary and fossil lineages is independent of elapsed time. I argue that these discrepancies are resolved by realistic accounts of environmental and evolutionary changes. First, the physical and biotic environment varies on all time-scales, leading to an indefinite increase in environmental variance over time. Secondly, the intensity and direction of natural selection are also likely to fluctuate over time, leading to an indefinite increase in phenotypic variance in any given evolving lineage. Finally, detailed long-term studies of selection in natural populations demonstrate that selection often changes in direction. I conclude that the traditional gradualist scheme of weak selection acting on polygenic variation should be supplemented by the view that adaptation is often based on oligogenic variation exposed to commonplace, strong, fluctuating natural selection.

Identification and dynamics of a beneficial mutation in a long-term evolution experiment with Escherichia coli

BACKGROUND

Twelve populations of E. coli were serially propagated for 20,000 generations in a glucose-supplemented minimal medium in order to study the dynamics of evolution. We sought to find and characterize one of the beneficial mutations responsible for the adaptation and other phenotypic changes, including increased cell size, in one of these populations.

RESULTS

We used transposon-tagging followed by P1-transduction into the ancestor, screening for increased cell size and fitness, co-transduction analysis, and DNA sequencing. We identified a 1-bp insertion in the BoxG1 region located upstream of glmUS, an operon involved in cell-wall biosynthesis. When transduced into the ancestor, this mutation increased competitive fitness by about 5%. This mutation spread through its population of origin between 500 and 1500 generations. Mutations in this region were not found in the other 11 evolving populations, even after 20,000 generations.

CONCLUSION

The 1-bp insertion in the BoxG1 region near glmUS was demonstrably beneficial in the environment in which it arose. The absence of similar mutations in the other evolved populations suggests that they substituted other mutations that rendered this particular mutation unimportant. These results show the unpredictability of adaptive evolution, whereas parallel substitutions at other loci in these same populations reveal the predictability.

E Unibus Plurum: genomic analysis of an experimentally evolved polymorphism in Escherichia coli

Microbial populations founded by a single clone and propagated under resource limitation can become polymorphic. We sought to elucidate genetic mechanisms whereby a polymorphism evolved in Escherichia coli under glucose limitation and persisted because of cross-feeding among multiple adaptive clones. Apart from a 29 kb deletion in the dominant clone, no large-scale genomic changes distinguished evolved clones from their common ancestor. Using transcriptional profiling on co-evolved clones cultured separately under glucose-limitation we identified 180 genes significantly altered in expression relative to the common ancestor grown under similar conditions. Ninety of these were similarly expressed in all clones, and many of the genes affected (e.g., mglBAC, mglD, and lamB) are in operons coordinately regulated by CRP and/or rpoS. While the remaining significant expression differences were clone-specific, 93% were exhibited by the majority clone, many of which are controlled by global regulators, CRP and CpxR. When transcriptional profiling was performed on adaptive clones cultured together, many expression differences that distinguished the majority clone cultured in isolation were absent, suggesting that CpxR may be activated by overflow metabolites removed by cross-feeding strains in co-culture. Relative to their common ancestor, shared expression differences among adaptive clones were partly attributable to early-arising shared mutations in the trans-acting global regulator, rpoS, and the cis-acting regulator, mglO. Gene expression differences that distinguished clones may in part be explained by mutations in trans-acting regulators malT and glpK, and in cis-acting sequences of acs. In the founder, a cis-regulatory mutation in acs (acetyl CoA synthetase) and a structural mutation in glpR (glycerol-3-phosphate repressor) likely favored evolution of specialists that thrive on overflow metabolites. Later-arising mutations that led to specialization emphasize the importance of compensatory rather than gain-of-function mutations in this system. Taken together, these findings underscore the importance of regulatory change, founder genotype, and the biotic environment in the adaptive evolution of microbes.

Experimental evolution of asexual species has shown that multiple genotypes can arise from a single ancestor and stably coexist (e unibus plurum). Although facilitated by environmental heterogeneity, this phenomenon also occurs in simple, homogeneous environments provisioned with a single limiting nutrient. We sought to discover genetic mechanisms that enabled an E. coli population founded by a single clone to become an interacting community composed of multiple clones. The founder of this population contained mutations that impair regulation of acetate and glycerol metabolism and likely favored the evolution of cross-feeding. Adaptive clones share cis- and trans-regulatory mutations shown elsewhere to enhance fitness under glucose limitation. Certain mutations that distinguish adaptive clones and underlie evolution of specialists were compensatory rather than gain-of-function, and all that we detected resulted in gene expression changes rather than protein structure changes. Evolved clones exhibited both common and clone-specific gene expression changes relative to their common ancestor; the pattern of gene expression in the dominant clone cultured alone differed from the pattern observed when it was cultured with variants feeding on its overflow metabolites. These findings illuminate the roles played by founder genotype, differential gene regulation, and the biotic environment in the adaptive evolution of bacteria.

Avida

Avida¹ is a software platform for experiments with self-replicating and evolving computer programs. It provides detailed control over experimental settings and protocols, a large array of measurement tools, and sophisticated methods to analyze and post-process experimental data. This chapter explains the general principles on which Avida is built, its main components and their interactions, and gives an overview of some prior research.

Perceiving molecular evolution processes in Escherichia coli by comprehensive metabolite and gene expression profiling

Transcript and metabolite abundance changes were analyzed in evolved and ancestor strains of Escherichia coli in three different evolutionary conditions

Background

Evolutionary changes that are due to different environmental conditions can be examined based on the various molecular aspects that constitute a cell, namely transcript, protein, or metabolite abundance. We analyzed changes in transcript and metabolite abundance in evolved and ancestor strains in three different evolutionary conditions - excess nutrient adaptation, prolonged stationary phase adaptation, and adaptation because of environmental shift - in two different strains of bacterium Escherichia coli K-12 (MG1655 and DH10B).

Results

Metabolite profiling of 84 identified metabolites revealed that most of the metabolites involved in the tricarboxylic acid cycle and nucleotide metabolism were altered in both of the excess nutrient evolved lines. Gene expression profiling using whole genome microarray with 4,288 open reading frames revealed over-representation of the transport functional category in all evolved lines. Excess nutrient adapted lines were found to exhibit greater degrees of positive correlation, indicating parallelism between ancestor and evolved lines, when compared with prolonged stationary phase adapted lines. Gene-metabolite correlation network analysis revealed over-representation of membrane-associated functional categories. Proteome analysis revealed the major role played by outer membrane proteins in adaptive evolution. GltB, LamB and YaeT proteins in excess nutrient lines, and FepA, CirA, OmpC and OmpA in prolonged stationary phase lines were found to be differentially over-expressed.

Conclusion

In summary, we report the vital involvement of energy metabolism and membrane-associated functional categories in all of the evolutionary conditions examined in this study within the context of transcript, outer membrane protein, and metabolite levels. These initial data obtained may help to enhance our understanding of the evolutionary process from a systems biology perspective.

The spread of a beneficial mutation in experimental bacterial populations: the influence of the environment and genotype on the fixation of rpoS mutations

The spread of beneficial mutations through populations is at the core of evolutionary change. A long-standing hindrance to understanding mutational sweeps was that beneficial mutations have been slow to be identified, even in commonly studied experimental populations. The lack of information on what constitutes a beneficial mutation has led to many uncertainties about the frequency, fitness benefit and fixation of beneficial mutations. A more complete picture is currently emerging for a limited set of identified mutations in bacterial populations. In turn, this will allow quantitation of several features of mutational sweeps. Most importantly, the 'benefit' of beneficial mutations can now be explained in terms of physiological function and how variations in the environment change the selectability of mutations. Here, the sweep of rpoS mutations in Escherichia coli, in both experimental and natural populations, is described in detail. These studies reveal the subtleties of physiology and regulation that strongly influence the benefit of a mutation and explain differences in sweeps between strains and between various environments.

Bacterial physiology, regulation and mutational adaptation in a chemostat environment

The chemostat was devised over 50 years ago and rapidly adopted for studies of bacterial physiology and mutation. Despite the long history and earlier analyses, the complexity of events in continuous cultures is only now beginning to be resolved. The application of techniques for following regulatory and mutational changes and the identification of mutated genes in chemostat populations has provided new insights into bacterial behaviour. Inoculation of bacteria into a chemostat culture results in a population competing for a limiting amount of a particular resource. Any utilizable carbon source or ion can be a limiting nutrient and bacteria respond to limitation through a regulated nutrient-specific hunger response. In addition to transcriptional responses to nutrient limitation, a second regulatory influence in a chemostat culture is the reduced growth rate fixed by the dilution rate in individual experiments. Sub-maximal growth rates and hunger result in regulation involving sigma factors and alarmones like cAMP and ppGpp. Reduced growth rate also results in increased mutation frequencies. The combination of a strongly selective environment (where mutants able to compete for limiting nutrient have a major fitness advantage) and elevated mutation rates (both endogenous and through the secondary enrichment of mutators) results in a population that changes rapidly and persistently over many generations. Contrary to common belief, the chemostat environment is never in "steady state" with fixed bacterial characteristics usable for clean comparisons of physiological or regulatory states. Adding to the complexity, chemostat populations do not simply exhibit a succession of mutational sweeps leading to a dominant winner clone. Instead, within 100 generations large populations become heterogeneous and evolving bacteria adopt alternative, parallel fitness strategies. Transport physiology, metabolism and respiration, as well as growth yields, are highly diverse in chemostat-evolved bacteria. The rich assortment of changes in an evolving chemostat provides an excellent experimental system for understanding bacterial evolution. The adaptive radiation or divergence of populations into a collection of individuals with alternative solutions to the challenge of chemostat existence provides an ideal model system for testing evolutionary and ecological theories on adaptive radiations and the generation of bacterial diversity.

Comparative genomic hybridization and physiological characterization of environmental isolates indicate that significant (eco-)physiological properties are highly conserved in the species Escherichia coli

Escherichia coli, the common inhabitant of the mammalian intestine, exhibits considerable intraspecies genomic variation, which has been suggested to reflect adaptation to different ecological niches. Also, regulatory trade-offs, e.g. between catabolic versatility and stress protection, are thought to result in significant physiological differences between strains. For these reasons, the relevance of experimental observations made for 'domesticated' E. coli strains with regard to the behaviour of this species in its natural environments is often questioned and doubts are frequently raised on the status of E. coli as a defined species. The variability of important (eco-)physiological functions, such as carbon substrate uptake and breakdown capabilities, as well as stress defence mechanisms, in the genomes of commensal and pathogenic E. coli strains were therefore investigated. Furthermore, (eco-)physiological properties of environmental strains were compared to standard laboratory strain K-12 MG1655. Catabolic, stress protection, and carbon- and energy source transport operons showed a very low intraspecies variability in 57 commensal and pathogenic E. coli. Environmental isolates adapted to glucose-limited growth in a similar way as E. coli MG1655, namely by increasing their catabolic flexibility and by inducing high-affinity substrate uptake systems. The results obtained indicate that significant (eco-)physiological properties are highly conserved in the natural population of E. coli. This questions the proposed dominant role of horizontal gene transfer for niche adaptation.

Constraints on microbial metabolism drive evolutionary diversification in homogeneous environments

Understanding the evolution of microbial diversity is an important and current problem in evolutionary ecology. In this paper, we investigated the role of two established biochemical trade-offs in microbial diversification using a model that connects ecological and evolutionary processes with fundamental aspects of biochemistry. The trade-offs that we investigated are as follows:(1) a trade-off between the rate and affinity of substrate transport; and (2) a trade-off between the rate and yield of ATP production. Our model shows that these biochemical trade-offs can drive evolutionary diversification under the simplest possible ecological conditions: a homogeneous environment containing a single limiting resource. We argue that the results of a number of microbial selection experiments are consistent with the predictions of our model.

Metabolic characterization of Escherichia coli strains adapted to growth on lactate

In comparison with intensive studies of genetic mechanisms related to biological evolutionary systems, much less analysis has been conducted on metabolic network responses to adaptive evolution that are directly associated with evolved metabolic phenotypes. Metabolic mechanisms involved in laboratory evolution of Escherichia coli on gluconeogenic carbon sources, such as lactate, were studied based on intracellular flux states determined from 13C tracer experiments and 13C-constrained flux analysis. At the end point of laboratory evolution, strains exhibited a more than doubling of the average growth rate and a 50% increase in the average biomass yield. Despite different evolutionary trajectories among parallel evolved populations, most improvements were obtained within the first 250 generations of evolution and were generally characterized by a significant increase in pathway capacity. Partitioning between gluconeogenic and pyruvate catabolic flux at the pyruvate node remained almost unchanged, while flux distributions around the key metabolites phosphoenolpyruvate, oxaloacetate, and acetyl-coenzyme A were relatively flexible over the course of evolution on lactate to meet energetic and anabolic demands during rapid growth on this gluconeogenic carbon substrate. There were no clear qualitative correlations between most transcriptional expression and metabolic flux changes, suggesting complex regulatory mechanisms at multiple levels of genetics and molecular biology. Moreover, higher fitness gains for cell growth on both evolutionary and alternative carbon sources were found for strains that adaptively evolved on gluconeogenic carbon sources compared to those that evolved on glucose. These results provide a novel systematic view of the mechanisms underlying microbial adaptation to growth on a gluconeogenic substrate.

Analysis of genetic systems using experimental evolution and whole-genome sequencing

Whole-genome sequencing of bacteria evolving in the laboratory promises to reveal the complex network of mutations that underlie adaptation.

The application of whole-genome sequencing to the study of microbial evolution promises to reveal the complex functional networks of mutations that underlie adaptation. A recent study of parallel evolution in populations of Escherichia coli shows how adaptation involves both functional changes to specific proteins as well as global changes in regulation.

Divergence and redundancy of transport and metabolic rate-yield strategies in a single Escherichia coli population

The energetic efficiency of nutrient uptake and conversion into biomass is a key factor in the ecological behavior of microorganisms. The constraints shaping the metabolic rate-yield trade-off in bacteria are not well understood. To examine whether metabolic rate-yield settings and physiological strategies evolve toward a particular optimum in a constant environment, we studied multiple Escherichia coli isolates evolving in a glucose-limited chemostat population. A major divergence in transport and metabolic strategies was observed, and the isolates included inefficient rate strategists (polluters or cheaters) and yield strategists (conservationists), as well as various hybrid rate-yield strategists and alternative ecotypes (dropouts). Sugar transport assays, strain comparisons based on metabolomics, and Biolog profiling revealed variance to the point of individuality within an evolving population. Only 68 of 177 metabolites assayed were not affected in 10 clonally related strains. The parallel enrichment of rate and yield strategists and the divergence in metabolic phylogenies indicate that bacteria do not converge on a particular rate-yield balance or unique evolutionary solutions. Redundancies in transport and metabolic pathways are proposed to have laid the framework for the multiplicity of bacterial adaptations.

High-resolution mutation mapping reveals parallel experimental evolution in yeast

Understanding the genetic basis of evolutionary adaptation is limited by our ability to efficiently identify the genomic locations of adaptive mutations. Here we describe a method that can quickly and precisely map the genetic basis of naturally and experimentally evolved complex traits using linkage analysis. A yeast strain that expresses the evolved trait is crossed to a distinct strain background and DNA from a large pool of progeny that express the trait of interest is hybridized to oligonucleotide microarrays that detect thousands of polymorphisms between the two strains. Adaptive mutations are detected by linkage to the polymorphisms from the evolved parent. We successfully tested our method by mapping five known genes to a precision of 0.2–24 kb (0.1–10 cM), and developed computer simulations to test the effect of different factors on mapping precision. We then applied this method to four yeast strains that had independently adapted to a fluctuating glucose–galactose environment. All four strains had acquired one or more missense mutations in GAL80, the repressor of the galactose utilization pathway. When transferred into the ancestral strain, the gal80 mutations conferred the fitness advantage that the evolved strains show in the transition from glucose to galactose. Our results show an example of parallel adaptation caused by mutations in the same gene.

An array hybridization method enables genetic mapping via linkage analysis; applied here this new method shows parallel adaptation to a fluctuating glucose-galactose environment evidenced by mutations in the GAL80 gene.

Global gene expression in Escherichia coli K-12 during short-term and long-term adaptation to glucose-limited continuous culture conditions

Microarray technology was used to study the cellular events that take place at the transcription level during short-term (physiological) and long-term (genetic) adaptation of the faecal indicator bacterium Escherichia coli K-12 to slow growth under limited nutrient supply. Short-term and long-term adaptation were assessed by comparing the mRNA levels isolated after 40 or 500 h of glucose-limited continuous culture at a dilution rate of 0.3 h(-1) with those from batch culture with glucose excess. A large number of genes encoding periplasmic binding proteins were upregulated, indicating that the cells are prepared for high-affinity uptake of all types of carbon sources during glucose-limited growth in continuous culture. All the genes belonging to the maltose (mal/lamB) and galactose (mgl/gal) operons were upregulated. A similar transcription pattern was observed for long-term cultures except that the expression factors were lower than in the short-term adaptation. The patterns of upregulation were confirmed by real-time RT-PCR. A switch from a fully operational citric acid cycle to the PEP-glyoxylate cycle was clearly observed in cells grown in glucose-limited continuous culture when compared to batch-grown cells and this was confirmed by transcriptome analysis. This transcriptome analysis confirms and extends the observations from previous proteome and catabolome studies in the authors' laboratory.

Adaptive radiation in microbial microcosms

It has often been argued that evolutionary diversification is the result of divergent natural selection for specialization on alternative resources. I provide a comprehensive review of experiments that examine the ecology and genetics of resource specialization and adaptive radiation in microbial microcosms. In these experiments, resource heterogeneity generates divergent selection for specialization on alternative resources. At a molecular level, the evolution of specialization is generally attributable to mutations that de-regulate the expression of existing biosynthetic and catabolic pathways. Trade-offs are associated with the evolution of resource specialization, but these trade-offs are often not the result of antagonistic pleiotropy. Replicate adaptive radiations result in the evolution of a similar assemblage of specialists, but the genetic basis of specialization differs in replicate radiations. The implications of microbial selection experiments for evolutionary theory are discussed and future directions of research are proposed.

Nutrient-scavenging stress response in an Escherichia coli strain lacking the phosphoenolpyruvate: carbohydrate phosphotransferase system, as explored by gene expression profile analysis

The physiological role of the phosphoenolpyruvate:carbohydrate phosphotransferase system (PTS) has been studied in Escherichia coli. It has been shown that it directly or indirectly regulates the activity of most catabolic genes involved in carbohydrate transport. Accordingly, strains lacking PTS have pleiotropic phenotypes and are impaired in their capacity to grow on glucose and other PTS sugars. We have previously reported the characterization of a mutant harboring a pts operon deletion (PB11) which, as expected, showed a severe reduction of its growth capacity when incubated on glucose as carbon source, as compared to that of the isogenic wild-type strain. These observations corroborate that PTS is the main determinant of the capacity to grow on glucose and confirm the existence of other systems that allow glucose utilization although at a reduced level. To explore the physiological state and the metabolic pathways involved in glucose utilization in a pts(-) background, we analyzed the global transcriptional response of the PB11 mutant when growing in minimal medium with glucose as carbon source. Genome-wide transcriptional analysis using microarrays revealed that, under this condition, expression of several genes related to carbon transport and metabolism was upregulated, as well as that of genes encoding transporters for certain nucleotides, nitrogen, phosphorus and sulfur sources. In addition, upregulation of rpoS and several genes transcribed by this sigma subunit was detected. These results indicate that the reduced capacity of glucose utilization present in the PB11 strain induces a general nutrient-scavenging response and this behavior is not dependent on a functional PTS. This condition is responsible of the utilization of secondary carbon sources in the presence of glucose.

Towards multidimensional genome annotation

Our information about the gene content of organisms continues to grow as more genomes are sequenced and gene products are characterized. Sequence-based annotation efforts have led to a list of cellular components, which can be thought of as a one-dimensional annotation. With growing information about component interactions, facilitated by the advancement of various high-throughput technologies, systemic, or two-dimensional, annotations can be generated. Knowledge about the physical arrangement of chromosomes will lead to a three-dimensional spatial annotation of the genome and a fourth dimension of annotation will arise from the study of changes in genome sequences that occur during adaptive evolution. Here we discuss all four levels of genome annotation, with specific emphasis on two-dimensional annotation methods.

Changes in the metabolome of Saccharomyces cerevisiae associated with evolution in aerobic glucose-limited chemostats

The effect of culture age on intra- and extracellular metabolite levels as well as on in vitro determined specific activities of enzymes of central carbon metabolism was investigated during evolution for over 90 generations of Saccharomyces cerevisiae CEN.PK 113-7D in an aerobic glucose/ethanol-limited chemostat at a specific dilution rate of 0.052 h(-1). It was found that the fluxes of consumed (O2, glucose/ethanol) and secreted compounds (CO2) did not change significantly during the entire cultivation period. However, morphological changes were observed, leading to an increased cellular surface area. During 90 generations of chemostat growth not only the residual glucose concentration decreased, also the intracellular concentrations of trehalose, glycolytic intermediates, TCA cycle intermediates and amino acids were found to have decreased with a factor 5-10. The only exception was glyoxylate which showed a fivefold increase in concentration. In addition to this the specific activities of most glycolytic enzymes also decreased by a factor 5-10 during long-term cultivation. Exceptions to this were hexokinase, phosphofructokinase, pyruvate kinase and 6-phosphogluconate dehydrogenase of which the activities remained unchanged. Furthermore, the concentrations of the adenylate nucleotides as well as the energy charge of the cells did not change in a significant manner. Surprisingly, the specific activities of glucose-6-phosphate dehydrogenase (G6PDH), malate synthase (MS) and isocitrate lyase (ICL) increased significantly during 90 generations of chemostat cultivation. These changes seem to indicate a pattern where metabolic overcapacities (for reversible reactions) and storage pools (trehalose, high levels of amino acids and excess protein in enzymes) are lost during the evolution period. The driving force is proposed to be a growth advantage in the absence of these metabolic overcapacities.

Homeostatic adjustment and metabolic remodeling in glucose-limited yeast cultures

We studied the physiological response to glucose limitation in batch and steady-state (chemostat) cultures of Saccharomyces cerevisiae by following global patterns of gene expression. Glucose-limited batch cultures of yeast go through two sequential exponential growth phases, beginning with a largely fermentative phase, followed by an essentially completely aerobic use of residual glucose and evolved ethanol. Judging from the patterns of gene expression, the state of the cells growing at steady state in glucose-limited chemostats corresponds most closely with the state of cells in batch cultures just before they undergo this "diauxic shift." Essentially the same pattern was found between chemostats having a fivefold difference in steady-state growth rate (the lower rate approximating that of the second phase respiratory growth rate in batch cultures). Although in both cases the cells in the chemostat consumed most of the glucose, in neither case did they seem to be metabolizing it primarily through respiration. Although there was some indication of a modest oxidative stress response, the chemostat cultures did not exhibit the massive environmental stress response associated with starvation that also is observed, at least in part, during the diauxic shift in batch cultures. We conclude that despite the theoretical possibility of a switch to fully aerobic metabolism of glucose in the chemostat under conditions of glucose scarcity, homeostatic mechanisms are able to carry out metabolic adjustment as if fermentation of the glucose is the preferred option until the glucose is entirely depleted. These results suggest that some aspect of actual starvation, possibly a component of the stress response, may be required for triggering the metabolic remodeling associated with the diauxic shift.

The role of isocitrate lyase and the glyoxylate cycle in Escherichia coli growing under glucose limitation

Escherichia coli changes its metabolism in response to environmental circumstances, and metabolic adaptations are evident in hungry bacteria growing slowly in glucose-limited chemostats. The role of isocitrate lyase (AceA) was examined in E. coli growing under glucose limitation. AceA activity was elevated in a strain-dependent manner in the commonly used E. coli K-12 laboratory strains MG1655 and MC4100, but an aceA disruption surprisingly increased fitness under glucose limitation in both strains. However, in bacteria adapted to limiting glucose in long-term chemostats, mutations outside aceA changed its role from a negative to a positive influence. These results suggest that a recently proposed pathway of central metabolism involving the glyoxylate cycle enzymes is redundant in wild-type bacteria, but may take on a beneficial role after context adaptation. Interestingly, the aceA gene sequence did not alter during prolonged selection, so mutations in unidentified genes changed the metabolic context of unaltered AceA from a negative to a positive influence in bacteria highly adapted to limiting glucose.