Genome-wide mutational diversity in an evolving population of Escherichia coli.

The level of genetic variation in a population is the result of a dynamic tension between evolutionary forces. Mutations create variation, certain frequency-dependent interactions may preserve diversity, and natural selection purges variation. New sequencing technologies offer unprecedented opportunities to discover and characterize the diversity present in evolving microbial populations on a whole-genome scale. By sequencing mixed-population samples, we have identified single-nucleotide polymorphisms (SNPs) present at various points in the history of an Escherichia coli population that has evolved for almost 20 years from a founding clone. With 50-fold genome coverage, we were able to catch beneficial mutations as they swept to fixation, discover contending beneficial alleles that were eliminated by clonal interference, and detect other minor variants possibly adapted to a new ecological niche. Additionally, there was a dramatic increase in genetic diversity late in the experiment after a mutator phenotype evolved. Still finer-resolution details of the structure of genetic variation and how it changes over time in microbial evolution experiments will enable new applications and quantitative tests of population genetic theory.

Historical and contingent factors affect re-evolution of a complex feature lost during mass extinction in communities of digital organisms.

Re-evolution of complex biological features following the extinction of taxa bearing them remains one of evolution's most interesting phenomena, but is not amenable to study in fossil taxa. We used communities of digital organisms (computer programs that self-replicate, mutate and evolve), subjected to periods of low resource availability, to study the evolution, loss and re-evolution of a complex computational trait, the function EQU (bit-wise logical equals). We focused our analysis on cases where the pre-extinction EQU clade had surviving descendents at the end of the extinction episode. To see if these clades retained the capacity to re-evolve EQU, we seeded one set of multiple subreplicate 'replay' populations using the most abundant survivor of the pre-extinction EQU clade, and another set with the actual end-extinction ancestor of the organism in which EQU re-evolved following the extinction episode. Our results demonstrate that stochastic, historical, genomic and ecological factors can lead to constraints on further adaptation, and facilitate or hinder re-evolution of a complex feature.

Epistasis between new mutations and genetic background and a test of genetic canalization.

The importance for fitness of epistatic interactions among mutations is poorly known, yet epistasis can exert important effects on the dynamics of evolving populations. We showed previously that epistatic interactions are common between pairs of random insertion mutations in the bacterium Escherichia coli. In this paper, we examine interactions between these mutations and other mutations by transducing each of twelve insertion mutations into two genetic backgrounds, one ancestral and the other having evolved in, and adapted to, a defined laboratory environment for 10,000 generations. To assess the effect of the mutation on fitness, we allowed each mutant to compete against its unmutated counterpart in that same environment. Overall, there was a strong positive correlation between the mutational effects on the two genetic backgrounds. Nonetheless, three of the twelve mutations had significantly different effects on the two backgrounds, indicating epistasis. There was no significant tendency for the mutations to be less harmful on the derived background. Thus, there is no evidence supporting the hypothesis that the derived bacteria had adapted, in part, by becoming buffered against the harmful effects of mutations.

Contribution of individual random mutations to genotype-by-environment interactions in Escherichia coli.

Numerous studies have shown genotype-by-environment (GxE) interactions for traits related to organismal fitness. However, the genetic architecture of the interaction is usually unknown because these studies used genotypes that differ from one another by many unknown mutations. These mutations were also present as standing variation in populations and hence had been subject to prior selection. Based on such studies, it is therefore impossible to say what fraction of new, random mutations contributes to GxE interactions. In this study, we measured the fitness in four environments of 26 genotypes of Escherichia coli, each containing a single random insertion mutation. Fitness was measured relative to their common progenitor, which had evolved on glucose at 37 degrees C for the preceding 10,000 generations. The four assay environments differed in limiting resource and temperature (glucose, 28 degrees C; maltose, 28 degrees C; glucose, 37 degrees C; and maltose, 37 degrees C). A highly significant interaction between mutation and resource was found. In contrast, there was no interaction involving temperature. The resource interaction reflected much higher among mutation variation for fitness in maltose than in glucose. At least 11 mutations (42%) contributed to this GxE interaction through their differential fitness effects across resources. Beneficial mutations are generally thought to be rare but, surprisingly, at least three mutations (12%) significantly improved fitness in maltose, a resource novel to the progenitor. More generally, our findings demonstrate that GxE interactions can be quite common, even for genotypes that differ by only one mutation and in environments differing by only a single factor.

Evolution of digital organisms at high mutation rates leads to survival of the flattest.

