Transition from positive to neutral in mutation fixation along with continuing rising fitness in thermal adaptive evolution

Evolution of pH-sensitive transcription termination during adaptation to repeated long-term starvation

Fluctuating environments that consist of regular cycles of co-occurring stress are a common challenge faced by cellular populations. For a population to thrive in constantly changing conditions, an ability to coordinate a rapid cellular response is essential. Here, we identify a mutation conferring an arginine-to-histidine (Arg to His) substitution in the transcription terminator Rho. The rho R109H mutation frequently arose in E. coli populations experimentally evolved under repeated long-term starvation conditions, during which feast and famine result in drastic environmental pH fluctuations. Metagenomic sequencing revealed that populations containing the rho mutation also possess putative loss-of-function mutations in ydcI, which encodes a recently characterized transcription factor associated with pH homeostasis. Genetic reconstructions of these mutations show that the rho allele confers a plastic alkaline-induced reduction of Rho function that, when found in tandem with a ΔydcI allele, leads to intracellular alkalinization and genetic assimilation of Rho mutant function. We further identify Arg to His substitutions at analogous sites in rho alleles from species originating from fluctuating alkaline environments. Our results suggest that Arg to His substitutions in global regulators of gene expression can serve to rapidly coordinate complex responses through pH sensing and shed light on how cellular populations across the tree of life use environmental cues to coordinate rapid responses to complex, fluctuating environments.

Enhanced Cell Wall and Cell Membrane Activity Promotes Heat Adaptation of Enterococcus faecium

Enterococcus faecium (E. faecium) is widely used in foods and is known as a probiotic to treat or prevent diarrhea in pets and livestock. However, the poor resistance of E. faecium to high temperature processing procedures limits its use. Strain domestication is a low-cost and effective method to obtain high-temperature-resistant strains. In this study, heat treatment was performed from 45 °C to 70 °C and the temperature was gradually increased by 5 °C every 3 days. After domestication, the survival rates of the high temperature adaptation strain RS047-wl under 65 °C water bath for 40 min was 11.5 times higher than WT RS047. Moreover, the saturated fatty acid (SFA) contents in cell membrane and the cell volume significantly increased in the RS047-wl. The combined transcriptomic, metabolomic, and proteomics analysis results showed a significant enhancement of cell wall and membrane synthesis ability in the RS047-wl. In conclusion, one of the main factors contributing to the improved high temperature resistance of RS047-wl was its enhanced ability to synthesize cell wall and membrane, which helped maintain normal cell morphology. Developing a high-temperature-resistant strain and understanding its mechanism enables it to adapt to high temperatures. This lays the groundwork for its future development and application.

The Potassium Efflux System Kef: Bacterial Protection against Toxic Electrophilic Compounds

Kef couples the potassium efflux with proton influx in gram-negative bacteria. The resulting acidification of the cytosol efficiently prevents the killing of the bacteria by reactive electrophilic compounds. While other degradation pathways for electrophiles exist, Kef is a short-term response that is crucial for survival. It requires tight regulation since its activation comes with the burden of disturbed homeostasis. Electrophiles, entering the cell, react spontaneously or catalytically with glutathione, which is present at high concentrations in the cytosol. The resulting glutathione conjugates bind to the cytosolic regulatory domain of Kef and trigger activation while the binding of glutathione keeps the system closed. Furthermore, nucleotides can bind to this domain for stabilization or inhibition. The binding of an additional ancillary subunit, called KefF or KefG, to the cytosolic domain is required for full activation. The regulatory domain is termed K⁺ transport-nucleotide binding (KTN) or regulator of potassium conductance (RCK) domain, and it is also found in potassium uptake systems or channels in other oligomeric arrangements. Bacterial RosB-like transporters and K⁺ efflux antiporters (KEA) of plants are homologs of Kef but fulfill different functions. In summary, Kef provides an interesting and well-studied example of a highly regulated bacterial transport system.

Whole-genome sequencing analysis of two heat-evolved Escherichia coli strains

BACKGROUND

High temperatures cause a suite of problems for cells, including protein unfolding and aggregation; increased membrane fluidity; and changes in DNA supercoiling, RNA stability, transcription and translation. Consequently, enhanced thermotolerance can evolve through an unknown number of genetic mechanisms even in the simple model bacterium Escherichia coli. To date, each E. coli study exploring this question resulted in a different set of mutations. To understand the changes that can arise when an organism evolves to grow at higher temperatures, we sequenced and analyzed two previously described E. coli strains, BM28 and BM28 ΔlysU, that have been laboratory adapted to the highest E. coli growth temperature reported to date.

RESULTS

We found three large deletions in the BM28 and BM28 ΔlysU strains of 123, 15 and 8.5 kb in length and an expansion of IS10 elements. We found that BM28 and BM28 ΔlysU have considerably different genomes, suggesting that the BM28 culture that gave rise to BM28 and BM28 ΔlysU was a mixed population of genetically different cells. Consistent with published findings of high GroESL expression in BM28, we found that BM28 inexplicitly carries the groESL bearing plasmid pOF39 that was maintained simply by high-temperature selection pressure. We identified over 200 smaller insertions, deletions, single nucleotide polymorphisms and other mutations, including changes in master regulators such as the RNA polymerase and the transcriptional termination factor Rho. Importantly, this genome analysis demonstrates that the commonly cited findings that LysU plays a crucial role in thermotolerance and that GroESL hyper-expression is brought about by chromosomal mutations are based on a previous misinterpretation of the genotype of BM28.

CONCLUSIONS

This whole-genome sequencing study describes genetically distinct mechanisms of thermotolerance evolution from those found in other heat-evolved E. coli strains. Studying adaptive laboratory evolution to heat in simple model organisms is important in the context of climate change. It is important to better understand genetic mechanisms of enhancing thermotolerance in bacteria and other organisms, both in terms of optimizing laboratory evolution methods for various organisms and in terms of potential genetic engineering of organisms most at risk or most important to our societies and ecosystems.

Latent functional diversity may accelerate microbial community responses to temperature fluctuations

How complex microbial communities respond to climatic fluctuations remains an open question. Due to their relatively short generation times and high functional diversity, microbial populations harbor great potential to respond as a community through a combination of strain-level phenotypic plasticity, adaptation, and species sorting. However, the relative importance of these mechanisms remains unclear. We conducted a laboratory experiment to investigate the degree to which bacterial communities can respond to changes in environmental temperature through a combination of phenotypic plasticity and species sorting alone. We grew replicate soil communities from a single location at six temperatures between 4°C and 50°C. We found that phylogenetically and functionally distinct communities emerge at each of these temperatures, with K-strategist taxa favored under cooler conditions and r-strategist taxa under warmer conditions. We show that this dynamic emergence of distinct communities across a wide range of temperatures (in essence, community-level adaptation) is driven by the resuscitation of latent functional diversity: the parent community harbors multiple strains pre-adapted to different temperatures that are able to ‘switch on’ at their preferred temperature without immigration or adaptation. Our findings suggest that microbial community function in nature is likely to respond rapidly to climatic temperature fluctuations through shifts in species composition by resuscitation of latent functional diversity.

