Adaptive Evolution That Requires Multiple Spontaneous Mutations. I. Mutations Involving an Insertion Sequence

Escherichia coli K12 strain chi 342LD requires two mutations in the bgl (beta-glucosidase) operon, bglR0----bglR+ and excision of IS103 from within bglF, in order to utilize salicin. In growing cells the two mutations occur at rates of 4 x 10(-8) per cell division and less than 2 x 10(-12) per cell division, respectively. In 2-3-week-old colonies on MacConkey salicin plates the double mutants occur at frequencies of 10(-8) per cell, yet the rate of an unselected mutation, resistance to valine, is unaffected. The two mutations occur sequentially. Colonies that are 8-12 days old contain from 1% to about 10% IS103 excision mutants, from which the Sal+ secondary bglR0----bglR+ mutants arise. It is shown that the excision mutants are not advantageous within colonies; thus, they must result from a burst of independent excisions late in the life of the colony. Excision of IS103 occurs only on medium containing salicin, despite the fact that the excision itself confers no detectable selective advantage and serves only to create the potential for a secondary selectively advantageous mutation.     

Activation of the bgl operon by adaptive mutation

In growing Escherichia coli K12 cells, the cryptic bgl operon is activated 98% of the time by insertions of IS1 or IS5 into the control region, designated bglR. The activated bgl operon permits utilization of the beta-glucoside sugar arbutin as a sole carbon and energy source. The bgl operon is also activated by late-occurring mutations during prolonged selection on arbutin. The late-occurring mutations that occurred during prolonged carbon starvation in the presence of arbutin were "adaptive mutations" because they were specific to the presence of arbutin, and they did not occur during prolonged starvation in the absence of arbutin. The spectrum of late-arising mutations differed from that of early-arising, growth-dependent mutations in that 20% of the late-arising mutants resulted from mutations at the hns locus. This provides the first direct evidence for adaptive mutagenesis mediated by the insertion of IS elements. Because no special genetic background is required to select Bgl+ mutants, this affords the opportunity to study IS-element-mediated adaptive mutagenesis in a variety of genetic backgrounds, including the backgrounds of natural isolates of E. coli.      

The structure of linkage disequilibrium around a selective sweep

The fixation of advantageous mutations by natural selection has a profound impact on patterns of linked neutral variation. While it has long been appreciated that such selective sweeps influence the frequency spectrum of nearby polymorphism, it has only recently become clear that they also have dramatic effects on local linkage disequilibrium. By extending previous results on the relationship between genealogical structure and linkage disequilibrium, I obtain simple expressions for the influence of a selective sweep on patterns of allelic association. I show that sweeps can increase, decrease, or even eliminate linkage disequilibrium (LD) entirely depending on the relative position of the selected and neutral loci. I also show the importance of the age of the neutral mutations in predicting their degree of association and describe the consequences of such results for the interpretation of empirical data. In particular, I demonstrate that while selective sweeps can eliminate LD, they generate patterns of genetic variation very different from those expected from recombination hotspots.     

Characterization of cycA mutants of Escherichia coli. An assay for measuring in vivo mutation rates

Quantitative assessment of the spontaneous or induced genomic mutation rate, a fundamental evolutionary parameter, usually requires the use of well-characterized mutant selection systems. Although there is a great number of genetic selection schemes available in Escherichia coli, the selection of D-cycloserine resistant mutants is shown here to be particularly useful to yield a general view of mutation rates and spectra. The combination of a well-defined experimental protocol with the Ma-Sandri-Sarkar maximum likelihood method of fluctuation analysis results in reproducible data, adequate for statistical comparisons. The straightforward procedure is based on a simple phenotype-genotype relationship, and detects mutations in the single-copy, chromosomal cycA gene, involved in the uptake of D-cycloserine. In contrast to the widely used rifampicin resistance assay, the procedure selects mutations which are neutral in respect of cell growth. No specific genetic background is needed, and practically the entire mutation spectrum (base substitutions, frameshifts, deletions, insertions) can simultaneously be measured. A systematic analysis of cycA mutations revealed a spontaneous mutation rate of 6.54 x 10(-8) in E. coli K-12 MG1655. The mutation spectrum was dominated by point mutations (base substitutions, frameshifts), spread over the entire gene. IS insertions, caused by IS1, IS2, IS3, IS4, IS5 and IS150, represented 24% of the mutations.     

The signature of positive selection at randomly chosen loci

