Assortative mating for fitness and the evolution of recombination

To understand selection on recombination, we need to consider how linkage disequilibria develop and how recombination alters these disequilibria. Any factor that affects the development of disequilibria, including nonrandom mating, can potentially change selection on recombination. Assortative mating is known to affect linkage disequilibria but its effects on the evolution of recombination have not been previously studied. Given that assortative mating for fitness can arise indirectly via a number of biologically realistic scenarios, it is plausible that weak assortative mating occurs across a diverse set of taxa. Using a modifier model, we examine how assortative mating for fitness affects the evolution of recombination under two evolutionary scenarios: selective sweeps and mutation-selection balance. We find there is no net effect of assortative mating during a selective sweep. In contrast, assortative mating could have a large effect on recombination when deleterious alleles are maintained at mutation-selection balance but only if assortative mating is sufficiently strong. Upon considering reasonable values for the number of loci affecting fitness components, the strength of selection, and the mutation rate, we conclude that the correlation in fitness between mates is unlikely to be sufficiently high for assortative mating to affect the evolution of recombination in most species.     

Effects of Interference Between Selected Loci on the Mutation Load,Inbreeding Depression,and Heterosis

A classical prediction from single-locus models is that inbreeding increases the efficiency of selection against partially recessive deleterious alleles (purging), thereby decreasing the mutation load and level of inbreeding depression. However, previous multilocus simulation studies found that increasing the rate of self-fertilization of individuals may not lead to purging and argued that selective interference among loci causes this effect. In this article, I derive simple analytical approximations for the mutation load and inbreeding depression, taking into account the effects of interference between pairs of loci. I consider two classical scenarios of nonrandomly mating populations: a single population undergoing partial selfing and a subdivided population with limited dispersal. In the first case, correlations in homozygosity between loci tend to reduce mean fitness and increase inbreeding depression. These effects are stronger when deleterious alleles are more recessive, but only weakly depend on the strength of selection against deleterious alleles and on recombination rates. In subdivided populations, interference increases inbreeding depression within demes, but decreases heterosis between demes. Comparisons with multilocus, individual-based simulations show that these analytical approximations are accurate as long as the effects of interference stay moderate, but fail for high deleterious mutation rates and low dominance coefficients of deleterious alleles.     

A resolution of the mutation load paradox in humans

Current information on the rate of mutation and the fraction of sites in the genome that are subject to selection suggests that each human has received, on average, at least two new harmful mutations from its parents. These mutations were subsequently removed by natural selection through reduced survival or fertility. It has been argued that the mutation load, the proportional reduction in population mean fitness relative to the fitness of an idealized mutation-free individual, allows a theoretical prediction of the proportion of individuals in the population that fail to reproduce as a consequence of these harmful mutations. Application of this theory to humans implies that at least 88% of individuals should fail to reproduce and that each female would need to have more than 16 offspring to maintain population size. This prediction is clearly at odds with the low reproductive excess of human populations. Here, we derive expressions for the fraction of individuals that fail to reproduce as a consequence of recurrent deleterious mutation () for a model in which selection occurs via differences in relative fitness, such as would occur through competition between individuals. We show that is much smaller than the value predicted by comparing fitness to that of a mutation-free genotype. Under the relative fitness model, we show that depends jointly on U and the selective effects of new deleterious mutations and that a species could tolerate 10's or even 100's of new deleterious mutations per genome each generation.     

Mutation-selection balance in mixed mating populations

An approximation to the average number of deleterious mutations per gamete, Q, is derived from a model allowing selection on both zygotes and male gametes. Progeny are produced by either outcrossing or self-fertilization with fixed probabilities. The genetic model is a standard in evolutionary biology: mutations occur at unlinked loci, have equivalent effects, and combine multiplicatively to determine fitness. The approximation developed here treats individual mutation counts with a generalized Poisson model conditioned on the distribution of selfing histories in the population. The approximation is accurate across the range of parameter sets considered and provides both analytical insights and greatly increased computational speed. Model predictions are discussed in relation to several outstanding problems, including the estimation of the genomic deleterious mutation rates (U), the generality of "selective interference" among loci, and the consequences of gametic selection for the joint distribution of inbreeding depression and mating system across species. Finally, conflicting results from previous analytical treatments of mutation-selection balance are resolved to assumptions about the life-cycle and the initial fate of mutations.     

Beneficial mutations, hitchhiking and the evolution of mutation rates in sexual populations

