The Effect of Adaptive Mutagenesis on Genetic Variation at a Linked, Neutral Locus

Based on recent studies in single-celled organisms, it has been argued that a fitness benefit associated with a mutation will increase the probability of that mutation occurring. This increase is independent of mutation rates at other loci and is called adaptive mutagenesis. We modeled the effect of adaptive mutagenesis on populations of haploid organisms with adaptive mutation rates ranging from 0 to 1 X 10(-5). Allele frequencies at the selected locus and a neutral linked locus were tracked. We also observed the amount of linkage disequilibrium during the selective sweep and the final heterozygosity after the sweep. The presence of adaptive mutagenesis increases the number of genetic backgrounds carrying the new fitter allele, making the outcomes more representative of the population before the selection. Therefore, more neutral genetic variation is preserved in simulations with adaptive mutagenesis than in those without it due to hitchhiking. Since adaptive mutagenesis is time-dependent, it can generate mutants when other mechanisms of mutation cannot. In addition, adaptive mutagenesis has the potential to confound both phylogeny construction and the detection of natural selection from patterns of nucleotide variation.     

Hitchhiking effect of a beneficial mutation spreading in a subdivided population

A central problem in population genetics is to detect and analyze positive natural selection by which beneficial mutations are driven to fixation. The hitchhiking effect of a rapidly spreading beneficial mutation, which results in local removal of standing genetic variation, allows such an analysis using DNA sequence polymorphism. However, the current mathematical theory that predicts the pattern of genetic hitchhiking relies on the assumption that a beneficial mutation increases to a high frequency in a single random-mating population, which is certainly violated in reality. Individuals in natural populations are distributed over a geographic space. The spread of a beneficial allele can be delayed by limited migration of individuals over the space and its hitchhiking effect can also be affected. To study this effect of geographic structure on genetic hitchhiking, we analyze a simple model of directional selection in a subdivided population. In contrast to previous studies on hitchhiking in subdivided populations, we mainly investigate the range of sufficiently high migration rates that would homogenize genetic variation at neutral loci. We provide a heuristic mathematical analysis that describes how the genealogical structure at a neutral locus linked to the locus under selection is expected to change in a population divided into two demes. Our results indicate that the overall strength of genetic hitchhiking--the degree to which expected heterozygosity decreases--is diminished by population subdivision, mainly because opportunity for the breakdown of hitchhiking by recombination increases as the spread of the beneficial mutation across demes is delayed when migration rate is much smaller than the strength of selection. Furthermore, the amount of genetic variation after a selective sweep is expected to be unequal over demes: a greater reduction in expected heterozygosity occurs in the subpopulation from which the beneficial mutation originates than in its neighboring subpopulations. This raises a possibility of detecting a "hidden" geographic structure of population by carefully analyzing the pattern of a selective sweep.     

Evolution of herbicide resistance in weeds: vertically transmitted fungal endophytes as genetic entities

The appearance of heritable resistance to herbicides in weeds is an evolutionary process driven by human selection. Assuming that spontaneous and random mutations originate herbicide resistance genes, which are selected by selection pressure imposed by herbicides, is the simplest model to understand how this phenomenon appears and increases in weed populations. However, the rate of herbicide resistance evolution is not only determined by the amount of genetic variation within the populations and the selection pressure exerted by herbicides, but also by factors related to genetics, biology and ecology of weeds. The inheritance of the resistance genes, the mating patterns of the populations, the relative fitness of susceptible and resistant phenotypes and gene flow processes also control the mentioned rate. Many cool season grasses are often infected by fungal symbiotic endophytes (Neotyphodium spp.). These organisms modify the physiology, ecology and reproductive biology of their hosts, conferring greater tolerance to biotic and abiotic stresses, greater competitive ability and the capacity of reducing ecosystem biodiversity. In this work, we present new empirical data and propose new theoretical support on how these microbial symbionts can modulate the evolution of herbicide resistance in weeds. Fungal endophytes are vertically transmitted, and may act as genetic entities altering the evolution of herbicide resistance by reducing herbicide efficacy (delaying effect on evolution). In addition, indirect evidence suggests that fungal endophytes might reduce the fitness penalty associated with the newly arisen resistant phenotypes. The importance and dynamic of these opposite effects is discussed.     

DNA Analysis of Herbarium Specimens of the Grass Weed Alopecurus myosuroides Reveals Herbicide Resistance Pre-Dated Herbicides