In Drosophila and humans, there are accumulating examples of loci with a significant excess of high-frequency-derived alleles or high levels of linkage disequilibrium, relative to a neutral model of a random-mating population of constant size. These are features expected after a recent selective sweep. Their prevalence suggests that positive directional selection may be widespread in both species. However, as I show here, these features do not persist long after the sweep ends: The high-frequency alleles drift to fixation and no longer contribute to polymorphism, while linkage disequilibrium is broken down by recombination. As a result, loci chosen without independent evidence of recent selection are not expected to exhibit either of these features, even if they have been affected by numerous sweeps in their genealogical history. How then can we explain the patterns in the data? One possibility is population structure, with unequal sampling from different subpopulations. Alternatively, positive selection may not operate as is commonly modeled. In particular, the rate of fixation of advantageous mutations may have increased in the recent past.     

The role of single-mutant intermediates in the generation of trpAB double revertants during prolonged selection.

Selection-induced mutations are nonrandom mutations that occur as specific, direct responses to environmental challenges and that occur more often when they are selectively advantageous than when they are selectively neutral. One of the most puzzling examples of selection-induced mutations involved the simultaneous reversions of two mutations, one in trpA and the other in trpB, at rates that were several orders of magnitude greater than would have been predicted if the two mutations had occurred as independent events (B. G. Hall, Proc. Natl. Acad. Sci. USA 88:5882-5886, 1991). Here I examine the possibility that the double mutations might be accounted for by sequential mutations with intervening growth.     

Conditions for mutation-order speciation

Two models for speciation via selection have been proposed. In the well-known model of ‘ecological speciation’, divergent natural selection between environments drives the evolution of reproductive isolation. In a second ‘mutation-order’ model, different, incompatible mutations (alleles) fix in different populations adapting to the same selective pressure. How to demonstrate mutation-order speciation has been unclear, although it has been argued that it can be ruled out when gene flow occurs because the same, most advantageous allele will fix in all populations. However, quantitative examination of the interaction of factors influencing the likelihood of mutation-order speciation is lacking. We used simulation models to study how gene flow, hybrid incompatibility, selective advantage, timing of origination of new mutations and an initial period of allopatric differentiation affect population divergence via the mutation-order process. We find that at least some population divergence can occur under a reasonably wide range of conditions, even with moderate gene flow. However, strong divergence (e.g. fixation of different alleles in different populations) requires very low gene flow, and is promoted when (i) incompatible mutations have similar fitness advantages, (ii) less fit mutations arise slightly earlier in evolutionary time than more fit alternatives, and (iii) allopatric divergence occurs prior to secondary contact.     

On the Specificity of Adaptive Mutations

Adaptive mutations are mutations that occur in nondividing or slowly dividing cells during prolonged nonlethal selection, and that appear to be specific to the challenge of the selection in the sense that the only mutations that arise are those that provide a growth advantage to the cell. The issue of the specificity has been controversial because it violates our most basic assumptions about the randomness of mutations with respect to their effect on the cell. Although a variety of experiments in several systems in both bacteria and yeast have claimed to demonstrate that specificity, those experiments have been subjected to a variety of technical criticisms suggesting that the specificity may not be real. Here I use the ebg system to provide evidence that when selection is applied to one specific nucleotide site within a gene, mutation occurs at that site but not at an alternative and equally mutable site within the same gene.     

'Haldane's Sieve' in a metapopulation: sifting through plant reproductive polymorphisms

An important result of population genetics is that advantageous mutations will be fixed by selection in a population with a greater probability if they are dominant rather than recessive. This selective filter on new variants entering a population, termed 'Haldane's Sieve', has hitherto been invoked to account for the greater role of dominant than completely recessive mutations in adaptive evolution. Here, we suggest that a process similar to Haldane's Sieve will act on migrants into subpopulations of a metapopulation, and that the repeated action of Haldane's Sieve on alleles maintained by frequency-dependent selection, such as those responsible for many plant reproductive polymorphisms, is expected to bias their frequency distribution in favour of dominant alleles. The genetic and phenotypic signatures left by these processes might provide additional indirect support for the contentious idea that metapopulation dynamics have had an important role in shaping the ecology and evolution of some plant species.     

Selection, adaptation, and bacterial operons

