Low Genetic Quality Alters Key Dimensions of the Mutational Spectrum

Variation in mutation, recombination, and transposition rates in Drosophila melanogaster and Drosophila simulans

The rates of mutation, recombination, and transposition are core parameters in models of evolution. They impact genetic diversity, responses to ongoing selection, and levels of genetic load. However, even for key evolutionary model species such as Drosophila melanogaster and Drosophila simulans, few estimates of these parameters are available, and we have little idea of how rates vary between individuals, sexes, or populations. Knowledge of this variation is fundamental for parameterizing models of genome evolution. Here, we provide direct estimates of mutation, recombination, and transposition rates and their variation in a West African and a European population of D. melanogaster and a European population of D. simulans Across 89 flies, we observe 58 single-nucleotide mutations, 286 crossovers, and 89 transposable element (TE) insertions. Compared to the European D. melanogaster, we find the West African population has a lower mutation rate (1.67 × 10-9 site-1 gen-1 vs. 4.86 × 10-9 site-1 gen-1) and a lower transposition rate (8.99 × 10-5 copy-1 gen-1 vs. 23.36 × 10-5 copy-1 gen-1), but a higher recombination rate (3.44 cM/Mb vs. 2.06 cM/Mb). The European D. simulans population has a similar mutation rate to European D. melanogaster, but a significantly higher recombination rate and a lower, but not significantly different, transposition rate. Overall, we find paternal-derived mutations are more frequent than maternal ones in both species. Our study quantifies the variation in rates of mutation, recombination, and transposition among different populations and sexes, and our direct estimates of these parameters in D. melanogaster and D. simulans will benefit future studies in population and evolutionary genetics.

Transcriptional and mutational signatures of the Drosophila ageing germline

Ageing is a complex biological process that is accompanied by changes in gene expression and mutational load. In many species, including humans, older fathers pass on more paternally derived de novo mutations; however, the cellular basis and cell types driving this pattern are still unclear. To explore the root causes of this phenomenon, we performed single-cell RNA sequencing on testes from young and old male Drosophila and genomic sequencing (DNA sequencing) on somatic tissues from the same flies. We found that early germ cells from old and young flies enter spermatogenesis with similar mutational loads but older flies are less able to remove mutations during spermatogenesis. Mutations in old cells may also increase during spermatogenesis. Our data reveal that old and young flies have distinct mutational biases. Many classes of genes show increased postmeiotic expression in the germlines of older flies. Late spermatogenesis-biased genes have higher dN/dS (ratio of non-synonymous to synonymous substitutions) than early spermatogenesis-biased genes, supporting the hypothesis that late spermatogenesis is a source of evolutionary innovation. Surprisingly, genes biased in young germ cells show higher dN/dS than genes biased in old germ cells. Our results provide new insights into the role of the germline in de novo mutation.

Salt stress alters the spectrum of de novo mutation available to selection during experimental adaptation of Chlamydomonas reinhardtii

The genetic basis of adaptation is driven by both selection and the spectrum of available mutations. Given that the rate of mutation is not uniformly distributed across the genome and varies depending on the environment, understanding the signatures of selection across the genome is aided by first establishing what the expectations of genetic change are from mutation. To determine the interaction between salt stress, selection, and mutation across the genome, we compared mutations observed in a selection experiment for salt tolerance in Chlamydomonas reinhardtii to those observed in mutation accumulation (MA) experiments with and without salt exposure. MA lines evolved under salt stress had a single-nucleotide mutation rate of 1.1 × 10 - 9 $1.1 \times 10^{-9}$ , similar to that of MA lines under standard conditions ( 9.6 × 10 - 10 $9.6 \times 10^{-10}$ ). However, we found that salt stress led to an increased rate of indel mutations, but that many of these mutations were removed under selection. Finally, lines adapted to salt also showed excess clustering of mutations in the genome and the co-expression network, suggesting a role for positive selection in retaining mutations in particular compartments of the genome during the evolution of salt tolerance. Our study shows that characterizing mutation rates and spectra expected under stress helps disentangle the effects of environment and selection during adaptation.

Is the mutation rate lower in genomic regions of stronger selective constraints?

A study of the plant Arabidopsis thaliana detected lower mutation rates in genomic regions where mutations are more likely to be deleterious, challenging the principle that mutagenesis is blind to its consequence. To examine the generality of this finding, we analyze large mutational data from baker's yeast and humans. The yeast data do not exhibit this trend, whereas the human data show an opposite trend that disappears upon the control of potential confounders. We find that the Arabidopsis study identified substantially more mutations than reported in the original data-generating studies and expected from Arabidopsis' mutation rate. These extra mutations are enriched in polynucleotide tracts and have relatively low sequencing qualities so are likely sequencing errors. Furthermore, the polynucleotide “mutations” can produce the purported mutational trend in Arabidopsis. Together, our results do not support lower mutagenesis of genomic regions of stronger selective constraints in the plant, fungal, and animal models examined.