Acetyl-CoA carboxylase (ACCase) alleles carrying one point mutation that confers resistance to herbicides have been identified in arable grass weed populations where resistance has evolved under the selective pressure of herbicides. In an effort to determine whether herbicide resistance evolves from newly arisen mutations or from standing genetic variation in weed populations, we used herbarium specimens of the grass weed Alopecurus myosuroides to seek mutant ACCase alleles carrying an isoleucine-to-leucine substitution at codon 1781 that endows herbicide resistance. These specimens had been collected between 1788 and 1975, i.e., prior to the commercial release of herbicides inhibiting ACCase. Among the 734 specimens investigated, 685 yielded DNA suitable for PCR. Genotyping the ACCase locus using the derived Cleaved Amplified Polymorphic Sequence (dCAPS) technique identified one heterozygous mutant specimen that had been collected in 1888. Occurrence of a mutant codon encoding a leucine residue at codon 1781 at the heterozygous state was confirmed in this specimen by sequencing, clearly demonstrating that resistance to herbicides can pre-date herbicides in weeds. We conclude that point mutations endowing resistance to herbicides without having associated deleterious pleiotropic effects can be present in weed populations as part of their standing genetic variation, in frequencies higher than the mutation frequency, thereby facilitating their subsequent selection by herbicide applications.     

Three Signatures of Adaptive Polymorphism Exemplified by Malaria-Associated Genes

Malaria has been one of the strongest selective pressures on our species. Many of the best-characterized cases of adaptive evolution in humans are in genes tied to malaria resistance. However, the complex evolutionary patterns at these genes are poorly captured by standard scans for nonneutral evolution. Here, we present three new statistical tests for selection based on population genetic patterns that are observed more than once among key malaria resistance loci. We assess these tests using forward-time evolutionary simulations and apply them to global whole-genome sequencing data from humans, and thus we show that they are effective at distinguishing selection from neutrality. Each test captures a distinct evolutionary pattern, here called Divergent Haplotypes, Repeated Shifts, and Arrested Sweeps, associated with a particular period of human prehistory. We clarify the selective signatures at known malaria-relevant genes and identify additional genes showing similar adaptive evolutionary patterns. Among our top outliers, we see a particular enrichment for genes involved in erythropoiesis and for genes previously associated with malaria resistance, consistent with a major role for malaria in shaping these patterns of genetic diversity. Polymorphisms at these genes are likely to impact resistance to malaria infection and contribute to ongoing host–parasite coevolutionary dynamics.     

Soft sweeps: molecular population genetics of adaptation from standing genetic variation

A population can adapt to a rapid environmental change or habitat expansion in two ways. It may adapt either through new beneficial mutations that subsequently sweep through the population or by using alleles from the standing genetic variation. We use diffusion theory to calculate the probabilities for selective adaptations and find a large increase in the fixation probability for weak substitutions, if alleles originate from the standing genetic variation. We then determine the parameter regions where each scenario-standing variation vs. new mutations-is more likely. Adaptations from the standing genetic variation are favored if either the selective advantage is weak or the selection coefficient and the mutation rate are both high. Finally, we analyze the probability of "soft sweeps," where multiple copies of the selected allele contribute to a substitution, and discuss the consequences for the footprint of selection on linked neutral variation. We find that soft sweeps with weaker selective footprints are likely under both scenarios if the mutation rate and/or the selection coefficient is high.     

Soft Selective Sweeps in Complex Demographic Scenarios

Adaptation from de novo mutation can produce so-called soft selective sweeps, where adaptive alleles of independent mutational origin sweep through the population at the same time. Population genetic theory predicts that such soft sweeps should be likely if the product of the population size and the mutation rate toward the adaptive allele is sufficiently large, such that multiple adaptive mutations can establish before one has reached fixation; however, it remains unclear how demographic processes affect the probability of observing soft sweeps. Here we extend the theory of soft selective sweeps to realistic demographic scenarios that allow for changes in population size over time. We first show that population bottlenecks can lead to the removal of all but one adaptive lineage from an initially soft selective sweep. The parameter regime under which such “hardening” of soft selective sweeps is likely is determined by a simple heuristic condition. We further develop a generalized analytical framework, based on an extension of the coalescent process, for calculating the probability of soft sweeps under arbitrary demographic scenarios. Two important limits emerge within this analytical framework: In the limit where population-size fluctuations are fast compared to the duration of the sweep, the likelihood of soft sweeps is determined by the harmonic mean of the variance effective population size estimated over the duration of the sweep; in the opposing slow fluctuation limit, the likelihood of soft sweeps is determined by the instantaneous variance effective population size at the onset of the sweep. We show that as a consequence of this finding the probability of observing soft sweeps becomes a function of the strength of selection. Specifically, in species with sharply fluctuating population size, strong selection is more likely to produce soft sweeps than weak selection. Our results highlight the importance of accurate demographic estimates over short evolutionary timescales for understanding the population genetics of adaptation from de novo mutation.     