Global coordination of the mutation and growth rates across the genetic and nutritional variety in Escherichia coli

Fitness and mutability are the primary traits of living organisms for adaptation and evolution. However, their quantitative linkage remained largely deficient. Whether there is any general relationship between the two features and how genetic and environmental variables influence them remained unclear and were addressed here. The mutation and growth rates of an assortment of Escherichia coli strain collections, including the wild-type strains and the genetically disturbed strains of either reduced genomes or deletion of the genes involved in the DNA replication fidelity, were evaluated in various media. The contribution of media to the mutation and growth rates was differentiated depending on the types of genetic disturbance. Nevertheless, the negative correlation between the mutation and growth rates was observed across the genotypes and was common in all media. It indicated the comprehensive association of the correlated mutation and growth rates with the genetic and medium variation. Multiple linear regression and support vector machine successfully predicted the mutation and growth rates and the categories of genotypes and media, respectively. Taken together, the study provided a quantitative dataset linking the mutation and growth rates, genotype, and medium and presented a simple and successful example of predicting bacterial growth and mutability by data-driven approaches.

Emergence of Resistant Escherichia coli Mutants in Microfluidic On-Chip Antibiotic Gradients

Spatiotemporal structures and heterogeneities are common in natural habitats, yet their role in the evolution of antibiotic resistance is still to be uncovered. We applied a microfluidic gradient generator device to study the emergence of resistant bacteria in spatial ciprofloxacin gradients. We observed biofilm formation in regions with sub-inhibitory concentrations of antibiotics, which quickly expanded into the high antibiotic regions. In the absence of an explicit structure of the habitat, this multicellular formation led to a spatial structure of the population with local competition and limited migration. Therefore, such structures can function as amplifiers of selection and aid the spread of beneficial mutations. We found that the physical environment itself induces stress-related mutations that later prove beneficial when cells are exposed to antibiotics. This shift in function suggests that exaptation occurs in such experimental scenarios. The above two processes pave the way for the subsequent emergence of highly resistant specific mutations.

Experimental Evolution Expands the Breadth of Adaptation to an Environmental Gradient Correlated With Genome Reduction

Whether and how adaptive evolution adjusts the breadth of adaptation in coordination with the genome are essential issues for connecting evolution with ecology. To address these questions, experimental evolution in five Escherichia coli strains carrying either the wild-type genome or a reduced genome was performed in a defined minimal medium (C0). The ancestral and evolved populations were subsequently subjected to fitness and chemical niche analyses across an environmental gradient with 29 combinations of eight chemical components of the minimal medium. The results showed that adaptation was achieved not only specific to the evolutionary condition (C0), but also generally, to the environmental gradient; that is, the breadth of adaptation to the eight chemical niches was expanded. The magnitudes of the adaptive improvement and the breadth increase were both correlated with genome reduction and were highly significant in two out of eight niches (i.e., glucose and sulfate). The direct adaptation-induced correlated adaptation to the environmental gradient was determined by only a few genome mutations. An additive increase in fitness associated with the stepwise fixation of mutations was consistently observed in the reduced genomes. In summary, this preliminary survey demonstrated that evolution finely tuned the breadth of adaptation correlated with genome reduction.

Primordial mimicry induces morphological change in Escherichia coli

The morphology of primitive cells has been the subject of extensive research. A spherical form was commonly presumed in prebiotic studies but lacked experimental evidence in living cells. Whether and how the shape of living cells changed are unclear. Here we exposed the rod-shaped bacterium Escherichia coli to a resource utilization regime mimicking a primordial environment. Oleate was given as an easy-to-use model prebiotic nutrient, as fatty acid vesicles were likely present on the prebiotic Earth and might have been used as an energy resource. Six evolutionary lineages were generated under glucose-free but oleic acid vesicle (OAV)-rich conditions. Intriguingly, fitness increase was commonly associated with the morphological change from rod to sphere and the decreases in both the size and the area-to-volume ratio of the cell. The changed cell shape was conserved in either OAVs or glucose, regardless of the trade-offs in carbon utilization and protein abundance. Highly differentiated mutations present in the genome revealed two distinct strategies of adaption to OAV-rich conditions, i.e., either directly targeting the cell wall or not. The change in cell morphology of Escherichia coli for adapting to fatty acid availability supports the assumption of the primitive spherical form.

Lu et al. investigate the evolution of the shape of living cells by generating six experimental lineages of the rod-shaped E. coli under glucose-free conditions in the presence of oleic acid mimicking a primordial environment. The authors show that the morphological changes from rod to sphere accompanied fitness increases and adaptation amongst fatty acid availability supports the assumption of a primitive spherical form.

Emerging Adaptive Strategies Under Temperature Fluctuations in a Laboratory Evolution Experiment of Escherichia Coli

Generalists and specialists are types of strategies individuals can employ that can evolve in fluctuating environments depending on the extremity and periodicity of the fluctuation. To evaluate whether the evolution of specialists or generalists occurs under environmental fluctuation regimes with different levels of periodicity, 24 populations of Escherichia coli underwent laboratory evolution with temperatures alternating between 15 and 43°C in three fluctuation regimes: two periodic regimes dependent on culture's cell density and one random (non-periodic) regime with no such dependency, serving as a control. To investigate contingencies on the genetic background, we seeded our experiment with two different strains. After the experiment, growth rate measurements at the two temperatures showed that the evolution of specialists was favored in the random regime, while generalists were favored in the periodic regimes. Whole genome sequencing demonstrated that several gene mutations were selected in parallel in the evolving populations with some dependency on the starting genetic background. Given the genes mutated, we hypothesized that the driving force behind the observed adaptations is the restoration of the internal physiology of the starting strains' unstressed states at 37°C, which may be a means of improving fitness in the new environments. Phenotypic array measurements supported our hypothesis by demonstrating a tendency of the phenotypic response of the evolved strains to move closer to the starting strains' response at the optimum of 37°C, especially for strains classified as generalists.

Adaptation of Bacillus thuringiensis to Plant Colonization Affects Differentiation and Toxicity