Intraspecific Variation in the Rates of Mutations Causing Structural Variation in Daphnia magna

Mutations that cause structural variation are important sources of genetic variation upon which other evolutionary forces can act, however, they are difficult to observe and therefore few direct estimates of their rate and spectrum are available. Understanding mutation rate evolution, however, requires adding to the limited number of species for which direct estimates are available, quantifying levels of intraspecific variation in mutation rates, and assessing whether rate estimates co-vary across types of mutation. Here, we report structural variation-causing mutation rates (svcMRs) for six categories of mutations (short insertions and deletions, long deletions and duplications, and deletions and duplications at copy number variable sites) from nine genotypes of Daphnia magna collected from three populations in Finland, Germany, and Israel using a mutation accumulation approach. Based on whole-genome sequence data and validated using simulations, we find svcMRs are high (two orders of magnitude higher than base substitution mutation rates measured in the same lineages), highly variable among populations, and uncorrelated across categories of mutation. Furthermore, to assess the impact of scvMRs on the genome, we calculated rates while adjusting for the lengths of events and ran simulations to determine if the mutations occur in genic regions more or less frequently than expected by chance. Our results pose a challenge to most prevailing theories aimed at explaining the evolution of the mutation rate, underscoring the importance of obtaining additional mutation rate estimates in more genotypes, for more types of mutation, in more species, in order to improve our future understanding of mutation rates, their variation, and their evolution.

Mutability of mononucleotide repeats, not oxidative stress, explains the discrepancy between laboratory-accumulated mutations and the natural allele-frequency spectrum in C. elegans

Important clues about natural selection can be gleaned from discrepancies between the properties of segregating genetic variants and of mutations accumulated experimentally under minimal selection, provided the mutational process is the same in the laboratory as in nature. The base-substitution spectrum differs between C. elegans laboratory mutation accumulation (MA) experiments and the standing site-frequency spectrum, which has been argued to be in part owing to increased oxidative stress in the laboratory environment. Using genome sequence data from C. elegans MA lines carrying a mutation (mev-1) that increases the cellular titer of reactive oxygen species (ROS), leading to increased oxidative stress, we find the base-substitution spectrum is similar between mev-1, its wild-type progenitor (N2), and another set of MA lines derived from a different wild strain (PB306). Conversely, the rate of short insertions is greater in mev-1, consistent with studies in other organisms in which environmental stress increased the rate of insertion–deletion mutations. Further, the mutational properties of mononucleotide repeats in all strains are different from those of nonmononucleotide sequence, both for indels and base-substitutions, and whereas the nonmononucleotide spectra are fairly similar between MA lines and wild isolates, the mononucleotide spectra are very different, with a greater frequency of A:T → T:A transversions and an increased proportion of ±1-bp indels. The discrepancy in mutational spectra between laboratory MA experiments and natural variation is likely owing to a consistent (but unknown) effect of the laboratory environment that manifests itself via different modes of mutability and/or repair at mononucleotide loci.

The spontaneous mutation rate of Drosophila pseudoobscura

The spontaneous mutation rate is a very variable trait that is subject to drift, selection and is sometimes highly plastic. Consequently, its variation between close species, or even between populations from the same species, can be very large. Here, I estimated the spontaneous mutation rate of Drosophila pseudoobscura and Drosophila persimilis crosses to explore the mutation rate variation within the Drosophila genus. All mutation rate estimations in Drosophila varied fourfold, probably explained by the sensitivity of the mutation rate to environmental and experimental conditions. Moreover, I found a very high mutation rate in the hybrid cross between D. pseudoobscura and D. persimilis, in agreement with known elevated mutation rate in hybrids. This mutation rate increase can be explained by heterozygosity and fitness decrease effects in hybrids.

Stability across the Whole Nuclear Genome in the Presence and Absence of DNA Mismatch Repair

