Search

The distribution of fitness effects (DFE) of new mutations plays a central role in evolutionary biology. Estimates of the DFE from experimental Mutation Accumulation (MA) lines are compromised by the complete linkage disequilibrium (LD) between mutations in different lines. To reduce LD, we constructed two sets of recombinant inbred lines from a cross of two C. elegans MA lines. One set of lines ("RIAILs") was intercrossed for ten generations prior to ten generations of selfing; the second set of lines ("RILs") omitted the intercrossing. Residual LD in the RIAILs is much less than in the RILs, which affects the inferred DFE 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 LD are pooled into haplotypes, the inferred DFE becomes increasingly negative-skewed and leptokurtic. We conclude that the conventional wisdom - most mutations have effects near zero, 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 DFE that infer a substantial fraction of mutations with positive effects are likely confounded by LD., (© The Author(s) 2024. Published by Oxford University Press on behalf of The Genetics Society of America. All rights reserved. For commercial re-use, please contact reprints@oup.com for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site—for further information please contact journals.permissions@oup.com.)

With within-species genetic diversity estimates that span the gambit of that seen across the entirety of animals, the Caenorhabditis genus of nematodes holds unique potential to provide insights into how population size and reproductive strategies influence gene and genome organization and evolution. Our study focuses on Caenorhabditis brenneri , currently known as one of the most genetically diverse nematodes within its genus and metazoan phyla. Here, we present a high-quality gapless genome assembly and annotation for C. brenneri , revealing a common nematode chromosome arrangement characterized by gene-dense central regions and repeat rich peripheral parts. Comparison of C. brenneri with other nematodes from the 'Elegans' group revealed conserved macrosynteny but a lack of microsynteny, characterized by frequent rearrangements and low correlation iof orthogroup sizes, indicative of high rates of gene turnover. We also assessed genome organization within corresponding syntenic blocks in selfing and outcrossing species, affirming that selfing species predominantly experience loss of both genes and intergenic DNA. Comparison of gene structures revealed strikingly small number of shared introns across species, yet consistent distributions of intron number and length, regardless of population size or reproductive mode, suggesting that their evolutionary dynamics are primarily reflective of functional constraints. Our study provides valuable insights into genome evolution and expands the nematode genome resources with the highly genetically diverse C. brenneri , facilitating research into various aspects of nematode biology and evolutionary processes.

Because of pleiotropy, mutations affect the expression and inheritance of multiple traits and, together with selection, are expected to shape standing genetic covariances between traits and eventual phenotypic divergence between populations. It is therefore important to find if the M matrix, describing mutational variances of each trait and covariances between traits, varies between genotypes. We here estimate the M matrix for six locomotion behavior traits in lines of two genotypes of the nematode Caenorhabditis elegans that accumulated mutations in a nearly neutral manner for 250 generations. We find significant mutational variance along at least one phenotypic dimension of the M matrices, but neither their size nor their orientation had detectable differences between genotypes. The number of generations of mutation accumulation, or the number of MA lines measured, was likely insufficient to sample enough mutations and detect potentially small differences between the two M matrices. We then tested if the M matrices were similar to one G matrix describing the standing genetic (co)variances of a population derived by the hybridization of several genotypes, including the two measured for M, and domesticated to a lab-defined environment for 140 generations. We found that the M and G were different because the genetic covariances caused by mutational pleiotropy in the two genotypes are smaller than those caused by linkage disequilibrium in the lab population. We further show that M matrices differed in their alignment with the lab population G matrix. If generalized to other founder genotypes of the lab population, these observations indicate that selection does not shape the evolution of the M matrix for locomotion behavior in the short-term of a few tens to hundreds of generations and suggests that the hybridization of C. elegans genotypes allows selection on new phenotypic dimensions of locomotion behavior., Competing Interests: Conflicts of interest None declared., (© The Author(s) 2022. Published by Oxford University Press on behalf of the Genetics Society of America.)

