Stochastic offspring distributions amplify selection bias in mutation accumulation experiments

Mutation accumulation (MA) experiments play an important role in understanding evolution. For microbial populations, such experiments often involve periods of population growth, such that a single individual can make a visible colony, followed by severe bottlenecks. Previous work has quantified the effect of positive and negative selection on MA experiments, demonstrating for example that with 20 generations of growth between bottlenecks, big-benefit mutations can be over-represented by a factor of five or more (Wahl and Agashe, 2022). This previous work assumed a deterministic model for population growth. We now develop a fully stochastic model, including realistic offspring distributions that incorporate genetic drift and allow for the loss of rare lineages. We demonstrate that when stochastic offspring distributions are considered, selection bias is even stronger than previously predicted. We describe several analytical and numerical methods that offer an accurate correction for the effects of selection on the observed distribution of fitness effects, describe the practical considerations in implementing each method, and demonstrate the use of this correction on simulated MA data.

The paradox of predictability provides a bridge between micro- and macroevolution

The relationship between the evolutionary dynamics observed in contemporary populations (microevolution) and evolution on timescales of millions of years (macroevolution) has been a topic of considerable debate. Historically, this debate centers on inconsistencies between microevolutionary processes and macroevolutionary patterns. Here, we characterize a striking exception: emerging evidence indicates that standing variation in contemporary populations and macroevolutionary rates of phenotypic divergence is often positively correlated. This apparent consistency between micro- and macroevolution is paradoxical because it contradicts our previous understanding of phenotypic evolution and is so far unexplained. Here, we explore the prospects for bridging evolutionary timescales through an examination of this "paradox of predictability." We begin by explaining why the divergence-variance correlation is a paradox, followed by data analysis to show that the correlation is a general phenomenon across a broad range of temporal scales, from a few generations to tens of millions of years. Then we review complementary approaches from quantitative genetics, comparative morphology, evo-devo, and paleontology to argue that they can help to address the paradox from the shared vantage point of recent work on evolvability. In conclusion, we recommend a methodological orientation that combines different kinds of short-term and long-term data using multiple analytical frameworks in an interdisciplinary research program. Such a program will increase our general understanding of how evolution works within and across timescales.

Reproductive and germ-cell mutagenic effects of poly-and perfluoroalkyl substances (PFAS) to Caenorhabditis elegans after multigenerational exposure

Per- and polyfluoroalkyl substances (PFAS) are a class of globally ubiquitous persistent organic pollutants (POPs). The developmental and reproductive toxicity of PFAS have attracted considerable attention. However, the influence of PFAS exposure on genomic stability of germ cells remains unexplored. In this study, we evaluated long-term reproductive toxicity of environmentally relevant levels of four long-chain PFAS compounds: perfluorooctanoic acid (PFOA, C8), perfluorononanoic acid (PFNA, C9), perfluorodecanoic acid (PFDA, C10), and perfluorooctanesulfonic acid (PFOS, C8), and examined their germ-cell mutagenicity in Caenorhabditis elegans. Our findings reveal that multigenerational exposure to PFAS exhibited minor impacts on development and reproduction of worms. Among all tested PFAS, PFNA significantly increased mutation frequencies of progeny by preferentially inducing T:A → C:G substitutions and small indels within repetitive regions. Further analysis of mutation spectra uncovered elevated frequencies of microhomology-mediated deletions and large deletions in PFOA-treated worms, indicating its potential activity in eliciting DNA double-strand breaks. This study provides the first comparative analysis of the genome-wide mutational profile of PFAS compounds, underscoring the importance of assessing germ-cell mutagenic actions of long-chain PFAS.

What can we infer about mutation calling by using time-series mutation accumulation data and a Bayesian Mutation Finder?

Accurate estimates of mutation rates derived from genome-wide mutation accumulation (MA) data are fundamental to understanding basic evolutionary processes. The rapidly improving high-throughput sequencing technologies provide unprecedented opportunities to identify single nucleotide mutations across genomes. However, such MA derived data are often difficult to analyze and the performance of the available methods of analysis is not well understood. In this study, we used the existing Bayesian Genotype Caller adapted for MA data that we refer to as Bayesian Mutation Finder (BMF) for identifying single nucleotide mutations while considering the characteristics of the data. We compared the performance of BMF with the widely used Genome Analysis Toolkit (GATK) by applying these two methods to time-series MA data as well as simulated data. The time-series data were obtained by propagating Daphnia pulex over an average of 188 generations and performing whole-genome sequencing of 14 MA lines across three time points. The results indicate that BMF enables more accurate identification of single nucleotide mutations than GATK especially when applied to the empirical data. Furthermore, BMF involves the use of fewer parameters and is more computationally efficient than GATK. Both BMF and GATK found surprisingly many candidate mutations that were not confirmed at later time points. We systematically infer causes of the unconfirmed candidate mutations, introduce a framework for estimating mutation rates based on genome-wide candidate mutations confirmed by subsequent sequencing, and provide an improved mutation rate estimate for D. pulex.

A simple, cost-effective, and efficient method for screening CRISPR/Cas9 mutants in plants

The clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9)-mediated genome editing system is widely used for targeted mutagenesis in a growing number of plant species. To streamline the screening process for mutants, especially those generated from low-efficiency editing events, there is a need for a rapid, cost-effective, and efficient method. Although several screening methods have been developed to process initial samples, these methods often tend to be time-consuming, expensive, or inefficient when dealing with larger sample sizes. Here we describe a simple, rapid, low-cost, and sensitive screening method for screening CRISPR/Cas9 mutants called PCR-Bsl I-associated analysis (PCR-BAA). This method requires only standard PCR and Bsl I restriction enzyme digestion, as well as agarose gel electrophoresis analysis. This method is particularly well suited for the efficient screening of mutants from larger populations of transformants. The simplicity, low cost, and high sensitivity of the PCR-BAA method make it particularly suitable for rapid screening of CRISPR/Cas9-induced mutants, especially those from low-efficiency editing events.

Reconciling theories of dominance with the relative rates of adaptive substitution on sex chromosomes and autosomes

The dominance of beneficial mutations is a key evolutionary parameter affecting the rate and genetic basis of adaptation, yet it is notoriously difficult to estimate. A leading method to infer it is to compare the relative rates of adaptive substitution for X-linked and autosomal genes, which-according to a classic model by Charlesworth et al. (1987)-is a simple function of the dominance of new beneficial mutations. Recent evidence that rates of adaptive substitution are faster for X-linked genes implies, accordingly, that beneficial mutations are usually recessive. However, this conclusion is incompatible with leading theories of dominance, which predict that beneficial mutations tend to be dominant or overdominant with respect to fitness. To address this incompatibility, we use Fisher's geometric model to predict the distribution of fitness effects of new mutations and the relative rates of positively selected substitution on the X and autosomes. Previous predictions of faster-X theory emerge as a special case of our model in which the phenotypic effects of mutations are small relative to the distance to the phenotypic optimum. But as mutational effects become large relative to the optimum, we observe an elevated tempo of positively selected substitutions on the X relative to the autosomes across a broader range of dominance conditions, including those predicted by theories of dominance. Our results imply that, contrary to previous models, dominant and overdominant beneficial mutations can plausibly generate patterns of faster-X adaptation. We discuss resulting implications for genomic studies of adaptation and inferences of dominance.

Direct inference of the distribution of fitness effects of spontaneous mutations from recombinant inbred Caenorhabditis elegans mutation accumulation lines

