Adaptation and heterogeneity of Escherichia coli MC1000 growing in complex environments

Bacteria evolve macroscopic multicellularity by the genetic assimilation of phenotypically plastic cell clustering

The evolutionary transition from unicellularity to multicellularity was a key innovation in the history of life. Experimental evolution is an important tool to study the formation of undifferentiated cellular clusters, the likely first step of this transition. Although multicellularity first evolved in bacteria, previous experimental evolution research has primarily used eukaryotes. Moreover, it focuses on mutationally driven (and not environmentally induced) phenotypes. Here we show that both Gram-negative and Gram-positive bacteria exhibit phenotypically plastic (i.e., environmentally induced) cell clustering. Under high salinity, they form elongated clusters of ~ 2 cm. However, under habitual salinity, the clusters disintegrate and grow planktonically. We used experimental evolution with Escherichia coli to show that such clustering can be assimilated genetically: the evolved bacteria inherently grow as macroscopic multicellular clusters, even without environmental induction. Highly parallel mutations in genes linked to cell wall assembly formed the genomic basis of assimilated multicellularity. While the wildtype also showed cell shape plasticity across high versus low salinity, it was either assimilated or reversed after evolution. Interestingly, a single mutation could genetically assimilate multicellularity by modulating plasticity at multiple levels of organization. Taken together, we show that phenotypic plasticity can prime bacteria for evolving undifferentiated macroscopic multicellularity.

Habitat heterogeneity induced by pyrogenic organic matter in wildfire-perturbed soils mediates bacterial community assembly processes

Although pyrogenic organic matter (PyOM) generated during wildfires plays a critical role in post-fire ecosystem recovery, the specific mechanisms by which PyOM controls soil microbial community assembly after wildfire perturbation remain largely uncharacterized. Herein we characterized the effect of PyOM on soil bacterial communities at two independent wildfire-perturbed forest sites. We observed that α-diversity of bacterial communities was the highest in wildfire-perturbed soils and that bacterial communities gradually changed along a sequence of unburnt soil → burnt soil → PyOM. The microbial communities reconstructed from unburnt soil and PyOM resembled the real bacterial communities in wildfire-perturbed soils in their α-diversity and community structure. Bacterial specialists in PyOM and soils clustered in phylogenetic coherent lineages with intra-lineage pH-niche conservatism and inter-lineage pH-niche divergence. Our results suggest that PyOM mediates bacterial community assembly in wildfire-perturbed soils by a combination of environmental selection and dispersal of phylogenetic coherent specialists with habitat preference in the heterogeneous microhabitats of burnt soils with distinct PyOM patches.

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.

Escherichia coli cultures maintain stable subpopulation structure during long-term evolution

How genetic variation is generated and maintained remains a central question in evolutionary biology. When presented with a complex environment, microbes can take advantage of genetic variation to exploit new niches. Here we present a massively parallel experiment where WT and repair-deficient (∆mutL) Escherichia coli populations have evolved over 3 y in a spatially heterogeneous and nutritionally complex environment. Metagenomic sequencing revealed that these initially isogenic populations evolved and maintained stable subpopulation structure in just 10 mL of medium for up to 10,000 generations, consisting of up to five major haplotypes with many minor haplotypes. We characterized the genomic, transcriptomic, exometabolomic, and phenotypic differences between clonal isolates, revealing subpopulation structure driven primarily by spatial segregation followed by differential utilization of nutrients. In addition to genes regulating the import and catabolism of nutrients, major polymorphisms of note included insertion elements transposing into fimE (regulator of the type I fimbriae) and upstream of hns (global regulator of environmental-change and stress-response genes), both known to regulate biofilm formation. Interestingly, these genes have also been identified as critical to colonization in uropathogenic E. coli infections. Our findings illustrate the complexity that can arise and persist even in small cultures, raising the possibility that infections may often be promoted by an evolving and complex pathogen population.

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.

Adaptation of Escherichia coli to long-term batch culture in various rich media

Experimental evolution studies have characterized the genetic strategies microbes utilize to adapt to their environments, mainly focusing on how microbes adapt to constant and/or defined environments. Using a system that incubates Escherichia coli in different complex media in long-term batch culture, we have focused on how heterogeneity and environment affects adaptive landscapes. In this system, there is no passaging of cells, and therefore genetic diversity is lost only through negative selection, without the experimentally-imposed bottlenecking common in other platforms. In contrast with other experimental evolution systems, because of cycling of nutrients and waste products, this is a heterogeneous environment, where selective pressures change over time, similar to natural environments. We determined that incubation in each environment leads to different adaptations by observing the growth advantage in stationary phase (GASP) phenotype. Re-sequencing whole genomes of populations identified both mutant alleles in a conserved set of genes and differences in evolutionary trajectories between environments. Reconstructing identified mutations in the parental strain background confirmed the adaptive advantage of some alleles, but also identified a surprising number of neutral or even deleterious mutations. This result indicates that complex epistatic interactions may be under positive selection within these heterogeneous environments.