The Bacillus cereus group (Bacillus cereus sensu lato) has a diverse ecology, including various species that are vertebrate or invertebrate pathogens. Few isolates from the B. cereus group have however been demonstrated to benefit plant growth. Therefore, it is crucial to explore how bacterial development and pathogenesis evolve during plant colonization. Herein, we investigated Bacillus thuringiensis (Cry-) adaptation to the colonization of Arabidopsis thaliana roots and monitored changes in cellular differentiation in experimentally evolved isolates. Isolates from two populations displayed improved iterative ecesis on roots and increased virulence against insect larvae. Molecular dissection and recreation of a causative mutation revealed the importance of a nonsense mutation in the rho transcription terminator gene. Transcriptome analysis revealed how Rho impacts various B. thuringiensis genes involved in carbohydrate metabolism and virulence. Our work suggests that evolved multicellular aggregates have a fitness advantage over single cells when colonizing plants, creating a trade-off between swimming and multicellularity in evolved lineages, in addition to unrelated alterations in pathogenicity. IMPORTANCE Biologicals-based plant protection relies on the use of safe microbial strains. During application of biologicals to the rhizosphere, microbes adapt to the niche, including genetic mutations shaping the physiology of the cells. Here, the experimental evolution of Bacillus thuringiensis lacking the insecticide crystal toxins was examined on the plant root to reveal how adaptation shapes the differentiation of this bacterium. Interestingly, evolution of certain lineages led to increased hemolysis and insect larva pathogenesis in B. thuringiensis driven by transcriptional rewiring. Further, our detailed study reveals how inactivation of the transcription termination protein Rho promotes aggregation on the plant root in addition to altered differentiation and pathogenesis in B. thuringiensis.

A role for bacterial experimental evolution in coral bleaching mitigation?

Coral reefs are rapidly declining because of widespread mass coral bleaching causing extensive coral mortality. Elevated seawater temperatures are the main drivers of coral bleaching, and climate change is increasing the frequency and severity of destructive marine heatwaves. Efforts to enhance coral thermal bleaching tolerance can be targeted at the coral host or at coral-associated microorganisms (e.g., dinoflagellate endosymbionts and bacteria). The literature on experimental evolution of bacteria suggests that it has value as a tool to increase coral climate resilience. We provide a workflow on how to experimentally evolve coral-associated bacteria to confer thermal tolerance to coral hosts and emphasize the value of implementing this approach in coral reef conservation and restoration efforts.

Thermal adaptation of acetic acid bacteria for practical high-temperature vinegar fermentation

Thermotolerant microorganisms are useful for high-temperature fermentation. Several thermally adapted strains were previously obtained from Acetobacter pasteurianus in a nutrient-rich culture medium, while these adapted strains could not grow well at high temperature in the nutrient-poor practical culture medium, "rice moromi." In this study, A. pasteurianus K-1034 originally capable of performing acetic acid fermentation in rice moromi was thermally adapted by experimental evolution using a "pseudo" rice moromi culture. The adapted strains thus obtained were confirmed to grow well in such the nutrient-poor media in flask or jar-fermentor culture up to 40 or 39 °C; the mutation sites of the strains were also determined. The high-temperature fermentation ability was also shown to be comparable with a low-nutrient adapted strain previously obtained. Using the practical fermentation system, "Acetofermenter," acetic acid production was compared in the moromi culture; the results showed that the adapted strains efficiently perform practical vinegar production under high-temperature conditions.

Identification and Characterization of Base-Substitution Mutations in the Macronuclear Genome of the Ciliate Tetrahymena thermophila

Polyploidy can provide adaptive advantages and drive evolution. Amitotic division of the polyploid macronucleus (MAC) in ciliates acts as a nonsexual genetic mechanism to enhance adaptation to stress conditions and thus provides a unique model to investigate the evolutionary role of polyploidy. Mutation is the primary source of the variation responsible for evolution and adaptation; however, to date, de novo mutations that occur in ciliate MAC genomes during these processes have not been characterized and their biological impacts are undefined. Here, we carried out long-term evolution experiments to directly explore de novo MAC mutations and their molecular features in the model ciliate, Tetrahymena thermophila. A simple but effective method was established to detect base-substitution mutations in evolving populations whereas filtering out most of the false positive base-substitutions caused by repetitive sequences and the programmed genome rearrangements. The detected mutations were rigorously validated using the MassARRAY system. Validated mutations showed a strong G/C→A/T bias, consistent with observations in other species. Moreover, a progressive increase in growth rate of the evolving populations suggested that some of these mutations might be responsible for cell fitness. The established mutation identification and validation methods will be an invaluable resource to make ciliates an important model system to study the role of polyploidy in evolution.

Drivers of Mating Type Composition in Tetrahymena thermophila

Sex offers advantages even in primarily asexual species. Some ciliates appear to utilize such reproductive strategy with many mating types. However, the factors determining the composition of mating types in the unicellular ciliate Tetrahymena thermophila are poorly understood, and this is further complicated by non-Mendelian determination of mating type in the offspring. We therefore developed a novel population genetics model to predict how various factors influence the dynamics of mating type composition, including natural selection. The model predicted either the coexistence of all seven mating types or fixation of a single mating type in a population, depending on parameter combinations, irrespective of natural selection. To understand what factor(s) may be more influential and to test the validity of theoretical prediction, five replicate populations were maintained in laboratory such that several factors could be controlled or measured. Whole-genome sequencing was used to identify newly arising mutations and determine mating type composition. Strikingly, all populations were found to be driven by strong selection on newly arising beneficial mutations to fixation of their carrying mating types, and the trajectories of speed to fixation agreed well with our theoretical predictions. This study illustrates the evolutionary strategies that T. thermophila can utilize to optimize population fitness.

Increased time sampling in an evolve-and-resequence experiment with outcrossing Saccharomyces cerevisiae reveals multiple paths of adaptive change

"Evolve and resequence" (E&R) studies combine experimental evolution and whole-genome sequencing to interrogate the genetics underlying adaptation. Due to ease of handling, E&R work with asexual organisms such as bacteria can employ optimized experimental design, with large experiments and many generations of selection. By contrast, E&R experiments with sexually reproducing organisms are more difficult to implement, and design parameters vary dramatically among studies. Thus, efforts have been made to assess how these differences, such as number of independent replicates, or size of experimental populations, impact inference. We add to this work by investigating the role of time sampling-the number of discrete time points sequence data are collected from evolving populations. Using data from an E&R experiment with outcrossing Saccharomyces cerevisiae in which populations were sequenced 17 times over ~540 generations, we address the following questions: (a) Do more time points improve the ability to identify candidate regions underlying selection? And (b) does high-resolution sampling provide unique insight into evolutionary processes driving adaptation? We find that while time sampling does not improve the ability to identify candidate regions, high-resolution sampling does provide valuable opportunities to characterize evolutionary dynamics. Increased time sampling reveals three distinct trajectories for adaptive alleles: one consistent with classic population genetic theory (i.e., models assuming constant selection coefficients), and two where trajectories suggest more context-dependent responses (i.e., models involving dynamic selection coefficients). We conclude that while time sampling has limited impact on candidate region identification, sampling eight or more time points has clear benefits for studying complex evolutionary dynamics.

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

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

The Single-Stranded RNA Bacteriophage Qβ Adapts Rapidly to High Temperatures: An Evolution Experiment

