Harnessing recombination to speed adaptive evolution in Escherichia coli

Evolutionary engineering typically involves asexual propagation of a strain to improve a desired phenotype. However, asexual populations suffer from extensive clonal interference, a phenomenon where distinct lineages of beneficial clones compete and are often lost from the population given sufficient time. Improved adaptive mutants can likely be generated by genetic exchange between lineages, thereby reducing clonal interference. We present a system that allows continuous in situ recombination by using an Esherichia coli F-based conjugation system lacking surface exclusion. Evolution experiments revealed that Hfr-mediated recombination significantly speeds adaptation in certain circumstances. These results show that our system is stable, effective, and suitable for use in evolutionary engineering applications.     

Evolutionary dynamics of Candida albicans during in vitro evolution

While mechanisms of resistance to major antifungal agents have been characterized in Candida albicans, little is known about the evolutionary trajectories during the emergence of drug resistance. Here, we examined the evolutionary dynamics of C. albicans that evolved in vitro in the presence or absence of fluconazole using the visualizing evolution in real-time (VERT) method, a novel experimental approach that facilitates the systematic isolation of adaptive mutants that arise in the population. We found an increase in the frequency of adaptive events in the presence of fluconazole compared to the no-drug controls. Analysis of the evolutionary dynamics revealed that mutations that led to increased drug resistance appeared frequently and that mutants with increased levels of resistance arose in independent lineages. Interestingly, most adaptive mutants with increased fitness in the presence of the drug did not exhibit a significant fitness decrease in the absence of the drug, supporting the idea that rapid resistance can arise from mutations in strains maintained in the population prior to exposure to the drug.     

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.     

Fitness evolution and the rise of mutator alleles in experimental Escherichia coli populations

We studied the evolution of high mutation rates and the evolution of fitness in three experimental populations of Escherichia coli adapting to a glucose-limited environment. We identified the mutations responsible for the high mutation rates and show that their rate of substitution in all three populations was too rapid to be accounted for simply by genetic drift. In two of the populations, large gains in fitness relative to the ancestor occurred as the mutator alleles rose to fixation, strongly supporting the conclusion that mutator alleles fixed by hitchhiking with beneficial mutations at other loci. In one population, no significant gain in fitness relative to the ancestor occurred in the population as a whole while the mutator allele rose to fixation, but a substantial and significant gain in fitness occurred in the mutator subpopulation as the mutator neared fixation. The spread of the mutator allele from rarity to fixation took >1000 generations in each population. We show that simultaneous adaptive gains in both the mutator and wild-type subpopulations (clonal interference) retarded the mutator fixation in at least one of the populations. We found little evidence that the evolution of high mutation rates accelerated adaptation in these populations.     

Clonal interference and the periodic selection of new beneficial mutations in Escherichia coli

The conventional model of adaptation in asexual populations implies sequential fixation of new beneficial mutations via rare selective sweeps that purge all variation and preserve the clonal genotype. However, in large populations multiple beneficial mutations may co-occur, causing competition among them, a phenomenon called "clonal interference." Clonal interference is thus expected to lead to longer fixation times and larger fitness effects of mutations that ultimately become fixed, as well as to a genetically more diverse population. Here, we study the significance of clonal interference in populations consisting of mixtures of differently marked wild-type and mutator strains of Escherichia coli that adapt to a minimal-glucose environment for 400 generations. We monitored marker frequencies during evolution and measured the competitive fitness of random clones from each marker state after evolution. The results demonstrate the presence of multiple beneficial mutations in these populations and slower and more erratic invasion of mutants than expected by the conventional model, showing the signature of clonal interference. We found that a consequence of clonal interference is that fitness estimates derived from invasion trajectories were less than half the magnitude of direct estimates from competition experiments, thus revealing fundamental problems with this fitness measure. These results force a reevaluation of the conventional model of periodic selection for asexual microbes.     

大肠杆菌乙酸耐受性菌株的构建及其耐受机制研究进展