The distribution of fitness effects of new mutations plays a central role in evolutionary biology. Estimates of the distribution of fitness effect from experimental mutation accumulation lines are compromised by the complete linkage disequilibrium between mutations in different lines. To reduce the linkage disequilibrium, we constructed 2 sets of recombinant inbred lines from a cross of 2 Caenorhabditis elegans mutation accumulation lines. One set of lines ("RIAILs") was intercrossed for 10 generations prior to 10 generations of selfing; the second set of lines ("RILs") omitted the intercrossing. Residual linkage disequilibrium in the RIAILs is much less than in the RILs, which affects the inferred distribution of fitness effect when the sets of lines are analyzed separately. The best-fit model estimated from all lines (RIAILs + RILs) infers a large fraction of mutations with positive effects (∼40%); models that constrain mutations to have negative effects fit much worse. The conclusion is the same using only the RILs. For the RIAILs, however, models that constrain mutations to have negative effects fit nearly as well as models that allow positive effects. When mutations in high linkage disequilibrium are pooled into haplotypes, the inferred distribution of fitness effect becomes increasingly negative-skewed and leptokurtic. We conclude that the conventional wisdom-most mutations have effects near 0, a handful of mutations have effects that are substantially negative, and mutations with positive effects are very rare-is likely correct, and that unless it can be shown otherwise, estimates of the distribution of fitness effect that infer a substantial fraction of mutations with positive effects are likely confounded by linkage disequilibrium.

Effects of spontaneous mutations on survival and reproduction of Drosophila serrata infected with Drosophila C virus

The impact of selection on host immune function genes has been widely documented. However, it remains essentially unknown how mutation influences the quantitative immune traits that selection acts on. Applying a classical mutation accumulation (MA) experimental design in Drosophila serrata, we found the mutational variation in susceptibility (median time of death, LT50) to Drosophila C virus (DCV) was of similar magnitude to that reported for intrinsic survival traits. Mean LT50 did not change as mutations accumulated, suggesting no directional bias in mutational effects. Maintenance of genetic variance in immune function is hypothesized to be influenced by pleiotropic effects on immunity and other traits that contribute to fitness. To investigate this, we assayed female reproductive output for a subset of MA lines with relatively long or short survival times under DCV infection. Longer survival time tended to be associated with lower reproductive output, suggesting that mutations affecting susceptibility to DCV had pleiotropic effects on investment in reproductive fitness. Further studies are needed to uncover the general patterns of mutational effect on immune responses and other fitness traits, and to determine how selection might typically act on new mutations via their direct and pleiotropic effects.

Accumulating waves of random mutations before fixation

Mutations provide variation for evolution to emerge. A quantitative analysis of how mutations arising in single individuals expand and possibly fixate in a population is essential for studying evolutionary processes. While it is intuitive to expect that a continuous influx of mutations will lead to a continuous flow of mutations fixating in a stable constant population, joint fixation of multiple mutations occur frequently in stochastic simulations even under neutral selection. We quantitatively measure and analyze the distribution of joint fixation events of neutral mutations in constant populations and discussed the connection with previous results. We propose a new concept, the mutation "waves," where multiple mutations reach given frequencies simultaneously. We show that all but the lowest frequencies of the variant allele frequency distribution are dominated by single mutation "waves," which approximately follow an exponential distribution in terms of size. Consequently, large swaths of empty frequencies are observed in the variant allele frequency distributions, with a few frequencies having numbers of mutations far in excess of the expected average values over multiple realizations. We quantify the amount of time each frequency is empty of mutations and further show that the discrete mutation waves average out to a continuous distribution named as the wave frequency distribution, the shape of which is predictable based on few model parameters.

Global and local genomic features together modulate the spontaneous single nucleotide mutation rate

Spontaneous mutations are evolutionary engines as they generate variants for the evolutionary downstream processes that give rise to speciation and adaptation. Single nucleotide mutations (SNM) are the most abundant type of mutations among them. Here, we perform a meta-analysis to quantify the influence of selected global genomic parameters (genome size, genomic GC content, genomic repeat fraction, number of coding genes, gene count, and strand bias in prokaryotes) and local genomic features (local GC content, repeat content, CpG content and the number of SNM at CpG islands) on spontaneous SNM rates across the tree of life (prokaryotes, unicellular eukaryotes, multicellular eukaryotes) using wild-type sequence data in two different taxon classification systems. We find that the spontaneous SNM rates in our data are correlated with many genomic features in prokaryotes and unicellular eukaryotes irrespective of their sample sizes. On the other hand, only the number of coding genes was correlated with the spontaneous SNM rates in multicellular eukaryotes primarily contributed by vertebrates data. Considering local features, we notice that local GC content and CpG content significantly were correlated with the spontaneous SNM rates in the unicellular eukaryotes, while local repeat fraction is an important feature in prokaryotes and certain specific uni- and multi-cellular eukaryotes. Such predictive features of the spontaneous SNM rates often support non-linear models as the best fit compared to the linear model. We also observe that the strand asymmetry in prokaryotes plays an important role in determining the spontaneous SNM rates but the SNM spectrum does not.

Promoters Constrain Evolution of Expression Levels of Essential Genes in Escherichia coli

Variability in expression levels in response to random genomic mutations varies among genes, influencing both the facilitation and constraint of phenotypic evolution in organisms. Despite its importance, both the underlying mechanisms and evolutionary origins of this variability remain largely unknown due to the mixed contributions of cis- and trans-acting elements. To address this issue, we focused on the mutational variability of cis-acting elements, that is, promoter regions, in Escherichia coli. Random mutations were introduced into the natural and synthetic promoters to generate mutant promoter libraries. By comparing the variance in promoter activity of these mutant libraries, we found no significant difference in mutational variability in promoter activity between promoter groups, suggesting the absence of a signature of natural selection for mutational robustness. In contrast, the promoters controlling essential genes exhibited a remarkable bias in mutational variability, with mutants displaying higher activities than the wild types being relatively rare compared to those with lower activities. Our evolutionary simulation on a rugged fitness landscape provided a rationale for this vulnerability. These findings suggest that past selection created nonuniform mutational variability in promoters biased toward lower activities of random mutants, which now constrains the future evolution of downstream essential genes toward higher expression levels.

Irc20 modulates LOH frequency and distribution in S. cerevisiae

Loss of Heterozygosity (LOH) due to mitotic recombination is frequently associated with the development of various cancers (e.g. retinoblastoma). LOH is also an important source of genetic diversity, especially in organisms where meiosis is infrequent. Irc20 is a putative helicase, and E3 ubiquitin ligase involved in DNA double-strand break repair pathway. We analyzed genome-wide LOH events, gross chromosomal changes, small insertion-deletions and single nucleotide mutations in eleven S. cerevisiae mutation accumulation lines of irc20∆, which underwent 50 mitotic bottlenecks. LOH enhancement in irc20∆ was small (1.6 fold), but statistically significant as compared to the wild type. Short (≤ 1 kb) and long (> 10 kb) LOH tracts were significantly enhanced in irc20∆. Both interstitial and terminal LOH events were also significantly enhanced in irc20∆ compared to the wild type. LOH events in irc20∆ were more telomere proximal and away from centromeres compared to the wild type. Gross chromosomal changes, single nucleotide mutations and in-dels were comparable between irc20∆ and wild type. Locus based and genome-wide analysis of meiotic recombination showed that meiotic crossover frequencies are not altered in irc20∆. These results suggest Irc20 primarily regulates mitotic recombination and does not affect meiotic crossovers. Our results suggest that the IRC20 gene is important for regulating LOH frequency and distribution.

Little impact of new mutations on mammalian trait variation