Caenorhabditis elegans strains with the heat-sensitive mortal germline phenotype become progressively sterile over the course of a few tens of generations when maintained at temperatures near the upper range of C. elegans' tolerance. Mortal germline is transgenerationally heritable, and proximately under epigenetic control. Previous studies have suggested that mortal germline presents a relatively large mutational target and that mortal germline is not uncommon in natural populations of C. elegans. The mortal germline phenotype is not monolithic. Some strains exhibit a strong mortal germline phenotype, in which individuals invariably become sterile over a few generations, whereas other strains show a weaker (less penetrant) phenotype in which the onset of sterility is slower and more stochastic. We present results in which we (1) quantify the rate of mutation to the mortal germline phenotype and (2) quantify the frequency of mortal germline in a collection of 95 wild isolates. Over the course of ∼16,000 meioses, we detected one mutation to a strong mortal germline phenotype, resulting in a point estimate of the mutation rate UMrt≈ 6×10-5/genome/generation. We detected no mutations to a weak mortal germline phenotype. Six out of 95 wild isolates have a strong mortal germline phenotype, and although quantification of the weak mortal germline phenotype is inexact, the weak mortal germline phenotype is not rare in nature. We estimate a strength of selection against mutations conferring the strong mortal germline phenotype s¯≈0.1%, similar to selection against mutations affecting competitive fitness. The appreciable frequency of weak mortal germline variants in nature combined with the low mutation rate suggests that mortal germline may be maintained by balancing selection., (© The Author(s) 2022. Published by Oxford University Press on behalf of Genetics Society of America.)

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.)

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

Mechanical properties such as substrate stiffness are a ubiquitous feature of a cell's environment. Many types of animal cells exhibit canonical phenotypic plasticity when grown on substrates of differing stiffness, in vitro and in vivo. Whether such plasticity is a multivariate optimum due to hundreds of millions of years of animal evolution, or instead is a compromise between conflicting selective demands, is unknown. We addressed these questions by means of experimental evolution of populations of mouse fibroblasts propagated for approximately 90 cell generations on soft or stiff substrates. The ancestral cells grow twice as fast on stiff substrate as on soft substrate and exhibit the canonical phenotypic plasticity. Soft-selected lines derived from a genetically diverse ancestral population increased growth rate on soft substrate to the ancestral level on stiff substrate and evolved the same multivariate phenotype. The pattern of plasticity in the soft-selected lines was opposite of the ancestral pattern, suggesting that reverse plasticity underlies the observed rapid evolution. Conversely, growth rate and phenotypes did not change in selected lines derived from clonal cells. Overall, our results suggest that the changes were the result of genetic evolution and not phenotypic plasticity per se. Whole-transcriptome analysis revealed consistent differentiation between ancestral and soft-selected populations, and that both emergent phenotypes and gene expression tended to revert in the soft-selected lines. However, the selected populations appear to have achieved the same phenotypic outcome by means of at least two distinct transcriptional architectures related to mechanotransduction and proliferation., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Metabolic disorders have a large heritable component, and have increased markedly in human populations over the past few generations. Genome-wide association studies of metabolic traits typically find a substantial unexplained fraction of total heritability, suggesting an important role of spontaneous mutation. An alternative explanation is that epigenetic effects contribute significantly to the heritable variation. Here, we report a study designed to quantify the cumulative effects of spontaneous mutation on adenosine metabolism in the nematode Caenorhabditis elegans, including both the activity and concentration of two metabolic enzymes and the standing pools of their associated metabolites. The only prior studies on the effects of mutation on metabolic enzyme activity, in Drosophila melanogaster, found that total enzyme activity presents a mutational target similar to that of morphological and life-history traits. However, those studies were not designed to account for short-term heritable effects. We find that the short-term heritable variance for most traits is of similar magnitude as the variance among MA lines. This result suggests that the potential heritable effects of epigenetic variation in metabolic disease warrant additional scrutiny., (© 2020 The Authors. Evolution © 2020 The Society for the Study of Evolution.)

Environmental dependence of mutation in microbes is well-known, but most experiments have investigated contexts in which growth rate is greatly reduced below optimum. A new experiment shows mutational variability extends to contexts in which growth is near optimum., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

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.