Weak selection and protein evolution

The "nearly neutral" theory of molecular evolution proposes that many features of genomes arise from the interaction of three weak evolutionary forces: mutation, genetic drift, and natural selection acting at its limit of efficacy. Such forces generally have little impact on allele frequencies within populations from generation to generation but can have substantial effects on long-term evolution. The evolutionary dynamics of weakly selected mutations are highly sensitive to population size, and near neutrality was initially proposed as an adjustment to the neutral theory to account for general patterns in available protein and DNA variation data. Here, we review the motivation for the nearly neutral theory, discuss the structure of the model and its predictions, and evaluate current empirical support for interactions among weak evolutionary forces in protein evolution. Near neutrality may be a prevalent mode of evolution across a range of functional categories of mutations and taxa. However, multiple evolutionary mechanisms (including adaptive evolution, linked selection, changes in fitness-effect distributions, and weak selection) can often explain the same patterns of genome variation. Strong parameter sensitivity remains a limitation of the nearly neutral model, and we discuss concave fitness functions as a plausible underlying basis for weak selection.     

Inferring Adaptive Codon Preference to Understand Sources of Selection Shaping Codon Usage Bias

Alternative synonymous codons are often used at unequal frequencies. Classically, studies of such codon usage bias (CUB) attempted to separate the impact of neutral from selective forces by assuming that deviations from a predicted neutral equilibrium capture selection. However, GC-biased gene conversion (gBGC) can also cause deviation from a neutral null. Alternatively, selection has been inferred from CUB in highly expressed genes, but the accuracy of this approach has not been extensively tested, and gBGC can interfere with such extrapolations (e.g., if expression and gene conversion rates covary). It is therefore critical to examine deviations from a mutational null in a species with no gBGC. To achieve this goal, we implement such an analysis in the highly AT rich genome of Dictyostelium discoideum, where we find no evidence of gBGC. We infer neutral CUB under mutational equilibrium to quantify “adaptive codon preference,” a nontautologous genome wide quantitative measure of the relative selection strength driving CUB. We observe signatures of purifying selection consistent with selection favoring adaptive codon preference. Preferred codons are not GC rich, underscoring the independence from gBGC. Expression-associated “preference” largely matches adaptive codon preference but does not wholly capture the influence of selection shaping patterns across all genes, suggesting selective constraints associated specifically with high expression. We observe patterns consistent with effects on mRNA translation and stability shaping adaptive codon preference. Thus, our approach to quantifying adaptive codon preference provides a framework for inferring the sources of selection that shape CUB across different contexts within the genome.     

From bad to good: Fitness reversals and the ascent of deleterious mutations

Deleterious mutations are considered a major impediment to adaptation, and there are straightforward expectations for the rate at which they accumulate as a function of population size and mutation rate. In a simulation model of an evolving population of asexually replicating RNA molecules, initially deleterious mutations accumulated at rates nearly equal to that of initially beneficial mutations, without impeding evolutionary progress. As the mutation rate was increased within a moderate range, deleterious mutation accumulation and mean fitness improvement both increased. The fixation rates were higher than predicted by many population-genetic models. This seemingly paradoxical result was resolved in part by the observation that, during the time to fixation, the selection coefficient (s) of initially deleterious mutations reversed to confer a selective advantage. Significantly, more than half of the fixations of initially deleterious mutations involved fitness reversals. These fitness reversals had a substantial effect on the total fitness of the genome and thus contributed to its success in the population. Despite the relative importance of fitness reversals, however, the probabilities of fixation for both initially beneficial and initially deleterious mutations were exceedingly small (on the order of 10⁻⁵ of all mutations).     

The rate of fitness-valley crossing in sexual populations

Biological traits result in part from interactions between different genetic loci. This can lead to sign epistasis, in which a beneficial adaptation involves a combination of individually deleterious or neutral mutations; in this case, a population must cross a "fitness valley" to adapt. Recombination can assist this process by combining mutations from different individuals or retard it by breaking up the adaptive combination. Here, we analyze the simplest fitness valley, in which an adaptation requires one mutation at each of two loci to provide a fitness benefit. We present a theoretical analysis of the effect of recombination on the valley-crossing process across the full spectrum of possible parameter regimes. We find that low recombination rates can speed up valley crossing relative to the asexual case, while higher recombination rates slow down valley crossing, with the transition between the two regimes occurring when the recombination rate between the loci is approximately equal to the selective advantage provided by the adaptation. In large populations, if the recombination rate is high and selection against single mutants is substantial, the time to cross the valley grows exponentially with population size, effectively meaning that the population cannot acquire the adaptation. Recombination at the optimal (low) rate can reduce the valley-crossing time by up to several orders of magnitude relative to that in an asexual population.     