乙酸是微生物发酵生产常见的副产物,也可作为碳源存在于木质纤维素水解液等非粮原料发酵培养基中。培养基中含有高浓度的乙酸/乙酸盐时会抑制细胞生长、降低生物量,影响目标产品的产量和产率。研究乙酸耐受性机制,改进菌株的乙酸耐受性,构建具有高乙酸耐受性工程菌株,对于以乙酸为碳源或利用含乙酸的原料进行高附加值产品发酵生产具有重要意义。本文综述了通过代谢工程、实验室适应性进化、全局转录机器工程和基于CRISPR可追踪基因组工程等方法构建大肠杆菌乙酸耐受性菌株的研究进展,进一步从乙酸同化代谢、氨基酸依赖型代谢、离子转运系统调节和细胞膜成分修饰等4个方面阐述了大肠杆菌乙酸耐受性菌株的耐受性应答机制,总结了大肠杆菌乙酸耐受菌株的生产应用,展望了提高大肠杆菌乙酸耐受方法和大肠杆菌乙酸耐受机制的研究方向。     

Cellular heterogeneity in DNA alkylation repair increases population genetic plasticity

DNA repair mechanisms fulfil a dual role, as they are essential for cell survival and genome maintenance. Here, we studied how cells regulate the interplay between DNA repair and mutation. We focused on the adaptive response that increases the resistance of Escherichia coli cells to DNA alkylation damage. Combination of single-molecule imaging and microfluidic-based single-cell microscopy showed that noise in the gene activation timing of the master regulator Ada is accurately propagated to generate a distinct subpopulation of cells in which all proteins of the adaptive response are essentially absent. Whereas genetic deletion of these proteins causes extreme sensitivity to alkylation stress, a temporary lack of expression is tolerated and increases genetic plasticity of the whole population. We demonstrated this by monitoring the dynamics of nascent DNA mismatches during alkylation stress as well as the frequency of fixed mutations that are generated by the distinct subpopulations of the adaptive response. We propose that stochastic modulation of DNA repair capacity by the adaptive response creates a viable hypermutable subpopulation of cells that acts as a source of genetic diversity in a clonal population.     

Continuing coevolution of virus and defective interfering particles and of viral genome sequences during undiluted passages: virus mutants exhibiting nearly complete resistance to formerly dominant defective interfering particles

We quantitatively analyzed the interference interactions between defective interfering (DI) particles and mutants of cloned vesicular stomatitis virus passaged undiluted hundreds of times in BHK-21 cells. DI particles which predominated at different times in these serial passages always interfered most strongly (and very efficiently) with virus isolated a number of passages before the isolation of the DI particles. Virus isolated at the same passage level as the predominant DI particles usually exhibited severalfold resistance to these DI particles. Virus mutants (Sdi- mutants) isolated during subsequent passages always showed increasing resistance to these DI particles, followed by decreasing resistance as new DI particles arose to predominate and exert their own selective pressures on the virus mutant population. It appears that such coevolution of virus and DI particle populations proceeds indefinitely through multiple cycles of selection of virus mutants resistant to a certain DI particle (or DI particle class), followed by mutants resistant to a newly predominant DI particle, etc. At the peak of resistance, virus mutants were isolated which were essentially completely resistant to a particular DI particle; i.e., they were several hundred thousand-fold resistant, and they formed plaques of normal size and numbers in the presence of extremely high multiplicities of the DI particle. However, they were sensitive to interference by other DI particles. Recurring population interactions of this kind can promote rapid virus evolution. Complete sequencing of the N (nucleocapsid) and NS (polymerase associated) genes of numerous Sdi- mutants collected at passage intervals showed very few changes in the NS protein, but the N gene gradually accumulated a series of stable nucleotide and amino acid substitutions, some of which correlated with extensive changes in the Sdi- phenotype. Likewise, the 5' termini (and their complementary plus-strand 3' termini) continued to accumulate extensive base substitutions which were strikingly confined to the first 47 nucleotides. We also observed addition and deletion mutations in noncoding regions of the viral genome at a level suggesting that they probably occur at a high frequency throughout the genome, but usually with lethal or debilitating consequences when they occur in coding regions.     

Optimal Strategy for Competence Differentiation in Bacteria