New mutations provide the source of all genetic variation but their impact on trait variation remains poorly understood. A new study published in PLOS Biology addresses this question, finding that new mutations exert only weak effects on some traits in mice.

Design principles of multi-map variation in biological systems

Complexity in biology is often described using a multi-map hierarchical architecture, where the genotype, representing the encoded information, is mapped to the functional level, known as the phenotype, which is then connected to a latent phenotype we refer to as fitness. This underlying architecture governs the processes driving evolution. Furthermore, natural selection, along with other neutral forces, can, in turn, modify these maps. At each level, variation is observed. Here, I propose the need to establish principles that can aid in understanding the transformation of variation within this multi-map architecture. Specifically, I will introduce three, related to the presence of modulators, constraints, and the modular channeling of variation. By comprehending these design principles in various biological systems, we can gain better insights into the mechanisms underlying these maps and how they ultimately contribute to evolutionary dynamics.

Mutations in yeast are deleterious on average regardless of the degree of adaptation to the testing environment

The role of spontaneous mutations in evolution depends on the distribution of their effects on fitness. Despite a general consensus that new mutations are deleterious on average, a handful of mutation accumulation experiments in diverse organisms instead suggest that beneficial and deleterious mutations can have comparable fitness impacts, i.e. the product of their respective rates and effects can be roughly equal. We currently lack a general framework for predicting when such a pattern will occur. One idea is that beneficial mutations will be more evident in genotypes that are not well adapted to the testing environment. We tested this prediction experimentally in the laboratory yeast Saccharomyces cerevisiae by allowing nine replicate populations to adapt to novel environments with complex sets of stressors. After >1000 asexual generations interspersed with 41 rounds of sexual reproduction, we assessed the mean effect of induced mutations on yeast growth in both the environment to which they had been adapting and the alternative novel environment. The mutations were deleterious on average, with the severity depending on the testing environment. However, we found no evidence that the adaptive match between genotype and environment is predictive of mutational fitness effects.

Mutation Rate and Effective Population Size of the Model Cooperative Bacterium Myxococcus xanthus

Intrinsic rates of genetic mutation have diverged greatly across taxa and exhibit statistical associations with several other parameters and features. These include effective population size (Ne), genome size, and gametic multicellularity, with the latter being associated with both increased mutation rates and decreased effective population sizes. However, data sufficient to test for possible relationships between microbial multicellularity and mutation rate (µ) are lacking. Here, we report estimates of two key population-genetic parameters, Ne and µ, for Myxococcus xanthus, a bacterial model organism for the study of aggregative multicellular development, predation, and social swarming. To estimate µ, we conducted an ∼400-day mutation accumulation experiment with 46 lineages subjected to regular single colony bottlenecks prior to clonal regrowth. Upon conclusion, we sequenced one clonal-isolate genome per lineage. Given collective evolution for 85,323 generations across all lines, we calculate a per base-pair mutation rate of ∼5.5 × 10-10 per site per generation, one of the highest mutation rates among free-living eubacteria. Given our estimate of µ, we derived Ne at ∼107 from neutral diversity at four-fold degenerate sites across two dozen M. xanthus natural isolates. This estimate is below average for eubacteria and strengthens an already clear negative correlation between µ and Ne in prokaryotes. The higher and lower than average mutation rate and Ne for M. xanthus, respectively, amplify the question of whether any features of its multicellular life cycle-such as group-size reduction during fruiting-body development-or its highly structured spatial distribution have significantly influenced how these parameters have evolved.

Epistasis between mutator alleles contributes to germline mutation spectrum variability in laboratory mice

Maintaining germline genome integrity is essential and enormously complex. Although many proteins are involved in DNA replication, proofreading, and repair, mutator alleles have largely eluded detection in mammals. DNA replication and repair proteins often recognize sequence motifs or excise lesions at specific nucleotides. Thus, we might expect that the spectrum of de novo mutations - the frequencies of C>T, A>G, etc. - will differ between genomes that harbor either a mutator or wild-type allele. Previously, we used quantitative trait locus mapping to discover candidate mutator alleles in the DNA repair gene Mutyh that increased the C>A germline mutation rate in a family of inbred mice known as the BXDs (Sasani et al., 2022, Ashbrook et al., 2021). In this study we developed a new method to detect alleles associated with mutation spectrum variation and applied it to mutation data from the BXDs. We discovered an additional C>A mutator locus on chromosome 6 that overlaps Ogg1, a DNA glycosylase involved in the same base-excision repair network as Mutyh (David et al., 2007). Its effect depends on the presence of a mutator allele near Mutyh, and BXDs with mutator alleles at both loci have greater numbers of C>A mutations than those with mutator alleles at either locus alone. Our new methods for analyzing mutation spectra reveal evidence of epistasis between germline mutator alleles and may be applicable to mutation data from humans and other model organisms.

Gametic selection favours polyandry and selfing

Competition among pollen or sperm (gametic selection) can cause evolution. Mating systems shape the intensity of gametic selection by determining the competitors involved, which can in turn cause the mating system itself to evolve. We model the bidirectional relationship between gametic selection and mating systems, focusing on variation in female mating frequency (monandry-polyandry) and self-fertilisation (selfing-outcrossing). First, we find that monandry and selfing both reduce the efficiency of gametic selection in removing deleterious alleles. This means that selfing can increase mutation load, in contrast to cases without gametic selection where selfing purges deleterious mutations and decreases mutation load. Second, we explore how mating systems evolve via their effect on gametic selection. By manipulating gametic selection, polyandry can evolve to increase the fitness of the offspring produced. However, this indirect advantage of post-copulatory sexual selection is weak and is likely to be overwhelmed by any direct fitness effects of mating systems. Nevertheless, gametic selection can be potentially decisive for selfing evolution because it significantly reduces inbreeding depression, which favours selfing. Thus, the presence of gametic selection could be a key factor driving selfing evolution.

Genetic drift promotes and recombination hinders speciation on holey fitness landscapes

Dobzhansky and Muller proposed a general mechanism through which microevolution, the substitution of alleles within populations, can cause the evolution of reproductive isolation between populations and, therefore, macroevolution. As allopatric populations diverge, many combinations of alleles differing between them have not been tested by natural selection and may thus be incompatible. Such genetic incompatibilities often cause low fitness in hybrids between species. Furthermore, the number of incompatibilities grows with the genetic distance between diverging populations. However, what determines the rate and pattern of accumulation of incompatibilities remains unclear. We investigate this question by simulating evolution on holey fitness landscapes on which genetic incompatibilities can be identified unambiguously. We find that genetic incompatibilities accumulate more slowly among genetically robust populations and identify two determinants of the accumulation rate: recombination rate and population size. In large populations with abundant genetic variation, recombination selects for increased genetic robustness and, consequently, incompatibilities accumulate more slowly. In small populations, genetic drift interferes with this process and promotes the accumulation of genetic incompatibilities. Our results suggest a novel mechanism by which genetic drift promotes and recombination hinders speciation.

Sex ratio affects sexual selection against mutant alleles in a locus-specific way

Higher male:female operational sex ratio (OSR) is often assumed to lead to stronger sexual selection on males. Yet, this premise has been directly tested by very few studies, with mixed outcomes. We investigated how OSR affects the strength of sexual selection against two deleterious alleles, a natural ebony mutant and a transgenic GFP insertion, in Drosophila melanogaster. To this end, we estimated the relative paternity share of homozygous mutant males competing against wild-type males under different OSRs (1:2, 1:1, 2:1). We also manipulated the mating pool density (18, 36, or 54 individuals) and assessed paternity over three consecutive days, during which the nature of sexual interaction changed. The strength of sexual selection against the ebony mutant increased with OSR, became weaker after the first day, and was little affected by density. In contrast, sexual selection against the GFP transgene was markedly affected by density: at the highest density, it increased with OSR, but at lower densities, it was strongest at 1:1 OSR, remaining strong throughout the experiment. Thus, while OSR can strongly affect the strength of sexual selection against "bad genes," it does not necessarily increase monotonically with male:female OSR. Furthermore, the pattern of relationship between OSR and the strength of sexual selection can be locus-specific, likely reflecting the specific phenotypic effects of the mutation.