We describe the contribution of DNA mismatch repair (MMR) to the stability of the eukaryotic nuclear genome as determined by whole-genome sequencing. To date, wild-type nuclear genome mutation rates are known for over 40 eukaryotic species, while measurements in mismatch repair-defective organisms are fewer in number and are concentrated on Saccharomyces cerevisiae and human tumors. Well-studied organisms include Drosophila melanogaster and Mus musculus, while less genetically tractable species include great apes and long-lived trees. A variety of techniques have been developed to gather mutation rates, either per generation or per cell division. Generational rates are described through whole-organism mutation accumulation experiments and through offspring–parent sequencing, or they have been identified by descent. Rates per somatic cell division have been estimated from cell line mutation accumulation experiments, from systemic variant allele frequencies, and from widely spaced samples with known cell divisions per unit of tissue growth. The latter methods are also used to estimate generational mutation rates for large organisms that lack dedicated germlines, such as trees and hyphal fungi. Mechanistic studies involving genetic manipulation of MMR genes prior to mutation rate determination are thus far confined to yeast, Arabidopsis thaliana, Caenorhabditis elegans, and one chicken cell line. A great deal of work in wild-type organisms has begun to establish a sound baseline, but far more work is needed to uncover the variety of MMR across eukaryotes. Nonetheless, the few MMR studies reported to date indicate that MMR contributes 100-fold or more to genome stability, and they have uncovered insights that would have been impossible to obtain using reporter gene assays.

High and Highly Variable Spontaneous Mutation Rates in Daphnia

The rate and spectrum of spontaneous mutations are critical parameters in basic and applied biology because they dictate the pace and character of genetic variation introduced into populations, which is a prerequisite for evolution. We use a mutation–accumulation approach to estimate mutation parameters from whole-genome sequence data from multiple genotypes from multiple populations of Daphnia magna, an ecological and evolutionary model system. We report extremely high base substitution mutation rates (µ-n,bs = 8.96 × 10⁻⁹/bp/generation [95% CI: 6.66–11.97 × 10⁻⁹/bp/generation] in the nuclear genome and µ-m,bs = 8.7 × 10⁻⁷/bp/generation [95% CI: 4.40–15.12 × 10⁻⁷/bp/generation] in the mtDNA), the highest of any eukaryote examined using this approach. Levels of intraspecific variation based on the range of estimates from the nine genotypes collected from three populations (Finland, Germany, and Israel) span 1 and 3 orders of magnitude, respectively, resulting in up to a ∼300-fold difference in rates among genomic partitions within the same lineage. In contrast, mutation spectra exhibit very consistent patterns across genotypes and populations, suggesting the mechanisms underlying the mutational process may be similar, even when the rates at which they occur differ. We discuss the implications of high levels of intraspecific variation in rates, the importance of estimating gene conversion rates using a mutation–accumulation approach, and the interacting factors influencing the evolution of mutation parameters. Our findings deepen our knowledge about mutation and provide both challenges to and support for current theories aimed at explaining the evolution of the mutation rate, as a trait, across taxa.

Measuring mutagenicity in ecotoxicology: A case study of Cd exposure in Chironomus riparius

Existing mutagenicity tests for metazoans lack the direct observation of enhanced germline mutation rates after exposure to anthropogenic substances, therefore being inefficient. Cadmium (Cd) is a metal described as a mutagen in mammalian cells and listed as a group 1 carcinogenic and mutagenic substance. But Cd mutagenesis mechanism is not yet clear. Therefore, in the present study, we propose a method coupling short-term mutation accumulation (MA) lines with subsequent whole genome sequencing (WGS) and a dedicated data analysis pipeline to investigate if chronic Cd exposure on Chironomus riparius can alter the rate at which de novo point mutations appear. Results show that Cd exposure did not affect the basal germline mutation rate nor the mutational spectrum in C. riparius, thereby arguing that exposed organisms might experience a range of other toxic effects before any mutagenic effect may occur. We show that it is possible to establish a practical and easily implemented pipeline to rapidly detect germ cell mutagens in a metazoan test organism. Furthermore, our data implicate that it is questionable to transfer mutagenicity assessments based on in vitro methods to complex metazoans.

Approximating the Coalescent Under Facultative Sex

Genome studies of facultative sexual species, which can either reproduce sexually or asexually, are providing insight into the evolutionary consequences of mixed reproductive modes. It is currently unclear to what extent the evolutionary history of facultative sexuals' genomes can be approximated by the standard coalescent, and if a coalescent effective population size Ne exists. Here, I determine if and when these approximations can be made. When sex is frequent (occurring at a frequency much greater than 1/N per reproduction per generation, for N the actual population size), the underlying genealogy can be approximated by the standard coalescent, with a coalescent Ne≈N. When sex is very rare (at frequency much lower than 1/N), approximations for the pairwise coalescent time can be obtained, which is strongly influenced by the frequencies of sex and mitotic gene conversion, rather than N. However, these terms do not translate into a coalescent Ne. These results are used to discuss the best sampling strategies for investigating the evolutionary history of facultative sexual species.

The population genetics of ploidy change in unicellular fungi