Bacteria are especially useful as systems to study the molecular basis of adaptive evolution. Selection for novel metabolic capabilities has allowed us to study the evolutionary potential of organisms and has shown that there are three major "strategies" for the evolution of new metabolic functions. (i) Regulatory mutations may allow a gene to be expressed under unusual conditions. If the product of that gene is already active toward a novel resource, then a regulatory mutation alone may confer a new metabolic capability. (ii) Structural gene mutations may alter the catalytic properties of enzymes so that they can act on novel substrates. These structural gene mutations may dramatically improve catalytic capabilities, and in some cases they can confer entirely new capabilities upon enzymes. In most cases both regulatory and structural gene mutations are required for the effective evolution of new metabolic functions. (iii) Operons that are normally silent, or cryptic, may be activated by either point mutations or by the action of mobile genetic elements. When activated, these operons can provide entirely new pathways for the metabolism of novel resources. Selection can also play a role in modulating the probability that a particular adaptive mutation will occur. In this paper I present evidence that a specific adaptive mutation, reversion of the metB1 mutation, occurs 60 to 80 times more frequently during prolonged selection on plates under conditions where the members of the population are not growing than it does in growing cells under nonselective conditions. This selective condition, methionine starvation, does not increase the frequency of other mutations unrelated to methionine biosynthesis.(ABSTRACT TRUNCATED AT 250 WORDS)     

Physical Analysis of Tn10- and Is10-Promoted Transpositions and Rearrangements

We have investigated by Southern blot hybridization the rate of IS10 transposition and other Tn10/IS10-promoted rearrangements in Escherichia coli and Salmonella strains bearing single chromosomal insertions of Tn10 or a related Tn10 derivative. We present evidence for three primary conclusions. First, the rate of IS10 transposition is approximately 10(-4) per cell per bacterial generation when overnight cultures are grown and plated on minimal media and is at least ten times more frequent than any other Tn10/IS10-promoted DNA alteration. Second, all of the chromosomal rearrangements observed can be accounted for by two previously characterized Tn10-promoted rearrangements: deletion/inversions and deletions. Together these rearrangements occur at about 10% the rate of IS10 transposition. Third, the data suggest that intramolecular Tn10-promoted rearrangements preferentially use nearby target sites, while the target sites for IS10 transposition events are scattered randomly around the chromosome.     

Transposon-mediated adaptive and directed mutations and their potential evolutionary benefits

Transposons, mobile genetic elements that can hop from one chromosomal location to another, are known to be both beneficial and deleterious to the cell that bears them. Their value in accelerating evolutionary adaptation is well recognized. We herein summarize published research dealing with these elements and then move on to review our own research efforts which focus on a small transposon that can induce mutations under the control of host factors in a process that phenotypically and mechanistically conforms to the definition of 'directed mutation'. Directed mutations occur at higher frequencies when they are beneficial, being induced by the stress condition that they relieve. Here, we review evidence for transposon-mediated directed mutation in Escherichia coli. Deletion mutants in the crp gene can not grow on glycerol (Glp(-)); however, these cells mutate specifically to efficient glycerol utilization (Glp(+)) at rates that are greatly enhanced by the presence of glycerol or the loss of the glycerol repressor (GlpR). These rates are greatly depressed by glucose or by glpR overexpression. Of the four tandem GlpR-binding sites (O1-O4) in the control region of the glpFK operon, O4 (downstream) specifically controls glpFK expression while O1 (upstream) controls mutation rate. Mutation is due to insertion of the small transposon IS5 into a specific site just upstream of the glpFK promoter. Mutational control by the glycerol regulon repressor GlpR is independent of the selection and assay procedures, and IS5 insertion into other gene activation sites is unaffected by the presence of glycerol or the loss of GlpR. The results establish the principle of transposon-mediated directed mutation, identify a protein responsible for its regulation, and define essential aspects of the mechanism.     

Mutational order: a major stochastic process in evolution

Computer simulations in which selection acts on a quantitative character show that the randomness of mutations can contribute significantly to evolutionary divergence between populations. In different populations, different advantageous mutations occur, and are selected to fixation, so that the populations diverge even when they are initially identical, and are subject to identical selection. This stochastic process is distinct from random genetic drift. In some circumstances (large populations or strong selection, or both) mutational order can be greatly more important than random drift in bringing about divergence. It can generate a 'disconnection' between evolution at the phenotypic and genotypic levels, and can give rise to a rough 'molecular clock', albeit episodic, that is driven by selection. In the absence of selection, mutational order has little or no effect.     

Reaction norms with bifurcations shaped by evolution

Two versions of a model for the evolution of seasonal polyphenism investigate the evolution of reaction norm bifurcation and branching. The first version is without a specific submodel for morphological development and the second has an explicit developmental map. Version 1 is evolutionarily relatively unconstrained: (i) reaction norms are specified by matrices containing the probabilities of occurrence of environment-phenotype combinations, (ii) all conceivable reaction norm matrices are reachable through a sequence of mutations, and (iii) small as well as large mutational effects occur. This version is used to find the evolutionarily stable strategy favoured by the population ecology that is characterized by stabilizing viability selection with a cyclically fluctuating selection optimum. When the strength of selection is large and when the lag between initiation of development and selection on mature phenotype is not a multiple of half the period of the environmental cycle, a branching reaction norm evolves. In the second model version, branching reaction norms occur for certain parameter combinations of the developmental submodel, but the evolution of this pattern is often constrained. The evolutionary trajectory becomes trapped in a local selective optimum for the parameters of the developmental system. Substantial developmental noise evolves, but mutations that produce a selectively advantageous branching pattern do not occur from there.     

