Male-driven biased gene conversion governs the evolution of base composition in human alu repeats

Regional biases in substitution pattern are likely to be responsible for the large-scale variation in base composition observed in vertebrate genomes. However, the evolutionary forces responsible for these biases are still not clearly defined. In order to study the processes of mutation and fixation across the entire human genome, we analyzed patterns of substitution in Alu repeats since their insertion. We also studied patterns of human polymorphism within the repeats. There is a highly significant effect of recombination rate on the pattern of substitution, whereas no such effect is seen on the pattern of polymorphism. These results suggest that regional biases in substitution are caused by biased gene conversion, a process that increases the probability of fixation of mutations that increase GC content. Furthermore, the strongest correlate of substitution patterns is found to be male recombination rates rather than female or sex-averaged recombination rates. This indicates that in addition to sexual dimorphism in recombination rates, the sexes also differ in the relative rates of crossover and gene conversion.     

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.     

Emergent neutrality in adaptive asexual evolution

In nonrecombining genomes, genetic linkage can be an important evolutionary force. Linkage generates interference interactions, by which simultaneously occurring mutations affect each other's chance of fixation. Here, we develop a comprehensive model of adaptive evolution in linked genomes, which integrates interference interactions between multiple beneficial and deleterious mutations into a unified framework. By an approximate analytical solution, we predict the fixation rates of these mutations, as well as the probabilities of beneficial and deleterious alleles at fixed genomic sites. We find that interference interactions generate a regime of emergent neutrality: all genomic sites with selection coefficients smaller in magnitude than a characteristic threshold have nearly random fixed alleles, and both beneficial and deleterious mutations at these sites have nearly neutral fixation rates. We show that this dynamic limits not only the speed of adaptation, but also a population's degree of adaptation in its current environment. We apply the model to different scenarios: stationary adaptation in a time-dependent environment and approach to equilibrium in a fixed environment. In both cases, the analytical predictions are in good agreement with numerical simulations. Our results suggest that interference can severely compromise biological functions in an adapting population, which sets viability limits on adaptive evolution under linkage.     

Effects of genotype on rates of substitution during experimental evolution

Rates of molecular evolution may vary widely between populations, yet the causes of this variation are still incompletely understood. Genetic differences between populations may make an important contribution to variation in rates of evolution, owing to differences in fitness, population size, mutation rates, or in the distribution of fitness effects (DFEs) of available beneficial mutations. By whole genome sequencing of Escherichia coli populations experimentally evolved in the presence of a quinolone antibiotic, we found that rates of substitution varied by genotype, with evidence for a contribution from a genotype's starting fitness. Subsequent targeted sequencing showed that genotypes with high average substitution rates were more likely to undergo the simultaneous fixation of several mutations, consistent with theoretical models of multiple mutation dynamics. Moreover, patterns of substitution were indicative of epistatic relationships between known resistance mutations.     

The speed of evolution and maintenance of variation in asexual populations

BACKGROUND: The rate at which beneficial mutations accumulate determines how fast asexual populations evolve, but this is only partially understood. Some recent clonal-interference models suggest that evolution in large asexual populations is limited because smaller beneficial mutations are outcompeted by larger beneficial mutations that occur in different lineages within the same population. This analysis assumes that the important mutations fix one at a time; it ignores multiple beneficial mutations that occur in the lineage of an earlier beneficial mutation, before the first mutation in the series can fix. We focus on the effects of such multiple mutations. RESULTS: Our analysis predicts that the variation in fitness maintained by a continuously evolving population increases as the logarithm of the population size and logarithm of the mutation rate and thus yields a similar logarithmic increase in the speed of evolution. To test these predictions, we evolved asexual budding yeast in glucose-limited media at a range of population sizes and mutation rates. CONCLUSIONS: We find that their evolution is dominated by the accumulation of multiple mutations of moderate effect. Our results agree with our theoretical predictions and are inconsistent with the one-by-one fixation of mutants assumed by recent clonal-interference analysis.     

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

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

Population dynamics of GC-changing mutations in humans and great apes