Changes in ploidy are a significant type of genetic variation, describing the number of chromosome sets per cell. Ploidy evolves in natural populations, clinical populations, and lab experiments, particularly in fungi. Despite a long history of theoretical work on this topic, predicting how ploidy will evolve has proven difficult, as it is often unclear why one ploidy state outperforms another. Here, we review what is known about contemporary ploidy evolution in diverse fungal species through the lens of population genetics. As with typical genetic variants, ploidy evolution depends on the rate that new ploidy states arise by mutation, natural selection on alternative ploidy states, and random genetic drift. However, ploidy variation also has unique impacts on evolution, with the potential to alter chromosomal stability, the rate and patterns of point mutation, and the nature of selection on all loci in the genome. We discuss how ploidy evolution depends on these general and unique factors and highlight areas where additional experimental evidence is required to comprehensively explain the ploidy transitions observed in the field and the lab.

Estimation of the SNP Mutation Rate in Two Vegetatively Propagating Species of Duckweed

Mutation rate estimates for vegetatively reproducing organisms are rare, despite their frequent occurrence across the tree of life. Here we report mutation rate estimates in two vegetatively reproducing duckweed species, Lemna minor and Spirodela polyrhiza We use a modified approach to estimating mutation rates by taking into account the reduction in mutation detection power that occurs when new individuals are produced from multiple cell lineages. We estimate an extremely low per generation mutation rate in both species of duckweed and note that allelic coverage at de novo mutation sites is very skewed. We also find no substantial difference in mutation rate between mutation accumulation lines propagated under benign conditions and those grown under salt stress. Finally, we discuss the implications of interpreting mutation rate estimates in vegetatively propagating organisms.

Intergenerational Transfer of Ageing: Parental Age and Offspring Lifespan

The extent to which the age of parents at reproduction can affect offspring lifespan and other fitness-related traits is important in our understanding of the selective forces shaping life history evolution. In this article, the widely reported negative effects of parental age on offspring lifespan (the 'Lansing effect') is examined. Outlined herein are the potential routes whereby a Lansing effect can occur, whether effects might accumulate across multiple generations, and how the Lansing effect should be viewed as part of a broader framework, considering how parental age affects offspring fitness. The robustness of the evidence for a Lansing effect produced so far, potential confounding variables, and how the underlying mechanisms might best be unravelled through carefully designed experimental studies are discussed.

Pervasive Strong Selection at the Level of Codon Usage Bias in Drosophila melanogaster

Codon usage bias (CUB), where certain codons are used more frequently than expected by chance, is a ubiquitous phenomenon and occurs across the tree of life. The dominant paradigm is that the proportion of preferred codons is set by weak selection. While experimental changes in codon usage have at times shown large phenotypic effects in contrast to this paradigm, genome-wide population genetic estimates have supported the weak selection model. Here we use deep genomic population sequencing of two Drosophila melanogaster populations to measure selection on synonymous sites in a way that allowed us to estimate the prevalence of both weak and strong purifying selection. We find that selection in favor of preferred codons ranges from weak (|Nes| ∼ 1) to strong (|Nes| > 10), with strong selection acting on 10-20% of synonymous sites in preferred codons. While previous studies indicated that selection at synonymous sites could be strong, this is the first study to detect and quantify strong selection specifically at the level of CUB. Further, we find that CUB-associated polymorphism accounts for the majority of strong selection on synonymous sites, with secondary contributions of splicing (selection on alternatively spliced genes, splice junctions, and spliceosome-bound sites) and transcription factor binding. Our findings support a new model of CUB and indicate that the functional importance of CUB, as well as synonymous sites in general, have been underestimated.

Exploring background mutational processes to decipher cancer genetic heterogeneity

Much remains unknown about the progression and heterogeneity of mutational processes in different cancers and their diagnostic and clinical potential. A growing body of evidence supports mutation rate dependence on the local DNA sequence context for various types of mutations. We propose several tools for the analysis of cancer context-dependent mutations, which are implemented in an online computational framework MutaGene. The framework explores DNA context-dependent mutational patterns and underlying somatic cancer mutagenesis, analyzes mutational profiles of cancer samples, identifies the combinations of underlying mutagenic processes including those related to infidelity of DNA replication and repair machinery, and various other endogenous and exogenous mutagenic factors. As a result, the combination of mutagenic processes can be identified in any query sample with subsequent comparison to mutational profiles derived from malignant and benign samples. In addition, mutagen or cancer-specific mutational background models are applied to calculate expected DNA and protein site mutability to decouple relative contributions of mutagenesis and selection in carcinogenesis, thus elucidating the site-specific driving events in cancer. MutaGene is freely available at https://www.ncbi.nlm.nih.gov/projects/mutagene/.

