Anaerobically Grown Escherichia coli Has an Enhanced Mutation Rate and Distinct Mutational Spectra

Ecological and evolutionary mechanisms driving within-patient emergence of antimicrobial resistance

The ecological and evolutionary mechanisms of antimicrobial resistance (AMR) emergence within patients and how these vary across bacterial infections are poorly understood. Increasingly widespread use of pathogen genome sequencing in the clinic enables a deeper understanding of these processes. In this Review, we explore the clinical evidence to support four major mechanisms of within-patient AMR emergence in bacteria: spontaneous resistance mutations; in situ horizontal gene transfer of resistance genes; selection of pre-existing resistance; and immigration of resistant lineages. Within-patient AMR emergence occurs across a wide range of host niches and bacterial species, but the importance of each mechanism varies between bacterial species and infection sites within the body. We identify potential drivers of such differences and discuss how ecological and evolutionary analysis could be embedded within clinical trials of antimicrobials, which are powerful but underused tools for understanding why these mechanisms vary between pathogens, infections and individuals. Ultimately, improving understanding of how host niche, bacterial species and antibiotic mode of action combine to govern the ecological and evolutionary mechanism of AMR emergence in patients will enable more predictive and personalized diagnosis and antimicrobial therapies.

Mutational spectra are associated with bacterial niche

As observed in cancers, individual mutagens and defects in DNA repair create distinctive mutational signatures that combine to form context-specific spectra within cells. We reasoned that similar processes must occur in bacterial lineages, potentially allowing decomposition analysis to detect both disruption of DNA repair processes and exposure to niche-specific mutagens. Here we reconstruct mutational spectra for 84 clades from 31 diverse bacterial species and find distinct mutational patterns. We extract signatures driven by specific DNA repair defects using hypermutator lineages, and further deconvolute the spectra into multiple signatures operating within different clades. We show that these signatures are explained by both bacterial phylogeny and replication niche. By comparing mutational spectra of clades from different environmental and biological locations, we identify niche-associated mutational signatures, and then employ these signatures to infer the predominant replication niches for several clades where this was previously obscure. Our results show that mutational spectra may be associated with sites of bacterial replication when mutagen exposures differ, and can be used in these cases to infer transmission routes for established and emergent human bacterial pathogens.

The divergence of mutation rates and spectra across the Tree of Life

Owing to advances in genome sequencing, genome stability has become one of the most scrutinized cellular traits across the Tree of Life. Despite its centrality to all things biological, the mutation rate (per nucleotide site per generation) ranges over three orders of magnitude among species and several-fold within individual phylogenetic lineages. Within all major organismal groups, mutation rates scale negatively with the effective population size of a species and with the amount of functional DNA in the genome. This relationship is most parsimoniously explained by the drift-barrier hypothesis, which postulates that natural selection typically operates to reduce mutation rates until further improvement is thwarted by the power of random genetic drift. Despite this constraint, the molecular mechanisms underlying DNA replication fidelity and repair are free to wander, provided the performance of the entire system is maintained at the prevailing level. The evolutionary flexibility of the mutation rate bears on the resolution of several prior conundrums in phylogenetic and population-genetic analysis and raises challenges for future applications in these areas.

The influence of oxygen and oxidative stress on de novo acquisition of antibiotic resistance in E. coli and Lactobacillus lactis

BACKGROUND

Bacteria can acquire resistance through DNA mutations in response to exposure to sub-lethal concentrations of antibiotics. According to the radical-based theory, reactive oxygen species (ROS), a byproduct of the respiratory pathway, and oxidative stress caused by reactive metabolic byproducts, play a role in cell death as secondary killing mechanism. In this study we address the question whether ROS also affects development of resistance, in the conditions that the cells is not killed by the antibiotic.

RESULTS

To investigate whether oxygen and ROS affect de novo acquisition of antibiotic resistance, evolution of resistance due to exposure to non-lethal levels of antimicrobials was compared in E. coli wildtype and ΔoxyR strains under aerobic and anaerobic conditions. Since Lactococcus lactis (L. lactis) does not have an active electron transport chain (ETC) even in the presence of oxygen, and thus forms much less ROS, resistance development in L. lactis was used to distinguish between oxygen and ROS. The resistance acquisition in E. coli wildtype under aerobic and anaerobic conditions did not differ much. However, the aerobically grown ΔoxyR strain gained resistance faster than the wildtype or anaerobic ΔoxyR. Inducing an ETC by adding heme increased the rate at which L. lactis acquired resistance. Whole genome sequencing identified specific mutations involved in the acquisition of resistance. These mutations were specific for each antibiotic. The lexA mutation in ΔoxyR strain under aerobic conditions indicated that the SOS response was involved in resistance acquisition.

CONCLUSIONS

The concept of hormesis can explain the beneficial effects of low levels of ROS and reactive metabolic byproducts, while high levels are lethal. DNA repair and mutagenesis may therefore expedite development of resistance. Taken together, the results suggest that oxygen as such barely affects resistance development. Nevertheless, non-lethal levels of ROS stimulate de novo acquisition of antibiotic resistance.

Translating eco-evolutionary biology into therapy to tackle antibiotic resistance

Antibiotic resistance is currently one of the most important public health problems. The golden age of antibiotic discovery ended decades ago, and new approaches are urgently needed. Therefore, preserving the efficacy of the antibiotics currently in use and developing compounds and strategies that specifically target antibiotic-resistant pathogens is critical. The identification of robust trends of antibiotic resistance evolution and of its associated trade-offs, such as collateral sensitivity or fitness costs, is invaluable for the design of rational evolution-based, ecology-based treatment approaches. In this Review, we discuss these evolutionary trade-offs and how such knowledge can aid in informing combination or alternating antibiotic therapies against bacterial infections. In addition, we discuss how targeting bacterial metabolism can enhance drug activity and impair antibiotic resistance evolution. Finally, we explore how an improved understanding of the original physiological function of antibiotic resistance determinants, which have evolved to reach clinical resistance after a process of historical contingency, may help to tackle antibiotic resistance.

