Fitness-dependent mutation rates in finite populations

Mutation rate may be condition dependent, whereby individuals in poor condition, perhaps from high mutation load, have higher mutation rates than individuals in good condition. Agrawal (J. Evol. Biol.15, 2002, 1004) explored the basic properties of fitness-dependent mutation rate (FDMR) in infinite populations and reported some heuristic results for finite populations. The key parameter governing how infinite populations evolve under FDMR is the curvature (k) of the relationship between fitness and mutation rate. We extend Agrawal's analysis to finite populations and consider dominance and epistasis. In finite populations, the probability of long-term existence depends on k. In sexual populations, positive curvature leads to low equilibrium mutation rate, whereas negative curvature results in high mutation rate. In asexual populations, negative curvature results in rapid extinction via 'mutational meltdown', whereas positive curvature sometimes allows persistence. We speculate that fitness-dependent mutation rate may provide the conditions for genetic architecture to diverge between sexual and asexual taxa.     

Genic mutation rates in mammals: local similarity,chromosomal heterogeneity,and X-versus-autosome disparity

The reduction of mutation rates on the mammalian X chromosome relative to autosomes is most often explained in the literature as evidence of male-driven evolution. This hypothesis attributes lowered mutation rates on the X chromosome to the fact that this chromosome spends less time in the germline of males than in the germline of females. In contrast to this majority view, two articles argued that the patterns of mutation rates across chromosomes are inconsistent with male-driven evolution. One article reported a 40% reduction in synonymous substitution rates (Ks) for X-linked genes relative to autosomes in the mouse-rat lineage. The authors argued that this reduction is too dramatic to be explained by male-driven evolution and concluded that selection has systematically reduced mutation rate on the X chromosome to a level optimal for this male-hemizygous chromosome. More recently, a second article found that chromosomal mutation rates in both the human-mouse and mouse-rat lineages were so heterogeneous that the X chromosome was not an outlier. Here again, the authors argued that this is at odds with male-driven evolution and suggested that selection has modulated chromosomal mutation rates to locally optimal levels, thus extending the argument of the first mentioned article to include autosomes. Here, we reexamine these conclusions using mouse-rat and human-mouse coding-region data. We find a more modest reduction of Ks on the X chromosome, but our results contradict the finding that the X chromosome is not distinct from autosomes. Multiple statistical tests show that Ks rates on the X chromosome differ systematically from the autosomes in both lineages. We conclude that the moderate reduction of mutation rate on the X chromosome of both lineages is consistent with male-driven evolution; however, the large variance in mutation rates across chromosomes suggests that mutation rates are affected by additional factors besides male-driven evolution. Investigation of mutation rates by synteny reveals that synteny blocks, rather than entire chromosomes, might represent the unit of mutation rate variation.     

Mutation rates and the evolution of germline structure

Genome sequencing studies of de novo mutations in humans have revealed surprising incongruities in our understanding of human germline mutation. In particular, the mutation rate observed in modern humans is substantially lower than that estimated from calibration against the fossil record, and the paternal age effect in mutations transmitted to offspring is much weaker than expected from our long-standing model of spermatogenesis. I consider possible explanations for these discrepancies, including evolutionary changes in life-history parameters such as generation time and the age of puberty, a possible contribution from undetected post-zygotic mutations early in embryo development, and changes in cellular mutation processes at different stages of the germline. I suggest a revised model of stem-cell state transitions during spermatogenesis, in which ‘dark’ gonial stem cells play a more active role than hitherto envisaged, with a long cycle time undetected in experimental observations. More generally, I argue that the mutation rate and its evolution depend intimately on the structure of the germline in humans and other primates.This article is part of the themed issue ‘Dating species divergences using rocks and clocks''.     

Higher intensity of purifying selection on >90% of the human genes revealed by the intrinsic replacement mutation rates

For over 3 decades, the rate of replacement mutations has been assumed to be equal to, and estimated from, the rate of "strictly" neutral sequence divergence in noncoding regions and in silent-codon positions where mutations do not alter the amino acid encoded. This assumption is fundamental to estimating the fraction of harmful protein mutations and to identifying adaptive evolution at individual codons and proteins. We show that the assumption is not justifiable because a much larger fraction of codon positions is involved in hypermutable CpG dinucleotides as compared with the introns, leading to a higher expected replacement mutation rate per site in a vast majority of the genes. Consideration of this difference reveals a higher intensity of purifying natural selection than previously inferred in human genes. We also show that a much smaller number of genes are expected to be evolving with positive selection than that predicted using sequence divergence at intron and silent positions in the human genome. These patterns indicate the need for using new approaches for estimating rates of amino acid-altering mutations in order to find positively selected genes and codons in genomes that contain hypermutable CpG's.     