A phylogenetically diverse subset of bacterial species are naturally competent for transformation by DNA. Transformation entails recombination of genes between different lineages, representing a form of bacterial sex that increases standing genetic variation. We first assess whether homologous recombination by transformation is favored by evolution. Using stochastic population genetic computer simulations in which beneficial and deleterious mutations occur at many loci throughout the whole genome, we find that transformation can increase both the rate of adaptive evolution and the equilibrium level of fitness. Secondly, motivated by experimental observations of Bacillus subtilis, we assume that competence additionally entails a weak persister phenotype, i.e., the rates of birth and death are reduced for these cells. Consequently, persisters evolve more slowly than non-persisters. We show via simulation that strains which stochastically switch into and out of the competent phenotype are evolutionarily favored over strains that express only a single phenotype. Our model''s simplicity enables us to derive and numerically solve a system of finite- deterministic equations that describe the evolutionary dynamics. The observed tradeoff between the benefit of recombination and the cost of persistence may explain the previously mysterious observation that only a fractional subpopulation of B. subtilis cells express competence. More generally, this work demonstrates that population genetic forces can give rise to phenotypic diversity even in an unchanging and homogeneous environment.     

Population dynamics and the evolution of antifungal drug resistance in Candida albicans

Candida albicans is an important human fungal pathogen. Resistance to all major antifungal agents has been observed in clinical isolates of Candida spp. and is a major clinical challenge. The rise and expansion of drug-resistant mutants during exposure to antifungal agents occurs through a process of adaptive evolution, with potentially complex population dynamics. Understanding the population dynamics during the emergence of drug resistance is important for determining the fundamental principles of how fungal pathogens evolve for resistance. While few detailed reports that focus on the population dynamics of C. albicans currently exist, several important features on the population structure and adaptive landscape can be elucidated from existing evolutionary studies in in vivo and in vitro systems.     

Stress-driven in vivo selection of a functional mini-gene from a randomized DNA library expressing combinatorial peptides in Escherichia coli

A plasmid-borne randomized mini-gene library expressing a population of combinatorial 20-mer peptides with no bias toward any biological function was used as an initial source of genetic variance during stress-driven evolution of Escherichia coli. The transformed bacteria were evolved under multiple rounds of selective pressure imposed by nearly lethal concentrations of NiCl(2), AgNO(3), or K(2)TeO(3). At the final stage, clones conferring resistance to NiCl(2) were found to carry identical functional mini-genes, which conferred significant nickel tolerance on the host cells. Expression of the mini-gene markedly improved growth parameters of the evolved clones at subinhibitory concentrations of NiCl(2) while being slightly detrimental in the absence of the stress. This substantial increase in resistance, as compared with control cultures adapted in the absence of the mini-gene, is shown to be largely due to coadaptation with changes elsewhere in the E. coli genome. Clones resistant to AgNO(3) and K(2)TeO(3) harbored plasmid variants with an inactive mini-gene and with a deleted mini-gene operon, respectively. In those cases, an exploration of the mini-gene sequence space apparently was fruitless, and the developed toxicity tolerance was likely to be exclusively associated with acquired adaptive mutations. Overall, the results demonstrate a very natural outcome in which the mini-genes were expected to be either successfully integrated into the bacterial genetic network or rejected depending on their effect on host fitness. This approach is immediately useful as a laboratory model to study the dynamics of bacterial adaptive evolution at the molecular level and is especially relevant as a rapid method to study cellular response to recently acquired genetic material.     

Contrasting dynamics of a mutator allele in asexual populations of differing size

Mutators have been shown to hitchhike in asexual populations when the anticipated beneficial mutation supply rate of the mutator subpopulation, NU(b) (for subpopulation of size N and beneficial mutation rate U(b)) exceeds that of the wild-type subpopulation. Here, we examine the effect of total population size on mutator dynamics in asexual experimental populations of Saccharomyces cerevisiae. Although mutators quickly hitchhike to fixation in smaller populations, mutator fixation requires more and more time as population size increases; this observed delay in mutator hitchhiking is consistent with the expected effect of clonal interference. Interestingly, despite their higher beneficial mutation supply rate, mutators are supplanted by the wild type in very large populations. We postulate that this striking reversal in mutator dynamics is caused by an interaction between clonal interference, the fitness cost of the mutator allele, and infrequent large-effect beneficial mutations in our experimental populations. Our work thus identifies a potential set of circumstances under which mutator hitchhiking can be inhibited in natural asexual populations, despite recent theoretical predictions that such populations should have a net tendency to evolve ever-higher genomic mutation rates.     