Shifts in Mutation Bias Promote Mutators by Altering the Distribution of Fitness Effects

AbstractRecent experimental evidence demonstrates that shifts in mutational biases-for example, increases in transversion frequency-can change the distribution of fitness effects of mutations (DFE). In particular, reducing or reversing a prevailing bias can increase the probability that a de novo mutation is beneficial. It has also been shown that mutator bacteria are more likely to emerge if the beneficial mutations they generate have a larger effect size than observed in the wild type. Here, we connect these two results, demonstrating that mutator strains that reduce or reverse a prevailing bias have a positively shifted DFE, which in turn can dramatically increase their emergence probability. Since changes in mutation rate and bias are often coupled through the gain and loss of DNA repair enzymes, our results predict that the invasion of mutator strains will be facilitated by shifts in mutation bias that offer improved access to previously undersampled beneficial mutations.

Shifts in mutation spectra enhance access to beneficial mutations

Biased mutation spectra are pervasive, with wide variation in the magnitude of mutational biases that influence genome evolution and adaptation. How do such diverse biases evolve? Our experiments show that changing the mutation spectrum allows populations to sample previously undersampled mutational space, including beneficial mutations. The resulting shift in the distribution of fitness effects is advantageous: Beneficial mutation supply and beneficial pleiotropy both increase, while deleterious load reduces. More broadly, simulations indicate that reducing or reversing the direction of a long-term bias is always selectively favored. Such changes in mutation bias can occur easily via altered function of DNA repair genes. A phylogenetic analysis shows that these genes are repeatedly gained and lost in bacterial lineages, leading to frequent bias shifts in opposite directions. Thus, shifts in mutation spectra may evolve under selection and can directly alter the outcome of adaptive evolution by facilitating access to beneficial mutations.

Distribution of mutation rates challenges evolutionary predictability

Natural selection is commonly assumed to act on extensive standing genetic variation. Yet, accumulating evidence highlights the role of mutational processes creating this genetic variation: to become evolutionarily successful, adaptive mutants must not only reach fixation, but also emerge in the first place, i.e. have a high enough mutation rate. Here, we use numerical simulations to investigate how mutational biases impact our ability to observe rare mutational pathways in the laboratory and to predict outcomes in experimental evolution. We show that unevenness in the rates at which mutational pathways produce adaptive mutants means that most experimental studies lack power to directly observe the full range of adaptive mutations. Modelling mutation rates as a distribution, we show that a substantially larger target size ensures that a pathway mutates more commonly. Therefore, we predict that commonly mutated pathways are conserved between closely related species, but not rarely mutated pathways. This approach formalizes our proposal that most mutations have a lower mutation rate than the average mutation rate measured experimentally. We suggest that the extent of genetic variation is overestimated when based on the average mutation rate.

Targeted Screening for Spontaneous Insertion Mutations in a Lactic Acid Bacterium, Tetragenococcus halophilus

Studies on the microorganisms used in food production are of interest because microbial genotypes are reflected in food qualities such as taste, flavor, and yield. However, several microbes are nonmodel organisms, and their analysis is often limited by the lack of genetic tools. Tetragenococcus halophilus, a halophilic lactic acid bacterium used in soy sauce fermentation starter culture, is one such microorganism. The lack of DNA transformation techniques for T. halophilus makes gene complementation and disruption assays difficult. Here, we report that the endogenous insertion sequence ISTeha4, belonging to the IS4 family, is translocated at an extremely high frequency in T. halophilus and causes insertional mutations at various loci. We developed a method named targeting spontaneous insertional mutations in genomes (TIMING), which combines high-frequency insertional mutations and efficient PCR screening, enabling the isolation of gene mutants of interest from a library. The method provides a reverse genetics and strain improvement tool, does not require the introduction of exogenous DNA constructs, and enables the analysis of nonmodel microorganisms lacking DNA transformation techniques. Our results highlight the important role of insertion sequences as a source of spontaneous mutagenesis and genetic diversity in bacteria. IMPORTANCE Genetic and strain improvement tools to manipulate a gene of interest are required for the nontransformable lactic acid bacterium Tetragenococcus halophilus. Here, we demonstrate that an endogenous transposable element, ISTeha4, is transposed into the host genome at an extremely high frequency. A genotype-based and non-genetically engineered screening system was constructed to isolate knockout mutants using this transposable element. The method described enables a better understanding of the genotype-phenotype relationship and serves as a tool to develop food-grade-appropriate mutants of T. halophilus.

Intermittent antibiotic treatment of bacterial biofilms favors the rapid evolution of resistance

Bacterial antibiotic resistance is a global health concern of increasing importance and intensive study. Although biofilms are a common source of infections in clinical settings, little is known about the development of antibiotic resistance within biofilms. Here, we use experimental evolution to compare selection of resistance mutations in planktonic and biofilm Escherichia coli populations exposed to clinically relevant cycles of lethal treatment with the aminoglycoside amikacin. Consistently, mutations in sbmA, encoding an inner membrane peptide transporter, and fusA, encoding the essential elongation factor G, are rapidly selected in biofilms, but not in planktonic cells. This is due to a combination of enhanced mutation rate, increased adhesion capacity and protective biofilm-associated tolerance. These results show that the biofilm environment favors rapid evolution of resistance and provide new insights into the dynamic evolution of antibiotic resistance in biofilms.

Mutation Rate and Spectrum of the Silkworm in Normal and Temperature Stress Conditions