A simple method for estimating the intensity of purifying selection in protein-coding genes.

We propose a method by which the intensity of purifying selection on a functional protein-coding gene is estimated by using three aligned homologous sequences: a processed pseudogene (psi), a functional paralog from the same species (g), and a functional ortholog from a different species (o). For each such trio, we calculate the numbers of nucleotide substitutions along the branches leading to psi and g, i.e., K psi and K(g). If we assume that the mutation rates are the same in the genes and the pseudogenes and that mutations occurring in a pseudogene do not affect the fitness of the organism, we can show that the fraction of mutations that are selectively neutral, fg, is equal to the ratio K(g)/K psi. Since advantageous mutations occur only very rarely, such that they do not contribute significantly to the rate of molecular evolution, the fraction of deleterious mutations that are subject to purifying selection is 1-fg. Therefore, the K(g)/K psi ratio can be used directly to estimate the intensity of purifying selection, thereby isolating its effects on the rate of evolution from those of mutation. We compared the selection intensities of 12 orthologous protein-coding pairs from humans and murids. As expected, the fraction of mutations that are subject to purifying selection is strongest in the second codon position and weakest in the third. Interestingly, the mean fractions of effectively neutral mutations in the third codon position were only 41% and 42% for murids and humans, respectively, indicating that many synonymous mutations are subject to selective constraint. In several orthologous genes, we found that the intensity of purifying selection is very different between murid and human orthologous genes. There was no statistically significant difference in overall intensity of purifying selection between humans and murids. Thus, purifying selection does not seem to be an important factor contributing to the observed differences in the rates of evolution between these two taxa.     

Quasi-independence, fitness, and advantageousness

I argue that the idea of ‘quasi-independence’ [Lewontin, R. C. (1978). Adaptation. Scientific American, 239(3), 212–230] cannot be understood without attending to the distinction between fitness and advantageousness [Sober, E. (1993). Philosophy of biology. Boulder: Westview Press]. Natural selection increases the frequency of fitter traits, not necessarily of advantageous ones. A positive correlation between an advantageous trait and a disadvantageous one may prevent the advantageous trait from evolving. The quasi-independence criterion is aimed at specifying the conditions under which advantageous traits will evolve by natural selection in this type of situation. Contrary to what others have argued [Sterelny, K. (1992). Evolutionary explanations of human behavior. Australian Journal of Philosophy, 70(2), 156–172, and Sterelny, K., &; Griffiths, P. (1999). Sex and death. Chicago: University of Chicago Press], these conditions must involve a precise quantitative measure of (a) the extent to which advantageous traits are beneficial, and (b) the degree to which they are correlated with other traits. Driscoll (2004) [Driscoll, C. (2004). Can behaviors be adaptations? Philosophy of Science, 71, 16–35] recognizes the need for such a measure, but I argue that she does not provide the correct formulation. The account of quasi-independence that I offer clarifies this point.     

Somatic selection for and against cancer

In multicellular organisms, cells cooperate within a well-defined developmental program. Cancer is a breakdown of such cooperation: cells mutate to phenotypes of uncoordinated proliferation. We study basic principles of the architecture of solid tissues that influence the rate of cancer initiation. In particular, we explore how somatic selection acts to prevent or to promote cancer. Cells with mutations in oncogenes or tumor suppressor genes often have increased proliferation rates. Somatic selection increases their abundance and thus enhances the risk of cancer. Many potentially harmful mutations, however, increase the probability of triggering apoptosis and, hence, initially lead to cells with reduced net proliferation rates. Such cells are eliminated by somatic selection, which therefore also works to reduce the risk of cancer. We show that a tissue organization into small compartments avoids the rapid spread of mutations in oncogenes and tumor suppressor genes, but promotes genetic instability. In small compartments, genetic instability, which confers a selective disadvantage for the cell, can spread by random drift. If both deleterious and advantageous mutations participate in tumor initiation, then we find an intermediate optimum for the compartment size.     

Population,evolutionary and genomic consequences of interference selection