Darwinian evolution favours genotypes with high replication rates, a process called 'survival of the fittest'. However, knowing the replication rate of each individual genotype may not suffice to predict the eventual survivor, even in an asexual population. According to quasi-species theory, selection favours the cloud of genotypes, interconnected by mutation, whose average replication rate is highest. Here we confirm this prediction using digital organisms that self-replicate, mutate and evolve. Forty pairs of populations were derived from 40 different ancestors in identical selective environments, except that one of each pair experienced a 4-fold higher mutation rate. In 12 cases, the dominant genotype that evolved at the lower mutation rate achieved a replication rate >1.5-fold faster than its counterpart. We allowed each of these disparate pairs to compete across a range of mutation rates. In each case, as mutation rate was increased, the outcome of competition switched to favour the genotype with the lower replication rate. These genotypes, although they occupied lower fitness peaks, were located in flatter regions of the fitness surface and were therefore more robust with respect to mutations.

Evolution of thermal dependence of growth rate of Escherichia coli populations during 20,000 generations in a constant environment.

Twelve experimental populations of the bacterium Escherichia coli evolved for 20,000 generations in a defined medium at 37 degrees C. We measured their maximum growth rates across a broad range of temperatures and at several evolutionary time points to quantify the extent to which they became thermal specialists with diminished performance at other temperatures. We also sought to determine whether antagonistic pleiotropy (genetic trade-offs) or mutation accumulation (drift decay) was primarily responsible for any thermal specialization. Populations showed consistent improvement in growth rate at moderate temperatures (27-39 degrees C), but tended to have decreased growth rate at both low (20 degrees C) and high (41-42 degrees C) temperatures. Most loss occurred early in the experiment, when adaptation was most rapid. This dynamic is predicted by antagonistic pleiotropy but not by mutation accumulation. Several populations evolved high mutation rates due to defects in their DNA repair, but they did not subsequently undergo a greater decrease in growth rate at thermal extremes than populations that retained low mutation rates, contrary to the acceleration of decay predicted by mutation accumulation. Antagonistic pleiotropy therefore is more likely to be responsible for the evolution of thermal specialization observed in maximum growth rate.

Mechanisms causing rapid and parallel losses of ribose catabolism in evolving populations of Escherichia coli B.

Twelve populations of Escherichia coli B all lost D-ribose catabolic function during 2,000 generations of evolution in glucose minimal medium. We sought to identify the population genetic processes and molecular genetic events that caused these rapid and parallel losses. Seven independent Rbs(-) mutants were isolated, and their competitive fitnesses were measured relative to that of their Rbs(+) progenitor. These Rbs(-) mutants were all about 1 to 2% more fit than the progenitor. A fluctuation test revealed an unusually high rate, about 5 x 10(-5) per cell generation, of mutation from Rbs(+) to Rbs(-), which contributed to rapid fixation. At the molecular level, the loss of ribose catabolic function involved the deletion of part or all of the ribose operon (rbs genes). The physical extent of the deletion varied between mutants, but each deletion was associated with an IS150 element located immediately upstream of the rbs operon. The deletions apparently involved transposition into various locations within the rbs operon; recombination between the new IS150 copy and the one upstream of the rbs operon then led to the deletion of the intervening sequence. To confirm that the beneficial fitness effect was caused by deletion of the rbs operon (and not some undetected mutation elsewhere), we used P1 transduction to restore the functional rbs operon to two Rbs(-) mutants, and we constructed another Rbs(-) strain by gene replacement with a deletion not involving IS150. All three of these new constructs confirmed that Rbs(-) mutants have a competitive advantage relative to their Rbs(+) counterparts in glucose minimal medium. The rapid and parallel evolutionary losses of ribose catabolic function thus involved both (i) an unusually high mutation rate, such that Rbs(-) mutants appeared repeatedly in all populations, and (ii) a selective advantage in glucose minimal medium that drove these mutants to fixation.

Evolutionary adaptation to temperature. VIII. Effects of temperature on growth rate in natural isolates of Escherichia coli and Salmonella enterica from different thermal environments.

Are enteric bacteria specifically adapted to the thermal environment of their hosts? In particular, do the optimal temperatures and thermal niches of the bacterial flora reflect seasonal, geographic, or phylogenetic differences in their hosts' temperatures? We examined these questions by measuring the relationship between the temperature-dependent growth rates of enteric bacteria in a free-living ectothermic host. We sampled two species of enteric bacteria (Escherichia coli and Salmonella enterica) from three natural populations of slider turtles (Trachemys scripta elegans) seasonally over two years. Despite pronounced differences in turtle body temperatures at different seasons and in different locations, we found no evidence that the thermal growth profiles of these bacteria mirrored this variation. Optimal temperatures and maximal growth rates in rich medium were nearly the same for both bacterial species (35-36 degrees C, 2.5 doublings per hour). The thermal niche (defined as the range of temperatures over which 75% of maximal growth rate occurred) was slightly higher for E. coli (28.5-41.0 degrees C) than for S. enterica (27.7-39.8 degrees C), but the niche breadth was about the same for both. We also measured the thermal dependence of growth rate in these same bacterial species isolated from mammalian hosts. Both bacterial species had temperatures of maximal growth and thermal niches that were about 2 degrees C higher than those of their respective conspecifics sampled from turtles; niche breadths were not different. These data suggest that these bacterial species are thermal generalists that do not track fine-scale changes in their thermal environments. Even major differences in body temperatures, as great as those between ectothermic and endothermic hosts, may result in the evolution of rather modest changes in thermal properties.