Mutation rate is a crucial parameter in evolutionary genetics. However, the mutation rate of most species as well as the extent to which the environment can alter the genome of multicellular organisms remain poorly understood. Here, we used parents–progeny sequencing to investigate the mutation rate and spectrum of the domestic silkworm (Bombyx mori) among normal and two temperature stress conditions (32 °C and 0 °C). The rate of single-nucleotide mutations in the normal temperature rearing condition was 0.41 × 10−8 (95% confidence interval, 0.33 × 10−8–0.49 × 10−8) per site per generation, which was up to 1.5-fold higher than in four previously studied insects. Moreover, the mutation rates of the silkworm under the stresses are significantly higher than in normal conditions. Furthermore, the mutation rate varies less in gene regions under normal and temperature stresses. Together, these findings expand the known diversity of the mutation rate among eukaryotes but also have implications for evolutionary analysis that assumes a constant mutation rate among species and environments.

Biofilm antimicrobial susceptibility through an experimental evolutionary lens

Experimental evolution experiments in which bacterial populations are repeatedly exposed to an antimicrobial treatment, and examination of the genotype and phenotype of the resulting evolved bacteria, can help shed light on mechanisms behind reduced susceptibility. In this review we present an overview of why it is important to include biofilms in experimental evolution, which approaches are available to study experimental evolution in biofilms and what experimental evolution has taught us about tolerance and resistance in biofilms. Finally, we present an emerging consensus view on biofilm antimicrobial susceptibility supported by data obtained during experimental evolution studies.

A mitochondria-specific mutational signature of aging: increased rate of A > G substitutions on the heavy strand

The mutational spectrum of the mitochondrial DNA (mtDNA) does not resemble any of the known mutational signatures of the nuclear genome and variation in mtDNA mutational spectra between different organisms is still incomprehensible. Since mitochondria are responsible for aerobic respiration, it is expected that mtDNA mutational spectrum is affected by oxidative damage. Assuming that oxidative damage increases with age, we analyse mtDNA mutagenesis of different species in regards to their generation length. Analysing, (i) dozens of thousands of somatic mtDNA mutations in samples of different ages (ii) 70053 polymorphic synonymous mtDNA substitutions reconstructed in 424 mammalian species with different generation lengths and (iii) synonymous nucleotide content of 650 complete mitochondrial genomes of mammalian species we observed that the frequency of AH > GH substitutions (H: heavy strand notation) is twice bigger in species with high versus low generation length making their mtDNA more AH poor and GH rich. Considering that AH > GH substitutions are also sensitive to the time spent single-stranded (TSSS) during asynchronous mtDNA replication we demonstrated that AH > GH substitution rate is a function of both species-specific generation length and position-specific TSSS. We propose that AH > GH is a mitochondria-specific signature of oxidative damage associated with both aging and TSSS.

Salt stress alters the spectrum of de novo mutation available to selection during experimental adaptation of Chlamydomonas reinhardtii

The genetic basis of adaptation is driven by both selection and the spectrum of available mutations. Given that the rate of mutation is not uniformly distributed across the genome and varies depending on the environment, understanding the signatures of selection across the genome is aided by first establishing what the expectations of genetic change are from mutation. To determine the interaction between salt stress, selection, and mutation across the genome, we compared mutations observed in a selection experiment for salt tolerance in Chlamydomonas reinhardtii to those observed in mutation accumulation (MA) experiments with and without salt exposure. MA lines evolved under salt stress had a single-nucleotide mutation rate of 1.1 × 10 - 9 $1.1 \times 10^{-9}$ , similar to that of MA lines under standard conditions ( 9.6 × 10 - 10 $9.6 \times 10^{-10}$ ). However, we found that salt stress led to an increased rate of indel mutations, but that many of these mutations were removed under selection. Finally, lines adapted to salt also showed excess clustering of mutations in the genome and the co-expression network, suggesting a role for positive selection in retaining mutations in particular compartments of the genome during the evolution of salt tolerance. Our study shows that characterizing mutation rates and spectra expected under stress helps disentangle the effects of environment and selection during adaptation.

A Highly Glyphosate-Resistant EPSPS Mutant from Laboratory Evolution

EPSP synthase is the target enzyme of glyphosate herbicides. Due to the extensive use of glyphosate, it is very important to obtain EPSPS genes with high glyphosate resistance for the development of transgenic crops. GR79-EPSPS is a class I EPSP synthase with certain glyphosate resistance isolated from glyphosate-contaminated soil. After more than 1000 generations, a Y40I substitution was identified, and the enzyme had a nearly 1.8-fold decrease in Km [PEP] and a 1.7-fold increase in Ki[glyphosate] compared to the wild-type enzyme. Enzyme dynamics and molecular dynamics analysis showed that the substitution was near the hinge region of EPSPS, and the affinity of glyphosate binding to amino acid residues of the active site decreased due to Y40I substitution, resulting in an increase in glyphosate resistance. These results provide more evidence for the combination of directed evolution and rational design of protein engineering.

Prochlorococcus have low global mutation rate and small effective population size

Prochlorococcus are the most abundant free-living photosynthetic carbon-fixing organisms in the ocean. Prochlorococcus show small genome sizes, low genomic G+C content, reduced DNA repair gene pool and fast evolutionary rates, which are typical features of endosymbiotic bacteria. Nevertheless, their evolutionary mechanisms are believed to be different. Evolution of endosymbiotic bacteria is dominated by genetic drift owing to repeated population bottlenecks, whereas Prochlorococcus are postulated to have extremely large effective population sizes (N_e) and thus drift has rarely been considered. However, accurately extrapolating N_e requires measuring an unbiased global mutation rate through mutation accumulation, which is challenging for Prochlorococcus. Here, we managed this experiment over 1,065 days using Prochlorococcus marinus AS9601, sequenced genomes of 141 mutant lines and determined its mutation rate to be 3.50 × 10^-10 per site per generation. Extrapolating N_e additionally requires identifying population boundaries, which we defined using PopCOGenT and over 400 genomes related to AS9601. Accordingly, we calculated its N_e to be 1.68 × 10⁷, which is only reasonably greater than that of endosymbiotic bacteria but surprisingly smaller than that of many free-living bacteria extrapolated using the same approach. Our results therefore suggest that genetic drift is a key driver of Prochlorococcus evolution.