Natural selection acts in three ways on heritable variation for mutation rates. A modifier allele that increases the mutation rate is (i) disfavored due to association with deleterious mutations, but is also favored due to (ii) association with beneficial mutations and (iii) the reduced costs of lower fidelity replication. When a unique beneficial mutation arises and sweeps to fixation, genetic hitchhiking may cause a substantial change in the frequency of a modifier of mutation rate. In previous studies of the evolution of mutation rates in sexual populations, this effect has been underestimated. This article models the long-term effect of a series of such hitchhiking events and determines the resulting strength of indirect selection on the modifier. This is compared to the indirect selection due to deleterious mutations, when both types of mutations are randomly scattered over a given genetic map. Relative to an asexual population, increased levels of recombination reduce the effects of beneficial mutations more rapidly than those of deleterious mutations. However, the role of beneficial mutations in determining the evolutionarily stable mutation rate may still be significant if the function describing the cost of high-fidelity replication has a shallow gradient.     

Hitchhiking of Deleterious Alleles and the Cost of Adaptation in Partially Selfing Species

Self-fertilization is generally seen to be disadvantageous in the long term. It increases genetic drift, which subsequently reduces polymorphism and the efficiency of selection, which also challenges adaptation. However, high selfing rates can increase the fixation probability of recessive beneficial mutations, but existing theory has generally not accounted for the effect of linked sites. Here, we analyze a model for the fixation probability of deleterious mutants that hitchhike with selective sweeps in diploid, partially selfing populations. Approximate analytical solutions show that, conditional on the sweep not being lost by drift, higher inbreeding rates increase the fixation probability of the deleterious allele, due to the resulting reduction in polymorphism and effective recombination. When extending the analysis to consider a distribution of deleterious alleles, as well as the average fitness increase after a sweep, we find that beneficial alleles generally need to be more recessive than the previously assumed dominance threshold (h < 1/2) for selfing to be beneficial from one-locus theory. Our results highlight that recombination aiding the efficiency of selection on multiple loci amplifies the fitness benefits of outcrossing over selfing, compared to results obtained from one-locus theory. This effect additionally increases the parameter range under which obligate outcrossing is beneficial over partial selfing.     

Evolution of mutation rates in bacteria

Evolutionary success of bacteria relies on the constant fine-tuning of their mutation rates, which optimizes their adaptability to constantly changing environmental conditions. When adaptation is limited by the mutation supply rate, under some conditions, natural selection favours increased mutation rates by acting on allelic variation of the genetic systems that control fidelity of DNA replication and repair. Mutator alleles are carried to high frequency through hitchhiking with the adaptive mutations they generate. However, when fitness gain no longer counterbalances the fitness loss due to continuous generation of deleterious mutations, natural selection favours reduction of mutation rates. Selection and counter-selection of high mutation rates depends on many factors: the number of mutations required for adaptation, the strength of mutator alleles, bacterial population size, competition with other strains, migration, and spatial and temporal environmental heterogeneity. Such modulations of mutation rates may also play a role in the evolution of antibiotic resistance.     

Sexual selection on spontaneous mutations strengthens the between‐sex genetic correlation for fitness

A proposed benefit to sexual selection is that it promotes purging of deleterious mutations from populations. For this benefit to be realized, sexual selection, which is usually stronger on males, must purge mutations deleterious to both sexes. Here, we experimentally test the hypothesis that sexual selection on males purges deleterious mutations that affect both male and female fitness. We measured male and female fitness in two panels of spontaneous mutation‐accumulation lines of the fly, Drosophila serrata, each established from a common ancestor. One panel of mutation accumulation lines limited both natural and sexual selection (LS lines), whereas the other panel limited natural selection, but allowed sexual selection to operate (SS lines). Although mutation accumulation caused a significant reduction in male and female fitness in both the LS and SS lines, sexual selection had no detectable effect on the extent of the fitness reduction. Similarly, despite evidence of mutational variance for fitness in males and females of both treatments, sexual selection had no significant impact on the amount of mutational genetic variance for fitness. However, sexual selection did reshape the between‐sex correlation for fitness: significantly strengthening it in the SS lines. After 25 generations, the between‐sex correlation for fitness was positive but considerably less than one in the LS lines, suggesting that, although most mutations had sexually concordant fitness effects, sex‐limited, and/or sex‐biased mutations contributed substantially to the mutational variance. In the SS lines this correlation was strong and could not be distinguished from unity. Individual‐based simulations that mimick the experimental setup reveal two conditions that may drive our results: (1) a modest‐to‐large fraction of mutations have sex‐limited (or highly sex‐biased) fitness effects, and (2) the average fitness effect of sex‐limited mutations is larger than the average fitness effect of mutations that affect both sexes similarly.     

