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

Thermally adapted Escherichia coli keeps transcriptomic response during temperature upshift exposure

The heat shock response is a cellular protection mechanism against sudden temperature upshifts extensively studied in Escherichia coli. However, the effects of thermal evolution on this response remain largely unknown. In this study, we investigated the early and late physiological and transcriptional responses to temperature upshift in a thermotolerant strain under continuous culture conditions. Adaptive laboratory evolution was performed on a metabolically engineered E. coli strain (JU15), designed for D-lactic acid production, to enable cellular growth and fermentation of glucose at 45 °C in batch cultures. The resulting homofermentative strain, ECL45, successfully adapted to 45 °C in a glucose-mineral medium at pH 7 under non-aerated conditions. The thermal-adapted ECL45 retained the parental strain's high volumetric productivity and product/substrate yield. Genomic sequencing of ECL45 revealed eight mutations, including one in a non-coding region and six within the coding regions of genes associated with metabolic, transport, and regulatory functions. Transcriptomic analysis comparing the evolved strain with its parental counterpart under early and late temperature upshifts indicated that the adaptation involved a controlled stringent response. This mechanism likely contributes to the strain's ability to maintain growth capacity at high temperatures. KEY POINTS: • The temperature upshift response of a thermally adapted strain in continuous culture was studied for the first time. • Genomic analyses revealed the presence of a double point mutation in the spoT gene. • The thermally adapted strain maintained underexpression of the spoT gene at high temperatures. • Supplementation of 0.15 g/L of hydrolyzed protein favored thermal adaptation at 45 °C.

Limited Diversity of Thermal Adaptation to a Critical Temperature in Zymomonas mobilis: Evidence from Multiple-Parallel Laboratory Evolution Experiments

Laboratory evolution is an effective means of understanding microbial adaptation to the environment. We previously isolated four thermoadapted Zymomonas mobilis mutants, which showed a 2 °C rise in the critical high temperature (CHT), by performing multiple-parallel adaptation experiments. In the present study, the individual mutations in these mutants were intensively analyzed. Two mutations in each adapted mutant were found to primarily contribute to the increase in the upper temperature limit. RNA sequencing (RNA-seq) analysis revealed that the two mutations led to the upregulation of 79-185 genes and the downregulation of 242-311 genes. The findings from transcriptomic and physiological experiments suggest two common and primary mechanisms for thermal resistance: a decrease in the activity of diacylglycerol kinase, which may change the structure of lipopolysaccharide (LPS) probably to strengthen the membrane structure, and an increase in the expression of genes for GroEL/GroES or cell wall hydrolase to repair the protein or membrane damage that occurs at such critical temperatures. Additionally, transporters including efflux pumps may contribute to intracellular homeostasis by expelling toxic compounds such as ethanol and acetate or by maintaining the K+ concentration. The results of this study on four independently thermoadapted mutants led to the conclusion that the mutants have almost the same thermal adaptation strategies and thus their molecular diversity is limited.

Rho-dependent termination: a bacterial evolutionary capacitor for stress resistance

Since the Modern Synthesis, interest has grown in resolving the "black box" between genotype and phenotype. Contained within this black box are highly plastic RNA and proteins with global effects on chromosome integrity and gene expression that serve as evolutionary capacitors - elements that enable the accumulation and buffering of genetic variation in normal conditions and reveal hidden genetic variation when induced by environmental stress. Discussion of evolutionary capacitors has primarily focused on eukaryotic translation factors and chaperones, such as Hsp90 and PSI+ prion. However, due to the coupling of transcription and translation in prokaryotes, transcription factors can be equally impactful in the modulation of gene expression and phenotypes. In this review, we discuss the prokaryotic transcription terminator Rho and how mutagenesis and plasticity of Rho influence epistasis, evolvability, and adaptation to stress in bacteria. We discuss the effects of variation in Rho generated by nature, laboratory mutagenesis, and experimental evolution; and how this variation is constrained or encouraged by Rho's extensive network of protein interactors. Exploring Rho's role as an evolutionary capacitor, along with identifying additional elements that can serve this function, can significantly advance our understanding of how organisms adapt to thrive in diverse environments.

Simultaneous saccharification and fermentation for D-lactic acid production using a metabolically engineered Escherichia coli adapted to high temperature

BACKGROUND