Parthenogenomics: Insights on mutation rates and nucleotide diversity in parthenogenetic Panagrolaimus nematodes

Asexual reproduction is assumed to lead to the accumulation of deleterious mutations, and reduced heterozygosity due to the absence of recombination. Panagrolaimid nematode species display different modes of reproduction. Sexual reproduction with distinct males and females, asexual reproduction through parthenogenesis in the genus Panagrolaimus, and hermaphroditism in Propanagrolaimus. Here, we compared genomic features of free-living nematodes in populations and species isolated from geographically distant regions to study diversity, and genome-wide differentiation under different modes of reproduction. We firstly estimated genome-wide spontaneous mutation rates in a triploid parthenogenetic Panagrolaimus, and a diploid hermaphroditic Propanagrolaimus via long-term mutation accumulation lines. Secondly, we calculated population genetic parameters including nucleotide diversity, and fixation index (F ST) between populations of asexually and sexually reproducing nematodes. Thirdly, we used phylogenetic network methods on sexually and asexually reproducing Panagrolaimus populations to understand evolutionary relationships between them. The estimated mutation rate was slightly lower for the asexual population, as expected for taxa with this reproductive mode. Natural polyploid asexual populations revealed higher nucleotide diversity. Despite their common ancestor, a gene network revealed a high level of genetic differentiation among asexual populations. The elevated heterozygosity found in the triploid parthenogens could be explained by the third genome copy. Given their tendentially lower mutation rates it can be hypothesized that this is part of the mechanism to evade Muller's ratchet. Our findings in parthenogenetic triploid nematode populations seem to challenge common expectations of evolution under asexuality.

Progress in gene editing tools, implications and success in plants: a review

Genetic modifications are made through diverse mutagenesis techniques for crop improvement programs. Among these mutagenesis tools, the traditional methods involve chemical and radiation-induced mutagenesis, resulting in off-target and unintended mutations in the genome. However, recent advances have introduced site-directed nucleases (SDNs) for gene editing, significantly reducing off-target changes in the genome compared to induced mutagenesis and naturally occurring mutations in breeding populations. SDNs have revolutionized genetic engineering, enabling precise gene editing in recent decades. One widely used method, homology-directed repair (HDR), has been effective for accurate base substitution and gene alterations in some plant species. However, its application has been limited due to the inefficiency of HDR in plant cells and the prevalence of the error-prone repair pathway known as non-homologous end joining (NHEJ). The discovery of CRISPR-Cas has been a game-changer in this field. This system induces mutations by creating double-strand breaks (DSBs) in the genome and repairing them through associated repair pathways like NHEJ. As a result, the CRISPR-Cas system has been extensively used to transform plants for gene function analysis and to enhance desirable traits. Researchers have made significant progress in genetic engineering in recent years, particularly in understanding the CRISPR-Cas mechanism. This has led to various CRISPR-Cas variants, including CRISPR-Cas13, CRISPR interference, CRISPR activation, base editors, primes editors, and CRASPASE, a new CRISPR-Cas system for genetic engineering that cleaves proteins. Moreover, gene editing technologies like the prime editor and base editor approaches offer excellent opportunities for plant genome engineering. These cutting-edge tools have opened up new avenues for rapidly manipulating plant genomes. This review article provides a comprehensive overview of the current state of plant genetic engineering, focusing on recently developed tools for gene alteration and their potential applications in plant research.

Using Computational Simulations to Model Deleterious Variation and Genetic Load in Natural Populations