Nature of Deleterious Mutation Load in Drosophila

Much population genetics and evolution theory depends on knowledge of genomic mutation rates and distributions of mutation effects for fitness, but most information comes from a few mutation accumulation experiments in Drosophila in which replicated chromosomes are sheltered from natural selection by a balancer chromosome. I show here that data from these experiments imply the existence of a large class of minor viability mutations with approximately equivalent effects. However, analysis of the distribution of viabilities of chromosomes exposed to EMS mutagenesis reveals a qualitatively different distribution of effects lacking such a minor effects class. A possible explanation for this difference is that transposable element insertions, a common class of spontaneous mutation event in Drosophila, frequently generate minor viability effects. This explanation would imply that current estimates of deleterious mutation rates are not generally applicable in evolutionary models, as transposition rates vary widely. Alternatively, much of the apparent decline in viability under spontaneous mutation accumulation could have been nonmutational, perhaps due to selective improvement of balancer chromosomes. This explanation accords well with the data and implies a spontaneous mutation rate for viability two orders of magnitude lower than previously assumed, with most mutation load attributable to major effects.     

Investigating the Consequences of Interference between Multiple CD8+ T Cell Escape Mutations in Early HIV Infection

During early human immunodeficiency virus (HIV) infection multiple CD8⁺ T cell responses are elicited almost simultaneously. These responses exert strong selective pressures on different parts of HIV’s genome, and select for mutations that escape recognition and are thus beneficial to the virus. Some studies reveal that the later these escape mutations emerge, the more slowly they go to fixation. This pattern of escape rate decrease(ERD) can arise by distinct mechanisms. In particular, in large populations with high beneficial mutation rates interference among different escape strains –an effect that can emerge in evolution with asexual reproduction and results in delayed fixation times of beneficial mutations compared to sexual reproduction– could significantly impact the escape rates of mutations. In this paper, we investigated how interference between these concurrent escape mutations affects their escape rates in systems with multiple epitopes, and whether it could be a source of the ERD pattern. To address these issues, we developed a multilocus Wright-Fisher model of HIV dynamics with selection, mutation and recombination, serving as a null-model for interference. We also derived an interference-free null model assuming initial neutral evolution before immune response elicitation. We found that interference between several equally selectively advantageous mutations can generate the observed ERD pattern. We also found that the number of loci, as well as recombination rates substantially affect ERD. These effects can be explained by the underexponential decline of escape rates over time. Lastly, we found that the observed ERD pattern in HIV infected individuals is consistent with both independent, interference-free mutations as well as interference effects. Our results confirm that interference effects should be considered when analyzing HIV escape mutations. The challenge in estimating escape rates and mutation-associated selective coefficients posed by interference effects cannot simply be overcome by improved sampling frequencies or sizes. This problem is a consequence of the fundamental shortcomings of current estimation techniques under interference regimes. Hence, accounting for the stochastic nature of competition between mutations demands novel estimation methodologies based on the analysis of HIV strains, rather than mutation frequencies.     

Consequences of life history for inbreeding depression and mating system evolution in plants

Many plants are perennial, but most studies of inbreeding depression and mating system evolution focus on annuals. This paper extends a population genetic model of inbreeding depression due to recessive deleterious mutations to perennials. The model incorporates life history and mating system variation, and multiplicative selection across many genetic loci. In the absence of substantial mitotic mutation, perennials have higher mean fitness and lower, or even negative, inbreeding depression than annuals with the same mating system. As in annuals, self fertilization exposes deleterious recessive mutations to selection, increasing mean fitness and decreasing inbreeding depression. Including mitotic mutation decreases mean fitness while increasing inbreeding depression. Perenniality introduces a kind of selective sieve, such that strongly recessive mutations contribute disproportionately to mean fitness and inbreeding depression. In the presence of high mitotic mutation, this selective sieve may provide a mechanistic basis for high inbreeding depression observed in some long lived perennials. Without substantial mitotic mutation, it is difficult to reconcile genetically based models of inbreeding depression with the empirical generalization that perennials outcross while related annuals self fertilize.     

Hill–Robertson interference maintained by Red Queen dynamics favours the evolution of sex