Organismal fitness is relevant in many contexts in biology. The most meaningful experimental measure of fitness is competitive fitness, when two or more entities (e.g., genotypes) are allowed to compete directly. In theory, competitive fitness is simple to measure: an experimental population is initiated with the different types in known proportions and allowed to evolve under experimental conditions to a predefined endpoint. In practice, there are several obstacles to obtaining robust estimates of competitive fitness in multicellular organisms, the most pervasive of which is simply the time it takes to count many individuals of different types from many replicate populations. Methods by which counting can be automated in high throughput are desirable, but for automated methods to be useful, the bias and technical variance associated with the method must be (a) known, and (b) sufficiently small relative to other sources of bias and variance to make the effort worthwhile. The nematode Caenorhabditis elegans is an important model organism, and the fitness effects of genotype and environmental conditions are often of interest. We report a comparison of three experimental methods of quantifying competitive fitness, in which wild-type strains are competed against GFP-marked competitors under standard laboratory conditions. Population samples were split into three replicates and counted (1) "by eye" from a saved image, (2) from the same image using CellProfiler image analysis software, and (3) with a large particle flow cytometer (a "worm sorter"). From 720 replicate samples, neither the frequency of wild-type worms nor the among-sample variance differed significantly between the three methods. CellProfiler and the worm sorter provide at least a tenfold increase in sample handling speed with little (if any) bias or increase in variance., Competing Interests: The authors have declared that no competing interests exist.

A fundamental issue in evolutionary systems biology is understanding the relationship between the topological architecture of a biological network, such as a metabolic network, and the evolution of the network. The rate at which an element in a metabolic network accumulates genetic variation via new mutations depends on both the size of the mutational target it presents and its robustness to mutational perturbation. Quantifying the relationship between topological properties of network elements and the mutability of those elements will facilitate understanding the variation in and evolution of networks at the level of populations and higher taxa. We report an investigation into the relationship between two topological properties of 29 metabolites in the C. elegans metabolic network and the sensitivity of those metabolites to the cumulative effects of spontaneous mutation. The correlations between measures of network centrality and mutability are not statistically significant, but several trends point toward a weak positive association between network centrality and mutational sensitivity. There is a small but significant negative association between the mutational correlation of a pair of metabolites ( r M ) and the shortest path length between those metabolites. Positive association between the centrality of a metabolite and its mutational heritability is consistent with centrally-positioned metabolites presenting a larger mutational target than peripheral ones, and is inconsistent with centrality conferring mutational robustness, at least in toto . The weakness of the correlation between r M and the shortest path length between pairs of metabolites suggests that network locality is an important but not overwhelming factor governing mutational pleiotropy. These findings provide necessary background against which the effects of other evolutionary forces, most importantly natural selection, can be interpreted.

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%.

The hermaphroditic nematode Caenorhabditis elegans has been one of the primary model systems in biology since the 1970s, but only within the last two decades has this nematode also become a useful model for experimental evolution. Here, we outline the goals and major foci of experimental evolution with C. elegans and related species, such as C . briggsae and C. remanei , by discussing the principles of experimental design, and highlighting the strengths and limitations of Caenorhabditis as model systems. We then review three exemplars of Caenorhabditis experimental evolution studies, underlining representative evolution experiments that have addressed the: (1) maintenance of genetic variation; (2) role of natural selection during transitions from outcrossing to selfing, as well as the maintenance of mixed breeding modes during evolution; and (3) evolution of phenotypic plasticity and its role in adaptation to variable environments, including host-pathogen coevolution. We conclude by suggesting some future directions for which experimental evolution with Caenorhabditis would be particularly informative., (Copyright © 2017 Teotónio et al.)

Model organisms are important in many areas of chemical biology. In metabolomics, model organisms can provide excellent samples for methods development as well as the foundation of comparative phylometabolomics, which will become possible as metabolomics applications expand. Comparative studies of conserved and unique metabolic pathways will help in the annotation of metabolites as well as provide important new targets of investigation in biology and biomedicine. However, most chemical biologists are not familiar with genetics, which needs to be considered when choosing a model organism. In this review we summarize the strengths and weaknesses of several genetic systems, including natural isolates, recombinant inbred lines, and genetic mutations. We also discuss methods to detect targets of selection on the metabolome., (Copyright © 2016 The Authors. Published by Elsevier Ltd.. All rights reserved.)

A properly functioning organism must maintain metabolic homeostasis. Deleterious mutations degrade organismal function, presumably at least in part via effects on metabolic function. Here we present an initial investigation into the mutational structure of the Caenorhabditis elegans metabolome by means of a mutation accumulation experiment. We find that pool sizes of 29 metabolites vary greatly in their vulnerability to mutation, both in terms of the rate of accumulation of genetic variance (the mutational variance, VM) and the rate of change of the trait mean (the mutational bias, ΔM). Strikingly, some metabolites are much more vulnerable to mutation than any other trait previously studied in the same way. Although we cannot statistically assess the strength of mutational correlations between individual metabolites, principal component analysis provides strong evidence that some metabolite pools are genetically correlated, but also that there is substantial scope for independent evolution of different groups of metabolites. Averaged over mutation accumulation lines, PC3 is positively correlated with relative fitness, but a model in which metabolites are uncorrelated with fitness is nearly as good by Akaike's Information Criterion., (© 2016 The Author(s). Evolution © 2016 The Society for the Study of Evolution.)

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.