AbstractDeleterious genetic variation is abundant in wild populations, and understanding the ecological and conservation implications of such variation is an area of active research. Genomic methods are increasingly used to quantify the impacts of deleterious variation in natural populations; however, these approaches remain limited by an inability to accurately predict the selective and dominance effects of mutations. Computational simulations of deleterious variation offer a complementary tool that can help overcome these limitations, although such approaches have yet to be widely employed. In this perspective article, we aim to encourage ecological and conservation genomics researchers to adopt greater use of computational simulations to aid in deepening our understanding of deleterious variation in natural populations. We first provide an overview of the components of a simulation of deleterious variation, describing the key parameters involved in such models. Next, we discuss several approaches for validating simulation models. Finally, we compare and validate several recently proposed deleterious mutation models, demonstrating that models based on estimates of selection parameters from experimental systems are biased toward highly deleterious mutations. We describe a new model that is supported by multiple orthogonal lines of evidence and provide example scripts for implementing this model (https://github.com/ckyriazis/simulations_review).

Epistasis between mutator alleles contributes to germline mutation spectra variability in laboratory mice

Maintaining germline genome integrity is essential and enormously complex. Although many proteins are involved in DNA replication, proofreading, and repair [1], mutator alleles have largely eluded detection in mammals. DNA replication and repair proteins often recognize sequence motifs or excise lesions at specific nucleotides. Thus, we might expect that the spectrum of de novo mutations - the frequencies of C>T, A>G, etc. - will differ between genomes that harbor either a mutator or wild-type allele. Previously, we used quantitative trait locus mapping to discover candidate mutator alleles in the DNA repair gene Mutyh that increased the C>A germline mutation rate in a family of inbred mice known as the BXDs [2,3]. In this study we developed a new method to detect alleles associated with mutation spectrum variation and applied it to mutation data from the BXDs. We discovered an additional C>A mutator locus on chromosome 6 that overlaps Ogg1, a DNA glycosylase involved in the same base-excision repair network as Mutyh [4]. Its effect depended on the presence of a mutator allele near Mutyh, and BXDs with mutator alleles at both loci had greater numbers of C>A mutations than those with mutator alleles at either locus alone. Our new methods for analyzing mutation spectra reveal evidence of epistasis between germline mutator alleles and may be applicable to mutation data from humans and other model organisms.

Evidence of Site-Specific and Male-Biased Germline Mutation Rate in a Wild Songbird

Germline mutations are the ultimate source of genetic variation and the raw material for organismal evolution. Despite their significance, the frequency and genomic locations of mutations, as well as potential sex bias, are yet to be widely investigated in most species. To address these gaps, we conducted whole-genome sequencing of 12 great reed warblers (Acrocephalus arundinaceus) in a pedigree spanning 3 generations to identify single-nucleotide de novo mutations (DNMs) and estimate the germline mutation rate. We detected 82 DNMs within the pedigree, primarily enriched at CpG sites but otherwise randomly located along the chromosomes. Furthermore, we observed a pronounced sex bias in DNM occurrence, with male warblers exhibiting three times more mutations than females. After correction for false negatives and adjusting for callable sites, we obtained a mutation rate of 7.16 × 10-9 mutations per site per generation (m/s/g) for the autosomes and 5.10 × 10-9 m/s/g for the Z chromosome. To demonstrate the utility of species-specific mutation rates, we applied our autosomal mutation rate in models reconstructing the demographic history of the great reed warbler. We uncovered signs of drastic population size reductions predating the last glacial period (LGP) and reduced gene flow between western and eastern populations during the LGP. In conclusion, our results provide one of the few direct estimates of the mutation rate in wild songbirds and evidence for male-driven mutations in accordance with theoretical expectations.

Contribution of Spontaneous Mutations to Quantitative and Molecular Variation at the Highly Repetitive rDNA Locus in Yeast

The ribosomal DNA array in Saccharomyces cerevisiae consists of many tandem repeats whose copy number is believed to be functionally important but highly labile. Regulatory mechanisms have evolved to maintain copy number by directed mutation, but how spontaneous variation at this locus is generated and selected has not been well characterized. We applied a mutation accumulation approach to quantify the impacts of mutation and selection on this unique genomic feature across hundreds of mutant strains. We find that mutational variance for this trait is relatively high, and that unselected mutations elsewhere in the genome can disrupt copy number maintenance. In consequence, copy number generally declines gradually, consistent with a previously proposed model of rDNA maintenance where a downward mutational bias is normally compensated by mechanisms that increase copy number when it is low. This pattern holds across ploidy levels and strains in the standard lab environment but differs under some stressful conditions. We identify several alleles, gene categories, and genomic features that likely affect copy number, including aneuploidy for chromosome XII. Copy number change is associated with reduced growth in diploids, consistent with stabilizing selection. Levels of standing variation in copy number are well predicted by a balance between mutation and stabilizing selection, suggesting this trait is not subject to strong diversifying selection in the wild. The rate and spectrum of point mutations within the rDNA locus itself are distinct from the rest of the genome and predictive of polymorphism locations. Our findings help differentiate the roles of mutation and selection and indicate that spontaneous mutation patterns shape several aspects of ribosomal DNA evolution.

Genus-Wide Transcriptional Landscapes Reveal Correlated Gene Networks Underlying Microevolutionary Divergence in Diatoms

Marine microbes like diatoms make up the base of marine food webs and drive global nutrient cycles. Despite their key roles in ecology, biogeochemistry, and biotechnology, we have limited empirical data on how forces other than adaptation may drive diatom diversification, especially in the absence of environmental change. One key feature of diatom populations is frequent extreme reductions in population size, which can occur both in situ and ex situ as part of bloom-and-bust growth dynamics. This can drive divergence between closely related lineages, even in the absence of environmental differences. Here, we combine experimental evolution and transcriptome landscapes (t-scapes) to reveal repeated evolutionary divergence within several species of diatoms in a constant environment. We show that most of the transcriptional divergence can be captured on a reduced set of axes, and that repeatable evolution can occur along a single major axis of variation defined by core ortholog expression comprising common metabolic pathways. Previous work has associated specific transcriptional changes in gene networks with environmental factors. Here, we find that these same gene networks diverge in the absence of environmental change, suggesting these pathways may be central in generating phenotypic diversity as a result of both selective and random evolutionary forces. If this is the case, these genes and the functions they encode may represent universal axes of variation. Such axes that capture suites of interacting transcriptional changes during diversification improve our understanding of both global patterns in local adaptation and microdiversity, as well as evolutionary forces shaping algal cultivation.

Inference of the distribution of fitness effects of mutations is affected by single nucleotide polymorphism filtering methods, sample size and population structure

The distribution of fitness effects (DFE) of new mutations has been of interest to evolutionary biologists since the concept of mutations arose. Modern population genomic data enable us to quantify the DFE empirically, but few studies have examined how data processing, sample size and cryptic population structure might affect the accuracy of DFE inference. We used simulated and empirical data (from Arabidopsis lyrata) to show the effects of missing data filtering, sample size, number of single nucleotide polymorphisms (SNPs) and population structure on the accuracy and variance of DFE estimates. Our analyses focus on three filtering methods-downsampling, imputation and subsampling-with sample sizes of 4-100 individuals. We show that (1) the choice of missing-data treatment directly affects the estimated DFE, with downsampling performing better than imputation and subsampling; (2) the estimated DFE is less reliable in small samples (<8 individuals), and becomes unpredictable with too few SNPs (<5000, the sum of 0- and 4-fold SNPs); and (3) population structure may skew the inferred DFE towards more strongly deleterious mutations. We suggest that future studies should consider downsampling for small data sets, and use samples larger than 4 (ideally larger than 8) individuals, with more than 5000 SNPs in order to improve the robustness of DFE inference and enable comparative analyses.

The Contribution of Mutation to Variation in Temperature-Dependent Sprint Speed in Zebrafish, Danio rerio

AbstractThe contribution of new mutations to phenotypic variation and the consequences of this variation for individual fitness are fundamental concepts for understanding genetic variation and adaptation. Here, we investigated how mutation influenced variation in a complex trait in zebrafish, Danio rerio. Typical of many ecologically relevant traits in ectotherms, swimming speed in fish is temperature dependent, with evidence of adaptive evolution of thermal performance. We chemically induced novel germline point mutations in males and measured sprint speed in their sons at six temperatures (between 16°C and 34°C). Heterozygous mutational effects on speed were strongly positively correlated among temperatures, resulting in statistical support for only a single axis of mutational variation, reflecting temperature-independent variation in speed (faster-slower mode). These results suggest pleiotropic effects on speed across different temperatures; however, spurious correlations arise via linkage or heterogeneity in mutation number when mutations have consistent directional effects on each trait. Here, mutation did not change mean speed, indicating no directional bias in mutational effects. The results contribute to emerging evidence that mutations may predominantly have synergistic cross-environment effects, in contrast to conditionally neutral or antagonistic effects that underpin thermal adaptation. We discuss several aspects of experimental design that may affect resolution of mutations with nonsynergistic effects.

Life history and deleterious mutation rate coevolution

The cost of germline maintenance gives rise to a trade-off between lowering the deleterious mutation rate and investing in life history functions. Therefore, life history and the mutation rate coevolve, but this coevolution is not well understood. We develop a mathematical model to analyse the evolution of resource allocation traits, which simultaneously affect life history and the deleterious mutation rate. First, we show that the invasion fitness of such resource allocation traits can be approximated by the basic reproductive number of the least-loaded class; the expected lifetime production of offspring without deleterious mutations born to individuals without deleterious mutations. Second, we apply the model to investigate (i) the coevolution of reproductive effort and germline maintenance and (ii) the coevolution of age-at-maturity and germline maintenance. This analysis provides two resource allocation predictions when exposure to environmental mutagens is higher. First, selection favours higher allocation to germline maintenance, even if it comes at the expense of life history functions, and leads to a shift in allocation towards reproduction rather than survival. Second, life histories tend to be faster, characterised by individuals with shorter lifespans and smaller body sizes at maturity. Our results suggest that mutation accumulation via the cost of germline maintenance can be a major force shaping life-history traits.

Purging due to self-fertilization does not prevent accumulation of expansion load

As species expand their geographic ranges, colonizing populations face novel ecological conditions, such as new environments and limited mates, and suffer from evolutionary consequences of demographic change through bottlenecks and mutation load accumulation. Self-fertilization is often observed at species range edges and, in addition to countering the lack of mates, is hypothesized as an evolutionary advantage against load accumulation through increased homozygosity and purging. We study how selfing impacts the accumulation of genetic load during range expansion via purging and/or speed of colonization. Using simulations, we disentangle inbreeding effects due to demography versus due to selfing and find that selfers expand faster, but still accumulate load, regardless of mating system. The severity of variants contributing to this load, however, differs across mating system: higher selfing rates purge large-effect recessive variants leaving a burden of smaller-effect alleles. We compare these predictions to the mixed-mating plant Arabis alpina, using whole-genome sequences from refugial outcrossing populations versus expanded selfing populations. Empirical results indicate accumulation of expansion load along with evidence of purging in selfing populations, concordant with our simulations, suggesting that while purging is a benefit of selfing evolving during range expansions, it is not sufficient to prevent load accumulation due to range expansion.

Inferring the distributions of fitness effects and proportions of strongly deleterious mutations

The distribution of fitness effects is a key property in evolutionary genetics as it has implications for several evolutionary phenomena including the evolution of sex and mating systems, the rate of adaptive evolution, and the prevalence of deleterious mutations. Despite the distribution of fitness effects being extensively studied, the effects of strongly deleterious mutations are difficult to infer since such mutations are unlikely to be present in a sample of haplotypes, so genetic data may contain very little information about them. Recent work has attempted to correct for this issue by expanding the classic gamma-distributed model to explicitly account for strongly deleterious mutations. Here, we use simulations to investigate one such method, adding a parameter (plth) to capture the proportion of strongly deleterious mutations. We show that plth can improve the model fit when applied to individual species but underestimates the true proportion of strongly deleterious mutations. The parameter can also artificially maximize the likelihood when used to jointly infer a distribution of fitness effects from multiple species. As plth and related parameters are used in current inference algorithms, our results are relevant with respect to avoiding model artifacts and improving future tools for inferring the distribution of fitness effects.

Evolutionary codependency: insights into the mitonuclear interaction landscape from experimental and wild Caenorhabditis nematodes

Aided by new technologies, the upsurgence of research into mitochondrial genome biology during the past 15 years suggests that we have misunderstood, and perhaps dramatically underestimated, the ongoing biological and evolutionary significance of our long-time symbiotic partner. While we have begun to scratch the surface of several topics, many questions regarding the nature of mutation and selection in the mitochondrial genome, and the nature of its relationship to the nuclear genome, remain unanswered. Although best known for their contributions to studies of developmental and aging biology, Caenorhabditis nematodes are increasingly recognized as excellent model systems to advance understanding in these areas. We review recent discoveries with relevance to mitonuclear coevolution and conflict and offer several fertile areas for future work.

Experimental estimates of germline mutation rate in eukaryotes: a phylogenetic meta-analysis

Mutation is the ultimate source of all genetic variation, and over the last 10 years the ready availability of whole-genome sequencing has permitted direct estimation of mutation rate for many non-model species across the tree of life. In this meta-analysis, we make a comprehensive search of the literature for mutation rate estimates in eukaryotes, identifying 140 mutation accumulation (MA) and parent-offspring (PO) sequencing studies covering 134 species. Based on these data, we revisit differences in the single-nucleotide mutation (SNM) rate between different phylogenetic lineages and update the known relationships between mutation rate and generation time, genome size, and nucleotide diversity-while accounting for phylogenetic nonindependence. We do not find a significant difference between MA and PO in estimated mutation rates, but we confirm that mammal and plant lineages have higher mutation rates than arthropods and that unicellular eukaryotes have the lowest mutation rates. We find that mutation rates are higher in species with longer generation times and larger genome sizes, even when accounting for phylogenetic relationships. Moreover, although nucleotide diversity is positively correlated with mutation rate, the gradient of the relationship is significantly less than one (on a logarithmic scale), consistent with higher mutation rates in populations with smaller effective size. For the 29 species for which data are available, we find that indel mutation rates are positively correlated with nucleotide mutation rates and that short deletions are generally more common than short insertions. Nevertheless, despite recent progress, no estimates of either SNM or indel mutation rates are available for the majority of deeply branching eukaryotic lineages-or even for most animal phyla. Even among charismatic megafauna, experimental mutation rate estimates remain unknown for amphibia and scarce for reptiles and fish.

Population transcriptogenomics highlights impaired metabolism and small population sizes in tree frogs living in the Chernobyl Exclusion Zone

BACKGROUND

Individual functional modifications shape the ability of wildlife populations to cope with anthropogenic environmental changes. But instead of adaptive response, human-altered environments can generate a succession of deleterious functional changes leading to the extinction of the population. To study how persistent anthropogenic changes impacted local species' population status, we characterised population structure, genetic diversity and individual response of gene expression in the tree frog Hyla orientalis along a gradient of radioactive contamination around the Chernobyl nuclear power plant.

RESULTS

We detected lower effective population size in populations most exposed to ionizing radiation in the Chernobyl Exclusion Zone that is not compensated by migrations from surrounding areas. We also highlighted a decreased body condition of frogs living in the most contaminated area, a distinctive transcriptomics signature and stop-gained mutations in genes involved in energy metabolism. While the association with dose will remain correlational until further experiments, a body of evidence suggests the direct or indirect involvement of radiation exposure in these changes.

CONCLUSIONS

Despite ongoing migration and lower total dose rates absorbed than at the time of the accident, our results demonstrate that Hyla orientalis specimens living in the Chernobyl Exclusion Zone are still undergoing deleterious changes, emphasizing the long-term impacts of the nuclear disaster.

An empirical evaluation of the estimation of inbreeding depression from molecular markers under suboptimal conditions

Inbreeding depression (ID), the reduction in fitness due to inbreeding, is typically measured by the regression of the phenotypic values of individuals for a particular trait on their corresponding inbreeding coefficients (F). While genealogical records can provide these coefficients, they may be unavailable or incomplete, making molecular markers a useful alternative. The power to detect ID and its accuracy depend on the variation of F values of individuals, the sample sizes available, and the accuracy in the estimation of individual fitness traits and F values. In this study, we used Drosophila melanogaster to evaluate the effectiveness of molecular markers in estimating ID under suboptimal conditions. We generated two sets of 100 pairs of unrelated individuals from a large panmictic population and mated them for two generations to produce non-inbred and unrelated individuals (F = 0) and inbred individuals (full-sib progeny; F = 0.25). Using these expected genealogical F values, we calculated inbreeding depression for two fitness-related traits, pupae productivity and competitive fitness. We then sequenced the males from 17 non-inbred pairs and 17 inbred pairs to obtain their genomic inbreeding coefficients and estimate ID for the two traits. The scenario assumed was rather restrictive in terms of estimation of ID because: (1) the individuals belonged to the same generation of a large panmictic population, leading to low variation in individual F coefficients; (2) the sample sizes were small; and (3) the traits measured depended on both males and females while only males were sequenced. Despite the challenging conditions of our study, we found that molecular markers provided estimates of ID that were comparable to those obtained from simple pedigree estimations with larger sample sizes. The results therefore suggest that genomic measures of inbreeding are useful to provide estimates of inbreeding depression even under very challenging scenarios.

A null model for the distribution of fitness effects of mutations

The distribution of fitness effects (DFE) of new mutations is key to our understanding of many evolutionary processes. Theoreticians have developed several models to help understand the patterns seen in empirical DFEs. Many such models reproduce the broad patterns seen in empirical DFEs but these models often rely on structural assumptions that cannot be tested empirically. Here, we investigate how much of the underlying "microscopic" biological processes involved in the mapping of new mutations to fitness can be inferred from "macroscopic" observations of the DFE. We develop a null model by generating random genotype-to-fitness maps and show that the null DFE is that with the largest possible information entropy. We further show that, subject to one simple constraint, this null DFE is a Gompertz distribution. Finally, we illustrate how the predictions of this null DFE match empirically measured DFEs from several datasets, as well as DFEs simulated from Fisher's geometric model. This suggests that a match between models and empirical data is often not a very strong indication of the mechanisms underlying the mapping of mutation to fitness.

Low Spontaneous Mutation Rate in Complex Multicellular Eukaryotes with a Haploid-Diploid Life Cycle

The spontaneous mutation rate µ is a crucial parameter to understand evolution and biodiversity. Mutation rates are highly variable across species, suggesting that µ is susceptible to selection and drift and that species life cycle and life history may impact its evolution. In particular, asexual reproduction and haploid selection are expected to affect the mutation rate, but very little empirical data are available to test this expectation. Here, we sequence 30 genomes of a parent-offspring pedigree in the model brown alga Ectocarpus sp.7, and 137 genomes of an interspecific cross of the closely related brown alga Scytosiphon to have access to the spontaneous mutation rate of representative organisms of a complex multicellular eukaryotic lineage outside animals and plants, and to evaluate the potential impact of life cycle on the mutation rate. Brown algae alternate between a haploid and a diploid stage, both multicellular and free living, and utilize both sexual and asexual reproduction. They are, therefore, excellent models to empirically test expectations of the effect of asexual reproduction and haploid selection on mutation rate evolution. We estimate that Ectocarpus has a base substitution rate of µbs = 4.07 × 10-10 per site per generation, whereas the Scytosiphon interspecific cross had µbs = 1.22 × 10-9. Overall, our estimations suggest that these brown algae, despite being multicellular complex eukaryotes, have unusually low mutation rates. In Ectocarpus, effective population size (Ne) could not entirely explain the low µbs. We propose that the haploid-diploid life cycle, combined with extensive asexual reproduction, may be additional key drivers of the mutation rate in these organisms.

Attempts to limit sporulation in the probiotic strain Bacillus subtilis BG01-4TM through mutation accumulation and selection

The use of bacterial spores in probiotics over viable loads of bacteria has many advantages, including the durability of spores, which allows spore-based probiotics to effectively traverse the various biochemical barriers present in the gastrointestinal tract. However, the majority of spore-based probiotics developed currently aim to treat adults, and there is a litany of differences between the adult and infant intestinal systems, including the immaturity and low microbial species diversity observed within the intestines of infants. These differences are only further exacerbated in premature infants with necrotizing enterocolitis (NEC) and indicates that what may be appropriate for an adult or even a healthy full-term infant may not be suited for an unhealthy premature infant. Complications from using spore-based probiotics for premature infants with NEC may involve the spores remaining dormant and adhering to the intestinal epithelia, the out-competing of commensal bacteria by spores, and most importantly the innate antibiotic resistance of spores. Also, the ability of Bacillus subtilis to produce spores under duress may result in less B. subtilis perishing within the intestines and releasing membrane branched-chain fatty acids. The isolate B. subtilis BG01-4^TM is a proprietary strain developed by Vernx Biotechnology through accumulating mutations within the BG01-4^TM genome in a serial batch culture. Strain BG01-4^TM was provided as a non-spore-forming B. subtilis , but a positive sporulation status for BG01-4^TM was confirmed through in vitro testing and suggested that selection for the sporulation defective genes could occur within an environment that would select against sporulation. The durability of key sporulation genes was ratified in this study, as the ability of BG01-4^TM to produce spores was not eliminated by the attempts to select against sporulation genes in BG01-4^TM by the epigenetic factors of high glucose and low pH. However, a variation in the genes in isolate BG01-4-8 involved in the regulation of sporulation is believed to have occurred during the mutation selection from the parent strain BG01-4^TM. An alteration in selected sporulation regulation genes is expected to have occurred from BG01-4^TM to BG01-4-8, with BG01-4-8 producing spores within 24 h, ~48 h quicker than BG01-4^TM.

Contributions of mutation and selection to regulatory variation: lessons from the Saccharomyces cerevisiae TDH3 gene

Heritable variation in gene expression is common within and among species and contributes to phenotypic diversity. Mutations affecting either cis- or trans-regulatory sequences controlling gene expression give rise to variation in gene expression, and natural selection acting on this variation causes some regulatory variants to persist in a population for longer than others. To understand how mutation and selection interact to produce the patterns of regulatory variation we see within and among species, my colleagues and I have been systematically determining the effects of new mutations on expression of the TDH3 gene in Saccharomyces cerevisiae and comparing them to the effects of polymorphisms segregating within this species. We have also investigated the molecular mechanisms by which regulatory variants act. Over the past decade, this work has revealed properties of cis- and trans-regulatory mutations including their relative frequency, effects, dominance, pleiotropy and fitness consequences. Comparing these mutational effects to the effects of polymorphisms in natural populations, we have inferred selection acting on expression level, expression noise and phenotypic plasticity. Here, I summarize this body of work and synthesize its findings to make inferences not readily discernible from the individual studies alone. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Adaptation of a quantitative trait to a changing environment: New analytical insights on the asexual and infinitesimal sexual models

Predicting the adaptation of populations to a changing environment is crucial to assess the impact of human activities on biodiversity. Many theoretical studies have tackled this issue by modeling the evolution of quantitative traits subject to stabilizing selection around an optimal phenotype, whose value is shifted continuously through time. In this context, the population fate results from the equilibrium distribution of the trait, relative to the moving optimum. Such a distribution may vary with the shape of selection, the system of reproduction, the number of loci, the mutation kernel or their interactions. Here, we develop a methodology that provides quantitative measures of population maladaptation and potential of survival directly from the entire profile of the phenotypic distribution, without any a priori on its shape. We investigate two different systems of reproduction (asexual and infinitesimal sexual models of inheritance), with various forms of selection. In particular, we recover that fitness functions such that selection weakens away from the optimum lead to evolutionary tipping points, with an abrupt collapse of the population when the speed of environmental change is too high. Our unified framework allows deciphering the mechanisms that lead to this phenomenon. More generally, it allows discussing similarities and discrepancies between the two systems of reproduction, which are ultimately explained by different constraints on the evolution of the phenotypic variance. We demonstrate that the mean fitness in the population crucially depends on the shape of the selection function in the infinitesimal sexual model, in contrast with the asexual model. In the asexual model, we also investigate the effect of the mutation kernel and we show that kernels with higher kurtosis tend to reduce maladaptation and improve fitness, especially in fast changing environments.

Distribution of mutation rates challenges evolutionary predictability

Natural selection is commonly assumed to act on extensive standing genetic variation. Yet, accumulating evidence highlights the role of mutational processes creating this genetic variation: to become evolutionarily successful, adaptive mutants must not only reach fixation, but also emerge in the first place, i.e. have a high enough mutation rate. Here, we use numerical simulations to investigate how mutational biases impact our ability to observe rare mutational pathways in the laboratory and to predict outcomes in experimental evolution. We show that unevenness in the rates at which mutational pathways produce adaptive mutants means that most experimental studies lack power to directly observe the full range of adaptive mutations. Modelling mutation rates as a distribution, we show that a substantially larger target size ensures that a pathway mutates more commonly. Therefore, we predict that commonly mutated pathways are conserved between closely related species, but not rarely mutated pathways. This approach formalizes our proposal that most mutations have a lower mutation rate than the average mutation rate measured experimentally. We suggest that the extent of genetic variation is overestimated when based on the average mutation rate.

Silver nanoparticles elevate mutagenesis of eukaryotic genomes

Metal nanoparticles, especially silver, have been used in various medical scenarios, due to their excellent antimicrobial effects. Recent studies have shown that AgNPs do not exert mutagenic effects on target bacteria, but the degree to which they compromise eukaryotic genomes remains unclear. To study this, we evaluated the mutagenic effects of AgNPs on the fission yeast Schizosaccharomyces pombe ATCC-16979, of which ∼23% genes are homologous to human ones, at single-nucleotide resolution, and whole-genome scale by running 283 mutation accumulation lines for ∼260,000 cell divisions in total. We also explored the action and mutagenesis mechanisms using differential gene-expression analysis based on RNAseq. Upon AgNPs treatment, the genomic base-substitution mutation rate of S. pombe at four-fold degenerate sites increased by 3.46×, and small indels were prone to occur in genomic regions that are not simple sequence repeats. The G:C → T:A transversion rate was also significantly increased, likely mostly from oxidative damage. Thus, in addition to their antimicrobial potency, AgNPs might pose slight genotoxicity threats to eukaryotic and possibly human genomes, though at a low magnitude.

Similar mutation rates but different mutation spectra in moderate and extremely halophilic archaea

Archaea are a major part of Earth's microbiota and extremely diverse. Yet, we know very little about the process of mutation that drives such diversification. To expand beyond previous work with the moderate halophilic archaeal species Haloferax volcanii, we performed a mutation-accumulation experiment followed by whole-genome sequencing in the extremely halophilic archaeon Halobacterium salinarum. Although Hfx. volcanii and Hbt. salinarum have different salt requirements, both species have highly polyploid genomes and similar GC content. We accumulated mutations for an average of 1250 generations in 67 mutation accumulation lines of Hbt. salinarum, and revealed 84 single-base substitutions and 10 insertion-deletion mutations. The estimated base-substitution mutation rate of 3.99 × 10-10 per site per generation or 1.0 × 10-3 per genome per generation in Hbt. salinarum is similar to that reported for Hfx. volcanii (1.2 × 10-3 per genome per generation), but the genome-wide insertion-deletion rate and spectrum of mutations are somewhat dissimilar in these archaeal species. The spectra of spontaneous mutations were AT biased in both archaea, but they differed in significant ways that may be related to differences in the fidelity of DNA replication/repair mechanisms or a simple result of the different salt concentrations.

Very Low Rates of Spontaneous Gene Deletions and Gene Duplications in Dictyostelium discoideum

The study of spontaneous mutation rates has revealed a wide range of heritable point mutation rates across species, but there are comparatively few estimates for large-scale deletion and duplication rates. The handful of studies that have directly calculated spontaneous rates of deletion and duplication using mutation accumulation lines have estimated that genes are duplicated and deleted at orders of magnitude greater rates than the spontaneous point mutation rate. In our study, we tested whether spontaneous gene deletion and gene duplication rates are also high in Dictyostelium discoideum, a eukaryote with among the lowest point mutation rates (2.5 × 10-11 per site per generation) and an AT-rich genome (GC content of 22%). We calculated mutation rates of gene deletions and duplications using whole-genome sequencing data originating from a mutation accumulation experiment and determined the association between the copy number mutations and GC content. Overall, we estimated an average of 3.93 × 10-8 gene deletions and 1.18 × 10-8 gene duplications per gene per generation. While orders of magnitude greater than their point mutation rate, these rates are much lower compared to gene deletion and duplication rates estimated from mutation accumulation lines in other organisms (that are on the order of ~ 10-6 per gene/generation). The deletions and duplications were enriched in regions that were AT-rich even compared to the genomic background, in contrast to our expectations if low GC content was contributing to low mutation rates. The low deletion and duplication mutation rates in D. discoideum compared to other eukaryotes mirror their low point mutation rates, supporting previous work suggesting that this organism has high replication fidelity and effective molecular machinery to avoid the accumulation of mutations in their genome.

Rates and spectra of de novo structural mutations in Chlamydomonas reinhardtii

Genetic variation originates from several types of spontaneous mutation, including single-nucleotide substitutions, short insertions and deletions (indels), and larger structural changes. Structural mutations (SMs) drive genome evolution and are thought to play major roles in evolutionary adaptation, speciation, and genetic disease, including cancers. Sequencing of mutation accumulation (MA) lines has provided estimates of rates and spectra of single-nucleotide and indel mutations in many species, yet the rate of new SMs is largely unknown. Here, we use long-read sequencing to determine the full mutation spectrum in MA lines derived from two strains (CC-1952 and CC-2931) of the green alga Chlamydomonas reinhardtii The SM rate is highly variable between strains and between MA lines, and SMs represent a substantial proportion of all mutations in both strains (CC-1952 6%; CC-2931 12%). The SM spectra differ considerably between the two strains, with almost all inversions and translocations occurring in CC-2931 MA lines. This variation is associated with heterogeneity in the number and type of active transposable elements (TEs), which comprise major proportions of SMs in both strains (CC-1952 22%; CC-2931 38%). In CC-2931, a Crypton and a previously undescribed type of DNA element have caused 71% of chromosomal rearrangements, whereas in CC-1952, a Dualen LINE is associated with 87% of duplications. Other SMs, notably large duplications in CC-2931, are likely products of various double-strand break repair pathways. Our results show that diverse types of SMs occur at substantial rates, and support prominent roles for SMs and TEs in evolution.

Multigenerational downregulation of insulin/IGF-1 signaling in adulthood improves lineage survival, reproduction, and fitness in Caenorhabditis elegans supporting the developmental theory of ageing

Adulthood-only downregulation of insulin/IGF-1 signaling (IIS), an evolutionarily conserved pathway regulating resource allocation between somatic maintenance and reproduction, increases life span without fecundity cost in the nematode, Caenorhabditis elegans. However, long-term multigenerational effects of reduced IIS remain unexplored and are proposed to carry costs for offspring quality. To test this hypothesis, we ran a mutation accumulation (MA) experiment and downregulated IIS in half of the 400 MA lines by silencing daf-2 gene expression using RNA interference (RNAi) across 40 generations. Contrary to the prediction, adulthood-only daf-2 RNAi reduced extinction of MA lines both under UV-induced and spontaneous MA. Fitness of the surviving UV-induced MA lines was higher under daf-2 RNAi. Reduced IIS increased intergenerational F1 offspring fitness under UV stress but had no quantifiable transgenerational effects. Functional hrde-1 was required for the benefits of multigenerational daf-2 RNAi. Overall, we found net benefit to fitness from multigenerational reduction of IIS and the benefits became more apparent under stress. Because reduced daf-2 expression during development carries fitness costs, we suggest that our findings are best explained by the developmental theory of ageing, which maintains that the decline in the force of selection with age results in poorly regulated gene expression in adulthood.

Endosymbiotic selective pressure at the origin of eukaryotic cell biology

The dichotomy that separates prokaryotic from eukaryotic cells runs deep. The transition from pro- to eukaryote evolution is poorly understood due to a lack of reliable intermediate forms and definitions regarding the nature of the first host that could no longer be considered a prokaryote, the first eukaryotic common ancestor, FECA. The last eukaryotic common ancestor, LECA, was a complex cell that united all traits characterising eukaryotic biology including a mitochondrion. The role of the endosymbiotic organelle in this radical transition towards complex life forms is, however, sometimes questioned. In particular the discovery of the asgard archaea has stimulated discussions regarding the pre-endosymbiotic complexity of FECA. Here we review differences and similarities among models that view eukaryotic traits as isolated coincidental events in asgard archaeal evolution or, on the contrary, as a result of and in response to endosymbiosis. Inspecting eukaryotic traits from the perspective of the endosymbiont uncovers that eukaryotic cell biology can be explained as having evolved as a solution to housing a semi-autonomous organelle and why the addition of another endosymbiont, the plastid, added no extra compartments. Mitochondria provided the selective pressures for the origin (and continued maintenance) of eukaryotic cell complexity. Moreover, they also provided the energetic benefit throughout eukaryogenesis for evolving thousands of gene families unique to eukaryotes. Hence, a synthesis of the current data lets us conclude that traits such as the Golgi apparatus, the nucleus, autophagosomes, and meiosis and sex evolved as a response to the selective pressures an endosymbiont imposes.