Evolutionary Dynamics of Asexual Hypermutators Adapting to a Novel Environment

How microbes adapt to a novel environment is a central question in evolutionary biology. Although adaptive evolution must be fueled by beneficial mutations, whether higher mutation rates facilitate the rate of adaptive evolution remains unclear. To address this question, we cultured Escherichia coli hypermutating populations, in which a defective methyl-directed mismatch repair pathway causes a 140-fold increase in single-nucleotide mutation rates. In parallel with wild-type E. coli, populations were cultured in tubes containing Luria-Bertani broth, a complex medium known to promote the evolution of subpopulation structure. After 900 days of evolution, in three transfer schemes with different population-size bottlenecks, hypermutators always exhibited similar levels of improved fitness as controls. Fluctuation tests revealed that the mutation rates of hypermutator lines converged evolutionarily on those of wild-type populations, which may have contributed to the absence of fitness differences. Further genome-sequence analysis revealed that, although hypermutator populations have higher rates of genomic evolution, this largely reflects strong genetic linkage. Despite these linkage effects, the evolved population exhibits parallelism in fixed mutations, including those potentially related to biofilm formation, transcription regulation, and mutation-rate evolution. Together, these results are generally inconsistent with a hypothesized positive relationship between the mutation rate and the adaptive speed of evolution, and provide insight into how clonal adaptation occurs in novel environments.

Choice of β-Lactam Resistance Pathway Depends Critically on Initial Antibiotic Concentration

Antibiotic resistance trajectories with different final resistance may critically depend on the first mutation, due to epistatic interactions. Here, we study the effect of mutation bias and the concentration-dependent effects on fitness of two clinically important mutations in TEM-1 β-lactamase in initiating alternative trajectories to cefotaxime resistance. We show that at low cefotaxime concentrations, the R164S mutation (a mutation of arginine to serine at position 164), which confers relatively low resistance, is competitively superior to the G238S mutation, conferring higher resistance, thus highlighting a critical influence of antibiotic concentration on long-term resistance evolution.

The Physics-Biology continuum challenges darwinism: Evolution is directed by the homeostasis-dependent bidirectional relation between genome and phenotype

The physics-biology continuum relies on the fact that life emerged from prebiotic molecules. Here, I argue that life emerged from the coupling between nucleic acid and protein synthesis during which proteins (or proto-phenotypes) maintained the physicochemical parameter equilibria (or proto-homeostasis) in the proximity of their encoding nucleic acids (or proto-genomes). This protected the proto-genome physicochemical integrity (i.e., atomic composition) from environmental physicochemical constraints, and therefore increased the probability of reproducing the proto-genome without variation. From there, genomes evolved depending on the biological activities they generated in response to environmental fluctuations. Thus, a genome maintaining homeostasis (i.e., internal physicochemical parameter equilibria), despite and in response to environmental fluctuations, maintains its physicochemical integrity and has therefore a higher probability to be reproduced without variation. Consequently, descendants have a higher probability to share the same phenotype than their parents. Otherwise, the genome is modified during replication as a consequence of the imbalance of the internal physicochemical parameters it generates, until new mutation-deriving biological activities maintain homeostasis in offspring. In summary, evolution depends on feedforward and feedback loops between genome and phenotype, as the internal physicochemical conditions that a genome generates ─ through its derived phenotype in response to environmental fluctuations ─ in turn either guarantee its stability or direct its variation. Evolution may not be explained by the Darwinism-derived, unidirectional principle (random mutations-phenotypes-natural selection) but rather by the bidirectional relationship between genome and phenotype, in which the phenotype in interaction with the environment directs the evolution of the genome it derives from.

Maintenance of Metabolic Plasticity despite Relaxed Selection in a Long-Term Evolution Experiment with Escherichia coli

AbstractTraits that are unused in a given environment are subject to processes that tend to erode them, leading to reduced fitness in other environments. Although this general tendency is clear, we know much less about why some traits are lost while others are retained and about the roles of mutation and selection in generating different responses. We addressed these issues by examining populations of a facultative anaerobe, Escherichia coli, that have evolved for >30 years in the presence of oxygen, with relaxed selection for anaerobic growth and the associated metabolic plasticity. We asked whether evolution led to the loss, improvement, or maintenance of anaerobic growth, and we analyzed gene expression and mutational data sets to understand the outcomes. We identified genomic signatures of both positive and purifying selection on aerobic-specific genes, while anaerobic-specific genes showed clear evidence of relaxed selection. We also found parallel evolution at two interacting loci that regulate anaerobic growth. We competed the ancestor and evolved clones from each population in an anoxic environment, and we found that anaerobic fitness had not decayed, despite relaxed selection. In summary, relaxed selection does not necessarily reduce an organism's fitness in other environments. Instead, the genetic architecture of the traits under relaxed selection and their correlations with traits under positive and purifying selection may sometimes determine evolutionary outcomes.

Biodiversity-based development and evolution: the emerging research systems in model and non-model organisms

Evolutionary developmental biology, or Evo-Devo for short, has become an established field that, broadly speaking, seeks to understand how changes in development drive major transitions and innovation in organismal evolution. It does so via integrating the principles and methods of many subdisciplines of biology. Although we have gained unprecedented knowledge from the studies on model organisms in the past decades, many fundamental and crucially essential processes remain a mystery. Considering the tremendous biodiversity of our planet, the current model organisms seem insufficient for us to understand the evolutionary and physiological processes of life and its adaptation to exterior environments. The currently increasing genomic data and the recently available gene-editing tools make it possible to extend our studies to non-model organisms. In this review, we review the recent work on the regulatory signaling of developmental and regeneration processes, environmental adaptation, and evolutionary mechanisms using both the existing model animals such as zebrafish and Drosophila, and the emerging nonstandard model organisms including amphioxus, ascidian, ciliates, single-celled phytoplankton, and marine nematode. In addition, the challenging questions and new directions in these systems are outlined as well.