Background: Cellular structures such as the nucleus, Golgi, centrioles, and spindle show remarkable diversity between species, but the mechanisms that produce these variations in cell biology are not known., Results: Here we investigate the mechanisms that contribute to variations in morphology and dynamics of the mitotic spindle, which orchestrates chromosome segregation in all Eukaryotes and positions the division plane in many organisms. We use high-throughput imaging of the first division in nematodes to demonstrate that the measured effects of spontaneous mutations, combined with stabilizing selection on cell size, are sufficient to quantitatively explain both the levels of within-species variation in the spindle and its diversity over ∼100 million years of evolution. Furthermore, our finding of extensive within-species variation for the spindle demonstrates that there is not just one "wild-type" form, rather that cellular structures can exhibit a surprisingly broad diversity of naturally occurring behaviors., Conclusions: Our results argue that natural selection acts predominantly on cell size and indirectly influences the spindle through the scaling of the spindle with cell size. Previous studies have shown that the spindle also scales with cell size during early development. Thus, the scaling of the spindle with cell size controls its variation over both ontogeny and phylogeny., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

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.)

We explored the relationship between relaxed selection, oxidative stress, and spontaneous mutation in a set of mutation-accumulation (MA) lines of the nematode Caenorhabditis elegans and in their common ancestor. We measured steady-state levels of free radicals and oxidatively damaged guanosine nucleosides in the somatic tissues of five MA lines for which nuclear genome base substitution and GC-TA transversion frequencies are known. The two markers of oxidative stress are highly correlated and are elevated in the MA lines relative to the ancestor; point estimates of the per-generation rate of mutational decay (ΔM) of these measures of oxidative stress are similar to those reported for fitness-related traits. Conversely, there is no significant relationship between either marker of oxidative stress and the per-generation frequencies of base substitution or GC-TA transversion. Although these results provide no direct evidence for a causative relationship between oxidative damage and base substitution mutations, to the extent that oxidative damage may be weakly mutagenic in the germline, the case for condition-dependent mutation is advanced.

Mutation rate often increases with environmental temperature, but establishing causality is complicated. Asymmetry between physiological stress and deviation from the optimal temperature means that temperature and stress are often confounded. We allowed mutations to accumulate in two species of Caenorhabditis for approximately 100 generations at 18°C and for approximately 165 generations at 26°C; 26°C is stressful for Caenorhabditis elegans but not for Caenorhabditis briggsae. We report mutation rates at a set of microsatellite loci and estimates of the per-generation decay of fitness (ΔM(w)), the genomic mutation rate for fitness (U) and the average effect of a new mutation (E[a]), assayed at both temperatures. In C. elegans, the microsatellite mutation rate is significantly greater at 26°C than at 18°C whereas in C. briggsae there is only a slight, non-significant increase in mutation rate at 26°C, consistent with stress-dependent mutation in C. elegans. The fitness data from both species qualitatively reinforce the microsatellite results. The fitness results of C. elegans are potentially complicated by selection but also suggest temperature-dependent mutation; the difference between the two species suggests that physiological stress plays a significant role in the mutational process.

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.

Variation among lineages in the mutation process has the potential to impact diverse biological processes ranging from susceptibilities to genetic disease to the mode and tempo of molecular evolution. The combination of high-throughput DNA sequencing (HTS) with mutation-accumulation (MA) experiments has provided a powerful approach to genome-wide mutation analysis, though insights into mutational variation have been limited by the vast evolutionary distances among the few species analyzed. We performed a HTS analysis of MA lines derived from four Caenorhabditis nematode natural genotypes: C. elegans N2 and PB306 and C. briggsae HK104 and PB800. Total mutation rates did not differ among the four sets of MA lines. A mutational bias toward G:C→A:T transitions and G:C→T:A transversions was observed in all four sets of MA lines. Chromosome-specific rates were mostly stable, though there was some evidence for a slightly elevated X chromosome mutation rate in PB306. Rates were homogeneous among functional coding sequence types and across autosomal cores, arms, and tips. Mutation spectra were similar among the four MA line sets but differed significantly when compared with patterns of natural base-substitution polymorphism for 13/14 comparisons performed. Our findings show that base-substitution mutation processes in these closely related animal lineages are mostly stable but differ from natural polymorphism patterns in these two species.