Although it is well established theoretically that selective interference among mutations (Hill–Robertson interference) favours meiotic recombination, genomewide mean rates of mutation and strengths of selection appear too low to support this as the mechanism favouring recombination in nature. A possible solution to this discrepancy between theory and observation is that selection is at least intermittently very strong due to the antagonistic coevolution between a host and its parasites. The Red Queen theory posits that such coevolution generates fitness epistasis among loci, which generates negative linkage disequilibrium among beneficial mutations, which in turn favours recombination. This theory has received only limited support. However, Red Queen dynamics without epistasis may provide the ecological conditions that maintain strong and frequent selective interference in finite populations that indirectly selects for recombination. This hypothesis is developed here through the simulation of Red Queen dynamics. This approach required the development of a method to calculate the exact frequencies of multilocus haplotypes after recombination. Simulations show that recombination is favoured by the moderately weak selection of many loci involved in the interaction between a host and its parasites, which results in substitution rates that are compatible with empirical estimates. The model also reproduces the previously reported rapid increase in the rate of outcrossing in Caenorhabditis elegans coevolving with a bacterial pathogen.     

Dynamic mutation-selection balance as an evolutionary attractor

The vast majority of mutations are deleterious and are eliminated by purifying selection. Yet in finite asexual populations, purifying selection cannot completely prevent the accumulation of deleterious mutations due to Muller's ratchet: once lost by stochastic drift, the most-fit class of genotypes is lost forever. If deleterious mutations are weakly selected, Muller's ratchet can lead to a rapid degradation of population fitness. Evidently, the long-term stability of an asexual population requires an influx of beneficial mutations that continuously compensate for the accumulation of the weakly deleterious ones. Hence any stable evolutionary state of a population in a static environment must involve a dynamic mutation-selection balance, where accumulation of deleterious mutations is on average offset by the influx of beneficial mutations. We argue that such a state can exist for any population size N and mutation rate U and calculate the fraction of beneficial mutations, ε, that maintains the balanced state. We find that a surprisingly low ε suffices to achieve stability, even in small populations in the face of high mutation rates and weak selection, maintaining a well-adapted population in spite of Muller's ratchet. This may explain the maintenance of mitochondria and other asexual genomes.     

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.     

The Effects of Background and Interference Selection on Patterns of Genetic Variation in Subdivided Populations

It is well known that most new mutations that affect fitness exert deleterious effects and that natural populations are often composed of subpopulations (demes) connected by gene flow. To gain a better understanding of the joint effects of purifying selection and population structure, we focus on a scenario where an ancestral population splits into multiple demes and study neutral diversity patterns in regions linked to selected sites. In the background selection regime of strong selection, we first derive analytic equations for pairwise coalescent times and F_ST as a function of time after the ancestral population splits into two demes and then construct a flexible coalescent simulator that can generate samples under complex models such as those involving multiple demes or nonconservative migration. We have carried out extensive forward simulations to show that the new methods can accurately predict diversity patterns both in the nonequilibrium phase following the split of the ancestral population and in the equilibrium between mutation, migration, drift, and selection. In the interference selection regime of many tightly linked selected sites, forward simulations provide evidence that neutral diversity patterns obtained from both the nonequilibrium and equilibrium phases may be virtually indistinguishable for models that have identical variance in fitness, but are nonetheless different with respect to the number of selected sites and the strength of purifying selection. This equivalence in neutral diversity patterns suggests that data collected from subdivided populations may have limited power for differentiating among the selective pressures to which closely linked selected sites are subject.     

The clarifying role of time series data in the population genetics of HIV

HIV can evolve remarkably quickly in response to antiretroviral therapies and the immune system. This evolution stymies treatment effectiveness and prevents the development of an HIV vaccine. Consequently, there has been a great interest in using population genetics to disentangle the forces that govern the HIV adaptive landscape (selection, drift, mutation, and recombination). Traditional population genetics approaches look at the current state of genetic variation and infer the processes that can generate it. However, because HIV evolves rapidly, we can also sample populations repeatedly over time and watch evolution in action. In this paper, we demonstrate how time series data can bound evolutionary parameters in a way that complements and informs traditional population genetic approaches. Specifically, we focus on our recent paper (Feder et al., 2016, eLife), in which we show that, as improved HIV drugs have led to fewer patients failing therapy due to resistance evolution, less genetic diversity has been maintained following the fixation of drug resistance mutations. Because soft sweeps of multiple drug resistance mutations spreading simultaneously have been previously documented in response to the less effective HIV therapies used early in the epidemic, we interpret the maintenance of post-sweep diversity in response to poor therapies as further evidence of soft sweeps and therefore a high population mutation rate (θ) in these intra-patient HIV populations. Because improved drugs resulted in rarer resistance evolution accompanied by lower post-sweep diversity, we suggest that both observations can be explained by decreased population mutation rates and a resultant transition to hard selective sweeps. A recent paper (Harris et al., 2018, PLOS Genetics) proposed an alternative interpretation: Diversity maintenance following drug resistance evolution in response to poor therapies may have been driven by recombination during slow, hard selective sweeps of single mutations. Then, if better drugs have led to faster hard selective sweeps of resistance, recombination will have less time to rescue diversity during the sweep, recapitulating the decrease in post-sweep diversity as drugs have improved. In this paper, we use time series data to show that drug resistance evolution during ineffective treatment is very fast, providing new evidence that soft sweeps drove early HIV treatment failure.     