Single-stranded (ss)RNA viruses are thought to evolve rapidly due to an inherently high mutation rate. However, it remains unclear how ssRNA viruses adapt to novel environments and/or how many and what types of substitutions are needed to facilitate this evolution. In this study, we followed the adaptation of the ssRNA bacteriophage Qβ using thermally adapted Escherichia coli as a host, which can efficiently grow at temperatures between 37.2 and 45.3 °C. This made it possible to evaluate Qβ adaptation to the highest known temperature that supports growth, 45.3 °C. We found that Qβ was capable of replication at this temperature; within 114 days (~1260 generations), we detected more than 34 novel point mutations in the genome of the thermally adapted Qβ population, representing 0.8% of the total Qβ genome. In addition, we returned the 45.3 °C-adapted Qβ populations to 37.2 °C and passaged them for 8 days (~124 generations). We found that the reverse-adapted Qβ population showed little to no decrease in fitness. These results indicate that Qβ can evolve in response to increasing temperatures in a short period of time with the accumulation of point mutations.

Changes in Cell Membrane Fatty Acid Composition of Streptococcus thermophilus in Response to Gradually Increasing Heat Temperature

In this study, a method of heat adaptation was implemented in an attempt to increase the upper thermal threshold of two Streptococcus thermophilus found in South Korea and identified the alterations in membrane fatty acid composition to adaptive response to heat. In order to develop heat tolerant lactic acid bacteria, heat treatment was continuously applied to bacteria by increasing temperature from 60°C until the point that no surviving cell was detected. Our results indicated significant increase in heat tolerance of heat-adapted strains compared to the wild type (WT) strains. In particular, the survival ratio of basically low heat-tolerant strain increased even more. In addition, the strains with improved heat tolerance acquired cross protection, which improved their survival ratio in acid, bile salts and osmotic conditions. A relation between heat tolerance and membrane fatty acid composition was identified. As a result of heat adaptation, the ratio of unsaturated to saturated fatty acids (UFA/SFA) and C18:1 relative concentration were decreased. C6:0 in only heatadapted strains and C22:0 in only the naturally high heat tolerant strain were detected. These results support the hypothesis, that the consequent increase of SFA ratio is a cellular response to environmental stresses such as high temperatures, and it is able to protect the cells from acid, bile salts and osmotic conditions via cross protection. This study demonstrated that the increase in heat tolerance can be utilized as a mean to improve bacterial tolerance against various environmental stresses.

Understanding metabolic adaptation by using bacterial laboratory evolution and trans-omics analysis

Many diseases such as metabolic syndrome, cancer, inflammatory diseases, and pathological phenomena can be understood as an adaptive reconstitution of the metabolic state (metabolic adaptation). One of the effective approaches to reveal the property of metabolic networks is using model organisms such as microorganisms that are easier to experiment with than higher organisms. Using the laboratory evolution approach, we can elucidate the evolutionary dynamics in various stress environments, which provide us an understanding of the metabolic adaptation. In addition, the integration of omics data and phenotypic data enables us to clarify the genetic and phenotypic alterations during adaptation. In this review, we describe our recent studies on bacterial laboratory evolution and the omics approach to clarify the stress adaptation process. We have also obtained high-dimensional phenotypic data using our automated culture system. By combining these genomic and transcriptomic data within high-throughput phenotypic data, we can discuss the complex trans-omics network of metabolic adaptation.

Experimental Challenges for Reduced Genomes: The Cell Model Escherichia coli

Genome reduction, as a top-down approach to obtain the minimal genetic information essential for a living organism, has been conducted with bacterial cells for decades. The most popular and well-studied cell models for genome reduction are Escherichia coli strains. As the previous literature intensively introduced the genetic construction and application of the genome-reduced Escherichia coli strains, the present review focuses the design principles and compares the reduced genome collections from the specific viewpoint of growth, which represents a fundamental property of living cells and is an important feature for their biotechnological application. For the extended simplification of the genomic sequences, the approach of experimental evolution and concern for medium optimization are newly proposed. The combination of the current techniques of genomic construction and the newly proposed methodologies could allow us to acquire growing Escherichia coli cells carrying the extensively reduced genome and to address the question of what the minimal genome essential for life is.

Community-level respiration of prokaryotic microbes may rise with global warming

Understanding how the metabolic rates of prokaryotes respond to temperature is fundamental to our understanding of how ecosystem functioning will be altered by climate change, as these micro-organisms are major contributors to global carbon efflux. Ecological metabolic theory suggests that species living at higher temperatures evolve higher growth rates than those in cooler niches due to thermodynamic constraints. Here, using a global prokaryotic dataset, we find that maximal growth rate at thermal optimum increases with temperature for mesophiles (temperature optima [Formula: see text]C), but not thermophiles ([Formula: see text]C). Furthermore, short-term (within-day) thermal responses of prokaryotic metabolic rates are typically more sensitive to warming than those of eukaryotes. Because climatic warming will mostly impact ecosystems in the mesophilic temperature range, we conclude that as microbial communities adapt to higher temperatures, their metabolic rates and therefore, biomass-specific CO[Formula: see text] production, will inevitably rise. Using a mathematical model, we illustrate the potential global impacts of these findings.

Capacity for survival in global warming: Adaptation of mesophiles to the temperature upper limit

The Intergovernmental Panel on Climate Change recommends keeping the increase in temperature to less than a two-degree increase by the end of the century, but the direct impact of global warming on ecosystems including microbes has not been investigated. Here we performed thermal adaptation of two species and three strains of mesophilic microbes for improvement of the survival upper limit of temperature, and the improvement was evaluated by a newly developed method. To understand the limitation and variation of thermal adaptation, experiments with mutators and by multiple cultures were performed. The results of experiments including genome sequencing and analysis of the characteristics of mutants suggest that these microbes bear a genomic potential to endure a 2–3°C rise in temperature but possess a limited variation of strategies for thermal adaptation.

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.

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

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

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

Transcription termination factor Rho and microbial phenotypic heterogeneity

Populations of genetically identical microorganisms exhibit high degree of cell-to-cell phenotypic diversity even when grown in uniform environmental conditions. Heterogeneity is a genetically determined trait, which ensures bacterial adaptation and survival in the ever changing environmental conditions. Fluctuations in gene expression (noise) at the level of transcription initiation largely contribute to cell-to-cell variability within population. Not surprisingly, the analyses of the mechanisms driving phenotypic heterogeneity are mainly focused on the activity of promoters and transcriptional factors. Less attention is currently given to a role of intrinsic and factor-dependent transcription terminators. Here, we discuss recent advances in understanding the regulatory role of the multi-functional transcription termination factor Rho, the major inhibitor of pervasive transcription in bacteria and the emerging global regulator of gene expression. We propose that termination activity of Rho might be among the mechanisms by which cells manage the intensity of transcriptional noise, thus affecting population heterogeneity.

Coordinated Changes in Mutation and Growth Rates Induced by Genome Reduction