Environmental heterogeneity enhances clonal interference

Clonal interference (CI) is a phenomenon that may be important in several asexual microbes. It occurs when population sizes are large and mutation rates to new beneficial alleles are of significant magnitude. Here we explore the role of gene flow and spatial heterogeneity in selection strength in the adaptation of asexuals. We consider a subdivided population of individuals that are adapting, through new beneficial mutations, and that migrate between different patches. The fitness effect of each mutation depends on the patch and all mutations considered are assumed to be unconditionally beneficial. We find that spatial variation in selection pressure affects the rate of adaptive evolution and its qualitative effects depend on the level of gene flow. In particular, we find that both low migration and high levels of heterogeneity lead to enhanced CI. In contrast, for high levels of migration the rate of fixation of adaptive mutations is higher when environmental heterogeneity is present. In addition, we observe that the level of fitness variation is higher and simultaneous fixation of multiple mutations tends to occur in the regime of low migration rates and high heterogeneity.     

Cloning, mutagenesis, and characterization of the microalga Parietochloris incisa acetohydroxyacid synthase, and its possible use as an endogenous selection marker

Parietochloris incisa is an oleaginous fresh water green microalga that accumulates an unusually high content of the valuable long-chain polyunsaturated fatty acid (LC-PUFA) arachidonic acid within triacylglycerols in cytoplasmic lipid bodies. Here, we describe cloning and mutagenesis of the P. incisa acetohydroxyacid synthase (PiAHAS) gene for use as an herbicide resistance selection marker for transformation. Use of an endogenous gene circumvents the risks and regulatory difficulties of cultivating antibiotic-resistant organisms. AHAS is present in plants and microorganisms where it catalyzes the first essential step in the synthesis of branched-chain amino acids. It is the target enzyme of the herbicide sulfometuron methyl (SMM), which effectively inhibits growth of bacteria and plants. Several point mutations of AHAS are known to confer herbicide resistance. We cloned the cDNA that encodes PiAHAS and introduced a W605S point mutation (PimAHAS). Catalytic activity and herbicide resistance of the wild-type and mutant proteins were characterized in the AHAS-deficient E. coli, BUM1 strain. Cloned PiAHAS wild-type and mutant genes complemented AHAS-deficient bacterial growth. Furthermore, bacteria expressing the mutant PiAHAS exhibited high resistance to SMM. Purified PiAHAS wild-type and mutant proteins were assayed for enzymatic activity and herbicide resistance. The W605S mutation was shown to cause a twofold decrease in enzymatic activity and in affinity for the Pyruvate substrate. However, the mutant exhibited 7 orders of magnitude higher resistance to the SMM herbicide than that of the wild type.     

The distribution of beneficial and fixed mutation fitness effects close to an optimum

The distribution of the selection coefficients of beneficial mutations is pivotal to the study of the adaptive process, both at the organismal level (theories of adaptation) and at the gene level (molecular evolution). A now famous result of extreme value theory states that this distribution is an exponential, at least when considering a well-adapted wild type. However, this prediction could be inaccurate under selection for an optimum (because fitness effect distributions have a finite right tail in this case). In this article, we derive the distribution of beneficial mutation effects under a general model of stabilizing selection, with arbitrary selective and mutational covariance between a finite set of traits. We assume a well-adapted wild type, thus taking advantage of the robustness of tail behaviors, as in extreme value theory. We show that, under these general conditions, both beneficial mutation effects and fixed effects (mutations escaping drift loss) are beta distributed. In both cases, the parameters have explicit biological meaning and are empirically measurable; their variation through time can also be predicted. We retrieve the classic exponential distribution as a subcase of the beta when there are a moderate to large number of weakly correlated traits under selection. In this case too, we provide an explicit biological interpretation of the parameters of the distribution. We show by simulations that these conclusions are fairly robust to a lower adaptation of the wild type and discuss the relevance of our findings in the context of adaptation theories and experimental evolution.     

The classical hitchhiking model with continuous mutational pressure and purifying selection

