Search

The distribution of fitness effects (DFE) for new mutations is one of the most theoretically important but difficult to estimate properties in population genetics. A crucial challenge to inferring the DFE from natural genetic variation is the sensitivity of the site frequency spectrum to factors like population size change, population substructure, genome structure, and nonrandom mating. Although inference methods aim to control for population size changes, the influence of nonrandom mating remains incompletely understood, despite being a common feature of many species. We report the DFE estimated from 326 genomes of Caenorhabditis elegans, a nematode roundworm with a high rate of self-fertilization. We evaluate the robustness of DFE inferences using simulated data that mimics the genomic structure and reproductive life history of C. elegans. Our observations demonstrate how the combined influence of self-fertilization, genome structure, and natural selection on linked sites can conspire to compromise estimates of the DFE from extant polymorphisms with existing methods. These factors together tend to bias inferences toward weakly deleterious mutations, making it challenging to have full confidence in the inferred DFE of new mutations as deduced from standing genetic variation in species like C. elegans. Improved methods for inferring the DFE are needed to appropriately handle strong linked selection and selfing. These results highlight the importance of understanding the combined effects of processes that can bias our interpretations of evolution in natural populations., (© The Author(s) 2021. Published by Oxford University Press on behalf of Genetics Society of America. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

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.

Androdioecious Caenorhabditis have a high frequency of self-compatible hermaphrodites and a low frequency of males. The effects of mutations on male fitness are of interest for two reasons. First, when males are rare, selection on male-specific mutations is less efficient than in hermaphrodites. Second, males may present a larger mutational target than hermaphrodites because of the different ways in which fitness accrues in the two sexes. We report the first estimates of male-specific mutational effects in an androdioecious organism. The rate of male-specific inviable or sterile mutations is ⩽5 × 10 -4 /generation, below the rate at which males would be lost solely due to those kinds of mutations. The rate of mutational decay of male competitive fitness is ~ 0.17%/generation; that of hermaphrodite competitive fitness is ~ 0.11%/generation. The point estimate of ~ 1.5X faster rate of mutational decay of male fitness is nearly identical to the same ratio in Drosophila. Estimates of mutational variance (V M ) for male mating success and competitive fitness are not significantly different from zero, whereas V M for hermaphrodite competitive fitness is similar to that of non-competitive fitness. Two independent estimates of the average selection coefficient against mutations affecting hermaphrodite competitive fitness agree to within two-fold, 0.33-0.5%.

Different types of phenotypic traits consistently exhibit different levels of genetic variation in natural populations. There are two potential explanations: Either mutation produces genetic variation at different rates or natural selection removes or promotes genetic variation at different rates. Whether mutation or selection is of greater general importance is a longstanding unresolved question in evolutionary genetics. We report mutational variances (VM) for 19 traits related to the first mitotic cell division in Caenorhabditis elegans and compare them to the standing genetic variances (VG) for the same suite of traits in a worldwide collection C. elegans Two robust conclusions emerge. First, the mutational process is highly repeatable: The correlation between VM in two independent sets of mutation accumulation lines is ∼0.9. Second, VM for a trait is a good predictor of VG for that trait: The correlation between VM and VG is ∼0.9. This result is predicted for a population at mutation-selection balance; it is not predicted if balancing selection plays a primary role in maintaining genetic variation., (Copyright © 2016 by the Genetics Society of America.)

Although the effects of deleterious alleles often are predicted to be greater in stressful environments, there is no theoretical basis for this prediction and the empirical evidence is mixed. Here we characterized the effects of three types of abiotic stress (thermal, oxidative and hyperosmotic) on two sets of nematode (Caenorhabditis elegans) mutation accumulation (MA) lines that differ by threefold in fitness. We compared the survival and egg-to-adult viability between environments (benign and stressful) and between fitness categories (high-fitness MA, low-fitness MA). If the environment and mutation load have synergistic effects on trait means, then the difference between the high and low-fitness MA lines should be larger in stressful environments. Although the stress treatments consistently decreased survival and/or viability, we did not detect significant interactions between fitness categories and environment types. In contrast, we did find consistent evidence for synergistic effects on (micro)environmental variation. The lack of signal in trait means likely reflects the very low starting fitness of some low-fitness MA lines, the potential for cross-stress responses and the context dependence of mutational effects. In addition, the large increases in the environmental variance in the stressful environments may have masked small changes in trait means. These results do not provide evidence for synergism between mutation and stress.

Understanding the genetic basis of susceptibility to pathogens is an important goal of medicine and of evolutionary biology. A key first step toward understanding the genetics and evolution of any phenotypic trait is characterizing the role of mutation. However, the rate at which mutation introduces genetic variance for pathogen susceptibility in any organism is essentially unknown. Here, we quantify the per-generation input of genetic variance by mutation (VM) for susceptibility of Caenorhabditis elegans to the pathogenic bacterium Pseudomonas aeruginosa (defined as the median time of death, LT50). VM for LT50 is slightly less than VM for a variety of life-history and morphological traits in this strain of C. elegans, but is well within the range of reported values in a variety of organisms. Mean LT50 did not change significantly over 250 generations of mutation accumulation. Comparison of VM to the standing genetic variance (VG) implies a strength of selection against new mutations of a few tenths of a percent. These results suggest that the substantial standing genetic variation for susceptibility of C. elegans to P. aeruginosa can be explained by polygenic mutation coupled with purifying selection., (© 2014 The Author(s). Evolution © 2014 The Society for the Study of Evolution.)

Evidence is accumulating that individuals in poor physiologic condition may accumulate mutational damage faster than individuals in good condition. If poor condition results from pre-existing deleterious mutations, the result is "fitness-dependent mutation rate," which has interesting theoretical implications. Here we report a study in which 10 mutation accumulation (MA) lines of the nematode Caenorhabditis elegans that had previously accumulated mutations for 250 generations under relaxed selection were expanded into sets of "second-order" MA lines and allowed to accumulate mutations for an additional 150 generations. The 10 lines were chosen on the basis of the relative change in fitness over the first 250 generations of MA, five high-fitness lines and five low-fitness lines. On average, the mutational properties (per-generation change in mean relative fitness, mutational variance, and Bateman-Mukai estimates of genomic mutation rate and average mutational effect) of the high-fitness and low-fitness did not differ significantly, and averaged over all lines, the point estimates were extremely close to those of the first-order MA experiment after 200 generations of MA. However, several nonsignificant trends indicate that low-fitness lines may in fact be more likely to suffer mutational damage than high-fitness lines.

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

Evolutionary theory predicts that the strength of natural selection to reduce the mutation rate should be stronger in self-fertilizing than in outcrossing taxa. However, the relative efficacy of selection on mutation rate relative to the many other factors influencing the evolution of any species is poorly understood. To address this question, we allowed mutations to accumulate for ∼100 generations in several sets of "mutation accumulation" (MA) lines in three species of gonochoristic (dieocious) Caenorhabditis (C. remanei, C. brenneri, C. sp. 5) as well as in a dioecious strain of the historically self-fertile hermaprohodite C. elegans. In every case, the rate of mutational decay is substantially greater in the gonochoristic taxa than in C. elegans (∼4× greater on average). Residual heterozygosity in the ancestral controls of these MA lines introduces some complications in interpreting the results, but circumstantial evidence suggests the results are not primarily due to inbreeding depression resulting from residual segregating variation. The results suggest that natural selection operates to optimize the mutation rate in Caenorhabditis and that the strength (or efficiency) of selection differs consistently on the basis of mating system, as predicted by theory. However, context-dependent environmental and/or synergistic epistasis could also explain the results., (© 2010 The Author(s). Evolution© 2010 The Society for the Study of Evolution.)

Genetic and developmental architecture may bias the mutationally available phenotypic spectrum. Although such asymmetries in the introduction of variation may influence possible evolutionary trajectories, we lack quantitative characterization of biases in mutationally inducible phenotypic variation, their genotype-dependence, and their underlying molecular and developmental causes. Here we quantify the mutationally accessible phenotypic spectrum of the vulval developmental system using mutation accumulation (MA) lines derived from four wild isolates of the nematodes Caenorhabditis elegans and C. briggsae. The results confirm that on average, spontaneous mutations degrade developmental precision, with MA lines showing a low, yet consistently increased, proportion of developmental defects and variants. This result indicates strong purifying selection acting to maintain an invariant vulval phenotype. Both developmental system and genotype significantly bias the spectrum of mutationally inducible phenotypic variants. First, irrespective of genotype, there is a developmental bias, such that certain phenotypic variants are commonly induced by MA, while others are very rarely or never induced. Second, we found that both the degree and spectrum of mutationally accessible phenotypic variation are genotype-dependent. Overall, C. briggsae MA lines exhibited a two-fold higher decline in precision than the C. elegans MA lines. Moreover, the propensity to generate specific developmental variants depended on the genetic background. We show that such genotype-specific developmental biases are likely due to cryptic quantitative variation in activities of underlying molecular cascades. This analysis allowed us to identify the mutationally most sensitive elements of the vulval developmental system, which may indicate axes of potential evolutionary variation. Consistent with this scenario, we found that evolutionary trends in the vulval system concern the phenotypic characters that are most easily affected by mutation. This study provides an empirical assessment of developmental bias and the evolution of mutationally accessible phenotypes and supports the notion that such bias may influence the directions of evolutionary change., Competing Interests: The authors have declared that no competing interests exist.

Background: Biological phenotypes are described as "canalized" if they are robust to minor variation of environment and/or genetic background. The existence of a robust phenotype logically implies that some underlying mechanism must be variable, in the sense of "able to vary", in order to compensate for variation in the environment and/or genetic effects. Several lines of evidence lead to the conclusion that deleterious mutations predictably render morphological, developmental, and life-history traits more sensitive to small random environmental perturbations--that is, mutations de-canalize the phenotype., Methodology/principal Findings: Using conventional dye-swap microarray methodology, we compared transcript abundance in a sample of >7,000 genes between four mutation accumulation (MA) lines of the nematode Caenorhabditis elegans and the common (unmutated) ancestor. There was significantly less environmental variance in the MA lines than in the ancestor, both among replicates of the same gene and among genes., Conclusions/significance: Deleterious mutations consistently decrease the within-line component of variance in transcript abundance, which is straightforwardly interpreted as reducing the sensitivity of gene expression to small random variation in the environment. This finding is consistent with the idea that underlying variability in gene expression might be mechanistically responsible for phenotypic robustness.

The evolutionary mechanisms maintaining genetic variation in life span, particularly post-reproductive life span, are poorly understood. We characterized the effects of spontaneous mutations on life span in the rhabditid nematodes Caenorhabditis elegans and C. briggsae and standing genetic variance for life span and correlation of life span with fitness in C. briggsae. Mutations decreased mean life span, a signature of directional selection. Mutational correlations between life span and fitness were consistently positive. The average selection coefficient against new mutations in C. briggsae was approximately 2% when homozygous. The pattern of phylogeographic variation in life span is inconsistent with global mutation-selection balance (MSB), but MSB appears to hold at the local level. Standing genetic correlations in C. briggsae reflect mutational correlations at a local scale but not at a broad phylogeographic level. At the local scale, results are broadly consistent with predictions of the "mutation accumulation" hypothesis for the evolution of aging.

The genetic variation present in a species depends on the interplay between mutation, population size, and natural selection. At mutation-(purifying) selection balance (MSB) in a large population, the standing genetic variance for a trait (VG) is predicted to be proportional to the mutational variance for the trait (VM); VM is proportional to the mutation rate for the trait. The ratio VM/VG predicts the average strength of selection (S) against a new mutation. Here we compare VM and VG for lifetime reproductive success (approximately fitness) and body volume in two species of self-fertilizing rhabditid nematodes, Caenorhabditis briggsae and C. elegans, which the evidence suggests have different mutation rates. Averaged over traits, species, and populations within species, the relationship between VG and VM is quite stable, consistent with the hypothesis that differences among groups in standing variance can be explained by differences in mutational input. The average (homozygous) selection coefficient inferred from VM/VG is a few percent, smaller than typical direct estimates from mutation accumulation (MA) experiments. With one exception, the variance present in a worldwide sample of these species is similar to the variance present within a sample from a single locale. These results are consistent with specieswide MSB and uniform purifying selection, but genetic draft (hitchhiking) is a plausible alternative possibility.

Understanding the evolutionary processes responsible for shaping genetic variation within and between species requires separating the effects of mutation and selection. Differences between the patterns of genetic variation observed in nature and when mutations are allowed to accumulate in the relative absence of selection can reveal biases imposed by selection. We characterize the genetic variation at dinucleotide microsatellite repeats in four sets of 250-generation mutation accumulation (MA) lines, two in the species Caenorhabditis briggsae and two in Caenorhabditis elegans, and compare the mutational variation with the standing variation in those species. We also compare the mutational properties of microsatellites with the cumulative effects of mutations on fitness in the same lines. Integrated over the whole genome, we infer that the mutation rate of C. briggsae is about twice that of C. elegans, consistent with the cumulative mutational effects on fitness. The mutational spectrum (ratio of insertions to deletions) differs between repeat types and, in some cases, between species. The per-locus mutation rate is significantly positively correlated with the standing genetic variation at the same locus in both species, providing justification for the common practice of using the standing genetic variance as a surrogate for the mutation rate.

The evolution of canalization, the robustness of the phenotype to environmental or genetic perturbation, has attracted considerable recent interest. A key step toward understanding the evolution of any phenotype is characterizing the rate at which mutation introduces genetic variation for the trait (the mutational variance, V(M)) and the average directional effects of mutations on the trait mean (DeltaM). In this study, the mutational parameters for canalization of productivity and body volume are quantified in two sets of mutation accumulation lines of nematodes in the genus Caenorhabditis and are compared with the mutational parameters for the traits themselves. Four results emerge: (1) spontaneous mutations consistently decanalize the phenotype; (2) the mutational parameters for decanalization, V(M) (quantified as mutational heritability) and DeltaM, are of the same order of magnitude as the same parameters for the traits themselves; (3) the mutational parameters for canalization are roughly correlated with the parameters for the traits themselves across taxa; and (4) there is no evidence that residual segregating overdominant loci contribute to the decay of canalization. These results suggest that canalization is readily evolvable and that any evolutionary factor that causes mutations to accumulate will, on average, decanalize the phenotype.

A basic knowledge about mutation rates is central to our understanding of a myriad of evolutionary phenomena, including the maintenance of sex and rates of molecular evolution. Although there is substantial evidence that mutation rates vary among taxa, relatively little is known about the factors that underlie this variation at an empirical level, particularly in multicellular eukaryotes. Here we integrate several disparate lines of theoretical and empirical inquiry into a unified framework to guide future studies that are aimed at understanding why and how mutation rates evolve in multicellular species.

Mutational bias is a potentially important agent of evolution, but it is difficult to disentangle the effects of mutation from those of natural selection. Mutation-accumulation experiments, in which mutations are allowed to accumulate at very small population size, thus minimizing the efficiency of natural selection, are the best way to separate the effects of mutation from those of selection. Body size varies greatly among species of nematode in the family rhabditidae; mutational biases are both a potential cause and a consequence of that variation. We report data on the cumulative effects of mutations that affect body size in three species of rhabditid nematode that vary fivefold in adult size. Results are very consistent with previous studies of mutations underlying fitness in the same strains: two strains of Caenorhabditis briggsae decline in body size about twice as fast as two strains of C. elegans, with a concomitant higher point estimate of the genomic mutation rate; the confamilial Oscheius myriophila is intermediate. There is an overall mutational bias, such that mutations reduce size on average, but the bias appears consistent between species. The genetic correlation between mutations that affect size and those underlying fitness is large and positive, on average.

It is often assumed that the mutation rate is an evolutionarily optimized property of a taxon. The relevant mutation rate is for mutations that affect fitness, U, but the strength of selection on the mutation rate depends on the average effect of a mutation. Determination of U is complicated by the possibility that mutational effects depend on the particular environmental context in which the organism exists. It has been suggested that the effects of deleterious mutations are typically magnified in stressful environments, but most studies confound genotype with environment, so it is unclear to what extent environmental specificity of mutations is specific to a particular starting genotype. We report a study designed to separate effects of species, genotype, and environment on the degradation of fitness resulting from new mutations. Mutations accumulated for >200 generations at 20 degrees in two strains of two species of nematodes that differ in thermal sensitivity. Caenorhabditis briggsae and C. elegans have similar demography at 20 degrees, but C. elegans suffers markedly reduced fitness at 25 degrees. We find little evidence that mutational properties differ depending on environmental conditions and mutational correlations between environments are close to those expected if effects were identical in both environments.

Deleterious mutations are of fundamental importance to all aspects of organismal biology. Evolutionary geneticists have expended tremendous effort to estimate the genome-wide rate of mutation and the effects of new mutations on fitness, but the degree to which genomic mutational properties vary within and between taxa is largely unknown, particularly in multicellular organisms. Beginning with two highly inbred strains from each of three species in the nematode family Rhabditidae (Caenorhabditis briggsae, Caenorhabditis elegans, and Oscheius myriophila), we allowed mutations to accumulate in the relative absence of natural selection for 200 generations. We document significant variation in the rate of decay of fitness because of new mutations between strains and between species. Estimates of the per-generation mutational decay of fitness were very consistent within strains between assays 100 generations apart. Rate of mutational decay in fitness was positively associated with genomic mutation rate and negatively associated with average mutational effect. These results provide unambiguous experimental evidence for substantial variation in genome-wide properties of mutation both within and between species and reinforce conclusions from previous experiments that the cumulative effects on fitness of new mutations can differ markedly among related taxa.