Weakly selected mutations are most likely to be physically clustered across genomes and, when sufficiently linked, they alter each others' fixation probability, a process we call interference selection (IS). Here we study population genetics and evolutionary consequences of IS on the selected mutations themselves and on adjacent selectively neutral variation. We show that IS reduces levels of polymorphism and increases low-frequency variants and linkage disequilibrium, in both selected and adjacent neutral mutations. IS can account for several well-documented patterns of variation and composition in genomic regions with low rates of crossing over in Drosophila. IS cannot be described simply as a reduction in the efficacy of selection and effective population size in standard models of selection and drift. Rather, IS can be better understood with models that incorporate a constant "traffic" of competing alleles. Our simulations also allow us to make genome-wide predictions that are specific to IS. We show that IS will be more severe at sites in the center of a region containing weakly selected mutations than at sites located close to the edge of the region. Drosophila melanogaster genomic data strongly support this prediction, with genes without introns showing significantly reduced codon bias in the center of coding regions. As expected, if introns relieve IS, genes with centrally located introns do not show reduced codon bias in the center of the coding region. We also show that reasonably small differences in the length of intermediate "neutral" sequences embedded in a region under selection increase the effectiveness of selection on the adjacent selected sequences. Hence, the presence and length of sequences such as introns or intergenic regions can be a trait subject to selection in recombining genomes. In support of this prediction, intron presence is positively correlated with a gene's codon bias in D. melanogaster. Finally, the study of temporal dynamics of IS after a change of recombination rate shows that nonequilibrium codon usage may be the norm rather than the exception.     

The reality and importance of founder speciation in evolution

A founder event occurs when a new population is established from a small number of individuals drawn from a large ancestral population. Mayr proposed that genetic drift in an isolated founder population could alter the selective forces in an epistatic system, an observation supported by recent studies. Carson argued that a period of relaxed selection could occur when a founder population is in an open ecological niche, allowing rapid population growth after the founder event. Selectable genetic variation can actually increase during this founder-flush phase due to recombination, enhanced survival of advantageous mutations, and the conversion of non-additive genetic variance into additive variance in an epistatic system, another empirically confirmed prediction. Templeton combined the theories of Mayr and Carson with population genetic models to predict the conditions under which founder events can contribute to speciation, and these predictions are strongly confirmed by the empirical literature. Much of the criticism of founder speciation is based upon equating founder speciation to an adaptive peak shift opposed by selection. However, Mayr, Carson and Templeton all modeled a positive interaction of selection and drift, and Templeton showed that founder speciation is incompatible with peak-shift conditions. Although rare, founder speciation can have a disproportionate importance in adaptive innovation and radiation, and examples are given to show that "rare" does not mean "unimportant" in evolution. Founder speciation also interacts with other speciation mechanisms such that a speciation event is not a one-dimensional process due to either selection alone or drift alone.     

Bayesian analysis suggests that most amino acid replacements in Drosophila are driven by positive selection

One of the principal goals of population genetics is to understand the processes by which genetic variation within species (polymorphism) becomes converted into genetic differences between species (divergence). In this transformation, selective neutrality, near neutrality, and positive selection may each play a role, differing from one gene to the next. Synonymous nucleotide sites are often used as a uniform standard of comparison across genes on the grounds that synonymous sites are subject to relatively weak selective constraints and so may, to a first approximation, be regarded as neutral. Synonymous sites are also interdigitated with nonsynonymous sites and so are affected equally by genomic context and demographic factors. Hence a comparison of levels of polymorphism and divergence between synonymous sites and amino acid replacement sites in a gene is potentially informative about the magnitude of selective forces associated with amino acid replacements. We have analyzed 56 genes in which polymorphism data from D. simulans are compared with divergence from a reference strain of D. melanogaster. The framework of the analysis is Bayesian and assumes that the distribution of selective effects (Malthusian fitnesses) is Gaussian with a mean that differs for each gene. In such a model, the average scaled selection intensity (gamma = N(e)s) of amino acid replacements eligible to become polymorphic or fixed is -7.31, and the standard deviation of selective effects within each locus is 6.79 (assuming homoscedasticity across loci). For newly arising mutations of this type that occur in autosomal or X-linked genes, the average proportion of beneficial mutations is 19.7%. Among the amino acid polymorphisms in the sample, the expected average proportion of beneficial mutations is 47.7%, and among amino acid replacements that become fixed the average proportion of beneficial mutations is 94.3%. The average scaled selection intensity of fixed mutations is +5.1. The presence of positive selection is pervasive with the single exception of kl-5, a Y-linked fertility gene. We find no evidence that a significant fraction of fixed amino acid replacements is neutral or nearly neutral or that positive selection drives amino acid replacements at only a subset of the loci. These results are model dependent and we discuss possible modifications of the model that might allow more neutral and nearly neutral amino acid replacements to be fixed.