Background selection in single genes may explain patterns of codon bias

Background selection involves the reduction in effective population size caused by the removal of recurrent deleterious mutations from a population. Previous work has examined this process for large genomic regions. Here we focus on the level of a single gene or small group of genes and investigate how the effects of background selection caused by nonsynonymous mutations are influenced by the lengths of coding sequences, the number and length of introns, intergenic distances, neighboring genes, mutation rate, and recombination rate. We generate our predictions from estimates of the distribution of the fitness effects of nonsynonymous mutations, obtained from DNA sequence diversity data in Drosophila. Results for genes in regions with typical frequencies of crossing over in Drosophila melanogaster suggest that background selection may influence the effective population sizes of different regions of the same gene, consistent with observed differences in codon usage bias along genes. It may also help to cause the observed effects of gene length and introns on codon usage. Gene conversion plays a crucial role in determining the sizes of these effects. The model overpredicts the effects of background selection with large groups of nonrecombining genes, because it ignores Hill-Robertson interference among the mutations involved.     

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).     

Sex‐antagonistic genes,XY recombination and feminized Y chromosomes

The canonical model of sex‐chromosome evolution predicts that sex‐antagonistic (SA) genes play an instrumental role in the arrest of XY recombination and ensuing Y chromosome degeneration. Although this model might account for the highly differentiated sex chromosomes of birds and mammals, it does not fit the situation of many lineages of fish, amphibians or nonavian reptiles, where sex chromosomes are maintained homomorphic through occasional XY recombination and/or high turnover rates. Such situations call for alternative explanatory frameworks. A crucial issue at stake is the effect of XY recombination on the dynamics of SA genes and deleterious mutations. Using individual‐based simulations, we show that a complete arrest of XY recombination actually benefits females, not males. Male fitness is maximized at different XY recombination rates depending on SA selection, but never at zero XY recombination. This should consistently favour some level of XY recombination, which in turn generates a recombination load at sex‐linked SA genes. Hill–Robertson interferences with deleterious mutations also impede the differentiation of sex‐linked SA genes, to the point that males may actually fix feminized phenotypes when SA selection and XY recombination are low. We argue that sex chromosomes might not be a good localization for SA genes, and sex conflicts seem better solved through the differential expression of autosomal genes.     

QUANTITATIVE VARIATION IN FINITE PARTHENOGENETIC POPULATIONS: WHAT STOPS MULLER'S RATCHET IN THE ABSENCE OF RECOMBINATION?

Finite parthenogenetic populations with high genomic mutation rates accumulate deleterious mutations if back mutations are rare. This mechanism, known as Muller's ratchet, can explain the rarity of parthenogenetic species among so called higher organisms. However, estimates of genomic mutation rates for deleterious alleles and their average effect in the diploid condition in Drosophila suggest that Muller's ratchet should eliminate parthenogenetic insect populations within several hundred generations, provided all mutations are unconditionally deleterious. This fact is inconsistent with the existence of obligatory parthenogenetic insect species. In this paper an analysis of the extent to which compensatory mutations can counter Muller's ratchet is presented. Compensatory mutations are defined as all mutations that compensate for the phenotypic effects of a deleterious mutation. In the case of quantitative traits under stabilizing selection, the rate of compensatory mutations is easily predicted. It is shown that there is a strong analogy between the Muller's ratchet model of Felsenstein (1974) and the quantitative genetic model considered here, except for the frequency of compensatory mutations. If the intensity of stabilizing selection is too small or the mutation rate too high, the optimal genotype becomes extinct and the population mean drifts from the optimum but still reaches a stationary distribution. This distance is essentially the same as predicted for sexually reproducing populations under the same circumstances. Hence, at least in the short run, compensatory mutations for quantitative characters are as effective as recombination in halting the decline of mean fitness otherwise caused by Muller's ratchet. However, it is questionable whether compensatory mutations can prevent Muller's ratchet in the long run because there might be a limit to the capacity of the genome to provide compensatory mutations without eliminating deleterious mutations at least during occasional episodes of sex.