Evolutionary adaptation to temperature. IX. Preadaptation to novel stressful environments of Escherichia coli adapted to high temperature.

Stressful environments may be considered as those that reduce fitness, sometimes due in part to the increased metabolic expenditure required to sustain life. Direct adaptation to a stressor is expected to increase fitness and reduce maintenance metabolism, with the latter leading to increased biomass production. In this study, we test the general hypothesis that such adaptation to one stressor can preadapt organisms to novel stressful environments. Six lines of Escherichia coli propagated for 2,000 generations at 41-42 degrees C (42 group), a stressful temperature, were compared to six control lines propagated for 2,000 generations at 37degrees C (37 group) and to the common ancestor of both groups. We assayed biovolume yield (a measure of growth efficiency) and competitive fitness in the 42 group's selective high temperature environment as well as five novel stressful environments-acid, alkali, ethanol, high osmolarity and peroxide. As previously reported, at high temperature the 42 group had both higher yield and fitness than the 37 group and ancestor. In the novel environments, the 42 group generally produced yields higher than the 37 group (and marginally higher than the ancestor), but we found no differences in competitive fitness among the 37 and 42 groups and the ancestor. We also found that the performance of lines within groups was not correlated across stressful environments for either yield or relative fitness. Because previous adaptation to one stressor did not improve our measure of Darwinian fitness in novel stressful environments, we conclude that the 42 group shows no useful pre-adaptation, or cross-tolerance, to these types of environments.

The population genetics of ecological specialization in evolving Escherichia coli populations.

When organisms adapt genetically to one environment, they may lose fitness in other environments. Two distinct population genetic processes can produce ecological specialization-mutation accumulation and antagonistic pleiotropy. In mutation accumulation, mutations become fixed by genetic drift in genes that are not maintained by selection; adaptation to one environment and loss of adaptation to another are caused by different mutations. Antagonistic pleiotropy arises from trade-offs, such that the same mutations that are beneficial in one environment are detrimental in another. In general, it is difficult to distinguish between these processes. We analysed the decay of unused catabolic functions in 12 lines of Escherichia coli propagated on glucose for 20,000 generations. During that time, several lines evolved high mutation rates. If mutation accumulation is important, their unused functions should decay more than the other lines, but no significant difference was observed. Moreover, most catabolic losses occurred early in the experiment when beneficial mutations were being rapidly fixed, a pattern predicted by antagonistic pleiotropy. Thus, antagonistic pleiotropy appears more important than mutation accumulation for the decay of unused catabolic functions in these populations.

Long-term experimental evolution in Escherichia coli. IX. Characterization of insertion sequence-mediated mutations and rearrangements.

As part of a long-term evolution experiment, two populations of Escherichia coli B adapted to a glucose minimal medium for 10,000 generations. In both populations, multiple IS-associated mutations arose that then went to fixation. We identify the affected genetic loci and characterize the molecular events that produced nine of these mutations. All nine were IS-mediated events, including simple insertions as well as recombination between homologous elements that generated inversions and deletions. Sequencing DNA adjacent to the insertions indicates that the affected genes are involved in central metabolism (knockouts of pykF and nadR), cell wall synthesis (adjacent to the promoter of pbpA-rodA), and ill-defined functions (knockouts of hokB-sokB and yfcU). These genes are candidates for manipulation and competition experiments to determine whether the mutations were beneficial or merely hitchhiked to fixation.

Developmental cheating in the social bacterium Myxococcus xanthus.

Cheating is a potential problem in any social system that depends on cooperation and in which actions that benefit a group are costly to individuals that perform them. Genetic mutants that fail to perform a group-beneficial function but that reap the benefits of belonging to the group should have a within-group selective advantage, provided that the mutants are not too common. Here we show that social cheating exists even among prokaryotes. The bacterium Myxococcus xanthus exhibits several social behaviours, including aggregation of cells into spore-producing fruiting bodies during starvation. We examined a number of M. xanthus genotypes that were defective for fruiting-body development, including several lines that evolved for 1,000 generations under asocial conditions and others carrying defined mutations in developmental pathways, to determine whether they behaved as cheaters when mixed with their developmentally proficient progenitor. Clones from several evolved lines and two defined mutants exhibited cheating during development, being overrepresented among resulting spores relative to their initial frequency in the mixture. The ease of finding anti-social behaviours suggests that cheaters may be common in natural populations of M. xanthus.