Temporal patterns of fruit fly (Drosophila) evolution revealed by mutation clocks

Drosophila melanogaster has been a canonical model organism to study genetics, development, behavior, physiology, evolution, and population genetics for nearly a century. Despite this emphasis and the completion of its nuclear genome sequence, the timing of major speciation events leading to the origin of this fruit fly remain elusive because of the paucity of extensive fossil records and biogeographic data. Use of molecular clocks as an alternative has been fraught with non-clock-like accumulation of nucleotide and amino-acid substitutions. Here we present a novel methodology in which genomic mutation distances are used to overcome these limitations and to make use of all available gene sequence data for constructing a fruit fly molecular time scale. Our analysis of 2977 pairwise sequence comparisons from 176 nuclear genes reveals a long-term fruit fly mutation clock ticking at a rate of 11.1 mutations per kilobase pair per Myr. Genomic mutation clock-based timings of the landmark speciation events leading to the evolution of D. melanogaster show that it shared most recent common ancestry 5.4 MYA with D. simulans, 12.6 MYA with D. erecta+D. orena, 12.8 MYA with D. yakuba+D. teisseri, 35.6 MYA with the takahashii subgroup, 41.3 MYA with the montium subgroup, 44.2 MYA with the ananassae subgroup, 54.9 MYA with the obscura group, 62.2 MYA with the willistoni group, and 62.9 MYA with the subgenus Drosophila. These and other estimates are compatible with those known from limited biogeographic and fossil records. The inferred temporal pattern of fruit fly evolution shows correspondence with the cooling patterns of paleoclimate changes and habitat fragmentation in the Cenozoic.     

Low abundance of Escherichia coli microsatellites is associated with an extremely low mutation rate

It is widely assumed that microsatellites are generated by replication slippage, a mutation process specific to repetitive DNA. Consistent with their high mutation rate, microsatellites are highly abundant in most eukaryotic genomes. In Escherichia coli, however, microsatellites are rare and short despite the fact that a high microsatellite mutation rate was described. We show that this high microsatellite instability depends on the presence of the F-plasmid. E. coli cells lacking the F-plasmid have extremely low microsatellite mutation rates. This result provides a possible explanation for the genome-wide low density of microsatellites in E. coli. Furthermore, we show that the F-plasmid induced microsatellite instability is independent of the mismatch repair pathway.     

The evolution of bacterial mutation rates under simultaneous selection by interspecific and social parasitism

Many bacterial populations harbour substantial numbers of hypermutable bacteria, in spite of hypermutation being associated with deleterious mutations. One reason for the persistence of hypermutators is the provision of novel mutations, enabling rapid adaptation to continually changing environments, for example coevolving virulent parasites. However, hypermutation also increases the rate at which intraspecific parasites (social cheats) are generated. Interspecific and intraspecific parasitism are therefore likely to impose conflicting selection pressure on mutation rate. Here, we combine theory and experiments to investigate how simultaneous selection from inter- and intraspecific parasitism affects the evolution of bacterial mutation rates in the plant-colonizing bacterium Pseudomonas fluorescens. Both our theoretical and experimental results suggest that phage presence increases and selection for public goods cooperation (the production of iron-scavenging siderophores) decreases selection for mutator bacteria. Moreover, phages imposed a much greater growth cost than social cheating, and when both selection pressures were imposed simultaneously, selection for cooperation did not affect mutation rate evolution. Given the ubiquity of infectious phages in the natural environment and clinical infections, our results suggest that phages are likely to be more important than social interactions in determining mutation rate evolution.     

The effects of spontaneous mutation on competitive fitness in Drosophila melanogaster

Abstract We have analysed the effect of 288 generations of mutation accumulation (MA) on chromosome II competitive fitness in 21 full‐sib lines of Drosophila melanogaster and in a large control population, all derived from the same isogenic base. The rate of mean log‐fitness decline and that of increase of the between‐line variance were consistent with a low rate (λ ≈ 0.03 per gamete and generation), and moderate average fitness effect [E(s) ≈ 0.1] of deleterious mutation. Subsequently, crosses were made between pairs of MA lines, and these were maintained with effective size on the order of a few tens. In these crosses, MA recombinant chromosomes quickly recovered to about the average fitness level of control chromosomes. Thus, deleterious mutations responsible for the fitness decline were efficiently selected against in relatively small populations, confirming that their effects were larger than a few percent.     

Evaluation of structural and evolutionary contributions to deleterious mutation prediction