Streptomycin and nalidixic acid elevate the spontaneous genome-wide mutation rate in Escherichia coli

Since antibiotic resistance is a growing public health problem worldwide, it is important to understand how antibiotics and spontaneous mutations cooperate and shape the genome-wide mutation rate and spectrum. Here, we quantitatively evaluate genome-wide mutational profiles of Escherichia coli after long-term subinhibitory exposure to a broad-spectrum (streptomycin) and a narrow-spectrum antibiotic (nalidixic acid), using a mutation accumulation design combined with whole-genome resequencing of replicate lines as a mutagenicity test. We determined that, while the genome-wide mutation rate is slightly higher in the streptomycin-treated lines compared to the control lines, there is a significant increase in the nalidixic acid-treated lines. Our findings suggest that both broad and narrow-spectrum antibiotics may elevate the mutation rates in E. coli, but mechanisms of action may affect the consequence, thus contribute to accelerating the rate of adaptation and conferring antibiotic resistance.

Warmer temperatures enhance beneficial mutation effects

Temperature determines the rates of all biochemical and biophysical processes, and is also believed to be a key driver of macroevolutionary patterns. It is suggested that physiological constraints at low temperatures may diminish the fitness advantages of otherwise beneficial mutations; by contrast, relatively high, benign, temperatures allow beneficial mutations to efficiently show their phenotypic effects. To experimentally test this "mutational effects" mechanism, we examined the fitness effects of mutations across a temperature gradient using bacterial genotypes from the early stage of a mutation accumulation experiment with Escherichia coli. While the incidence of beneficial mutations did not significantly change across environmental temperatures, the number of mutations that conferred strong beneficial fitness effects was greater at higher temperatures. The results therefore support the hypothesis that warmer temperatures increase the chance and magnitude of positive selection, with implications for explaining the geographic patterns in evolutionary rates and understanding contemporary evolution under global warming.

Short- and Long-Term Transcriptomic Responses of Escherichia coli to Biocides: a Systems Analysis

The mechanisms of the bacterial response to biocides are poorly understood, despite their broad application. To identify the genetic basis and pathways implicated in the biocide stress response, we exposed Escherichia coli populations to 10 ubiquitous biocides. By comparing the transcriptional responses between a short-term exposure (30 min) and a long-term exposure (8 to 12 h) to biocide stress, we established the common gene and pathway clusters that are implicated in general and biocide-specific stress responses. Our analysis revealed a temporal choreography, starting from the upregulation of chaperones to the subsequent repression of motility and chemotaxis pathways and the induction of an anaerobic pool of enzymes and biofilm regulators. A systematic analysis of the transcriptional data identified a zur-regulated gene cluster to be highly active in the stress response against sodium hypochlorite and peracetic acid, presenting a link between the biocide stress response and zinc homeostasis. Susceptibility assays with knockout mutants further validated our findings and provide clear targets for downstream investigation of the implicated mechanisms of action.IMPORTANCE Antiseptics and disinfectant products are of great importance to control and eliminate pathogens, especially in settings such as hospitals and the food industry. Such products are widely distributed and frequently poorly regulated. Occasional outbreaks have been associated with microbes resistant to such compounds, and researchers have indicated potential cross-resistance with antibiotics. Despite that, there are many gaps in knowledge about the bacterial stress response and the mechanisms of microbial resistance to antiseptics and disinfectants. We investigated the stress response of the bacterium Escherichia coli to 10 common disinfectant and antiseptic chemicals to shed light on the potential mechanisms of tolerance to such compounds.

Dynamics of genetic variation in transcription factors and its implications for the evolution of regulatory networks in Bacteria

The evolution of regulatory networks in Bacteria has largely been explained at macroevolutionary scales through lateral gene transfer and gene duplication. Transcription factors (TF) have been found to be less conserved across species than their target genes (TG). This would be expected if TFs accumulate mutations faster than TGs. This hypothesis is supported by several lab evolution studies which found TFs, especially global regulators, to be frequently mutated. Despite these studies, the contribution of point mutations in TFs to the evolution of regulatory network is poorly understood. We tested if TFs show greater genetic variation than their TGs using whole-genome sequencing data from a large collection of Escherichia coli isolates. TFs were less diverse than their TGs across natural isolates, with TFs of large regulons being more conserved. In contrast, TFs showed higher mutation frequency in adaptive laboratory evolution experiments. However, over long-term laboratory evolution spanning 60 000 generations, mutation frequency in TFs gradually declined after a rapid initial burst. Extrapolating the dynamics of genetic variation from long-term laboratory evolution to natural populations, we propose that point mutations, conferring large-scale gene expression changes, may drive the early stages of adaptation but gene regulation is subjected to stronger purifying selection post adaptation.

Prolonged culture in aerobic environments alters Escherichia coli H2 production capacity

Growing interest in renewable energy continues to motivate new work on microbial biohydrogen production and in particular utilizing Escherichia coli a well-studied, facultative anaerobe. Here we characterize, for the first time the H2 production rate and capacity, of E coli isolates from the 50 000th generation of the Long-Term Evolution Experiment. Under these reaction conditions, peak production rates near or above 5 mL per hour for 100 mL of lysogeny broth (LB media) was established for the ancestral strains and batch efficiencies between 0.15 and 0.22 mL H2 produced per 1 mL LB media were achieved. All 11 isolates studied, which had been aerobically cultured in minimal media since 1988, exhibited a decreased H2 production rate or capacity with many strains unable to grow under anaerobic conditions at all. The genomes of these strains have been sequenced and a preliminary analysis of the correlations between genotype and phenotype shows that mutations in gene ydjO are exclusively observed in the two isolates which produce H2, potentially suggesting a role for this gene in the maintenance of wild type metabolic pathways in the context of diverse mutational backgrounds. These results provide hints towards uncovering new genetic targets for the pursuit of bacterial strains with increased capacity for H2 production as well as a case study in speciation and the control of phenotypic switching.