Pervasive compensatory adaptation in Escherichia coli.

To investigate compensatory adaptation (CA), we used genotypes of Escherichia coli which were identical except for one or two deleterious mutations. We compared CA for (i) deleterious mutations with large versus small effects, (ii) genotypes carrying one versus two mutations, and (iii) pairs of deleterious mutations which interact in a multiplicative versus synergistic fashion. In all, we studied 14 different genotypes, plus a control strain which was not mutated. Most genotypes showed CA during 200 generations of experimental evolution, where we define CA as a fitness increase which is disproportionately large relative to that in evolving control lines, coupled with retention of the original deleterious mutation(s). We observed greater CA for mutations of large effect than for those of small effect, which can be explained by the greater benefit to recovery in severely handicapped genotypes given the dynamics of selection. The rates of CA were similar for double and single mutants whose initial fitnesses were approximately equal. CA was faster for synergistic than for multiplicative pairs, presumably because the marginal gain which results from CA for one of the component mutations is greater in that case. The most surprising result in our view, is that compensation should be so readily achieved in an organism which is haploid and has little genetic redundancy This finding suggests a degree of versatility in the E. coil genome which demands further study from both genetic and physiological perspectives.

Genome complexity, robustness and genetic interactions in digital organisms.

Digital organisms are computer programs that self-replicate, mutate and adapt by natural selection. They offer an opportunity to test generalizations about living systems that may extend beyond the organic life that biologists usually study. Here we have generated two classes of digital organism: simple programs selected solely for rapid replication, and complex programs selected to perform mathematical operations that accelerate replication through a set of defined 'metabolic' rewards. To examine the differences in their genetic architecture, we introduced millions of single and multiple mutations into each organism and measured the effects on the organism's fitness. The complex organisms are more robust than the simple ones with respect to the average effects of single mutations. Interactions among mutations are common and usually yield higher fitness than predicted from the component mutations assuming multiplicative effects; such interactions are especially important in the complex organisms. Frequent interactions among mutations have also been seen in bacteria, fungi and fruitflies. Our findings support the view that interactions are a general feature of genetic systems.

Mutation, recombination, and incipient speciation of bacteria in the laboratory.

Mutations in the DNA mismatch repair system increase mutation and recombination. They may thereby promote the genetic divergence that underlies speciation, after which the reacquisition of a functional repair system may sustain that divergence by creating a barrier to recombination. We tested several lines of Escherichia coli, derived from a common ancestor and evolved for 20,000 generations, for their recombination ability. Some lines, but not others, had become mismatch repair-defective mutators during experimental evolution, providing different opportunities for DNA sequence divergence. We knocked out the repair system in lines that had retained this function, and we restored function to those lines that had become defective. We then estimated recombination rates in various crosses between these repair-deficient and -proficient strains. The effect of the mismatch repair system on recombination was greatest in those lines that had evolved nonfunctional repair, indicating they had undergone more sequence divergence and, consequently, were more sensitive to the recombination-inhibiting effect of a functional repair system. These results demonstrate the establishment of an incipient genetic barrier between formerly identical lines, and they support a model in which the mismatch repair system can influence speciation dynamics through its simultaneous effects on mutation and recombination.

Genomic evolution during a 10,000-generation experiment with bacteria.

Molecular methods are used widely to measure genetic diversity within populations and determine relationships among species. However, it is difficult to observe genomic evolution in action because these dynamics are too slow in most organisms. To overcome this limitation, we sampled genomes from populations of Escherichia coli evolving in the laboratory for 10,000 generations. We analyzed the genomes for restriction fragment length polymorphisms (RFLP) using seven insertion sequences (IS) as probes; most polymorphisms detected by this approach reflect rearrangements (including transpositions) rather than point mutations. The evolving genomes became increasingly different from their ancestor over time. Moreover, tremendous diversity accumulated within each population, such that almost every individual had a different genetic fingerprint after 10,000 generations. As has been often suggested, but not previously shown by experiment, the rates of phenotypic and genomic change were discordant, both across replicate populations and over time within a population. Certain pivotal mutations were shared by all descendants in a population, and these are candidates for beneficial mutations, which are rare and difficult to find. More generally, these data show that the genome is highly dynamic even over a time scale that is, from an evolutionary perspective, very brief.

Diminishing returns from mutation supply rate in asexual populations.

Mutator genotypes with increased mutation rates may be especially important in microbial evolution if genetic adaptation is generally limited by the supply of mutations. In experimental populations of the bacterium Escherichia coli, the rate of evolutionary adaptation was proportional to the mutation supply rate only in particular circumstances of small or initially well-adapted populations. These experiments also demonstrate a "speed limit" on adaptive evolution in asexual populations, one that is independent of the mutation supply rate.