Variation in rates of molecular evolution has been attributed to numerous, interrelated causes, including metabolic rate, body size, and generation time. Speculation concerning the influence of metabolic rate on rates of evolution often invokes the putative mutagenic effects of oxidative stress. To isolate the effects of oxidative stress on the germline from the effects of metabolic rate, generation time, and other factors, we allowed mutations to accumulate under relaxed selection for 125 generations in two strains of the nematode Caenorhabditis elegans, the canonical wild-type strain (N2) and a mutant strain with elevated steady-state oxidative stress (mev-1). Contrary to our expectation, the mutational decline in fitness did not differ between N2 and mev-1. This result suggests that the mutagenic effects of oxidative stress in C. elegans are minor relative to the effects of other types of mutations, such as errors during DNA replication. However, mev-1 MA lines did go extinct more frequently than wild-type lines; some possible explanations for the difference in extinction rate are discussed.

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.

Mitochondrial DNA (mtDNA) mutations underlie a variety of human genetic disorders and are associated with the aging process. mtDNA polymorphisms are widely used in a variety of evolutionary applications. Although mtDNA mutation spectra are known to differ between distantly related model organisms, the extent to which mtDNA mutation processes vary between more closely related species and within species remains enigmatic. We analyzed mtDNA divergence in two sets of 250-generation Caenorhabditis briggsae mutation-accumulation (MA) lines, each derived from a different natural isolate progenitor: strain HK104 from Okayama, Japan, and strain PB800 from Ohio, United States. Both sets of C. briggsae MA lines accumulated numerous large heteroplasmic mtDNA deletions, whereas only one similar event was observed in a previous analysis of Caenorhabditis elegans MA line mtDNA. Homopolymer length change mutations were frequent in both sets of C. briggsae MA lines and occurred in both intergenic and protein-coding gene regions. The spectrum of C. briggsae mtDNA base substitution mutations differed from the spectrum previously observed in C. elegans. In C. briggsae, the HK104 MA lines experienced many different base substitution types, whereas the PB800 lines displayed only C:G --> T:A transitions, although the difference was not significant. Over half of the mtDNA base substitutions detected in the C. briggsae MA lines were in a heteroplasmic state, whereas all those previously characterized in C. elegans MA line mtDNA were fixed changes, indicating a narrower mtDNA bottleneck in C. elegans as compared with C. briggsae. Our results show that C. briggsae mtDNA is highly susceptible to large deletions and that the mitochondrial mutation process varies between Caenorhabditis nematode species.

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.

Knowledge of mutation processes is central to understanding virtually all evolutionary phenomena and the underlying nature of genetic disorders and cancers. However, the limitations of standard molecular mutation detection methods have historically precluded a genome-wide understanding of mutation rates and spectra in the nuclear genomes of multicellular organisms. We applied two high-throughput DNA sequencing technologies to identify and characterize hundreds of spontaneously arising base-substitution mutations in 10 Caenorhabditis elegans mutation-accumulation (MA)-line nuclear genomes. C. elegans mutation rate estimates were similar to previous calculations based on smaller numbers of mutations. Mutations were distributed uniformly within and among chromosomes and were not associated with recombination rate variation in the MA lines, suggesting that intragenomic variation in genetic hitchhiking and/or background selection are primarily responsible for the chromosomal distribution patterns of polymorphic nucleotides in C. elegans natural populations. A strong mutational bias from G/C to A/T nucleotides was detected in the MA lines, implicating oxidative DNA damage as a major endogenous mutagenic force in C. elegans. The observed mutational bias also suggests that the C. elegans nuclear genome cannot be at equilibrium because of mutation alone. Transversions dominate the spectrum of spontaneous mutations observed here, whereas transitions dominate patterns of allegedly neutral polymorphism in natural populations of C. elegans and many other animal species; this observation challenges the assumption that natural patterns of molecular variation in noncoding regions of the nuclear genome accurately reflect underlying mutation processes.

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.