Physicochemical Foundations of Life that Direct Evolution: Chance and Natural Selection are not Evolutionary Driving Forces

The current framework of evolutionary theory postulates that evolution relies on random mutations generating a diversity of phenotypes on which natural selection acts. This framework was established using a top-down approach as it originated from Darwinism, which is based on observations made of complex multicellular organisms and, then, modified to fit a DNA-centric view. In this article, it is argued that based on a bottom-up approach starting from the physicochemical properties of nucleic and amino acid polymers, we should reject the facts that (i) natural selection plays a dominant role in evolution and (ii) the probability of mutations is independent of the generated phenotype. It is shown that the adaptation of a phenotype to an environment does not correspond to organism fitness, but rather corresponds to maintaining the genome stability and integrity. In a stable environment, the phenotype maintains the stability of its originating genome and both (genome and phenotype) are reproduced identically. In an unstable environment (i.e., corresponding to variations in physicochemical parameters above a physiological range), the phenotype no longer maintains the stability of its originating genome, but instead influences its variations. Indeed, environment- and cellular-dependent physicochemical parameters define the probability of mutations in terms of frequency, nature, and location in a genome. Evolution is non-deterministic because it relies on probabilistic physicochemical rules, and evolution is driven by a bidirectional interplay between genome and phenotype in which the phenotype ensures the stability of its originating genome in a cellular and environmental physicochemical parameter-depending manner.

Irregularities in genetic variation and mutation rates with environmental stresses

The appearance of new mutations is determined by the equilibrium between DNA error formation and repair. In bacteria like Escherichia coli, stresses are thought shift this balance towards increased mutagenesis. Recent findings, however, suggest a very uneven relationship between stress and mutations. Only a subset of stressful environments increase the net rate of mutation and different forms of nutritional stress (such as oxygen, carbon or phosphorus limitations) result in markedly different mutation rates after similar reductions in growth rate. Moreover, different stresses result in altered mutational spectra, with some increasing transposition and others increasing indel formation. Single-base substitution rates are lower with some stresses than in unstressed bacteria. Indeed, changes to the mix of mutations with stress are more widespread than a marked increase in net mutation rate. Much remains to be learned on how environments have unique mutational signatures and why some stresses are more mutagenic than others. Even beyond stress-induced genetic variation, the fundamental unresolved question in the stress-mutation relationship is the adaptive value of different types of mutations and mutation rates; is transposition, for example, more advantageous under anaerobic conditions? It remains to be investigated whether stress-specific genetic variation impacts on evolvability differentially in distinct environments.

Mutation bias and GC content shape antimutator invasions

Mutators represent a successful strategy in rapidly adapting asexual populations, but theory predicts their eventual extinction due to their unsustainably large deleterious load. While antimutator invasions have been documented experimentally, important discrepancies among studies remain currently unexplained. Here we show that a largely neglected factor, the mutational idiosyncrasy displayed by different mutators, can play a major role in this process. Analysing phylogenetically diverse bacteria, we find marked and systematic differences in the protein-disruptive effects of mutations caused by different mutators in species with different GC compositions. Computer simulations show that these differences can account for order-of-magnitude changes in antimutator fitness for a realistic range of parameters. Overall, our results suggest that antimutator dynamics may be highly dependent on the specific genetic, ecological and evolutionary history of a given population. This context-dependency further complicates our understanding of mutators in clinical settings, as well as their role in shaping bacterial genome size and composition.

Rapid cell division of Staphylococcus aureus during colonization of the human nose

BACKGROUND

Staphylococcus aureus is an important opportunistic pathogen and a commensal bacterium, thriving in the nasal cavities of 20% of the human population. Little is known about the dynamics of asymptomatic colonization and the occasional transition to infectious disease.

RESULTS

In this study, we inferred that S. aureus cells replicate every one to three hours on average while colonizing the human nose, based on two independent lines of genomic evidence. First, we collected nasal swab samples from human subjects, extracted and sequenced metagenomic DNA, and analyzed the distribution of sequencing coverage along the staphylococcal chromosome. Calibration of this data by comparison to a laboratory culture enabled measuring S. aureus cell division rates in nasal samples. Second, we applied mutation accumulation experiments paired with genome sequencing to measure spontaneous mutation rates at a genome scale. Relating these mutation rates to annual evolutionary rates confirmed that nasal S. aureus continuously pass several thousand cell divisions per year when averaged over large, globally distributed populations and over many years, corresponding to generation times of less than two hours.

CONCLUSIONS

The cell division rates we determined were higher than the fastest documented rates during fulminant disease progression (in a mouse model of systemic infection) and much higher than those previously measured in expectorated sputum from cystic fibrosis patients. This paper supplies absolute in-vivo generation times for an important bacterial commensal, indicating that colonization of the human upper respiratory tract is characterized by a highly dynamic equilibrium between bacterial growth and removal.

Predicting mutational routes to new adaptive phenotypes

Predicting evolutionary change poses numerous challenges. Here we take advantage of the model bacterium Pseudomonas fluorescens in which the genotype-to-phenotype map determining evolution of the adaptive 'wrinkly spreader' (WS) type is known. We present mathematical descriptions of three necessary regulatory pathways and use these to predict both the rate at which each mutational route is used and the expected mutational targets. To test predictions, mutation rates and targets were determined for each pathway. Unanticipated mutational hotspots caused experimental observations to depart from predictions but additional data led to refined models. A mismatch was observed between the spectra of WS-causing mutations obtained with and without selection due to low fitness of previously undetected WS-causing mutations. Our findings contribute toward the development of mechanistic models for forecasting evolution, highlight current limitations, and draw attention to challenges in predicting locus-specific mutational biases and fitness effects.