Dynamics and genetic diversification of Escherichia coli during experimental adaptation to an anaerobic environment

Background

Many bacteria are facultative anaerobes, and can proliferate in both anoxic and oxic environments. Under anaerobic conditions, fermentation is the primary means of energy generation in contrast to respiration. Furthermore, the rates and spectra of spontaneous mutations that arise during anaerobic growth differ to those under aerobic growth. A long-term selection experiment was undertaken to investigate the genetic changes that underpin how the facultative anaerobe, Escherichia coli, adapts to anaerobic environments.

Methods

Twenty-one populations of E. coli REL4536, an aerobically evolved 10,000th generation descendent of the E. coli B strain, REL606, were established from a clonal ancestral culture. These were serially sub-cultured for 2,000 generations in a defined minimal glucose medium in strict aerobic and strict anaerobic environments, as well as in a treatment that fluctuated between the two environments. The competitive fitness of the evolving lineages was assessed at approximately 0, 1,000 and 2,000 generations, in both the environment of selection and the alternative environment. Whole genome re-sequencing was performed on random colonies from all lineages after 2,000-generations. Mutations were identified relative to the ancestral genome, and based on the extent of parallelism, traits that were likely to have contributed towards adaptation were inferred.

Results

There were increases in fitness relative to the ancestor among anaerobically evolved lineages when tested in the anaerobic environment, but no increases were found in the aerobic environment. For lineages that had evolved under the fluctuating regime, relative fitness increased significantly in the anaerobic environment, but did not increase in the aerobic environment. The aerobically-evolved lineages did not increase in fitness when tested in either the aerobic or anaerobic environments. The strictly anaerobic lineages adapted more rapidly to the anaerobic environment than did the fluctuating lineages. Two main strategies appeared to predominate during adaptation to the anaerobic environment: modification of energy generation pathways, and inactivation of non-essential functions. Fermentation pathways appeared to alter through selection for mutations in genes such as nadR, adhE, dcuS/R, and pflB. Mutations were frequently identified in genes for presumably dispensable functions such as toxin-antitoxin systems, prophages, virulence and amino acid transport. Adaptation of the fluctuating lineages to the anaerobic environments involved mutations affecting traits similar to those observed in the anaerobically evolved lineages.

Discussion

There appeared to be strong selective pressure for activities that conferred cell yield advantages during anaerobic growth, which include restoring activities that had previously been inactivated under long-term continuous aerobic evolution of the ancestor.

Escherichia coli populations adapt to complex, unpredictable fluctuations by minimizing trade-offs across environments

In nature, organisms are simultaneously exposed to multiple stresses (i.e. complex environments) that often fluctuate unpredictably. Although both these factors have been studied in isolation, the interaction of the two remains poorly explored. To address this issue, we selected laboratory populations of Escherichia coli under complex (i.e. stressful combinations of pH, H₂ O₂ and NaCl) unpredictably fluctuating environments for ~900 generations. We compared the growth rates and the corresponding trade-off patterns of these populations to those that were selected under constant values of the component stresses (i.e. pH, H₂ O₂ and NaCl) for the same duration. The fluctuation-selected populations had greater mean growth rate and lower variation for growth rate over all the selection environments experienced. However, whereas the populations selected under constant stresses experienced trade-offs in the environments other than those in which they were selected, the fluctuation-selected populations could bypass the across-environment trade-offs almost entirely. Interestingly, trade-offs were found between growth rates and carrying capacities. The results suggest that complexity and fluctuations can strongly affect the underlying trade-off structure in evolving populations.

Escherichia coli populations in unpredictably fluctuating environments evolve to face novel stresses through enhanced efflux activity