Genome size is determined during evolution, but it can also be altered by genetic engineering in laboratories. The systematic characterization of reduced genomes provides valuable insights into the cellular properties that are quantitatively described by the global parameters related to the dynamics of growth and mutation. In the present study, we analyzed a small collection of W3110 Escherichia coli derivatives containing either the wild-type genome or reduced genomes of various lengths to examine whether the mutation rate, a global parameter representing genomic plasticity, was affected by genome reduction. We found that the mutation rates of these cells increased with genome reduction. The correlation between genome length and mutation rate, which has been reported for the evolution of bacteria, was also identified, intriguingly, for genome reduction. Gene function enrichment analysis indicated that the deletion of many of the genes encoding membrane and transport proteins play a role in the mutation rate changes mediated by genome reduction. Furthermore, the increase in the mutation rate with genome reduction was highly associated with a decrease in the growth rate in a nutrition-dependent manner; thus, poorer media showed a larger change that was of higher significance. This negative correlation was strongly supported by experimental evidence that the serial transfer of the reduced genome improved the growth rate and reduced the mutation rate to a large extent. Taken together, the global parameters corresponding to the genome, growth, and mutation showed a coordinated relationship, which might be an essential working principle for balancing the cellular dynamics appropriate to the environment.

Genome reduction is a powerful approach for investigating the fundamental rules for living systems. Whether genetically disturbed genomes have any specific properties that are different from or similar to those of natively evolved genomes has been under investigation. In the present study, we found that Escherichia coli cells with reduced genomes showed accelerated nucleotide substitution errors (mutation rates), although these cells retained the normal DNA mismatch repair systems. Intriguingly, this finding of correlation between reduced genome size and a higher mutation rate was consistent with the reported evolution of mutation rates. Furthermore, the increased mutation rate was quantitatively associated with a decreased growth rate, indicating that the global parameters related to the genome, growth, and mutation, which represent the amount of genetic information, the efficiency of propagation, and the fidelity of replication, respectively, are dynamically coordinated.

Improvement of isopropanol tolerance of Escherichia coli using adaptive laboratory evolution and omics technologies

Isopropanol (IPA) is the secondary alcohol that can be dehydrated to yield propylene. To produce IPA using microorganisms, a significant issue is that the toxicity of IPA causes retardation or inhibition of cell growth, decreasing the yield. One possible strategy to overcome this problem is to improve IPA tolerance of production organisms. For the understanding of tolerance to IPA, we performed parallel adaptive laboratory evolution (ALE) of Escherichia coli under IPA stress. To identify the genotypic change during ALE, we performed genome re-sequencing analyses of obtained tolerant strains. To verify which mutations were contributed to IPA tolerance, we constructed the mutant strains and quantify the IPA tolerance of the constructed mutants. From these analyses, we found that five mutations (relA, marC, proQ, yfgO, and rraA) provided the increase of IPA tolerance. To understand the phenotypic change during ALE, we performed transcriptome analysis of tolerant strains. From transcriptome analysis, we found that expression levels of genes related to biosynthetic pathways of amino acids, iron ion homeostasis, and energy metabolisms were changed in the tolerant strains. Results from these experiments provide fundamental bases for designing IPA tolerant strains for industrial purposes.

Global relationships in fluctuation and response in adaptive evolution

Cells change their internal state to adapt to environmental changes, and evolve in response to the new conditions. The phenotype changes first via adaptation in response to environmental changes, and then through mutational changes in the genomic sequence, followed by selection in evolution. Here, we analysed simulated adaptive evolution using a simple cell model consisting of thousands of intracellular components, and found that the changes in their concentrations by adaptation are proportional to those by evolution across all the components, where the proportion coefficient between the two agreed well with the change in the growth rate of a cell. Furthermore, we demonstrate that the phenotypic variance in concentrations of cellular components due to (non-genetic) noise and to genomic alternations is proportional across all components. This implies that the specific phenotypes that are highly evolvable were already given by non-genetic fluctuations. These global relationships in cellular states were also supported by phenomenological theory based on steady reproduction and transcriptome analysis of laboratory evolution in Escherichia coli. These findings demonstrate that a possible evolutionary change in phenotypic state is highly restricted. Our results provide a basis for the development of a quantitative theory of plasticity and robustness in phenotypic evolution.

First-Step Mutations during Adaptation Restore the Expression of Hundreds of Genes

The temporal change of phenotypes during the adaptive process remains largely unexplored, as do the genetic changes that affect these phenotypic changes. Here we focused on three mutations that rose to high frequency in the early stages of adaptation within 12 Escherichia coli populations subjected to thermal stress (42 °C). All the mutations were in the rpoB gene, which encodes the RNA polymerase beta subunit. For each mutation, we measured the growth curves and gene expression (mRNAseq) of clones at 42 °C. We also compared growth and gene expression with their ancestor under unstressed (37 °C) and stressed conditions (42 °C). Each of the three mutations changed the expression of hundreds of genes and conferred large fitness advantages, apparently through the restoration of global gene expression from the stressed toward the prestressed state. These three mutations had a similar effect on gene expression as another single mutation in a distinct domain of the rpoB protein. Finally, we compared the phenotypic characteristics of one mutant, I572L, with two high-temperature adapted clones that have this mutation plus additional background mutations. The background mutations increased fitness, but they did not substantially change gene expression. We conclude that early mutations in a global transcriptional regulator cause extensive changes in gene expression, many of which are likely under positive selection for their effect in restoring the prestress physiology.

Bacterial transcriptome reorganization in thermal adaptive evolution

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Functional specialization in regulation and quality control in thermal adaptive evolution

Distinctive survival strategies, specialized in regulation and in quality control, were observed in thermal adaptive evolution with a laboratory Escherichia coli strain. The two specialists carried a single mutation either within rpoH or upstream of groESL, which led to the activated global regulation by sigma factor 32 or an increased amount of GroEL/ES chaperonins, respectively. Although both specialists succeeded in thermal adaptation, the common winner of the evolution was the specialist in quality control, that is, the strategy of chaperonin-mediated protein folding. To understand this evolutionary consequence, multilevel analyses of cellular status, for example, transcriptome, protein and growth fitness, were carried out. The specialist in quality control showed less change in transcriptional reorganization responding to temperature increase, which was consistent with the finding of that the two specialists showed the biased expression of molecular chaperones. Such repressed changes in gene expression seemed to be advantageous for long-term sustainability because a specific increase in chaperonins not only facilitated the folding of essential gene products but also saved cost in gene expression compared with the overall transcriptional increase induced by rpoH regulation. Functional specialization offered two strategies for successful thermal adaptation, whereas the evolutionary advantageous was more at the points of cost-saving in gene expression and the essentiality in protein folding.

Phenotypic convergence in bacterial adaptive evolution to ethanol stress

BACKGROUND

Bacterial cells have a remarkable ability to adapt to environmental changes, a phenomenon known as adaptive evolution. During adaptive evolution, phenotype and genotype dynamically changes; however, the relationship between these changes and associated constraints is yet to be fully elucidated.