The nucleotide composition of the genome is a balance between the origin and fixation rates of different mutations. For example, it is well-known that transitions occur more frequently than transversions, particularly at CpG sites. Differences in fixation rates of mutation types are less explored. Specifically, recombination-associated GC-biased gene conversion (gBGC) may differentially impact GC-changing mutations, due to differences in their genomic distributions and efficiency of mismatch repair mechanisms. Given that recombination evolves rapidly across species, we explore gBGC of different mutation types across human populations and great ape species. We report a stronger correlation between segregating GC frequency and recombination for transitions than for transversions. Notably, CpG transitions are most strongly affected by gBGC in humans and chimpanzees. We show that the overall strength of gBGC is generally correlated with effective population sizes in humans, with some notable exceptions, such as a stronger effect of gBGC on non-CpG transitions in populations of European descent. Furthermore, species of the Gorilla and Pongo genus have a greatly reduced gBGC effect on CpG sites. We also study the dependence of gBGC dynamics on flanking nucleotides and show that some mutation types evolve in opposition to the gBGC expectation, likely due to the hypermutability of specific nucleotide contexts. Our results highlight the importance of different gBGC dynamics experienced by GC-changing mutations and their impact on nucleotide composition evolution.     

Evidence for widespread convergent evolution around human microsatellites

Microsatellites are a major component of the human genome, and their evolution has been much studied. However, the evolution of microsatellite flanking sequences has received less attention, with reports of both high and low mutation rates and of a tendency for microsatellites to cluster. From the human genome we generated a database of many thousands of (AC)_n flanking sequences within which we searched for common characteristics. Sequences flanking microsatellites of similar length show remarkable levels of convergent evolution, indicating shared mutational biases. These biases extend 25–50 bases either side of the microsatellite and may therefore affect more than 30% of the entire genome. To explore the extent and absolute strength of these effects, we quantified the observed convergence. We also compared homologous human and chimpanzee loci to look for evidence of changes in mutation rate around microsatellites. Most models of DNA sequence evolution assume that mutations are independent and occur randomly. Allowances may be made for sites mutating at different rates and for general mutation biases such as the faster rate of transitions over transversions. Our analysis suggests that these models may be inadequate, in that proximity to even very short microsatellites may alter the rate and distribution of mutations that occur. The elevated local mutation rate combined with sequence convergence, both of which we find evidence for, also provide a possible resolution for the apparently contradictory inferences of mutation rates in microsatellite flanking sequences.     

Male phenotypes in a female framework: Evidence for degeneration in sperm produced by male snails from asexual lineages

How changes in selective regimes affect trait evolution is an important open biological question. We take advantage of naturally occurring and repeated transitions from sexual to asexual reproduction in a New Zealand freshwater snail species, Potamopyrgus antipodarum, to address how evolution in an asexual context—including the potential for relaxed selection on male‐specific traits—influences sperm morphology. The occasional production of male offspring by the otherwise all‐female asexual P. antipodarum lineages affords a unique and powerful opportunity to assess the fate of sperm traits in a context where males are exceedingly rare. These comparisons revealed that the sperm produced by ‘asexual’ males are markedly distinct from sexual counterparts. We also found that the asexual male sperm harboured markedly higher phenotypic variation and was much more likely to be morphologically abnormal. Together, these data suggest that transitions to asexual reproduction might be irreversible, at least in part because male function is likely to be compromised. These results are also consistent with a scenario where relaxed selection and/or mutation accumulation in the absence of sex translates into rapid trait degeneration.     

FISHER'S MODEL AND THE GENOMICS OF ADAPTATION: RESTRICTED PLEIOTROPY,HETEROGENOUS MUTATION,AND PARALLEL EVOLUTION