Mutational Landscape of Spontaneous Base Substitutions and Small Indels in Experimental Caenorhabditis elegans Populations of Differing Size

Experimental investigations into the rates and fitness effects of spontaneous mutations are fundamental to our understanding of the evolutionary process. To gain insights into the molecular and fitness consequences of spontaneous mutations, we conducted a mutation accumulation (MA) experiment at varying population sizes in the nematode Caenorhabditis elegans, evolving 35 lines in parallel for 409 generations at three population sizes (N = 1, 10, and 100 individuals). Here, we focus on nuclear SNPs and small insertion/deletions (indels) under minimal influence of selection, as well as their accrual rates in larger populations under greater selection efficacy. The spontaneous rates of base substitutions and small indels are 1.84 (95% C.I. ± 0.14) × 10-9 substitutions and 6.84 (95% C.I. ± 0.97) × 10-10 changes/site/generation, respectively. Small indels exhibit a deletion bias with deletions exceeding insertions by threefold. Notably, there was no correlation between the frequency of base substitutions, nonsynonymous substitutions, or small indels with population size. These results contrast with our previous analysis of mitochondrial DNA mutations and nuclear copy-number changes in these MA lines, and suggest that nuclear base substitutions and small indels are under less stringent purifying selection compared to the former mutational classes. A transition bias was observed in exons as was a near universal base substitution bias toward A/T. Strongly context-dependent base substitutions, where 5'-Ts and 3'-As increase the frequency of A/T → T/A transversions, especially at the boundaries of A or T homopolymeric runs, manifest as higher mutation rates in (i) introns and intergenic regions relative to exons, (ii) chromosomal cores vs. arms and tips, and (iii) germline-expressed genes.

Evolution of the Mutational Process under Relaxed Selection in Caenorhabditis elegans

The mutational process varies at many levels, from within genomes to among taxa. Many mechanisms have been linked to variation in mutation, but understanding of the evolution of the mutational process is rudimentary. Physiological condition is often implicated as a source of variation in microbial mutation rate and may contribute to mutation rate variation in multicellular organisms.Deleterious mutations are an ubiquitous source of variation in condition. We test the hypothesis that the mutational process depends on the underlying mutation load in two groups of Caenorhabditis elegans mutation accumulation (MA) lines that differ in their starting mutation loads. "First-order MA" (O1MA) lines maintained under minimal selection for ∼250 generations were divided into high-fitness and low-fitness groups and sets of "second-order MA" (O2MA) lines derived from each O1MA line were maintained for ∼150 additional generations. Genomes of 48 O2MA lines and their progenitors were sequenced. There is significant variation among O2MA lines in base-substitution rate (µbs), but no effect of initial fitness; the indel rate is greater in high-fitness O2MA lines. Overall, µbs is positively correlated with recombination and proximity to short tandem repeats and negatively correlated with 10 bp and 1 kb GC content. However, probability of mutation is sufficiently predicted by the three-nucleotide motif alone. Approximately 90% of the variance in standing nucleotide variation is explained by mutability. Total mutation rate increased in the O2MA lines, as predicted by the "drift barrier" model of mutation rate evolution. These data, combined with experimental estimates of fitness, suggest that epistasis is synergistic.

Yeast Spontaneous Mutation Rate and Spectrum Vary with Environment

Mutation is the ultimate genetic source of evolution and biodiversity, but to what extent the environment impacts mutation rate and spectrum is poorly understood. Past studies discovered mutagenesis induced by antibiotic treatment or starvation, but its relevance and importance to long-term evolution is unclear because these severe stressors typically halt cell growth and/or cause substantial cell deaths. Here, we quantify the mutation rate and spectrum in Saccharomyces cerevisiae by whole-genome sequencing following mutation accumulation in each of seven environments with relatively rapid cell growths and minimal cell deaths. We find the point mutation rate per generation to differ by 3.6-fold among the seven environments, generally increasing in environments with slower cell growths. This trend renders the mutation rate per year more constant than that per generation across environments, which has implications for neutral evolution and the molecular clock. Additionally, we find substantial among-environment variations in mutation spectrum, such as the transition to transversion ratio and AT mutational bias. Other main mutation types, including small insertion or deletion, segmental duplication or deletion, and chromosome gain or loss also tend to occur more frequently in environments where yeast grows more slowly. In contrast to these findings from the nuclear genome, the yeast mitochondrial mutation rate rises with the growth rate, consistent with the metabolic rate hypothesis. Together, these observations indicate that environmental changes, which are ubiquitous in nature, influence not only natural selection, but also the amount and type of mutations available to selection, and suggest that ignoring the latter impact, as is currently practiced, may mislead evolutionary inferences.