RESULTS

In this study, we analyzed phenotypic and genotypic changes in Escherichia coli cells during adaptive evolution to ethanol stress. Phenotypic changes were quantified by transcriptome and metabolome analyses and were similar among independently evolved ethanol tolerant populations, which indicate the existence of evolutionary constraints in the dynamics of adaptive evolution. Furthermore, the contribution of identified mutations in one of the tolerant strains was evaluated using site-directed mutagenesis. The result demonstrated that the introduction of all identified mutations cannot fully explain the observed tolerance in the tolerant strain.

CONCLUSIONS

The results demonstrated that the convergence of adaptive phenotypic changes and diverse genotypic changes, which suggested that the phenotype-genotype mapping is complex. The integration of transcriptome and genome data provides a quantitative understanding of evolutionary constraints.

Evolutionary Consequence of a Trade-Off between Growth and Maintenance along with Ribosomal Damages

Microorganisms in nature are constantly subjected to a limited availability of resources and experience repeated starvation and nutrition. Therefore, microbial life may evolve for both growth fitness and sustainability. By contrast, experimental evolution, as a powerful approach to investigate microbial evolutionary strategies, often targets the increased growth fitness in controlled, steady-state conditions. Here, we address evolutionary changes balanced between growth and maintenance while taking nutritional fluctuations into account. We performed a 290-day-long evolution experiment with a histidine-requiring Escherichia coli strain that encountered repeated histidine-rich and histidine-starved conditions. The cells that experienced seven rounds of starvation and re-feed grew more sustainably under prolonged starvation but dramatically lost growth fitness under rich conditions. The improved sustainability arose from the evolved capability to use a trace amount of histidine for cell propagation. The reduced growth rate was attributed to mutations genetically disturbing the translation machinery, that is, the ribosome, ultimately slowing protein translation. This study provides the experimental demonstration of slow growth accompanied by an enhanced affinity to resources as an evolutionary adaptation to oscillated environments and verifies that it is possible to evolve for reduced growth fitness. Growth economics favored for population increase under extreme resource limitations is most likely a common survival strategy adopted by natural microbes.

Molecular Clock of Neutral Mutations in a Fitness-Increasing Evolutionary Process

The molecular clock of neutral mutations, which represents linear mutation fixation over generations, is theoretically explained by genetic drift in fitness-steady evolution or hitchhiking in adaptive evolution. The present study is the first experimental demonstration for the molecular clock of neutral mutations in a fitness-increasing evolutionary process. The dynamics of genome mutation fixation in the thermal adaptive evolution of Escherichia coli were evaluated in a prolonged evolution experiment in duplicated lineages. The cells from the continuously fitness-increasing evolutionary process were subjected to genome sequencing and analyzed at both the population and single-colony levels. Although the dynamics of genome mutation fixation were complicated by the combination of the stochastic appearance of adaptive mutations and clonal interference, the mutation fixation in the population was simply linear over generations. Each genome in the population accumulated 1.6 synonymous and 3.1 non-synonymous neutral mutations, on average, by the spontaneous mutation accumulation rate, while only a single genome in the population occasionally acquired an adaptive mutation. The neutral mutations that preexisted on the single genome hitchhiked on the domination of the adaptive mutation. The successive fixation processes of the 128 mutations demonstrated that hitchhiking and not genetic drift were responsible for the coincidence of the spontaneous mutation accumulation rate in the genome with the fixation rate of neutral mutations in the population. The molecular clock of neutral mutations to the fitness-increasing evolution suggests that the numerous neutral mutations observed in molecular phylogenetic trees may not always have been fixed in fitness-steady evolution but in adaptive evolution.

Mutations that have little influence on biological function are referred to as neutral mutations and frequently appear in molecular phylogenetic analyses. The fixation of neutral mutations in populations has been attributed to genetic drift in fitness-steady evolutionary processes or hitchhiking in adaptive evolution. We examined the fitness-increasing evolution of Escherichia coli for thermal adaptation to observe the fixation dynamics of genome-wide mutations. In the adaptive evolution, all genomes in the population equally accumulated neutral mutations by replication errors. The infrequent occurrence of an adaptive mutation on one of the genomes by chance resulted in the fixation of the neutral mutations that had pre-accumulated in the same genome by hitchhiking. Via successive hitchhiking events, the neutral mutations were fixed in the population linearly over generations at the same rate as the spontaneous mutation accumulation rate in the genome. The molecular clock of neutral mutations thus functions even in adaptive evolution. The evolutionary period characterized by the accumulation of numerous neutral mutations observed in molecular phylogenetic trees may not be specific to neutral evolution but may occur in adaptive evolution as well.

Thermal and solvent stress cross-tolerance conferred to Corynebacterium glutamicum by adaptive laboratory evolution

Reinforcing microbial thermotolerance is a strategy to enable fermentation with flexible temperature settings and thereby to save cooling costs. Here, we report on adaptive laboratory evolution (ALE) of the amino acid-producing bacterium Corynebacterium glutamicum under thermal stress. After 65 days of serial passage of the transgenic strain GLY3, in which the glycolytic pathway is optimized for alanine production under oxygen deprivation, three strains adapted to supraoptimal temperatures were isolated, and all the mutations they acquired were identified by whole-genome resequencing. Of the 21 mutations common to the three strains, one large deletion and two missense mutations were found to promote growth of the parental strain under thermal stress. Additive effects on thermotolerance were observed among these mutations, and the combination of the deletion with the missense mutation on otsA, encoding a trehalose-6-phosphate synthase, allowed the parental strain to overcome the upper limit of growth temperature. Surprisingly, the three evolved strains acquired cross-tolerance for isobutanol, which turned out to be partly attributable to the genomic deletion associated with the enhanced thermotolerance. The deletion involved loss of two transgenes, pfk and pyk, encoding the glycolytic enzymes, in addition to six native genes, and elimination of the transgenes, but not the native genes, was shown to account for the positive effects on thermal and solvent stress tolerance, implying a link between energy-producing metabolism and bacterial stress tolerance. Overall, the present study provides evidence that ALE can be a powerful tool to refine the phenotype of C. glutamicum and to investigate the molecular bases of stress tolerance.

Characterization of a single mutation in TraQ in a strain of Escherichia coli partially resistant to Qβ infection