Genetic theories of adaptation generally overlook the genes in which beneficial substitutions occur, and the likely variation in their mutational effects. We investigate the consequences of heterogeneous mutational effects among loci on the genetics of adaptation. We use a generalization of Fisher's geometrical model, which assumes multivariate Gaussian stabilizing selection on multiple characters. In our model, mutation has a distinct variance–covariance matrix of phenotypic effects for each locus. Consequently, the distribution of selection coefficients s varies across loci. We assume each locus can only affect a limited number of independent linear combinations of phenotypic traits (restricted pleiotropy), which differ among loci, an effect we term “orientation heterogeneity.” Restricted pleiotropy can sharply reduce the overall proportion of beneficial mutations. Orientation heterogeneity has little impact on the shape of the genomic distribution, but can substantially increase the probability of parallel evolution (the repeated fixation of beneficial mutations at the same gene in independent populations), which is highest with low pleiotropy. We also consider variation in the degree of pleiotropy and in the mean s across loci. The latter impacts the genomic distribution of s, but has a much milder effect on parallel evolution. We discuss these results in the light of evolution experiments.     

Assessing the underlying pattern of human germline mutations: lessons from the factor IX gene.

Germline mutations cause or predispose to most disease. Hemophilia B is a useful model for studying the underlying pattern of recent germline mutations in humans because the observed pattern of mutation in factor IX more closely reflects the underlying pattern of mutation than the observed pattern for many other genes. In addition, it is possible to identify and correct for biases inherent in ascertaining only those mutations that cause hemophilia. Aspects of the pattern of germline mutation in the factor IX gene are becoming clear: 1) in the United States, two-thirds of mutations causing mild disease arose from three founders whereas almost all the mutations resulting in either moderate or severe disease arose independently, generally within the past 150 years; 2) direct estimates of the rates of mutation in humans indicate that transitions are more frequent than transversions, which in turn are more frequent than deletions and insertions; 3) transitions at CpG are elevated approximately 24-fold relative to transitions at non-CpG dinucleotides; 4) transversions at CpG are elevated approximately eightfold relative to transversions at non-CpG dinucleotides; 5) the sum total of the dinucleotide mutation rates produces a bias against G and C bases that would be sufficient to maintain the G+C content of the factor IX gene at its evolutionarily conserved level of 40%; and 6) the pattern of mutation is similar for Caucasians residing in the United States and for Asians residing in Asia. Two ideas emerge from this and from an analysis of the pattern of recent deleterious mutations compared with ancient neutral mutations that have been fixed during evolution into the factor IX gene. First, the bulk of germline mutations are likely to arise from endogenous processes rather than environmental mutagens. Second, the factor IX protein is composed mostly of two classes of amino acids: critical residues in which all single-base missense changes will disrupt protein function, and "spacer" residues in which the precise nature of the residue is unimportant but the peptide bond is necessary to keep the critical residues in register. More work is necessary to assess the veracity and generality of these ideas.     

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

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

Fixation probability of rare nonmutator and evolution of mutation rates

Although mutations drive the evolutionary process, the rates at which the mutations occur are themselves subject to evolutionary forces. Our purpose here is to understand the role of selection and random genetic drift in the evolution of mutation rates, and we address this question in asexual populations at mutation‐selection equilibrium neglecting selective sweeps. Using a multitype branching process, we calculate the fixation probability of a rare nonmutator in a large asexual population of mutators and find that a nonmutator is more likely to fix when the deleterious mutation rate of the mutator population is high. Compensatory mutations in the mutator population are found to decrease the fixation probability of a nonmutator when the selection coefficient is large. But, surprisingly, the fixation probability changes nonmonotonically with increasing compensatory mutation rate when the selection is mild. Using these results for the fixation probability and a drift‐barrier argument, we find a novel relationship between the mutation rates and the population size. We also discuss the time to fix the nonmutator in an adapted population of asexual mutators, and compare our results with experiments.     

Effects of pollen availability and the mutation bias on the fixation of mutations disabling the male specificity of self‐incompatibility