There is considerable understanding about how laboratory populations respond to predictable (constant or deteriorating environment) selection for single environmental variables such as temperature or pH. However, such insights may not apply when selection environments comprise multiple variables that fluctuate unpredictably, as is common in nature. To address this issue, we grew replicate laboratory populations of Escherichia coli in nutrient broth whose pH and concentrations of salt (NaCl) and hydrogen peroxide (H2 O2 ) were randomly changed daily. After ~170 generations, the fitness of the selected populations had not increased in any of the three selection environments. However, these selected populations had significantly greater fitness in four novel environments which have no known fitness-correlation with tolerance to pH, NaCl or H2 O2 . Interestingly, contrary to expectations, hypermutators did not evolve. Instead, the selected populations evolved an increased ability for energy-dependent efflux activity that might enable them to throw out toxins, including antibiotics, from the cell at a faster rate. This provides an alternate mechanism for how evolvability can evolve in bacteria and potentially lead to broad-spectrum antibiotic resistance, even in the absence of prior antibiotic exposure. Given that environmental variability is increasing in nature, this might have serious consequences for public health.

Acid-adapted strains of Escherichia coli K-12 obtained by experimental evolution

Enteric bacteria encounter a wide range of pHs throughout the human intestinal tract. We conducted experimental evolution of Escherichia coli K-12 to isolate clones with increased fitness during growth under acidic conditions (pH 4.5 to 4.8). Twenty-four independent populations of E. coli K-12 W3110 were evolved in LBK medium (10 g/liter tryptone, 5 g/liter yeast extract, 7.45 g/liter KCl) buffered with homopiperazine-N,N'-bis-2-(ethanosulfonic acid) and malate at pH 4.8. At generation 730, the pH was decreased to 4.6 with HCl. By 2,000 generations, all populations had achieved higher endpoint growth than the ancestor at pH 4.6 but not at pH 7.0. All evolving populations showed a progressive loss of activity of lysine decarboxylase (CadA), a major acid stress enzyme. This finding suggests a surprising association between acid adaptation and moderation of an acid stress response. At generation 2,000, eight clones were isolated from four populations, and their genomes were sequenced. Each clone showed between three and eight missense mutations, including one in a subunit of the RNA polymerase holoenzyme (rpoB, rpoC, or rpoD). Missense mutations were found in adiY, the activator of the acid-inducible arginine decarboxylase (adiA), and in gcvP (glycine decarboxylase), a possible acid stress component. For tests of fitness relative to that of the ancestor, lacZ::kan was transduced into each strain. All acid-evolved clones showed a high fitness advantage at pH 4.6. With the cytoplasmic pH depressed by benzoate (at external pH 6.5), acid-evolved clones showed decreased fitness; thus, there was no adaptation to cytoplasmic pH depression. At pH 9.0, acid-evolved clones showed no fitness advantage. Thus, our acid-evolved clones showed a fitness increase specific to low external pH.

Mutations in global regulators lead to metabolic selection during adaptation to complex environments

Adaptation to ecologically complex environments can provide insights into the evolutionary dynamics and functional constraints encountered by organisms during natural selection. Adaptation to a new environment with abundant and varied resources can be difficult to achieve by small incremental changes if many mutations are required to achieve even modest gains in fitness. Since changing complex environments are quite common in nature, we investigated how such an epistatic bottleneck can be avoided to allow rapid adaptation. We show that adaptive mutations arise repeatedly in independently evolved populations in the context of greatly increased genetic and phenotypic diversity. We go on to show that weak selection requiring substantial metabolic reprogramming can be readily achieved by mutations in the global response regulator arcA and the stress response regulator rpoS. We identified 46 unique single-nucleotide variants of arcA and 18 mutations in rpoS, nine of which resulted in stop codons or large deletions, suggesting that subtle modulations of ArcA function and knockouts of rpoS are largely responsible for the metabolic shifts leading to adaptation. These mutations allow a higher order metabolic selection that eliminates epistatic bottlenecks, which could occur when many changes would be required. Proteomic and carbohydrate analysis of adapting E. coli populations revealed an up-regulation of enzymes associated with the TCA cycle and amino acid metabolism, and an increase in the secretion of putrescine. The overall effect of adaptation across populations is to redirect and efficiently utilize uptake and catabolism of abundant amino acids. Concomitantly, there is a pronounced spread of more ecologically limited strains that results from specialization through metabolic erosion. Remarkably, the global regulators arcA and rpoS can provide a “one-step” mechanism of adaptation to a novel environment, which highlights the importance of global resource management as a powerful strategy to adaptation.