Bacteria and virulent bacteriophages are in a prey–predator relationship. Experimental models under simplified conditions with the presence of bacteria and bacteriophages have been used to elucidate the mechanisms that have enabled both prey and predator to coexist over long periods. In experimental coevolution conducted with Escherichia coli and the virulent RNA bacteriophage Qβ in serial transfer, both coexisted for at least for 54 days, during which time they continued to change genetically and phenotypically. By day 16, an E. coli strain partially resistant to Qβ appeared and caused an approximately 10⁴-fold decrease in Qβ amplification. Whole-genome analysis of this strain suggested that a single mutation in TraQ was responsible for the partially resistant phenotype. TraQ interacts with propilin, encoded by the traA gene and a precursor of pilin, which is a component of the F pilus. The present study was performed to elucidate the mechanism underlying the coexistence of E. coli and Qβ by investigating how a mutation in TraQ altered the physiological state of E. coli, and thus the amplification of Qβ. Overexpression of wild-type TraQ in the partially resistant E. coli strain resulted in recovery of both TraA protein content, including propilin and pilin, and Qβ amplification to levels comparable to those observed in the susceptible strain. Intriguingly, overexpression of the mutant TraQ in the partially resistant strains also increased the levels of TraA protein and Qβ amplification, but these increases were smaller than those observed in the wild-type strain or the partially resistant strain expressing wild-type TraQ. The results of this study represent an example of how E. coli can become partially resistant to RNA bacteriophage infection via changes in a protein involved in maturation of a receptor rather than in the receptor itself and of how E. coli can stably coexist with virulent RNA bacteriophages.

Directed evolution of cell size in Escherichia coli

Background

In bacteria, cell size affects chromosome replication, the assembly of division machinery, cell wall synthesis, membrane synthesis and ultimately growth rate. In addition, cell size can also be a target for Darwinian evolution for protection from predators. This strong coupling of cell size and growth, however, could lead to the introduction of growth defects after size evolution. An important question remains: can bacterial cell size change and/or evolve without imposing a growth burden?

Results

The directed evolution of particular cell sizes, without a growth burden, was tested with a laboratory Escherichia coli strain. Cells of defined size ranges were collected by a cell sorter and were subsequently cultured. This selection-propagation cycle was repeated, and significant changes in cell size were detected within 400 generations. In addition, the width of the size distribution was altered. The changes in cell size were unaccompanied by a growth burden. Whole genome sequencing revealed that only a few mutations in genes related to membrane synthesis conferred the size evolution.

Conclusions

In conclusion, bacterial cell size could evolve, through a few mutations, without growth reduction. The size evolution without growth reduction suggests a rapid evolutionary change to diverse cell sizes in bacterial survival strategies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0257-1) contains supplementary material, which is available to authorized users.

Nutrient-dependent growth defects and mutability of mutators in Escherichia coli

So-called mutators emerge when mismatch repair and proofreading mechanisms are defective. Mutators not only accelerate the accumulation of mutations that are beneficial for adaptation but also cause a large number of deleterious mutations that are disadvantageous for cell growth. However, such growth defects may be compensated by nutrient availability. How the growth burden is associated with high mutability in relation to nutritional variation is an intriguing question. To address this question, we constructed a variety of Escherichia coli mutator strains through combinatorial deletions of mismatch repair and proofreading genes and quantitatively evaluated their growth and mutation rates under different nutritional conditions. Growth defects caused by high mutation rates were commonly observed in all mutators, and these defects were alleviated by nutrient supplementation in most mutators. In addition, the mutation rates of the mutators fluctuated greatly in response to nutritional conditions, in contrast to the nearly constant mutation rate of the wild-type strain under varying nutritional conditions. The results showed conditional growth defects and nutrition-sensitive mutability as general features of mutators. This study indicates the importance of modulating mutability in response to changing nutrient conditions to minimize the risk of extinction due to genetic load.

Contribution of silent mutations to thermal adaptation of RNA bacteriophage Qβ

Changes in protein function and other biological properties, such as RNA structure, are crucial for adaptation of organisms to novel or inhibitory environments. To investigate how mutations that do not alter amino acid sequence may be positively selected, we performed a thermal adaptation experiment using the single-stranded RNA bacteriophage Qβ in which the culture temperature was increased from 37.2°C to 41.2°C and finally to an inhibitory temperature of 43.6°C in a stepwise manner in three independent lines. Whole-genome analysis revealed 31 mutations, including 14 mutations that did not result in amino acid sequence alterations, in this thermal adaptation. Eight of the 31 mutations were observed in all three lines. Reconstruction and fitness analyses of Qβ strains containing only mutations observed in all three lines indicated that five mutations that did not result in amino acid sequence changes but increased the amplification ratio appeared in the course of adaptation to growth at 41.2°C. Moreover, these mutations provided a suitable genetic background for subsequent mutations, altering the fitness contribution from deleterious to beneficial. These results clearly showed that mutations that do not alter the amino acid sequence play important roles in adaptation of this single-stranded RNA virus to elevated temperature.

IMPORTANCE

Recent studies using whole-genome analysis technology suggested the importance of mutations that do not alter the amino acid sequence for adaptation of organisms to novel environmental conditions. It is necessary to investigate how these mutations may be positively selected and to determine to what degree such mutations that do not alter amino acid sequences contribute to adaptive evolution. Here, we report the roles of these silent mutations in thermal adaptation of RNA bacteriophage Qβ based on experimental evolution during which Qβ showed adaptation to growth at an inhibitory temperature. Intriguingly, four synonymous mutations and one mutation in the untranslated region that spread widely in the Qβ population during the adaptation process at moderately high temperature provided a suitable genetic background to alter the fitness contribution of subsequent mutations from deleterious to beneficial at a higher temperature.

Evolution of Escherichia coli to 42 °C and subsequent genetic engineering reveals adaptive mechanisms and novel mutations

Adaptive laboratory evolution (ALE) has emerged as a valuable method by which to investigate microbial adaptation to a desired environment. Here, we performed ALE to 42 °C of ten parallel populations of Escherichia coli K-12 MG1655 grown in glucose minimal media. Tightly controlled experimental conditions allowed selection based on exponential-phase growth rate, yielding strains that uniformly converged toward a similar phenotype along distinct genetic paths. Adapted strains possessed as few as 6 and as many as 55 mutations, and of the 144 genes that mutated in total, 14 arose independently across two or more strains. This mutational recurrence pointed to the key genetic targets underlying the evolved fitness increase. Genome engineering was used to introduce the novel ALE-acquired alleles in random combinations into the ancestral strain, and competition between these engineered strains reaffirmed the impact of the key mutations on the growth rate at 42 °C. Interestingly, most of the identified key gene targets differed significantly from those found in similar temperature adaptation studies, highlighting the sensitivity of genetic evolution to experimental conditions and ancestral genotype. Additionally, transcriptomic analysis of the ancestral and evolved strains revealed a general trend for restoration of the global expression state back toward preheat stressed levels. This restorative effect was previously documented following evolution to metabolic perturbations, and thus may represent a general feature of ALE experiments. The widespread evolved expression shifts were enabled by a comparatively scant number of regulatory mutations, providing a net fitness benefit but causing suboptimal expression levels for certain genes, such as those governing flagellar formation, which then became targets for additional ameliorating mutations. Overall, the results of this study provide insight into the adaptation process and yield lessons important for the future implementation of ALE as a tool for scientific research and engineering.

The power to detect quantitative trait loci using resequenced, experimentally evolved populations of diploid, sexual organisms