Detecting selective sweeps driven by strong positive selection and localizing the targets of selection in the genome play a major role in modern population genetics and genomics. Most of these analyses are based on the classical model of genetic hitchhiking proposed by Maynard Smith and Haigh (1974, Genetical Research, 23, 23). Here, we consider extensions of the classical two‐locus model. Introducing mutation at the strongly selected site, we analyze the conditions under which soft sweeps may arise. We identify a new parameter (the ratio of the beneficial mutation rate to the selection coefficient) that characterizes the occurrence of multiple‐origin soft sweeps. Furthermore, we quantify the hitchhiking effect when the polymorphism at the linked locus is not neutral but maintained in a mutation‐selection balance. In this case, we find a smaller relative reduction of heterozygosity at the linked site than for a neutral polymorphism. In our analysis, we use a semi‐deterministic approach; i.e., we analyze the frequency process of the beneficial allele in an infinitely large population when its frequency is above a certain threshold; however, for very small frequencies in the initial phase after the onset of selection we rely on diffusion theory.     

Reduced variation around drug-resistant dhfr alleles in African Plasmodium falciparum

We have measured microsatellite diversity at 26 markers around the dhfr gene in pyrimethamine-sensitive and -resistant parasites collected in southeast Africa. Through direct comparison with diversity on sensitive chromosomes we have found significant loss of diversity across a region of 70 kb around the most highly resistant allele which is evidence of a selective sweep attributable to selection through widespread use of pyrimethamine (in combination with sulfadoxine) as treatment for malaria. Retrospective analysis through four years of direct and continuous selection from use of sulfadoxine-pyrimethamine as first-line malaria treatment on a Plasmodium falciparum population in KwaZulu Natal, South Africa, has revealed how recombination significantly narrowed the margins of the selective sweep over time. A deterministic model incorporating selection coefficients measured during the same interval indicates that the transition was toward a state of recombination-selection equilibrium. We compared loss of diversity around the same resistance allele in two populations at either extreme of the range of entomological inoculation rates (EIRs), namely, under one infective bite per year in Mpumalanga, South Africa, and more than one per day in southern Tanzania. EIRs determine effective recombination rates and are expected to profoundly influence the dimensions of the selective sweep. Surprisingly, the dimensions were broadly consistent across both populations. We conclude that despite different recombination rates and contrasting drug selection histories in neighboring countries, the region-wide movement of resistant parasites has played a key role in the establishment of resistance in these populations and the dimensions of the selective sweep are dominated by the influence of high initial starting frequencies.     

Drought exposure leads to rapid acquisition and inheritance of herbicide resistance in the weed Alopecurus myosuroides

Globally, herbicide resistance in weeds poses a threat to food security. Resistance evolves rapidly through the co‐option of a suite of physiological mechanisms that evolved to allow plants to survive environmental stress. Consequently, we hypothesize that stress tolerance and herbicide resistance are functionally linked. We address two questions: (i) does exposure to stress in a parental generation promote the evolution of resistance in the offspring? (ii) Is such evolution mediated through non‐genetic mechanisms? We exposed individuals of a grass weed to drought, and tested whether this resulted in herbicide resistance in the first generation. In terms of both survival and dry mass, we find enhanced resistance to herbicide in the offspring of parents that had been exposed to drought. Our results suggest that exposure of weeds to drought can confer herbicide resistance in subsequent generations, and that the mechanism conferring heritability of herbicide resistance is non‐genetic.     

Costs of transgenic herbicide resistance introgressed from Brassica napus into weedy B. rapa

Wild relatives of genetically engineered crops can acquire transgenic traits such as herbicide resistance via spontaneous crop–wild hybridization. In agricultural weeds, resistance to herbicides is often a beneficial trait, but little is known about possible costs that could affect the persistence of this trait when herbicides are not used. We tested for costs associated with transgenic resistance to glufosinate when introgressed into weedy Brassica rapa . Crosses were made between transgenic B. napus and wild B. rapa from Denmark. F₁ progeny were backcrossed to B. rapa and BC₁ plants were selected for chromosome numbers similar to B. rapa . Further backcrossing resulted in a BC₂ generation that was hemizygous for herbicide resistance. We quantified the reproductive success of 457 BC₃ progeny representing six full-sib families raised in growth rooms (plants were pollinated by captive bumblebees). Pollen fertility and seed production of BC₃ plants were as great as those of B. rapa raised in the same growth rooms. Segregation for herbicide resistance in BC₃ plants was 1:1 overall, but the frequency of resistant progeny was lower than expected in one family and higher than expected in another. There were no significant differences between transgenic and nontransgenic plants in survival or the number of seeds per plant, indicating that costs associated with the transgene are probably negligible. Results from this growth-chamber study suggest that transgenic resistance to glufosinate is capable of introgressing into populations of B. rapa and persisting, even in the absence of selection due to herbicide application.