Old Trade, New Tricks: Insights into the Spontaneous Mutation Process from the Partnering of Classical Mutation Accumulation Experiments with High-Throughput Genomic Approaches

Mutations spawn genetic variation which, in turn, fuels evolution. Hence, experimental investigations into the rate and fitness effects of spontaneous mutations are central to the study of evolution. Mutation accumulation (MA) experiments have served as a cornerstone for furthering our understanding of spontaneous mutations for four decades. In the pregenomic era, phenotypic measurements of fitness-related traits in MA lines were used to indirectly estimate key mutational parameters, such as the genomic mutation rate, new mutational variance per generation, and the average fitness effect of mutations. Rapidly emerging next-generating sequencing technology has supplanted this phenotype-dependent approach, enabling direct empirical estimates of the mutation rate and a more nuanced understanding of the relative contributions of different classes of mutations to the standing genetic variation. Whole-genome sequencing of MA lines bears immense potential to provide a unified account of the evolutionary process at multiple levels-the genetic basis of variation, and the evolutionary dynamics of mutations under the forces of selection and drift. In this review, we have attempted to synthesize key insights into the spontaneous mutation process that are rapidly emerging from the partnering of classical MA experiments with high-throughput sequencing, with particular emphasis on the spontaneous rates and molecular properties of different mutational classes in nuclear and mitochondrial genomes of diverse taxa, the contribution of mutations to the evolution of gene expression, and the rate and stability of transgenerational epigenetic modifications. Future advances in sequencing technologies will enable greater species representation to further refine our understanding of mutational parameters and their functional consequences.

The impact of growth rate and environmental factors on mutation rates and spectra in Escherichia coli

Genetic variation in bacterial populations is remarkably sensitive to environmental influences, including simple, nutritional differences. Not only the rate but also the kind of mutational changes is biased by the nutritional state of bacteria. Here we investigate the mutational consequences of a universal variable for free-living bacteria, namely the growth rate. By controlling growth in chemostats, the rate and mix of mutations was investigated for populations of Escherichia coli subject to different specific growth rates. Both aerobic and anaerobic cultures were compared with see if growth rate is a factor in the commonest respiratory conditions for E. coli. We find mutation rates are raised markedly with decreasing growth rate. Base pair substitutions and 1-bp insertions and deletions increase with reduced growth rate, but less so in anaerobic cultures. Insertion sequence movements are particularly sensitive to growth rate, with IS2 being optimal at intermediate growth rates whereas IS1 and IS150 movements are highest at the slowest tested growth rate. A comprehensive comparison of growth rate effects, as well as six other environmental factors, provides the most complete picture yet of the range of mutational signatures in bacterial genetic variation.

Escherichia coli mutation rates and spectra with combinations of environmental limitations

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

Temperature responses of mutation rate and mutational spectrum in an Escherichia coli strain and the correlation with metabolic rate

Background

Temperature is a major determinant of spontaneous mutation, but the precise mode, and the underlying mechanisms, of the temperature influences remain less clear. Here we used a mutation accumulation approach combined with whole-genome sequencing to investigate the temperature dependence of spontaneous mutation in an Escherichia coli strain. Experiments were performed under aerobic conditions at 25, 28 and 37 °C, three temperatures that were non-stressful for the bacterium but caused significantly different bacterial growth rates.

Results

Mutation rate did not differ between 25 and 28 °C, but was higher at 37 °C. Detailed analyses of the molecular spectrum of mutations were performed; and a particularly interesting finding is that higher temperature led to a bias of mutation to coding, relative to noncoding, DNA. Furthermore, the temperature response of mutation rate was extremely similar to that of metabolic rate, consistent with an idea that metabolic rate predicts mutation rate.

Conclusions

Temperature affects mutation rate and the types of mutation supply, both being crucial for the opportunity of natural selection. Our results help understand how temperature drives evolutionary speed of organisms and thus the global patterns of biodiversity. This study also lend support to the metabolic theory of ecology for linking metabolic rate and molecular evolution rate.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1252-8) contains supplementary material, which is available to authorized users.

Insertion sequence elements-mediated structural variations in bacterial genomes

Mobile genetic elements (MGEs) impact the evolution and stability of their host genomes. Insertion sequence (IS) elements are the most common MGEs in bacterial genomes and play a crucial role in mediating large-scale variations in bacterial genomes. It is understood that IS elements and MGEs in general coexist in a dynamical equilibrium with their respective hosts. Current studies indicate that the spontaneous movement of IS elements does not follow a constant rate in different bacterial genomes. However, due to the paucity and sparsity of the data, these observations are yet to be conclusive. In this paper, we conducted a comparative analysis of the IS-mediated genome structural variations in ten mutation accumulation (MA) experiments across eight strains of five bacterial species containing IS elements, including four strains of the E. coli. We used GRASPER algorithm, a denovo structural variation (SV) identification algorithm designed to detect SVs involving repetitive sequences in the genome. We observed highly diverse rates of IS insertions and IS-mediated recombinations across different bacterial species as well as across different strains of the same bacterial species. We also observed different rates of the elements from the same IS family in different bacterial genomes, suggesting that the distinction in rates might not be due to the different composition of IS elements across bacterial genomes.

Evolution of Antibiotic Resistance in Biofilm and Planktonic Pseudomonas aeruginosa Populations Exposed to Subinhibitory Levels of Ciprofloxacin