Fitness and Genomic Consequences of Chronic Exposure to Low Levels of Copper and Nickel in Daphnia pulex Mutation Accumulation Lines

In at least some unicellular organisms, mutation rates are temporarily raised upon exposure to environmental stress, potentially contributing to the evolutionary response to stress. Whether this is true for multicellular organisms, however, has received little attention. This study investigated the effects of chronic mild stress, in the form of low-level copper and nickel exposure, on mutational processes in Daphnia pulex using a combination of mutation accumulation, whole genome sequencing and life-history assays. After over 100 generations of mutation accumulation, we found no effects of metal exposure on the rates of single nucleotide mutations and of loss of heterozygosity events, the two mutation classes that occurred in sufficient numbers to allow statistical analysis. Similarly, rates of decline in fitness, as measured by intrinsic rate of population increase and of body size at first reproduction, were negligibly affected by metal exposure. We can reject the possibility that Daphnia were insufficiently stressed to invoke genetic responses as we have previously shown rates of large-scale deletions and duplications are elevated under metal exposure in this experiment. Overall, the mutation accumulation lines did not significantly depart from initial values for phenotypic traits measured, indicating the lineage used was broadly mutationally robust. Taken together, these results indicate that the mutagenic effects of chronic low-level exposure to these metals are restricted to certain mutation classes and that fitness consequences are likely minor and therefore unlikely to be relevant in determining the evolutionary responses of populations exposed to these stressors.

Accelerated rates of large-scale mutations in the presence of copper and nickel

Mutation rate variation has been under intense investigation for decades. Despite these efforts, little is known about the extent to which environmental stressors accelerate mutation rates and influence the genetic load of populations. Moreover, most studies on stressors have focused on unicellular organisms and point mutations rather than large-scale deletions and duplications (copy number variations [CNVs]). We estimated mutation rates in Daphnia pulex exposed to low levels of environmental stressors as well as the effect of selection on de novo mutations. We conducted a mutation accumulation (MA) experiment in which selection was minimized, coupled with an experiment in which a population was propagated under competitive conditions in a benign environment. After an average of 103 generations of MA propagation, we sequenced 60 genomes and found significantly accelerated rates of deletions and duplications in MA lines exposed to ecologically relevant concentrations of metals. Whereas control lines had gene deletion and duplication rates comparable to other multicellular eukaryotes (1.8 × 10-6 per gene per generation), the presence of nickel and copper increased these rates fourfold. The realized mutation rate under selection was reduced to 0.4× that of control MA lines, providing evidence that CNVs contribute to mutational load. Our CNV breakpoint analysis revealed that nonhomologous recombination associated with regions of DNA fragility is the primary source of CNVs, plausibly linking metal-induced DNA strand breaks with higher CNV rates. Our findings suggest that environmental stress, in particular multiple stressors, can have profound effects on large-scale mutation rates and mutational load of multicellular organisms.

Old Trade, New Tricks: Insights into the Spontaneous Mutation Process from the Partnering of Classical Mutation Accumulation Experiments with High-Throughput Genomic Approaches

Mutations spawn genetic variation which, in turn, fuels evolution. Hence, experimental investigations into the rate and fitness effects of spontaneous mutations are central to the study of evolution. Mutation accumulation (MA) experiments have served as a cornerstone for furthering our understanding of spontaneous mutations for four decades. In the pregenomic era, phenotypic measurements of fitness-related traits in MA lines were used to indirectly estimate key mutational parameters, such as the genomic mutation rate, new mutational variance per generation, and the average fitness effect of mutations. Rapidly emerging next-generating sequencing technology has supplanted this phenotype-dependent approach, enabling direct empirical estimates of the mutation rate and a more nuanced understanding of the relative contributions of different classes of mutations to the standing genetic variation. Whole-genome sequencing of MA lines bears immense potential to provide a unified account of the evolutionary process at multiple levels-the genetic basis of variation, and the evolutionary dynamics of mutations under the forces of selection and drift. In this review, we have attempted to synthesize key insights into the spontaneous mutation process that are rapidly emerging from the partnering of classical MA experiments with high-throughput sequencing, with particular emphasis on the spontaneous rates and molecular properties of different mutational classes in nuclear and mitochondrial genomes of diverse taxa, the contribution of mutations to the evolution of gene expression, and the rate and stability of transgenerational epigenetic modifications. Future advances in sequencing technologies will enable greater species representation to further refine our understanding of mutational parameters and their functional consequences.

An experimental test of the mutation-selection balance model for the maintenance of genetic variance in fitness components