A novel approach for dissecting complex traits is to experimentally evolve laboratory populations under a controlled environment shift, resequence the resulting populations, and identify single nucleotide polymorphisms (SNPs) and/or genomic regions highly diverged in allele frequency. To better understand the power and localization ability of such an evolve and resequence (E&R) approach, we carried out forward-in-time population genetics simulations of 1 Mb genomic regions under a large combination of experimental conditions, then attempted to detect significantly diverged SNPs. Our analysis indicates that the ability to detect differentiation between populations is primarily affected by selection coefficient, population size, number of replicate populations, and number of founding haplotypes. We estimate that E&R studies can detect and localize causative sites with 80% success or greater when the number of founder haplotypes is over 500, experimental populations are replicated at least 25-fold, population size is at least 1,000 diploid individuals, and the selection coefficient on the locus of interest is at least 0.1. More achievable experimental designs (less replicated, fewer founder haplotypes, smaller effective population size, and smaller selection coefficients) can have power of greater than 50% to identify a handful of SNPs of which one is likely causative. Similarly, in cases where s ≥ 0.2, less demanding experimental designs can yield high power.

Genomic diversity of "deep ecotype" Alteromonas macleodii isolates: evidence for Pan-Mediterranean clonal frames

We have compared genomes of Alteromonas macleodii “deep ecotype” isolates from two deep Mediterranean sites and two surface samples from the Aegean and the English Channel. A total of nine different genomes were analyzed. They belong to five clonal frames (CFs) that differ among them by approximately 30,000 single-nucleotide polymorphisms (SNPs) over their core genomes. Two of the CFs contain three strains each with nearly identical genomes (∼100 SNPs over the core genome). One of the CFs had representatives that were isolated from samples taken more than 1,000 km away, 2,500 m deeper, and 5 years apart. These data mark the longest proven persistence of a CF in nature (outside of clinical settings). We have found evidence for frequent recombination events between or within CFs and even with the distantly related A. macleodii surface ecotype. The different CFs had different flexible genomic islands. They can be classified into two groups; one type is additive, that is, containing different numbers of gene cassettes, and is very variable in short time periods (they often varied even within a single CF). The other type was more stable and produced the complete replacement of a genomic fragment by another with different genes. Although this type was more conserved within each CF, we found examples of recombination among distantly related CFs including English Channel and Mediterranean isolates.

Microbial evolution in vivo and in silico: methods and applications

Microbial evolution has been extensively studied in the past fifty years, which has lead to seminal discoveries that have shaped our understanding of evolutionary forces and dynamics. It is only recently however, that transformative technologies and computational advances have enabled a larger in-scale and in-depth investigation of the genetic basis and mechanistic underpinnings of evolutionary adaptation. In this review we focus on the strengths and limitations of in vivo and in silico techniques for studying microbial evolution in the laboratory, and we discuss how these complementary approaches can be integrated in a unifying framework for elucidating microbial evolution.

Emerging tools for synthetic genome design

Synthetic biology is an emerging discipline for designing and synthesizing predictable, measurable, controllable, and transformable biological systems. These newly designed biological systems have great potential for the development of cheaper drugs, green fuels, biodegradable plastics, and targeted cancer therapies over the coming years. Fortunately, our ability to quickly and accurately engineer biological systems that behave predictably has been dramatically expanded by significant advances in DNA-sequencing, DNA-synthesis, and DNA-editing technologies. Here, we review emerging technologies and methodologies in the field of building designed biological systems, and we discuss their future perspectives.

Multilevel comparative analysis of the contributions of genome reduction and heat shock to the Escherichia coli transcriptome

Background

Both large deletions in genome and heat shock stress would lead to alterations in the gene expression profile; however, whether there is any potential linkage between these disturbances to the transcriptome have not been discovered. Here, the relationship between the genomic and environmental contributions to the transcriptome was analyzed by comparing the transcriptomes of the bacterium Escherichia coli (strain MG1655 and its extensive genomic deletion derivative, MDS42) grown in regular and transient heat shock conditions.

Results

The transcriptome analysis showed the following: (i) there was a reorganization of the transcriptome in accordance with preferred chromosomal periodicity upon genomic or heat shock perturbation; (ii) there was a considerable overlap between the perturbed regulatory networks and the categories enriched for differentially expressed genes (DEGs) following genome reduction and heat shock; (iii) the genes sensitive to genome reduction tended to be located close to genomic scars, and some were also highly responsive to heat shock; and (iv) the genomic and environmental contributions to the transcriptome displayed not only a positive correlation but also a negatively compensated relationship (i.e., antagonistic epistasis).

Conclusion

The contributions of genome reduction and heat shock to the Escherichia coli transcriptome were evaluated at multiple levels. The observations of overlapping perturbed networks, directional similarity in transcriptional changes, positive correlation and epistatic nature linked the two contributions and suggest somehow a crosstalk guiding transcriptional reorganization in response to both genetic and environmental disturbances in bacterium E. coli.

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.

Current understanding of the formation and adaptation of metabolic systems based on network theory

Formation and adaptation of metabolic networks has been a long-standing question in biology. With recent developments in biotechnology and bioinformatics, the understanding of metabolism is progressively becoming clearer from a network perspective. This review introduces the comprehensive metabolic world that has been revealed by a wide range of data analyses and theoretical studies; in particular, it illustrates the role of evolutionary events, such as gene duplication and horizontal gene transfer, and environmental factors, such as nutrient availability and growth conditions, in evolution of the metabolic network. Furthermore, the mathematical models for the formation and adaptation of metabolic networks have also been described, according to the current understanding from a perspective of metabolic networks. These recent findings are helpful in not only understanding the formation of metabolic networks and their adaptation, but also metabolic engineering.

Parallel genetic changes and nonparallel gene-environment interactions characterize the evolution of drug resistance in yeast

Beneficial mutations are required for adaptation to novel environments, yet the range of mutational pathways that are available to a population has been poorly characterized, particularly in eukaryotes. We assessed the genetic changes of the first mutations acquired during adaptation to a novel environment (exposure to the fungicide, nystatin) in 35 haploid lines of Saccharomyces cerevisiae. Through whole-genome resequencing we found that the genomic scope for adaptation was narrow; all adapted lines acquired a mutation in one of four late-acting genes in the ergosterol biosynthesis pathway, with very few other mutations found. Lines that acquired different ergosterol mutations in the same gene exhibited very similar tolerance to nystatin. All lines were found to have a cost relative to wild type in an unstressful environment; the level of this cost was also strongly correlated with the ergosterol gene bearing the mutation. Interestingly, we uncovered both positive and negative effects on tolerance to other harsh environments for mutations in the different ergosterol genes, indicating that these beneficial mutations have effects that differ in sign among environmental challenges. These results demonstrate that although the genomic target was narrow, different adaptive mutations can lead populations down different evolutionary pathways, with respect to their ability to tolerate (or succumb to) other environmental challenges.