Laboratory natural selection and artificial selection are vital tools for addressing specific questions about evolutionary patterns of variation. Laboratory natural selection can illuminate whether a putative selective agent is capable of generating long-term, sustained changes in individual traits and suites of traits. Artificial selection is the essential tool for understanding the general evolvability of traits and the extent to which genetic correlations constrain evolution. We review the contexts in which each type of experiment seems capable of offering key insights into important evolutionary issues. We also discuss theoretical and methodological considerations that play critical roles in designing selection experiments that are relevant to evolutionary patterns of trait variation. In particular, we focus on the critical role of selection intensity and the consequences of experiments with different intensities. While selection experiments are not practical in many cases, sophisticated selection experiments-designed with careful consideration of the theory of selection-should be taken beyond model organisms and used in well-chosen natural systems to understand natural patterns of variation.

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.

The exceptional diversity of insects is often attributed to the effects of specialized relationships between insects and their hosts. Parasite-host interactions are influenced by current natural selection and dispersal, in addition to historical effects that may include past selection, vicariance, and random genetic drift. Both current and historical events can lead to reduced fitness on some hosts. If trade-offs in fitness on alternate hosts are common, adaptation to one host can prevent adaptation to another, giving rise to genetic differentiation among host-associated lineages. Previous studies of Diaeretiella rapae (Hymenoptera: Aphidiidae), a parasitoid of aphids, have revealed additive genetic differences in performance between populations that parasitize different aphid host species. To determine whether D. rapae populations collected from different aphid hosts have diverged into genetically independent lineages, we constructed a haplotype network based on sequence variation in mitochondrial DNA (mtDNA). We used single strand conformation polymorphism (SSCP) analysis to examine 2041 base pairs of mtDNA and to identify nucleotide sequences of 42 unique SSCP haplotypes. We found no association between mtDNA haplotypes and host species in either the ancestral range (Europe, Mediterranean region, Middle East, Asia) or part of the introduced range (western North America). Haplotypes likely to be ancestral were geographically widespread and found on both hosts, suggesting that the ability to use both hosts evolved prior to the diversification of the mtDNA. Ongoing gene flow appears to prevent the formation of host races.

Explaining the repeated evolution of similar sets of traits under similar environmental conditions is an important issue in evolutionary biology. The extreme alternative classes of explanations for correlated suites of traits are optimal adaptation and genetic constraint resulting from pleiotropy. Adaptive explanations presume that individual traits are free to evolve to their local optima and that convergent evolution represents particularly adaptive combinations of traits. Alternatively, if pleiotropy is strong and difficult to break, strong selection on one or a few particularly important characters would be expected to result in consistent correlated evolution of associated traits. If pleiotropy is common, we predict that the pattern of divergence among populations will consistently reflect the within-population genetic architecture. To test the idea that the multivariate life-history phenotype is largely a byproduct of strong selection on body size, we imposed divergent artificial selection on size at maturity upon two populations of the cladoceran Daphnia pulicaria, chosen on the basis of their extreme divergence in body size. Overall, the trajectory of divergence between the two natural populations did not differ from that predicted by the genetic architecture within each population. However, the pattern of correlated responses suggested the presence of strong pleiotropic constraints only for adult body size and not for other life-history traits. One trait, offspring size, appears to have evolved in a way different from that expected from the within-population genetic architecture and may be under stabilizing selection.

Modern population genetics underwent a major paradigm shift during the last decade of the 20th century with the discovery that thousands of genes of known function and position in a genome can be analyzed simultaneously in a single individual. The impact of this technology on insect population genetics is potentially profound. Sampling distributions of genetic statistics can now be derived from many individual loci or among many segregating sites within a gene. Inferences regarding random mating, gene flow, effective population sizes, disequilibrium, and relatedness among populations can now be based on patterns of variation at many loci. More importantly, genome-wide sampling enables population geneticists to distinguish effects that act on the whole genome from those that act on individual loci or nucleotides. We introduce the term "population genomics" to describe the process of simultaneous sampling of numerous variable loci within a genome and the inference of locus-specific effects from the sample distributions. The four critical assumptions implicit in the population genomics approach are explained in detail. Studies adopting this paradigm are reviewed, and the steps necessary to complete a population genomics study are outlined.