Despite decades of research, the factors that maintain genetic variation for fitness are poorly understood. It is unclear what fraction of the variance in a typical fitness component can be explained by mutation-selection balance (MSB) and whether fitness components differ in this respect. In theory, the level of standing variance in fitness due to MSB can be predicted using the rate of fitness decline under mutation accumulation, and this prediction can be directly compared to the standing variance observed. This approach allows for controlled statistical tests of the sufficiency of the MSB model, and could be used to identify traits or populations where genetic variance is maintained by other factors. For example, some traits may be influenced by sexually antagonistic balancing selection, resulting in an excess of standing variance beyond that generated by deleterious mutations. We describe the underlying theory and use it to test the MSB model for three traits in Drosophila melanogaster We find evidence for differences among traits, with MSB being sufficient to explain genetic variance in larval viability but not male mating success or female fecundity. Our results are consistent with balancing selection on sexual fitness components, and demonstrate the feasibility of rigorous statistical tests of the MSB model.

Rates and Patterns of Mutation in Tandem Repetitive DNA in Six Independent Lineages of Chlamydomonas reinhardtii

The mutational patterns of large tandem arrays of short sequence repeats remain largely unknown, despite observations of their high levels of variation in sequence and genomic abundance within and between species. Many factors can influence the dynamics of tandem repeat evolution; however, their evolution has only been examined over a limited phylogenetic sample of taxa. Here, we use publicly available whole-genome sequencing data of 85 haploid mutation accumulation lines derived from six geographically diverse Chlamydomonas reinhardtii isolates to investigate genome-wide mutation rates and patterns in tandem repeats in this species. We find that tandem repeat composition differs among ancestral strains, both in genome-wide abundance and presence/absence of individual repeats. Estimated mutation rates (repeat copy number expansion and contraction) were high, averaging 4.3×10⁻⁴ per generation per single unit copy. Although orders of magnitude higher than other types of mutation previously reported in C. reinhardtii, these tandem repeat mutation rates were one order of magnitude lower than what has recently been found in Daphnia pulex, even after correcting for lower overall genome-wide satellite abundance in C. reinhardtii. Most high-abundance repeats were related to others by a single mutational step. Correlations of repeat copy number changes within genomes revealed clusters of closely related repeats that were strongly correlated positively or negatively, and similar patterns of correlation arose independently in two different mutation accumulation experiments. Together, these results paint a dynamic picture of tandem repeat evolution in this unicellular alga.

The genome-wide rate and spectrum of spontaneous mutations differ between haploid and diploid yeast

By altering the dynamics of DNA replication and repair, alternative ploidy states may experience different rates and types of new mutations, leading to divergent evolutionary outcomes. We report a direct comparison of the genome-wide spectrum of spontaneous mutations arising in haploids and diploids following a mutation-accumulation experiment in the budding yeast Saccharomyces cerevisiae Characterizing the number, types, locations, and effects of thousands of mutations revealed that haploids were more prone to single-nucleotide mutations (SNMs) and mitochondrial mutations, while larger structural changes were more common in diploids. Mutations were more likely to be detrimental in diploids, even after accounting for the large impact of structural changes, contrary to the prediction that mutations would have weaker effects, due to masking, in diploids. Haploidy is expected to reduce the opportunity for conservative DNA repair involving homologous chromosomes, increasing the insertion-deletion rate, but we found little support for this idea. Instead, haploids were more susceptible to SNMs in late-replicating genomic regions, resulting in a ploidy difference in the spectrum of substitutions. In diploids, we detect mutation rate variation among chromosomes in association with centromere location, a finding that is supported by published polymorphism data. Diploids are not simply doubled haploids; instead, our results predict that the spectrum of spontaneous mutations will substantially shape the dynamics of genome evolution in haploid and diploid populations.

Mutation accumulation in populations of varying size: large effect mutations cause most mutational decline in the rotifer Brachionus calyciflorus under UV-C radiation

Theory predicts that fitness decline via mutation accumulation will depend on population size, but there are only a few direct tests of this key idea. To gain a qualitative understanding of the fitness effect of new mutations, we performed a mutation accumulation experiment with the facultative sexual rotifer Brachionus calyciflorus at six different population sizes under UV-C radiation. Lifetime reproduction assays conducted after ten and sixteen UV-C radiations showed that while small populations lost fitness, fitness losses diminished rapidly with increasing population size. Populations kept as low as 10 individuals were able to maintain fitness close to the nonmutagenized populations throughout the experiment indicating that selection was able to remove the majority of large effect mutations in small populations. Although our results also seem to imply that small populations are effectively immune to mutational decay, we caution against this interpretation. Given sufficient time, populations of moderate to large size can experience declines in fitness from accumulating weakly deleterious mutations as demonstrated by fitness estimates from simulations and, tentatively, from a long-term experiment with populations of moderate size. There is mounting evidence to suggest that mutational distributions contain a heavier tail of large effects. Our results suggest that this is also true when the mutational spectrum is altered by UV radiation.