The evolution of self‐compatibility (SC) by the loss of self‐incompatibility (SI) is regarded as one of the most frequent transitions in flowering plants. SI systems are generally characterized by specific interactions between the male and female specificity genes encoded at the S‐locus. Recent empirical studies have revealed that the evolution of SC is often driven by male SC‐conferring mutations at the S‐locus rather than by female mutations. In this study, using a forward simulation model, we compared the fixation probabilities of male vs. female SC‐conferring mutations at the S‐locus. We explicitly considered the effects of pollen availability in the population and bias in the occurrence of SC‐conferring mutations on the male and female specificity genes. We found that male SC‐conferring mutations were indeed more likely to be fixed than were female SC‐conferring mutations in a wide range of parameters. This pattern was particularly strong when pollen availability was relatively high. Under such a condition, even if the occurrence of mutations was biased strongly towards the female specificity gene, male SC‐conferring mutations were much more often fixed. Our study demonstrates that fixation probabilities of those two types of mutation vary strongly depending on ecological and genetic conditions, although both types result in the same evolutionary consequence—the loss of SI.     

Biomechanical trade-offs bias rates of evolution in the feeding apparatus of fishes

Morphological diversification does not proceed evenly across the organism. Some body parts tend to evolve at higher rates than others, and these rate biases are often attributed to sexual and natural selection or to genetic constraints. We hypothesized that variation in the rates of morphological evolution among body parts could also be related to the performance consequences of the functional systems that make up the body. Specifically, we tested the widely held expectation that the rate of evolution for a trait is negatively correlated with the strength of biomechanical trade-offs to which it is exposed. We quantified the magnitude of trade-offs acting on the morphological components of three feeding-related functional systems in four radiations of teleost fishes. After accounting for differences in the rates of morphological evolution between radiations, we found that traits that contribute more to performance trade-offs tend to evolve more rapidly, contrary to the prediction. While ecological and genetic factors are known to have strong effects on rates of phenotypic evolution, this study highlights the role of the biomechanical architecture of functional systems in biasing the rates and direction of trait evolution.     

The effects of genotype, age, and social environment on male ornamentation, mating behavior, and attractiveness

The traits thought to advertise genetic quality are often highly susceptible to environmental variation and prone to change with age. These factors may either undermine or reinforce the potential for advertisement traits to signal quality depending on the magnitude of age-dependent expression, environmental variation, and genotype-age and genotype-environment interaction. Measurements of the magnitude of these effects are thus a necessary step toward assessing the implications of age dependence and environmental variability for the evolution of signals of quality. We conducted a longitudinal study of male guppies (Poecilia reticulata) from 22 full-sibling families. Each fish was assigned at maturity to one of three treatments in order to manipulate his allocation of resources to reproduction: a control in which the male was kept alone, a courtship-only treatment in which he could see and court a female across a clear partition, and a mating treatment in which he interacted freely with a female. We measured each male's size, ornamental color patterns, courtship, attractiveness to females, and mating success at three ages. Size was influenced by treatment and age-treatment interactions, indicating that courtship and mating may impose costs on growth. Tail size and color patterns were influenced by age but not by treatment, suggesting fixed age-dependent trajectories in these advertisement traits. By contrast, display rate and attempted sneak copulation rate differed among treatments but not among ages, suggesting greater plasticity of these behavioral traits. As a result of the different patterns of variation in ornamentation and behavior, male attractiveness and mating success responded to male age, treatment, and the interaction between age and treatment. Neither age nor treatment obscured the presence of genetic variation, and the genetic relationship between male ornamentation and attractiveness remained the same among treatments. Our findings suggest that neither age-dependent variation nor environmentally induced variation in reproductive effort is likely to undermine the reliability of male signaling.     

TESTING THE EFFECTS OF MATING SYSTEM VARIATION ON RATES OF MOLECULAR EVOLUTION IN PRIMATES