The opportunistic Gram-negative pathogen Pseudomonas aeruginosa, known for its intrinsic and acquired antibiotic resistance, has a notorious ability to form biofilms, which often facilitate chronic infections. The evolutionary paths to antibiotic resistance have mainly been investigated in planktonic cultures and are less studied in biofilms. We experimentally evolved P. aeruginosa PAO1 colony biofilms and stationary-phase planktonic cultures for seven passages in the presence of subinhibitory levels (0.1 mg/liter) of ciprofloxacin (CIP) and performed a genotypic (whole-bacterial population sequencing) and phenotypic assessment of the populations. We observed a higher proportion of CIP resistance in the CIP-evolved biofilm populations than in planktonic populations exposed to the same drug concentrations. However, the MICs of ciprofloxacin were lower in CIP-resistant isolates selected from the biofilm population than the MICs of CIP-resistant isolates from the planktonic cultures. We found common evolutionary trajectories between the different lineages, with mutations in known CIP resistance determinants as well as growth condition-dependent adaptations. We observed a general trend toward a reduction in type IV-pilus-dependent motility (twitching) in CIP-evolved populations and a loss of virulence-associated traits in the populations evolved in the absence of antibiotic. In conclusion, our data indicate that biofilms facilitate the development of low-level mutational resistance, probably due to the lower effective drug exposure than in planktonic cultures. These results provide a framework for the selection process of resistant variants and the evolutionary mechanisms involved under the two different growth conditions.

The distribution of bacterial doubling times in the wild

Generation time varies widely across organisms and is an important factor in the life cycle, life history and evolution of organisms. Although the doubling time (DT) has been estimated for many bacteria in the laboratory, it is nearly impossible to directly measure it in the natural environment. However, an estimate can be obtained by measuring the rate at which bacteria accumulate mutations per year in the wild and the rate at which they mutate per generation in the laboratory. If we assume the mutation rate per generation is the same in the wild and in the laboratory, and that all mutations in the wild are neutral, an assumption that we show is not very important, then an estimate of the DT can be obtained by dividing the latter by the former. We estimate the DT for five species of bacteria for which we have both an accumulation and a mutation rate estimate. We also infer the distribution of DTs across all bacteria from the distribution of the accumulation and mutation rates. Both analyses suggest that DTs for bacteria in the wild are substantially greater than those in the laboratory, that they vary by orders of magnitude between different species of bacteria and that a substantial fraction of bacteria double very slowly in the wild.

Genomic and Transcriptomic Insights into How Bacteria Withstand High Concentrations of Benzalkonium Chloride Biocides

Benzalkonium chlorides (BAC) are commonly used biocides in broad-spectrum disinfectant solutions. How microorganisms cope with BAC exposure remains poorly understood, despite its importance for disinfection and disinfectant-induced antibiotic resistance. To provide insights into these issues, we exposed two isolates of an opportunistic pathogen, Pseudomonas aeruginosa, to increasing concentrations of BAC. One isolate was preadapted to BAC, as it originated from a bioreactor fed with subinhibitory concentrations of BAC for 3 years, while the other originated from a bioreactor that received no BAC. Replicated populations of both isolates were able to survive high concentrations of BAC, up to 1,200 and 1,600 mg/liter for the non- and preadapted strains, respectively, exceeding typical application doses. Transcriptome sequencing (RNA-seq) analysis revealed upregulation of efflux pump genes and decreased expression of porins related to BAC transport as well as reduced growth rate. Increased expression of spermidine (a polycation) synthase genes and mutations in the pmrB (polymyxin resistance) gene, which cause a reduction in membrane negative charge, suggested that a major adaptation to exposure to the cationic surfactant BAC was to actively stabilize cell surface charge. Collectively, these results revealed that P. aeruginosa adapts to BAC exposure by a combination of mechanisms and provided genetic markers to monitor BAC-resistant organisms that may have applications in the practice of disinfection.IMPORTANCE BAC are widely used as biocides in disinfectant solutions, food-processing lines, domestic households, and health care facilities. Due to their wide use and mode of action, there has been rising concern that BAC may promote antibiotic resistance. Consistent with this idea, at least 40 outbreaks have been attributed to infection by disinfectant- and antibiotic-resistant pathogens such as P. aeruginosa However, the underlying molecular mechanisms that bacteria use to deal with BAC exposure remain poorly elucidated. Elucidating these mechanisms may be important for monitoring and limiting the spread of disinfectant-resistant pathogens. Using an integrated approach that combined genomics and transcriptomics with physiological characterization of BAC-adapted isolates, this study provided a comprehensive understanding of the BAC resistance mechanisms in P. aeruginosa Our findings also revealed potential genetic markers to detect and monitor the abundance of BAC-resistant pathogens across clinical or environmental settings. This work contributes new knowledge about high concentrations of benzalkonium chlorides disinfectants-resistance mechanisms at the whole-cell genomic and transcriptomic level.

Sources of spontaneous mutagenesis in bacteria

Mutations in an organism's genome can arise spontaneously, that is, in the absence of exogenous stress and prior to selection. Mutations are often neutral or deleterious to individual fitness but can also provide genetic diversity driving evolution. Mutagenesis in bacteria contributes to the already serious and growing problem of antibiotic resistance. However, the negative impacts of spontaneous mutagenesis on human health are not limited to bacterial antibiotic resistance. Spontaneous mutations also underlie tumorigenesis and evolution of drug resistance. To better understand the causes of genetic change and how they may be manipulated in order to curb antibiotic resistance or the development of cancer, we must acquire a mechanistic understanding of the major sources of mutagenesis. Bacterial systems are particularly well-suited to studying mutagenesis because of their fast growth rate and the panoply of available experimental tools, but efforts to understand mutagenic mechanisms can be complicated by the experimental system employed. Here, we review our current understanding of mutagenic mechanisms in bacteria and describe the methods used to study mutagenesis in bacterial systems.

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

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

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

Background

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

Methods

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

Results

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

Discussion

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