Group selection has historically been an important and controversial subject in evolutionary biology. There is now a compelling body of evidence, both theoretical and experimental, that group selection not only can be effective, but can be effective in situations when individual selection is not. However, experiments in which true population-level traits have been shown to evolve in response to group selection are currently limited to two species of flour beetle in the genus Tribolium and RNA viruses. Here we report the results of an experiment wherein we imposed group selection via differential extinction for increased and decreased population size at 6-week intervals, a true population-level trait, in the poeciliid fish Heterandria formosa. In contrast to most other group selection experiments, we observed no evolutionary response after six rounds of group selection in either the up- or down-selected lines. Populational heritability for population size was low, if not actually negative. Our results suggest that group selection via differential extinction may be effective only if population sizes are very small and/or migration rates are low.

Tradeoffs in performance or fitness across environments have important implications regarding the nature of evolutionary constraints. It remains controversial whether tradeoffs such as these reflect genetic correlations that are genuine evolutionary constraints. However, if such long-term genetic constraints do exist, they must be due to underlying pleiotropy such that alleles that confer high performance in one environment invariably confer low performance in another. The distribution of genetic correlations within and among populations can provide insight about the existence of such pleiotropic tradeoffs. The long-term association of certain teleost fish taxa with particular abiotic environments suggests that tradeoffs in performance across environments have constrained the geographic distribution of those taxa. Here we report the results of an experiment in which we artificially selected on acute heat- and cold-stress tolerance in two stocks of the poeciliid fish Heterandria formosa from source populations with different thermal histories. Unexpectedly, we observed no direct responses to selection. Under certain conditions, fish from the different source populations differed significantly in cold tolerance, but not in heat tolerance. The results suggest there are no strong pleiotropic tradeoffs between heat- and cold-stress tolerance in these populations.

Variation among loci in the distribution of allele frequencies among subpopulations is well known; how to tell when the variation exceeds that expected when all loci are subject to uniform evolutionary processes is not well known. If locus-specific effects are important, the ability to detect those effects should vary with the level of gene flow. Populations with low gene flow should exhibit greater variation among loci in Fst than populations with high gene flow, because gene flow acts to homogenize allele frequencies among subpopulations. Here I use Lewontin and Krakauer's k statistic to describe the variance among allozyme loci in 102 published data sets from fishes. As originally proposed, k >> 2 was considered evidence that the variation in Fst among loci is greater than expected from neutral evolution. Although that interpretation is invalid, large differences in k in different populations suggest that locus-specific forces may be important in shaping genetic diversity. In these data, k is not greater for populations with expected low levels of gene flow than for populations with expected high levels of gene flow. There is thus no evidence that locus-specific forces are of general importance in shaping the distribution of allele frequencies at enzyme loci among populations of fishes.

The phylogeography of the freshwater fish fauna of the southeastern United States has almost achieved paradigm status in evolutionary biology (Avise 1992), and the major geographic features responsible for shaping species distributions are well-characterized. Nevertheless, variation among species in distributions of allele or haplotype frequencies suggests that species-specific processes (e.g., migration) may also play a role in establishing those distributions. There has also been relatively little investigation into how population structure may differ among subregions in the Southeast, for example, on the Florida peninsula versus the U.S. mainland to the northwest and/or northeast. The geology of the peninsula is such that both physical and biotic fluctuations may have been (and still be) particularly important in establishing the population structure of freshwater taxa. This possibility leads to two interesting questions in population genetics. (1) Does gene flow in freshwater species of the region better approximate a one- or two-dimensional pattern? (2) Are populations on the peninsula farther from migration-genetic drift equilibrium than their counterparts on the mainland? These questions were addressed by examining the population strucuture of a livebearing fish, Heterandria formosa; several features of the biology of the species make it particularly likely that recent gene flow has been important in its evolution. I surveyed electrophoretic variation in 34 populations distributed throughout the species range. The phylogeographic patterns observed are in general concordance with those found in other species, although with some differences. A two-dimensional hypothesis of gene flow on the Florida peninsula better explains the data than does a one-dimensional one. There is no evidence that populations on the peninsula are farther from migration-drift equilibrium than those to the northwest. Populations in the northeast have lower genetic diversity than those to the south and west and show no isolation by distance; those results are consistent with a recent range expansion into the northeast, although smaller historical effective population sizes could also explain the pattern., (© 1998 The Society for the Study of Evolution.)