Deep sequencing of natural and experimental populations of Drosophila melanogaster reveals biases in the spectrum of new mutations

Mutations provide the raw material of evolution, and thus our ability to study evolution depends fundamentally on having precise measurements of mutational rates and patterns. We generate a data set for this purpose using (1) de novo mutations from mutation accumulation experiments and (2) extremely rare polymorphisms from natural populations. The first, mutation accumulation (MA) lines are the product of maintaining flies in tiny populations for many generations, therefore rendering natural selection ineffective and allowing new mutations to accrue in the genome. The second, rare genetic variation from natural populations allows the study of mutation because extremely rare polymorphisms are relatively unaffected by the filter of natural selection. We use both methods in Drosophila melanogaster, first generating our own novel data set of sequenced MA lines and performing a meta-analysis of all published MA mutations (∼2000 events) and then identifying a high quality set of ∼70,000 extremely rare (≤0.1%) polymorphisms that are fully validated with resequencing. We use these data sets to precisely measure mutational rates and patterns. Highlights of our results include: a high rate of multinucleotide mutation events at both short (∼5 bp) and long (∼1 kb) genomic distances, showing that mutation drives GC content lower in already GC-poor regions, and using our precise context-dependent mutation rates to predict long-term evolutionary patterns at synonymous sites. We also show that de novo mutations from independent MA experiments display similar patterns of single nucleotide mutation and well match the patterns of mutation found in natural populations.

Direct Determination of the Mutation Rate in the Bumblebee Reveals Evidence for Weak Recombination-Associated Mutation and an Approximate Rate Constancy in Insects

Accurate knowledge of the mutation rate provides a base line for inferring expected rates of evolution, for testing evolutionary hypotheses and for estimation of key parameters. Advances in sequencing technology now permit direct estimates of the mutation rate from sequencing of close relatives. Within insects there have been three prior such estimates, two in nonsocial insects (Drosophila: 2.8 × 10-9 per bp per haploid genome per generation; Heliconius: 2.9 × 10-9) and one in a social species, the honeybee (3.4 × 10-9). Might the honeybee's rate be ∼20% higher because it has an exceptionally high recombination rate and recombination may be directly or indirectly mutagenic? To address this possibility, we provide a direct estimate of the mutation rate in the bumblebee (Bombus terrestris), this being a close relative of the honeybee but with a much lower recombination rate. We confirm that the crossover rate of the bumblebee is indeed much lower than honeybees (8.7 cM/Mb vs. 37 cM/Mb). Importantly, we find no significant difference in the mutation rates: we estimate for bumblebees a rate of 3.6 × 10-9 per haploid genome per generation (95% confidence intervals 2.38 × 10-9 and 5.37 × 10-9) which is just 5% higher than the estimate that of honeybees. Both genomes have approximately one new mutation per haploid genome per generation. While we find evidence for a direct coupling between recombination and mutation (also seen in honeybees), the effect is so weak as to leave almost no footprint on any between-species differences. The similarity in mutation rates suggests an approximate constancy of the mutation rate in insects.

Hsp90 and Physiological Stress Are Linked to Autonomous Transposon Mobility and Heritable Genetic Change in Nematodes

Transposable elements (TEs) have been recognized as potentially powerful drivers of genomic evolutionary change, but factors affecting their mobility and regulation remain poorly understood. Chaperones such as Hsp90 buffer environmental perturbations by regulating protein conformation, but are also part of the PIWI-interacting RNA pathway, which regulates genomic instability arising from mobile TEs in the germline. Stress-induced mutagenesis from TE movement could thus arise from functional trade-offs in the dual roles of Hsp90. We examined the functional constraints of Hsp90 and its role as a regulator of TE mobility by exposing nematodes (Caenorhabditis elegans and Caenorhabditis briggsae) to environmental stress, with and without RNAi-induced silencing of Hsp90. TE excision frequency increased with environmental stress intensity at multiple loci in several strains of each species. These effects were compounded by RNAi-induced knockdown of Hsp90. Mutation frequencies at the unc-22 marker gene in the offspring of animals exposed to environmental stress and Hsp90 RNAi mirrored excision frequency in response to these treatments. Our results support a role for Hsp90 in the suppression of TE mobility, and demonstrate that that the Hsp90 regulatory pathway can be overwhelmed with moderate environmental stress. By compromising genomic stability in germline cells, environmentally induced mutations arising from TE mobility and insertion can have permanent and heritable effects on both the phenotype and the genotype of subsequent generations.