Post‐copulatory sexual selection has been proposed to drive the rapid evolution of reproductive proteins, and, more recently, to increase genome‐wide mutation rates. Comparisons of rates of molecular evolution between lineages with different levels of female multiple mating represent a promising, but under‐utilized, approach for testing the effects of sperm competition on sequence evolution. Here, I use comparisons between primate species with divergent mating systems to examine the effects of sperm competition on reproductive protein evolution, as well as on sex‐averaged mutation rates. Rates of nonsynonymous substitution are higher for testis‐specific genes along the chimpanzee lineage in comparison to the human lineage, consistent with expectations. However, the data reported here do not allow firm conclusions concerning the effects of mating system on genome‐wide mutation rates, with different results obtained from different species pairs. Ultimately, comparative studies encompassing a range of mating systems and other life history traits will be required to make broad generalizations concerning the genomic effects of sperm competition.     

Fixation biases affecting human SNPs

Under neutrality all classes of mutation have an equal probability of becoming fixed in a population. In this article, we describe our analysis of the frequency distributions of >5000 human SNPs and provide evident of biases in the process of fixation of certain classes of point mutation that are most likely to be attributable to biased gene conversion. The results indicate an increased fixation probability of mutations that result in the incorporation of a GC base pair. Furthermore, in transcribed regions this process exhibits strand asymmetry, and is biased towards preserving a G base on the coding strand. Biased gene conversion has the potential to explain both existence of isochores and the compositional asymmetry in mammalian transcribed regions.     

Molecular Clock of Neutral Mutations in a Fitness-Increasing Evolutionary Process

The molecular clock of neutral mutations, which represents linear mutation fixation over generations, is theoretically explained by genetic drift in fitness-steady evolution or hitchhiking in adaptive evolution. The present study is the first experimental demonstration for the molecular clock of neutral mutations in a fitness-increasing evolutionary process. The dynamics of genome mutation fixation in the thermal adaptive evolution of Escherichia coli were evaluated in a prolonged evolution experiment in duplicated lineages. The cells from the continuously fitness-increasing evolutionary process were subjected to genome sequencing and analyzed at both the population and single-colony levels. Although the dynamics of genome mutation fixation were complicated by the combination of the stochastic appearance of adaptive mutations and clonal interference, the mutation fixation in the population was simply linear over generations. Each genome in the population accumulated 1.6 synonymous and 3.1 non-synonymous neutral mutations, on average, by the spontaneous mutation accumulation rate, while only a single genome in the population occasionally acquired an adaptive mutation. The neutral mutations that preexisted on the single genome hitchhiked on the domination of the adaptive mutation. The successive fixation processes of the 128 mutations demonstrated that hitchhiking and not genetic drift were responsible for the coincidence of the spontaneous mutation accumulation rate in the genome with the fixation rate of neutral mutations in the population. The molecular clock of neutral mutations to the fitness-increasing evolution suggests that the numerous neutral mutations observed in molecular phylogenetic trees may not always have been fixed in fitness-steady evolution but in adaptive evolution.     

MUTATIONAL CONTRIBUTIONS TO GENETIC VARIANCE-COVARIANCE MATRICES: AN EXPERIMENTAL APPROACH USING INDUCED MUTATIONS IN ARABIDOPSIS THALIANA

Genetic potential for evolutionary change and covariational constraints are typically summarized as the genetic variance-covariance matrix G , and there is currently debate over the extent to which G remains effectively constant during the course of adaptive evolution. However, G provides only a temporally restricted view of constraints that ignores possible biases in how new mutations affect multivariate phenotypes. We used chemical mutagenesis to study the effect of mutations as summarized by the mutational covariance matrix, M , in Arabidopsis thaliana. By introducing mutations into three isogenic strains of A. thaliana, we were able to quantify M directly as the genetic variance-covariance matrix of mutagenized lines. Induced mutations generally did not alter the means of the six morphology and life-history traits we measured, but they did affect the levels of available genetic variation and the covariances among traits. However, these effects were not consistent among the three isogenic lines; that is, there were significant differences among the lines in both the number of mutations produced by ethyl-methane-sulfonate treatment and the M matrices they induced. The evolutionary implications of the dependence of M on the number of mutations, the particular genetic background, and the mutagenic sampling of loci in the genome are discussed in light of commonly applied models of multivariate evolution and the potential for the genetic architecture itself to change in ways that facilitate the coordinated evolution of complex phenotypes.