Methods for automated prediction of deleterious protein mutations have utilized both structural and evolutionary information but the relative contribution of these two factors remains unclear. To address this, we have used a variety of structural and evolutionary features to create simple deleterious mutation models that have been tested on both experimental mutagenesis and human allele data. We find that the most accurate predictions are obtained using a solvent-accessibility term, the C(beta) density, and a score derived from homologous sequences, SIFT. A classification tree using these two features has a cross-validated prediction error of 20.5% on an experimental mutagenesis test set when the prior probability for deleterious and neutral cases is equal, whereas this prediction error is 28.8% and 22.2% using either the C(beta) density or SIFT alone. The improvement imparted by structure increases when fewer homologs are available: when restricted to three homologs the prediction error improves from 26.9% using SIFT alone to 22.4% using SIFT and the C(beta) density, or 24.8% using SIFT and a noisy C(beta) density term approximating the inaccuracy of ab initio structures modeled by the Rosetta method. We conclude that methods for deleterious mutation prediction should include structural information when fewer than five to ten homologs are available, and that ab initio predicted structures may soon be useful in such cases when high-resolution structures are unavailable.     

A long-term evolutionary pressure on the amount of noncoding DNA

A significant part of eukaryotic noncoding DNA is viewed as the passive result of mutational processes, such as the proliferation of mobile elements. However, sequences lacking an immediate utility can nonetheless play a major role in the long-term evolvability of a lineage, for instance by promoting genomic rearrangements. They could thus be subject to an indirect selection. Yet, such a long-term effect is difficult to isolate either in vivo or in vitro. Here, by performing in silico experimental evolution, we demonstrate that, under low mutation rates, the indirect selection of variability promotes the accumulation of noncoding sequences: Even in the absence of self-replicating elements and mutational bias, noncoding sequences constituted an important fraction of the evolved genome because the indirectly selected genomes were those that were variable enough to discover beneficial mutations. On the other hand, high mutation rates lead to compact genomes, much like the viral ones, although no selective cost of genome size was applied: The indirectly selected genomes were those that were small enough for the genetic information to be reliably transmitted. Thus, the spontaneous evolution of the amount of noncoding DNA strongly depends on the mutation rate. Our results suggest the existence of an additional pressure on the amount of noncoding DNA, namely the indirect selection of an appropriate trade-off between the fidelity of the transmission of the genetic information and the exploration of the mutational neighborhood. Interestingly, this trade-off resulted robustly in the accumulation of noncoding DNA so that the best individual leaves one offspring without mutation (or only neutral ones) per generation.     

Characterizing the time dependency of human mitochondrial DNA mutation rate estimates

Previous research has established a discrepancy of nearly anorder of magnitude between pedigree-based and phylogeny-based(human vs. chimpanzee) estimates of the mitochondrial DNA (mtDNA)control region mutation rate. We characterize the time dependencyof the human mitochondrial hypervariable region one mutationrate by generating 14 new phylogeny-based mutation rate estimatesusing within-human comparisons and archaeological dates. Rateestimates based on population events between 15,000 and 50,000years ago are at least 2-fold lower than pedigree-based estimates.These within-human estimates are also higher than estimatesgenerated from phylogeny-based human–chimpanzee comparisons.Our new estimates establish a rapid decay in evolutionary mutationrate between approximately 2,500 and 50,000 years ago and aslow decay from 50,000 to 6 Ma. We then extend this analysisto the mtDNA-coding region. Our within-human coding region mutationrate estimates display a similar, though less rapid, time-dependentdecay. We explore the possibility that multiple hits explainthe discrepancy between pedigree-based and phylogeny-based mutationrates. We conclude that whereas nucleotide substitution modelsincorporating multiple hits do provide a possible explanationfor the discrepancy between pedigree-based and human–chimpanzeemutation rate estimates, they do not explain the rapid declineof within-human rate estimates. We propose that demographicprocesses such as serial bottlenecks prior to the Holocene couldexplain the difference between rates estimated before and after15,000 years ago. Our findings suggest that human mtDNA estimatesof dates of population and phylogenetic events should be adjustedin light of this time dependency of the mutation rate estimates.     

A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences

Summary Some simple formulae were obtained which enable us to estimate evolutionary distances in terms of the number of nucleotide substitutions (and, also, the evolutionary rates when the divergence times are known). In comparing a pair of nucleotide sequences, we distinguish two types of differences; if homologous sites are occupied by different nucleotide bases but both are purines or both pyrimidines, the difference is called type I (or transition type), while, if one of the two is a purine and the other is a pyrimidine, the difference is called type II (or transversion type). Letting P and Q be respectively the fractions of nucleotide sites showing type I and type II differences between two sequences compared, then the evolutionary distance per site is K = — (1/2) ln {(1 — 2P — Q) }. The evolutionary rate per year is then given by k = K/(2T), where T is the time since the divergence of the two sequences. If only the third codon positions are compared, the synonymous component of the evolutionary base substitutions per site is estimated by K'_S = — (1/2) ln (1 — 2P — Q). Also, formulae for standard errors were obtained. Some examples were worked out using reported globin sequences to show that synonymous substitutions occur at much higher rates than amino acid-altering substitutions in evolution.Contribution No. 1330 from the National Institute of Genetics, Mishima, 411 Japan     

Context-dependent mutation rates may cause spurious signatures of a fixation bias favoring higher GC-content in humans

Understanding the proximate and ultimate causes underlying the evolution of nucleotide composition in mammalian genomes is of fundamental interest to the study of molecular evolution. Comparative genomics studies have revealed that many more substitutions occur from G and C nucleotides to A and T nucleotides than the reverse, suggesting that mammalian genomes are not at equilibrium for base composition. Analysis of human polymorphism data suggests that mutations that increase GC-content tend to be at much higher frequencies than those that decrease or preserve GC-content when the ancestral allele is inferred via parsimony using the chimpanzee genome. These observations have been interpreted as evidence for a fixation bias in favor of G and C alleles due to either positive natural selection or biased gene conversion. Here, we test the robustness of this interpretation to violations of the parsimony assumption using a data set of 21,488 noncoding single nucleotide polymorphisms (SNPs) discovered by the National Institute of Environmental Health Sciences (NIEHS) SNPs project via direct resequencing of n = 95 individuals. Applying standard nonparametric and parametric population genetic approaches, we replicate the signatures of a fixation bias in favor of G and C alleles when the ancestral base is assumed to be the base found in the chimpanzee outgroup. However, upon taking into account the probability of misidentifying the ancestral state of each SNP using a context-dependent mutation model, the corrected distribution of SNP frequencies for GC-content increasing SNPs are nearly indistinguishable from the patterns observed for other types of mutations, suggesting that the signature of fixation bias is a spurious artifact of the parsimony assumption.     

Coevolution with phages does not influence the evolution of bacterial mutation rates in soil

Coevolution with phages drive the evolution of high bacterial mutation rates in vitro, but the relevance of this finding to natural populations is unclear. Here, we investigated how coevolution affects mutation rate evolution in soil, in the presence and absence of the rest of the natural microbial community. Although mutation rate on average increased threefold, neither coevolving phages nor the rest of natural community significantly affected mutation rates. Our results suggest that features of the soil over and above directly interacting organisms constrain the evolution of strong mutators, helping to explain their relatively low frequency compared with some laboratory and clinical settings.     

Direct determination of the influence of extreme temperature on transposition and structural mutation rates of Drosophila melanogaster mobile elements

Two sets of mutation accumulation lines, one reared at 28°C and the other at 24°C, were compared for their transposition and rearrangement rates of eleven transposable element families. The changes affecting mobile elements were analysed by the Southern technique and in situ hybridization. No differences were found between treated and control lines. The role of the host genotype in transposition control and the significance of structural mutations in transposable element dynamics are discussed.     

Disentangling the effects of mating systems and mutation rates on cytoplamic diversity in gynodioecious Silene nutans and dioecious Silene otites

Many flowering plant species exhibit a variety of distinct sexual morphs, the two mostcommon cases being the co-occurrence of females and males (dioecy) or the co-occurrence ofhermaphrodites and females (gynodioecy). In this study, we compared DNA sequencevariability of the three genomes (nuclear, mitochondrial and chloroplastic) of agynodioecious species, Silene nutans, with that of a closely related dioeciousspecies, Silene otites. In the light of theoretical models, we expect cytoplasmicdiversity to differ between the two species due to the selective dynamics that acts oncytoplasmic genomes in gynodioecious species: under an epidemic scenario, thegynodioecious species is expected to exhibit lower cytoplasmic diversity than thedioecious species, while the opposite is expected in the case of balancing selectionmaintaining sterility cytoplasms in the gynodioecious species. We found no differencebetween the species for nuclear gene diversity, but, for the cytoplasmic loci, thegynodioecious S. nutans had more haplotypes, and higher nucleotide diversity,than the dioecious relative, S. otites, even though the latter has a relativelyhigh rate of mitochondrial synonymous substitutions, and therefore presumably a highermutation rate. Therefore, as the mitochondrial mutation rate cannot account for the highercytoplasmic diversity found in S. nutans, our findings support the hypothesisthat gynodioecy in S. nutans has been maintained by balancing selection ratherthan by epidemic-like dynamics.