Changing environmental conditions are the norm in biology. However, understanding adaptation to complex environments presents many challenges. For example, adaptation to resource-rich environments can potentially have many successful evolutionary trajectories to increased fitness. Even in conditions of plenty, the utilization of numerous but novel resources can require multiple mutations before a benefit is accrued. We evolved two bacterial species isolated from the gut of healthy humans in two different, resource-rich media commonly used in the laboratory. We anticipated that under weak selection the population would evolve tremendous genetic diversity. Despite such a complex genetic background we were able to identify a strong degree of parallel evolution and using a combination of population proteomic and population genomic approaches we show that two global regulators, arcA and rpoS, are the principle targets of selection. Up-regulation of the different metabolic pathways that are controlled by these global regulators in combination with up-regulation of transporters that transport nutrients into the cell revealed increased use of the novel resources. Thus global regulators can provide a one-step model to shift metabolism efficiently and provide rapid a one-step reprogramming of the cell metabolic profile.

Plasticity predicts evolution in a marine alga

Under global change, populations have four possible responses: ‘migrate, acclimate, adapt or die’ (Gienapp et al. 2008 Climate change and evolution: disentangling environmental and genetic response. Mol. Ecol.17, 167–178. (doi:10.1111/j.1365-294X.2007.03413.x)). The challenge is to predict how much migration, acclimatization or adaptation populations are capable of. We have previously shown that populations from more variable environments are more plastic (Schaum et al. 2013 Variation in plastic responses of a globally distributed picoplankton species to ocean acidification. Nature3, 298–230. (doi:10.1038/nclimate1774)), and here we use experimental evolution with a marine microbe to learn that plastic responses predict the extent of adaptation in the face of elevated partial pressure of CO₂ (pCO₂). Specifically, plastic populations evolve more, and plastic responses in traits other than growth can predict changes in growth in a marine microbe. The relationship between plasticity and evolution is strongest when populations evolve in fluctuating environments, which favour the evolution and maintenance of plasticity. Strikingly, plasticity predicts the extent, but not direction of phenotypic evolution. The plastic response to elevated pCO₂ in green algae is to increase cell division rates, but the evolutionary response here is to decrease cell division rates over 400 generations until cells are dividing at the same rate their ancestors did in ambient CO₂. Slow-growing cells have higher mitochondrial potential and withstand further environmental change better than faster growing cells. Based on this, we hypothesize that slow growth is adaptive under CO₂ enrichment when associated with the production of higher quality daughter cells.

Differential stress resistance and metabolic traits underlie coexistence in a sympatrically evolved bacterial population

Following intermittent batch growth in Luria-Bertani (LB) broth for about 1000 generations, differentially evolved forms were found in a population of Escherichia coli cells. Studies on this population revealed the emergence of key polymorphisms, as evidenced by analysis of both whole genome sequences and transcription analysis. Here, we investigated the phenotypic nature of several key forms and found a remarkable (interactive) coexistence of forms which highlights the presence of different ecological roles pointing at a dichotomy in: (i) tolerance to environmental stresses and (ii) the capacity to utilize particular carbon sources such as galactose. Both forms differed from their common ancestor by different criteria. This apparent coexistence of two diverged forms points at the occurrence of niche partitioning as a consequence of dichotomous adaptive evolution. Remarkably, the two forms were shown to continue to coexist - in varying ratio's - in an experiment that cycled them through periods of nutrient feast (plentiful growth substrates) and famine (growth-restrictive - stress conditions). The results further indicated that the equilibrium of the coexistence was destroyed when one of the parameters was high tuned, jeopardizing the stability of the coexisting pair.

Comparative genomics and transcriptomics analysis of experimentally evolved Escherichia coli MC1000 in complex environments

It has recently become feasible to study the basis and nature of evolutionary changes in bacteria in an experimental setting using defined media. However, assessment of adaptive changes in complex environments has been scarce. In an effort to describe the responses in such environments, we unravel, in a comparative approach, the transcriptional and genetic profiles of 19 Escherichia coli strains that evolved in Luria Bertani medium under three different oxygen regimes over 1000 generations. A positive relationship between upregulation of gene expression and the number of mutations was observed, suggesting that a number of metabolic pathways were activated. Phenotypic polymorphisms were observed in parallel cultures, of which some were related with mutations at the regulatory level. Non-parallel responses were observed at the intrapopulational level, which is indicative of diversifying selection. Parallel responses encompassed transcriptome diversity, and their effects were directly affected by differing genomic backgrounds. A fluctuating selective force produced higher phenotypic diversity compared with constant forces. This study demonstrates how phenotypic innovations may depend on the relationship between genomic changes and local ecological conditions. Using both comparative genomics and transcriptomics approaches, the results help elucidating various adaptive responses in cultures in unexplored complex environments.