The dynamics of HIV-1 adaptation in early infection

Human immunodeficiency virus type 1 (HIV-1) undergoes a severe population bottleneck during sexual transmission and yet adapts extremely rapidly to the earliest immune responses. The bottleneck has been inferred to typically consist of a single genome, and typically eight amino acid mutations in viral proteins spread to fixation by the end of the early chronic phase of infection in response to selection by CD8(+) T cells. Stochastic simulation was used to examine the effects of the transmission bottleneck and of potential interference among spreading immune-escape mutations on the adaptive dynamics of the virus in early infection. If major viral population genetic parameters are assigned realistic values that permit rapid adaptive evolution, then a bottleneck of a single genome is not inconsistent with the observed pattern of adaptive fixations. One requirement is strong selection by CD8(+) T cells that decreases over time. Such selection may reduce effective population sizes at linked loci through genetic hitchhiking. However, this effect is predicted to be minor in early infection because the transmission bottleneck reduces the effective population size to such an extent that the resulting strong selection and weak mutation cause beneficial mutations to fix sequentially and thus avoid interference.     

THE EFFECTS OF POPULATION BOTTLENECKS ON CLONAL INTERFERENCE, AND THE ADAPTATION EFFECTIVE POPULATION SIZE

Clonal interference refers to the competition that arises in asexual populations when multiple beneficial mutations segregate simultaneously. A large body of theoretical and experimental work now addresses this issue. Although much of the experimental work is performed in populations that grow exponentially between periodic population bottlenecks, the theoretical work to date has addressed only populations of a constant size. We derive an analytical approximation for the rate of adaptation in the presence of both clonal interference and bottlenecks, and compare this prediction to the results of an individual-based simulation, showing excellent agreement in the parameter regime in which clonal interference prevails. We also derive an appropriate definition for the effective population size for adaptive evolution experiments in the presence of population bottlenecks. This "adaptation effective population size" allows for a good approximation of the expected rate of adaptation, either in the strong-selection weak-mutation regime, or when clonal interference comes into play. In the multiple mutation regime, when the product of the population size and mutation rate is extremely large, these results no longer hold.     

Deterministic and stochastic regimes of asexual evolution on rugged fitness landscapes

We study the adaptation dynamics of an initially maladapted asexual population with genotypes represented by binary sequences of length L. The population evolves in a maximally rugged fitness landscape with a large number of local optima. We find that whether the evolutionary trajectory is deterministic or stochastic depends on the effective mutational distance d(eff) up to which the population can spread in genotype space. For d(eff) = L, the deterministic quasi-species theory operates while for d(eff) < 1, the evolution is completely stochastic. Between these two limiting cases, the dynamics are described by a local quasi-species theory below a crossover time T(x) while above T(x) the population gets trapped at a local fitness peak and manages to find a better peak via either stochastic tunneling or double mutations. In the stochastic regime d(eff) < 1, we identify two subregimes associated with clonal interference and uphill adaptive walks, respectively. We argue that our findings are relevant to the interpretation of evolution experiments with microbial populations.     

SLOWLY SWITCHING BETWEEN ENVIRONMENTS FACILITATES REVERSE EVOLUTION IN SMALL POPULATIONS

Natural populations must constantly adapt to ever‐changing environmental conditions. A particularly interesting question is whether such adaptations can be reversed by returning the population to an ancestral environment. Such evolutionary reversals have been observed in both natural and laboratory populations. However, the factors that determine the reversibility of evolution are still under debate. The time scales of environmental change vary over a wide range, but little is known about how the rate of environmental change influences the reversibility of evolution. Here, we demonstrate computationally that slowly switching between environments increases the reversibility of evolution for small populations that are subject to only modest clonal interference. For small populations, slow switching reduces the mean number of mutations acquired in a new environment and also increases the probability of reverse evolution at each of these “genetic distances.” As the population size increases, slow switching no longer reduces the genetic distance, thus decreasing the evolutionary reversibility. We confirm this effect using both a phenomenological model of clonal interference and also a Wright–Fisher stochastic simulation that incorporates genetic diversity. Our results suggest that the rate of environmental change is a key determinant of the reversibility of evolution, and provides testable hypotheses for experimental evolution.