Escherichia coli JU15 is a metabolically engineered strain capable to metabolize C5 and C6 sugars with a high yield of D-lactic acid production at its optimal growth temperature (37 °C). The simultaneous saccharification and fermentation process allow to use lignocellulosic biomass as a cost-effective and high-yield strategy. However, this process requires microorganisms capable of growth at a temperature close to 50 °C, at which the activity of cellulolytic enzymes works efficiently.

RESULTS

The thermotolerant strain GT48 was generated by adaptive laboratory evolution in batch and chemostat cultures under temperature increments until 48 °C. The strain GT48 was able to grow and ferment glucose to D-lactate at 47 °C. It was found that a pH of 6.3 conciliated with GT48 growth and cellulase activity of a commercial cocktail. Hence, this pH was used for the SSF of a diluted acid-pretreated corn stover (DAPCS) at a solid load of 15% (w/w), 15 FPU/g-DAPCS, and 47 °C. Under such conditions, the strain GT48 exhibited remarkable performance, producing D-lactate at a level of 1.41, 1.42, and 1.48-fold higher in titer, productivity, and yield, respectively, compared to parental strain at 45 °C.

CONCLUSIONS

In general, our results show for the first time that a thermal-adapted strain of E. coli is capable of being used in the simultaneous saccharification and fermentation process without pre-saccharification stage at high temperatures.

Genomic evidence of Escherichia coli gut population diversity translocation in leukemia patients

Escherichia coli, a commensal species of the human gut, is an opportunistic pathogen that can reach extra-intestinal compartments, including the bloodstream and the bladder, among others. In non-immunosuppressed patients, purifying or neutral evolution of E. coli populations has been reported in the gut. Conversely, it has been suggested that when migrating to extra-intestinal compartments, E. coli genomes undergo diversifying selection as supported by strong evidence for adaptation. The level of genomic polymorphism and the size of the populations translocating from gut to extra-intestinal compartments is largely unknown. To gain insights into the pathophysiology of these translocations, we investigated the level of polymorphism and the evolutionary forces acting on the genomes of 77 E. coli isolated from various compartments in three immunosuppressed patients. Each patient had a unique strain, which was a mutator in one case. In all instances, we observed that translocation encompasses much of the genomic diversity present in the gut. The same signature of selection, whether purifying or diversifying, and as anticipated, neutral for mutator isolates, was observed in both the gut and bloodstream. Additionally, we found a limited number of non-specific mutations among compartments for non-mutator isolates. In all cases, urine isolates were dominated by neutral selection. These findings indicate that substantial proportions of populations are undergoing translocation and that they present a complex compartment-specific pattern of selection at the patient level.IMPORTANCEIt has been suggested that intra and extra-intestinal compartments differentially constrain the evolution of E. coli strains. Whether host particular conditions, such as immunosuppression, could affect the strain evolutionary trajectories remains understudied. We found that, in immunosuppressed patients, large fractions of E. coli gut populations are translocating with variable modifications of the signature of selection for commensal and pathogenic isolates according to the compartment and/or the patient. Such multiple site sampling should be performed in large cohorts of patients to gain a better understanding of E. coli extra-intestinal diseases.

Evolution of pH-sensitive transcription termination in Escherichia coli 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 Escherichia coli populations experimentally evolved under repeated long-term starvation conditions, during which the accumulation of metabolic waste followed by transfer into fresh media results in drastic environmental pH fluctuations associated with feast and famine. 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 plasticity via an alkaline-induced reduction of Rho function that, when found in tandem with a ΔydcI allele, leads to intracellular alkalization and genetic assimilation of Rho mutant function. We further identify Arg to His substitutions at analogous sites in rho alleles from species that regularly experience neutral to alkaline pH fluctuations in their environments. Our results suggest that Arg to His substitutions in Rho may serve to rapidly coordinate complex physiological responses through pH sensing and shed light on how cellular populations use environmental cues to coordinate rapid responses to complex, fluctuating environments.

Development of specialized devices for microbial experimental evolution

Experimental evolution of microbial cells provides valuable information on evolutionary dynamics, such as mutations that contribute to fitness gain under given selection pressures. Although experimental evolution is a promising tool in evolutionary biology and bioengineering, long-term culture experiments under multiple environmental conditions often impose an excessive workload on researchers. Therefore, the development of automated systems significantly contributes to the advancement of experimental evolutionary research. This review presents several specialized devices designed for experimental evolution studies, such as an automated system for high-throughput culture experiments, a culture device that generate a temperature gradient, and an automated ultraviolet (UV) irradiation culture device. The ongoing development of such specialized devices is poised to continually expand new frontiers in experimental evolution research.

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.