Direct estimation of per nucleotide and genomic deleterious mutation rates in Drosophila

Mitotic evolution of Plasmodium falciparum shows a stable core genome but recombination in antigen families

Malaria parasites elude eradication attempts both within the human host and across nations. At the individual level, parasites evade the host immune responses through antigenic variation. At the global level, parasites escape drug pressure through single nucleotide variants and gene copy amplification events conferring drug resistance. Despite their importance to global health, the rates at which these genomic alterations emerge have not been determined. We studied the complete genomes of different Plasmodium falciparum clones that had been propagated asexually over one year in the presence and absence of drug pressure. A combination of whole-genome microarray analysis and next-generation deep resequencing (totaling 14 terabases) revealed a stable core genome with only 38 novel single nucleotide variants appearing in seventeen evolved clones (avg. 5.4 per clone). In clones exposed to atovaquone, we found cytochrome b mutations as well as an amplification event encompassing the P. falciparum multidrug resistance associated protein (mrp1) on chromosome 1. We observed 18 large-scale (>1 kb on average) deletions of telomere-proximal regions encoding multigene families, involved in immune evasion (9.5×10⁻⁶ structural variants per base pair per generation). Six of these deletions were associated with chromosomal crossovers generated during mitosis. We found only minor differences in rates between genetically distinct strains and between parasites cultured in the presence or absence of drug. Using these derived mutation rates for P. falciparum (1.0–9.7×10⁻⁹ mutations per base pair per generation), we can now model the frequency at which drug or immune resistance alleles will emerge under a well-defined set of assumptions. Further, the detection of mitotic recombination events in var gene families illustrates how multigene families can arise and change over time in P. falciparum. These results will help improve our understanding of how P. falciparum evolves to evade control efforts within both the individual hosts and large populations.

Malaria is one of the six diseases that together are responsible for 90% of all infectious disease deaths throughout the world. The five species of Plasmodium that cause human malaria take over 655,000 lives each year. Parasites evade the immune response through antigenic variation and develop resistance to anti-malarial drugs through genetic changes in either the drug target or genes conferring resistance. We used whole-genome sequencing and microarray techniques to study evolution in P. falciparum parasites propagated in vitro for up to 180 generations. We determined the mutation rate and found that the core genome of a single clone is stable, while the subtelomeric regions are prone to acquire structural variants. These changes occur mainly in multigene families involved in immune evasion. Our findings indicate that the parasite specifically increases the sequence variability in multigene families through mitotic recombination. This high plasticity of the parasite genome suggests that multiple haplotypes will be present in a natural infection initiated by a single parasite.

Genomic background and generation time influence deleterious mutation rates in Daphnia

Understanding how genetic variation is generated and how selection shapes mutation rates over evolutionary time requires knowledge of the factors influencing mutation and its effects on quantitative traits. We explore the impact of two factors, genomic background and generation time, on deleterious mutation in Daphnia pulicaria, a cyclically parthenogenic aquatic microcrustacean, using parallel mutation-accumulation experiments. The deleterious mutational properties of life-history characters for individuals from two different populations, and for individuals maintained at two different generation times, were quantified and compared. Mutational properties varied between populations, especially for clutch size, suggesting that genomic background influences mutational properties for some characters. Generation time was found to have a greater effect on mutational properties, with higher per-generation deleterious mutation rates in lines with longer generation times. These results suggest that differences in genetic architecture among populations and species may be explained in part by demographic features that significantly influence generation time and therefore the rate of mutation.

Non-random mutation: the evolution of targeted hypermutation and hypomutation

A widely accepted tenet of evolutionary biology is that spontaneous mutations occur randomly with regard to their fitness effect. However, since the mutation rate varies along a genome and this variation can be subject to selection, organisms might evolve lower mutation rates at loci where mutations are most deleterious or increased rates where mutations are most needed. In fact, mechanisms of targeted hypermutation are known in organisms ranging from bacteria to humans. Here we review the main forces driving the evolution of local mutation rates and identify the main limiting factors. Both targeted hyper- and hypomutation can evolve, although the former is restricted to loci under very frequent positive selection and the latter is severely limited by genetic drift. Nevertheless, we show how an association of repair with transcription or chromatin-associated proteins could overcome the drift limit and lead to non-random hypomutation along the genome in most organisms.

The effect of spontaneous mutations on competitive ability

Understanding the impact of spontaneous mutations on fitness has many theoretical and practical applications in biology. Although mutational effects on individual morphological or life-history characters have been measured in several classic genetic model systems, there are few estimates of the rate of decline due to mutation for complex fitness traits. Here, we estimate the effects of mutation on competitive ability, an important complex fitness trait, in a model system for ecological and evolutionary genomics, Daphnia. Competition assays were performed to compare fitness between mutation-accumulation (MA) lines and control lines from eight different genotypes from two populations of Daphnia pulicaria after 30 and 65 generations of mutation accumulation. Our results show a fitness decline among MA lines relative to controls as expected, but highlight the influence of genomic background on this effect. In addition, in some assays, MA lines outperform controls providing insight into the frequency of beneficial mutations.

Population genomics of the Wolbachia endosymbiont in Drosophila melanogaster

Wolbachia are maternally inherited symbiotic bacteria, commonly found in arthropods, which are able to manipulate the reproduction of their host in order to maximise their transmission. The evolutionary history of endosymbionts like Wolbachia can be revealed by integrating information on infection status in natural populations with patterns of sequence variation in Wolbachia and host mitochondrial genomes. Here we use whole-genome resequencing data from 290 lines of Drosophila melanogaster from North America, Europe, and Africa to predict Wolbachia infection status, estimate relative cytoplasmic genome copy number, and reconstruct Wolbachia and mitochondrial genome sequences. Overall, 63% of Drosophila strains were predicted to be infected with Wolbachia by our in silico analysis pipeline, which shows 99% concordance with infection status determined by diagnostic PCR. Complete Wolbachia and mitochondrial genomes show congruent phylogenies, consistent with strict vertical transmission through the maternal cytoplasm and imperfect transmission of Wolbachia. Bayesian phylogenetic analysis reveals that the most recent common ancestor of all Wolbachia and mitochondrial genomes in D. melanogaster dates to around 8,000 years ago. We find evidence for a recent global replacement of ancestral Wolbachia and mtDNA lineages, but our data suggest that the derived wMel lineage arose several thousand years ago, not in the 20th century as previously proposed. Our data also provide evidence that this global replacement event is incomplete and is likely to be one of several similar incomplete replacement events that have occurred since the out-of-Africa migration that allowed D. melanogaster to colonize worldwide habitats. This study provides a complete genomic analysis of the evolutionary mode and temporal dynamics of the D. melanogaster–Wolbachia symbiosis, as well as important resources for further analyses of the impact of Wolbachia on host biology.

Host–microbe interactions play important roles in the physiology, development, and ecology of many organisms. Studying how hosts and their microbial symbionts evolve together over time is crucial for understanding the impact that microbes have on host biology. With the advent of high-throughput sequencing technologies, it is now possible to obtain complete genomic information for hosts and their associated microbes. Here we use whole-genome sequences from ∼300 strains of the fruitfly Drosophila melanogaster to reveal the evolutionary history of this model species and its intracellular bacterial symbiont Wolbachia. The major findings of this study are that Wolbachia in D. melanogaster is inherited strictly through the egg with no evidence of horizontal transfer from other species, that the genealogies of Wolbachia and mitochondrial genomes are virtually the same, and that both Wolbachia and mitochondrial genomes show evidence for a recent incomplete global replacement event, which has left remnant lineages in North America, Europe, and Africa. We also use the fact that Wolbachia and mitochondrial genomes have the same genealogy to estimate the rate of molecular evolution for Wolbachia, which allows us to put dates on key events in the history of this important host–microbe model system.

Properties of spontaneous mutational variance and covariance for wing size and shape in Drosophila melanogaster

We estimated mutational variance-covariance matrices, M, for wing shape and size in two genotypes of Drosophila melanogaster after 192 generations of mutation accumulation. We characterized 21 potentially independent aspects of wing shape and size using geometric morphometrics, and analyzed the data using a likelihood-based factor-analytic approach. We implement a previously unused analysis that describes those directions with the greatest difference in evolvability between pairs of matrices. There are significant mutational effects on 19 of 21 possible aspects of wing form, consistent with the high dimensionality of standing genetic variation for wing shape previously identified in D. melanogaster. Mutations have partially recessive effects, consistent with average dominance around 0.25. Sex-specific matrices are relatively similar, although male-specific matrices are slightly larger, as expected due to dosage compensation on the X chromosome. Genotype-specific matrices are quite different. Matrices may differ both because of sampling error based on small samples of mutations with large phenotypic effects, and because of the mutational properties of the genotypes. Genotypic differences are likely to be involved, as the two genotypes have different molecular mutation rates and properties.

Mutation Load: The Fitness of Individuals in Populations Where Deleterious Alleles Are Abundant

Many multicellular eukaryotes have reasonably high per-generation mutation rates. Consequently, most populations harbor an abundance of segregating deleterious alleles. These alleles, most of which are of small effect individually, collectively can reduce substantially the fitness of individuals relative to what it would be otherwise; this is mutation load. Mutation load can be lessened by any factor that causes more mutations to be removed per selective death, such as inbreeding, synergistic epistasis, population structure, or harsh environments. The ecological effects of load are not clear-cut because some conditions (such as selection early in life, sexual selection, reproductive compensation, and intraspecific competition) reduce the effects of load on population size and persistence, but other conditions (such as interspecific competition and load on resource use efficiency) can cause small amounts of load to have strong effects on the population, even extinction. We suggest a series of studies to improve our understanding of the effects of mutation load.

Male-biased fitness effects of spontaneous mutations in Drosophila melanogaster

In populations with males and females, sexual selection may often represent a major component of overall selection. Sexual selection could act to eliminate deleterious alleles in concert with other forms of selection, thereby improving the fitness of sexual populations. Alternatively, the divergent reproductive strategies of the sexes could promote the maintenance of sexually antagonistic variation, causing sexual populations to be less fit. The net impact of sexual selection on fitness is not well understood, due in part to limited data on the sex-specific effects of spontaneous mutations on total fitness. Using a set of mutation accumulation lines of Drosophila melanogaster, we found that mutations were deleterious in both sexes and had larger effects on fitness in males than in females. This pattern is expected to reduce the mutation load of sexual females and promote the maintenance of sexual reproduction.

Modeling the effect of changing selective pressures on polymorphism and divergence

The most common models of sequence evolution used to make inferences about adaptation rely on the assumption that selective pressures at a site remain constant through time. Instead, one might plausibly imagine that a change in the environment renders an allele beneficial and that when it fixes, the site is now constrained-until another change in the environment occurs that affects the selective pressures at that site. With this view in mind, we introduce a simple dynamic model for the evolution of coding regions, in which non-synonymous sites alternate between being fixed for the favored allele and being neutral with respect to other alleles. We use the pruning algorithm to derive closed forms for observable patterns of polymorphism and divergence in terms of the model parameters. Using our model, estimates of the fraction of beneficial substitutions α would remain similar to those obtained from existing approaches. In this framework, however, it becomes natural to ask how often adaptive substitutions originate from previously constrained or previously neutral sites, i.e., about the source of adaptive substitutions. We show that counts of coding sites that are both polymorphic in a sample from one species and divergent between two others carry information about this parameter. We also extend the basic model to include the effects of weakly deleterious mutations and discuss the importance of assumptions about the distribution of deleterious mutations among constrained non-synonymous sites. Finally, we derive a likelihood function for the parameters and apply it to a toy example, variation data for coding regions from chromosome 2 of the Drosophila melanogaster subgroup. This modeling work underscores how restrictive assumptions about adaptation have been to date, and how further work in this area will help to reveal unexplored and yet basic characteristics of adaptation.

Estimating divergence dates and substitution rates in the Drosophila phylogeny

An absolute timescale for evolution is essential if we are to associate evolutionary phenomena, such as adaptation or speciation, with potential causes, such as geological activity or climatic change. Timescales in most phylogenetic studies use geologically dated fossils or phylogeographic events as calibration points, but more recently, it has also become possible to use experimentally derived estimates of the mutation rate as a proxy for substitution rates. The large radiation of drosophilid taxa endemic to the Hawaiian islands has provided multiple calibration points for the Drosophila phylogeny, thanks to the "conveyor belt" process by which this archipelago forms and is colonized by species. However, published date estimates for key nodes in the Drosophila phylogeny vary widely, and many are based on simplistic models of colonization and coalescence or on estimates of island age that are not current. In this study, we use new sequence data from seven species of Hawaiian Drosophila to examine a range of explicit coalescent models and estimate substitution rates. We use these rates, along with a published experimentally determined mutation rate, to date key events in drosophilid evolution. Surprisingly, our estimate for the date for the most recent common ancestor of the genus Drosophila based on mutation rate (25-40 Ma) is closer to being compatible with independent fossil-derived dates (20-50 Ma) than are most of the Hawaiian-calibration models and also has smaller uncertainty. We find that Hawaiian-calibrated dates are extremely sensitive to model choice and give rise to point estimates that range between 26 and 192 Ma, depending on the details of the model. Potential problems with the Hawaiian calibration may arise from systematic variation in the molecular clock due to the long generation time of Hawaiian Drosophila compared with other Drosophila and/or uncertainty in linking island formation dates with colonization dates. As either source of error will bias estimates of divergence time, we suggest mutation rate estimates be used until better models are available.

Estimate of the spontaneous mutation rate in Chlamydomonas reinhardtii

The nature of spontaneous mutations, including their rate, distribution across the genome, and fitness consequences, is of central importance to biology. However, the low rate of mutation has made it difficult to study spontaneous mutagenesis, and few studies have directly addressed these questions. Here, we present a direct estimate of the mutation rate and a description of the properties of new spontaneous mutations in the unicellular green alga Chlamydomonas reinhardtii. We conducted a mutation accumulation experiment for ∼350 generations followed by whole-genome resequencing of two replicate lines. Our analysis identified a total of 14 mutations, including 5 short indels and 9 single base mutations, and no evidence of larger structural mutations. From this, we estimate a total mutation rate of 3.23 × 10(-10)/site/generation (95% C.I. 1.82 × 10(-10) to 5.23 × 10(-10)) and a single base mutation rate of 2.08 × 10(-10)/site/generation (95% C.I., 1.09 × 10(-10) to 3.74 × 10(-10)). We observed no mutations from A/T → G/C, suggesting a strong mutational bias toward A/T, although paradoxically, the GC content of the C. reinhardtii genome is very high. Our estimate is only the second direct estimate of the mutation rate from plants and among the lowest spontaneous base-substitution rates known in eukaryotes.

Whole genome sequencing of mutation accumulation lines reveals a low mutation rate in the social amoeba Dictyostelium discoideum

Spontaneous mutations play a central role in evolution. Despite their importance, mutation rates are some of the most elusive parameters to measure in evolutionary biology. The combination of mutation accumulation (MA) experiments and whole-genome sequencing now makes it possible to estimate mutation rates by directly observing new mutations at the molecular level across the whole genome. We performed an MA experiment with the social amoeba Dictyostelium discoideum and sequenced the genomes of three randomly chosen lines using high-throughput sequencing to estimate the spontaneous mutation rate in this model organism. The mitochondrial mutation rate of 6.76×10(-9), with a Poisson confidence interval of 4.1×10(-9) - 9.5×10(-9), per nucleotide per generation is slightly lower than estimates for other taxa. The mutation rate estimate for the nuclear DNA of 2.9×10(-11), with a Poisson confidence interval ranging from 7.4×10(-13) to 1.6×10(-10), is the lowest reported for any eukaryote. These results are consistent with low microsatellite mutation rates previously observed in D. discoideum and low levels of genetic variation observed in wild D. discoideum populations. In addition, D. discoideum has been shown to be quite resistant to DNA damage, which suggests an efficient DNA-repair mechanism that could be an adaptation to life in soil and frequent exposure to intracellular and extracellular mutagenic compounds. The social aspect of the life cycle of D. discoideum and a large portion of the genome under relaxed selection during vegetative growth could also select for a low mutation rate. This hypothesis is supported by a significantly lower mutation rate per cell division in multicellular eukaryotes compared with unicellular eukaryotes.

Revising the human mutation rate: implications for understanding human evolution

It is now possible to make direct measurements of the mutation rate in modern humans using next-generation sequencing. These measurements reveal a value that is approximately half of that previously derived from fossil calibration, and this has implications for our understanding of demographic events in human evolution and other aspects of population genetics. Here, we discuss the implications of a lower-than-expected mutation rate in relation to the timescale of human evolution.

Experimental evolution

Experimental evolution is the study of evolutionary processes occurring in experimental populations in response to conditions imposed by the experimenter. This research approach is increasingly used to study adaptation, estimate evolutionary parameters, and test diverse evolutionary hypotheses. Long applied in vaccine development, experimental evolution also finds new applications in biotechnology. Recent technological developments provide a path towards detailed understanding of the genomic and molecular basis of experimental evolutionary change, while new findings raise new questions that can be addressed with this approach. However, experimental evolution has important limitations, and the interpretation of results is subject to caveats resulting from small population sizes, limited timescales, the simplified nature of laboratory environments, and, in some cases, the potential to misinterpret the selective forces and other processes at work.

Phylogeny, integration and expression of sigma virus-like genes in Drosophila

The recent and surprising discovery of widespread NIRVs (non-retroviral integrated RNA viruses) has highlighted the importance of genomic interactions between non-retroviral RNA viruses and their eukaryotic hosts. Among the viruses with integrated representatives are the rhabdoviruses, a family of negative sense single-stranded RNA viruses. We identify sigma virus-like NIRVs of Drosophila spp. that represent unique cases where NIRVs are closely related to exogenous RNA viruses in a model host organism. We have used a combination of bioinformatics and laboratory methods to explore the evolution and expression of sigma virus-like NIRVs in Drosophila. Recent integrations in Drosophila provide a promising experimental system to study functionality of NIRVs. Moreover, the genomic architecture of recent NIRVs provides an unusual evolutionary window on the integration mechanism. For example, we found that a sigma virus-like polymerase associated protein (P) gene appears to have been integrated by template switching of the blastopia-like LTR retrotransposon. The sigma virus P-like NIRV is present in multiple retroelement fused open reading frames on the X and 3R chromosomes of Drosophila yakuba - the X-linked copy is transcribed to produce an RNA product in adult flies. We present the first account of sigma virus-like NIRVs and the first example of NIRV expression in a model animal system, and therefore provide a platform for further study of the possible functions of NIRVs in animal hosts.

Limits to the rate of adaptive substitution in sexual populations

In large populations, many beneficial mutations may be simultaneously available and may compete with one another, slowing adaptation. By finding the probability of fixation of a favorable allele in a simple model of a haploid sexual population, we find limits to the rate of adaptive substitution, Λ, that depend on simple parameter combinations. When variance in fitness is low and linkage is loose, the baseline rate of substitution is Λ₀ = 2NU , where N is the population size, U is the rate of beneficial mutations per genome, and is their mean selective advantage. Heritable variance v in log fitness due to unlinked loci reduces Λ by e⁻⁴(v) under polygamy and e⁻⁸ (v) under monogamy. With a linear genetic map of length R Morgans, interference is yet stronger. We use a scaling argument to show that the density of adaptive substitutions depends on s, N, U, and R only through the baseline density: Λ/R = F (Λ₀/R). Under the approximation that the interference due to different sweeps adds up, we show that Λ/R ~(Λ₀/R) / (1 +2Λ₉/R) , implying that interference prevents the rate of adaptive substitution from exceeding one per centimorgan per 200 generations. Simulations and numerical calculations confirm the scaling argument and confirm the additive approximation for Λ₀/R ~ 1; for higher Λ₀/R , the rate of adaptation grows above R/2, but only very slowly. We also consider the effect of sweeps on neutral diversity and show that, while even occasional sweeps can greatly reduce neutral diversity, this effect saturates as sweeps become more common-diversity can be maintained even in populations experiencing very strong interference. Our results indicate that for some organisms the rate of adaptive substitution may be primarily recombination-limited, depending only weakly on the mutation supply and the strength of selection.

Collapse

Deficiency screening for genomic regions with effects on environmental sensitivity of the sensory bristles of Drosophila melanogaster

Environmental canalization is defined as a reduction in the effect of external environmental perturbations on a phenotype, while phenotypic plasticity is defined as the production of different phenotypes in alternative environments. These terms describe different aspects of the same phenomenon, that is, the sensitivity of the phenotype to the environment. Genetic regulation of the environmental sensitivity has been a central topic in the field of evolutionary biology. In this study, we performed deficiency screening to detect genomic regions with effects on the environmental sensitivity of Drosophila melanogaster sensory bristles. We used a collection of isogenic deficiency strains established by the DrosDel Project for screening. We screened 423 genomic deficiencies that encompassed approximately 63.6% of the entire D. melanogaster genome. We identified 29 genomic deficiencies showing significant effects on environmental sensitivity, suggesting that multiple genomic regions may influence phenotypic variation. We also found significant correlations among the effects of deficiencies on environmental sensitivity for different bristle traits, suggesting that the same genetic mechanism can regulate environmental sensitivity of multiple traits. Current high-resolution mapping will facilitate the examination of individual candidate genes using mutations or RNAi approaches in future studies.

Estimating the strength of selective sweeps from deep population diversity data

Selective sweeps are typically associated with a local reduction of genetic diversity around the adaptive site. However, selective sweeps can also quickly carry neutral mutations to observable population frequencies if they arise early in a sweep and hitchhike with the adaptive allele. We show that the interplay between mutation and exponential amplification through hitchhiking results in a characteristic frequency spectrum of the resulting novel haplotype variation that depends only on the ratio of the mutation rate and the selection coefficient of the sweep. On the basis of this result, we develop an estimator for the selection coefficient driving a sweep. Since this estimator utilizes the novel variation arising from mutations during a sweep, it does not rely on preexisting variation and can also be applied to loci that lack recombination. Compared with standard approaches that infer selection coefficients from the size of dips in genetic diversity around the adaptive site, our estimator requires much shorter sequences but sampled at high population depth to capture low-frequency variants; given such data, it consistently outperforms standard approaches. We investigate analytically and numerically how the accuracy of our estimator is affected by the decay of the sweep pattern over time as a consequence of random genetic drift and discuss potential effects of recombination, soft sweeps, and demography. As an example for its use, we apply our estimator to deep sequencing data from human immunodeficiency virus populations.

Whole genome sequencing for quantifying germline mutation frequency in humans and model species: cautious optimism

Factors affecting the type and frequency of germline mutations in animals are of significant interest from health and toxicology perspectives. However, studies in this field have been limited by the use of markers with low detection power or uncertain relevance to phenotype. Whole genome sequencing (WGS) is now a potential option to directly determine germline mutation type and frequency in family groups at all loci simultaneously. Medical studies have already capitalized on WGS to identify novel mutations in human families for clinical purposes, such as identifying candidate genes contributing to inherited conditions. However, WGS has not yet been used in any studies of vertebrates that aim to quantify changes in germline mutation frequency as a result of environmental factors. WGS is a promising tool for detecting mutation induction, but it is currently limited by several technical challenges. Perhaps the most pressing issue is sequencing error rates that are currently high in comparison to the intergenerational mutation frequency. Different platforms and depths of coverage currently result in a range of 10-10(3) false positives for every true mutation. In addition, the cost of WGS is still relatively high, particularly when comparing mutation frequencies among treatment groups with even moderate sample sizes. Despite these challenges, WGS offers the potential for unprecedented insight into germline mutation processes. Refinement of available tools and emergence of new technologies may be able to provide the improved accuracy and reduced costs necessary to make WGS viable in germline mutation studies in the very near future. To streamline studies, researchers may use multiple family triads per treatment group and sequence a targeted (reduced) portion of each genome with high (20-40 ×) depth of coverage. We are optimistic about the application of WGS for quantifying germline mutations, but caution researchers regarding the resource-intensive nature of the work using existing technology.

Mutation accumulation in a selfing population: consequences of different mutation rates between selfers and outcrossers

Currently existing theories predict that because deleterious mutations accumulate at a higher rate, selfing populations suffer from more intense genetic degradation relative to outcrossing populations. This prediction may not always be true when we consider a potential difference in deleterious mutation rate between selfers and outcrossers. By analyzing the evolutionary stability of selfing and outcrossing in an infinite population, we found that the genome-wide deleterious mutation rate would be lower in selfing than in outcrossing organisms. When this difference in mutation rate was included in simulations, we found that in a small population, mutations accumulated more slowly under selfing rather than outcrossing. This result suggests that under frequent and intense bottlenecks, a selfing population may have a lower risk of genetic extinction than an outcrossing population.

Evidence for elevated mutation rates in low-quality genotypes

The deleterious mutation rate plays a key role in a number of important topics in biology, from mating system evolution to human health. Despite this broad significance, the nature and causes of variation in mutation rate are poorly understood, especially in multicellular organisms. We test whether genetic quality, the presence or absence of deleterious alleles, affects the mutation rate in Drosophila melanogaster by using a modified mutation accumulation approach. We find evidence that genotypes constructed to carry deleterious "treatment" alleles on one chromosome during mutation accumulation experience an elevated mutation rate on a different chromosome. Further, this elevation is correlated with the effect of the treatment alleles on phenotypic condition, measured as body mass. Treatment alleles that reduce mass by 10% cause a doubling in the rate of mutational decline. Our results show that mutation rates are sensitive to genetic stress, such that individuals with low-quality genotypes will produce offspring of even lower genetic quality, in a mutational positive feedback loop. This type of variation in mutation rate is expected to alter a variety of predictions based on mutation load theory and accelerate adaptation to new environments. Positive mutational feedback could affect human health by increasing the rate of germline mutation, and possibly somatic mutation, in individuals of poor health because of genetic or environmental stress.

The effect of variation in the effective population size on the rate of adaptive molecular evolution in eukaryotes

The role of adaptation is a fundamental question in molecular evolution. Theory predicts that species with large effective population sizes should undergo a higher rate of adaptive evolution than species with low effective population sizes if adaptation is limited by the supply of mutations. Previous analyses have appeared to support this conjecture because estimates of the proportion of nonsynonymous substitutions fixed by adaptive evolution, α, tend to be higher in species with large N_e. However, α is a function of both the number of advantageous and effectively neutral substitutions, either of which might depend on N_e. Here, we investigate the relationship between N_e and ω_a, the rate of adaptive evolution relative to the rate of neutral evolution, using nucleotide polymorphism and divergence data from 13 independent pairs of eukaryotic species. We find a highly significant positive correlation between ω_a and N_e. We also find some evidence that the rate of adaptive evolution varies between groups of organisms for a given N_e. The correlation between ω_a and N_e does not appear to be an artifact of demographic change or selection on synonymous codon use. Our results suggest that adaptation is to some extent limited by the supply of mutations and that at least some adaptation depends on newly occurring mutations rather than on standing genetic variation. Finally, we show that the proportion of nearly neutral nonadaptive substitutions declines with increasing N_e. The low rate of adaptive evolution and the high proportion of effectively neutral substitution in species with small N_e are expected to combine to make it difficult to detect adaptive molecular evolution in species with small N_e.

Variation in base-substitution mutation in experimental and natural lineages of Caenorhabditis nematodes

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.

Genome features of "Dark-fly", a Drosophila line reared long-term in a dark environment

Organisms are remarkably adapted to diverse environments by specialized metabolisms, morphology, or behaviors. To address the molecular mechanisms underlying environmental adaptation, we have utilized a Drosophila melanogaster line, termed “Dark-fly”, which has been maintained in constant dark conditions for 57 years (1400 generations). We found that Dark-fly exhibited higher fecundity in dark than in light conditions, indicating that Dark-fly possesses some traits advantageous in darkness. Using next-generation sequencing technology, we determined the whole genome sequence of Dark-fly and identified approximately 220,000 single nucleotide polymorphisms (SNPs) and 4,700 insertions or deletions (InDels) in the Dark-fly genome compared to the genome of the Oregon-R-S strain, a control strain. 1.8% of SNPs were classified as non-synonymous SNPs (nsSNPs: i.e., they alter the amino acid sequence of gene products). Among them, we detected 28 nonsense mutations (i.e., they produce a stop codon in the protein sequence) in the Dark-fly genome. These included genes encoding an olfactory receptor and a light receptor. We also searched runs of homozygosity (ROH) regions as putative regions selected during the population history, and found 21 ROH regions in the Dark-fly genome. We identified 241 genes carrying nsSNPs or InDels in the ROH regions. These include a cluster of alpha-esterase genes that are involved in detoxification processes. Furthermore, analysis of structural variants in the Dark-fly genome showed the deletion of a gene related to fatty acid metabolism. Our results revealed unique features of the Dark-fly genome and provided a list of potential candidate genes involved in environmental adaptation.

The role of background selection in shaping patterns of molecular evolution and variation: evidence from variability on the Drosophila X chromosome

In the putatively ancestral population of Drosophila melanogaster, the ratio of silent DNA sequence diversity for X-linked loci to that for autosomal loci is approximately one, instead of the expected "null" value of 3/4. One possible explanation is that background selection (the hitchhiking effect of deleterious mutations) is more effective on the autosomes than on the X chromosome, because of the lack of crossing over in male Drosophila. The expected effects of background selection on neutral variability at sites in the middle of an X chromosome or an autosomal arm were calculated for different models of chromosome organization and methods of approximation, using current estimates of the deleterious mutation rate and distributions of the fitness effects of deleterious mutations. The robustness of the results to different distributions of fitness effects, dominance coefficients, mutation rates, mapping functions, and chromosome size was investigated. The predicted ratio of X-linked to autosomal variability is relatively insensitive to these variables, except for the mutation rate and map length. Provided that the deleterious mutation rate per genome is sufficiently large, it seems likely that background selection can account for the observed X to autosome ratio of variability in the ancestral population of D. melanogaster. The fact that this ratio is much less than one in D. pseudoobscura is also consistent with the model's predictions, since this species has a high rate of crossing over. The results suggest that background selection may play a major role in shaping patterns of molecular evolution and variation.

Hsp90 is important for fecundity, longevity, and buffering of cryptic deleterious variation in wild fly populations

BACKGROUND

In the laboratory, the Drosophila melanogaster heat shock protein Hsp90 can buffer the phenotypic effects of genetic variation. Laboratory experiments either manipulate Hsp90 activity pharmacologically, or they induce mutations with strong effects in the gene Hsp83, the single-copy fly gene encoding Hsp90. It is unknown whether observations from such laboratory experiments are relevant in the wild.

RESULTS

We here study naturally occurring mutations in Hsp83, and their effects on fitness and phenotypic buffering in flies derived from wild populations. We examined more than 4500 flies from 42 Drosophila populations distributed world-wide for insertions or deletions of mobile DNA in or near the Hsp83 gene. The insertions we observed occur at low population frequencies, and reduce Hsp83 gene expression. In competition experiments, mutant flies performed much more poorly than wild-type flies. Mutant flies were also significantly less fecund and shorter-lived than wild-type flies, as well as less well buffered against cryptic deleterious variation, as we show through inbreeding experiments. Specifically, in Hsp83 mutant flies female fecundity dropped to much lower levels after inbreeding than in wild-type flies. At even slightly elevated temperatures, inbred mutant Hsp83 populations went extinct, whereas inbred wild-type populations persisted.

CONCLUSIONS

Our work shows that Hsp90, a regulator of the stress response and of signaling, helps buffer deleterious variation in fruit flies derived from wild population, and that its buffering role becomes even more important under heat stress.

Fitness of Arabidopsis thaliana mutation accumulation lines whose spontaneous mutations are known

Despite the fundamental importance of mutation to the evolutionary process, we have little knowledge of the direct consequences of specific spontaneous mutations to the fitness of the organism. Combining results of whole-genome sequencing with repeated field assays of survival and reproduction, we quantify the combined effects on fitness of spontaneous mutations identified in Arabidopsis thaliana. We demonstrate that the effects are beneficial, deleterious, or neutral depending on the environmental context. Some lines, bearing mutations disrupting known loci, differ strongly in fitness from the founder or premutation genotype. Those effects vary across environments, for example, a line with a major deletion spanning a transcription factor gene expressed lower fitness than the founder under most conditions but exceeded the founder's fitness in one environment. The large contribution of genotype by environment interaction (G × E) to mutation effects on fitness implies spatial and/or temporal variation in selection on new mutations and could contribute to the maintenance of standing genetic variation.

Selfish cellular networks and the evolution of complex organisms

Human gametogenesis takes years and involves many cellular divisions, particularly in males. Consequently, gametogenesis provides the opportunity to acquire multiple de novo mutations. A significant portion of these is likely to impact the cellular networks linking genes, proteins, RNA and metabolites, which constitute the functional units of cells. A wealth of literature shows that these individual cellular networks are complex, robust and evolvable. To some extent, they are able to monitor their own performance, and display sufficient autonomy to be termed "selfish". Their robustness is linked to quality control mechanisms which are embedded in and act upon the individual networks, thereby providing a basis for selection during gametogenesis. These selective processes are equally likely to affect cellular functions that are not gamete-specific, and the evolution of the most complex organisms, including man, is therefore likely to occur via two pathways: essential housekeeping functions would be regulated and evolve during gametogenesis within the parents before being transmitted to their progeny, while classical selection would operate on other traits of the organisms that shape their fitness with respect to the environment.

Orthologous genes identified by transcriptome sequencing in the spider genus Stegodyphus

Background

The evolution of sociality in spiders involves a transition from an outcrossing to a highly inbreeding mating system, a shift to a female biased sex ratio, and an increase in the reproductive skew among individuals. Taken together, these features are expected to result in a strong reduction in the effective population size. Such a decline in effective population size is expected to affect population genetic and molecular evolutionary processes, resulting in reduced genetic diversity and relaxed selective constraint across the genome. In the genus Stegodyphus, permanent sociality and regular inbreeding has evolved independently three times from periodic-social (outcrossing) ancestors. This genus is therefore an ideal model for comparative studies of the molecular evolutionary and population genetic consequences of the transition to a regularly inbreeding mating system. However, no genetic resources are available for this genus.

Results

We present the analysis of high throughput transcriptome sequencing of three Stegodyphus species. Two of these are periodic-social (Stegodyphus lineatus and S.tentoriicola) and one is permanently social (S. mimosarum). From non-normalized cDNA libraries, we obtained on average 7,000 putative uni-genes for each species. Three-way orthology, as predicted from reciprocal BLAST, identified 1,792 genes that could be used for cross-species comparison. Open reading frames (ORFs) could be deduced from 1,345 of the three-way alignments. Preliminary molecular analyses suggest a five- to ten-fold reduction in heterozygosity in the social S. mimosarum compared with the periodic-social species. Furthermore, an increased ratio of non-synonymous to synonymous polymorphisms in the social species indicated relaxed efficiency of selection. However, there was no sign of relaxed selection on the phylogenetic branch leading to S. mimosarum.

Conclusions

The 1,792 three-way ortholog genes identified in this study provide a unique resource for comparative studies of the eco-genomics, population genetics and molecular evolution of repeated evolution of inbreeding sociality within the Stegodyphus genus. Preliminary analyses support theoretical expectations of depleted heterozygosity and relaxed selection in the social inbreeding species. Relaxed selection could not be detected in the S. mimosarum lineage, suggesting that there has been a recent transition to sociality in this species.

Rates and fitness consequences of new mutations in humans

The human mutation rate per nucleotide site per generation (μ) can be estimated from data on mutation rates at loci causing Mendelian genetic disease, by comparing putatively neutrally evolving nucleotide sequences between humans and chimpanzees and by comparing the genome sequences of relatives. Direct estimates from genome sequencing of relatives suggest that μ is about 1.1 × 10(-8), which is about twofold lower than estimates based on the human-chimp divergence. This implies that an average of ~70 new mutations arise in the human diploid genome per generation. Most of these mutations are paternal in origin, but the male:female mutation rate ratio is currently uncertain and might vary even among individuals within a population. On the basis of a method proposed by Kondrashov and Crow, the genome-wide deleterious mutation rate (U) can be estimated from the product of the number of nucleotide sites in the genome, μ, and the mean selective constraint per site. Although the presence of many weakly selected mutations in human noncoding DNA makes this approach somewhat problematic, estimates are U ≈ 2.2 for the whole diploid genome per generation and 0.35 for mutations that change an amino acid of a protein-coding gene. A genome-wide deleterious mutation rate of 2.2 seems higher than humans could tolerate if natural selection is "hard," but could be tolerated if selection acts on relative fitness differences between individuals or if there is synergistic epistasis. I argue that in the foreseeable future, an accumulation of new deleterious mutations is unlikely to lead to a detectable decline in fitness of human populations.

Fitness effects of mutations in bacteria

Mutation is the primary source of variation in any organism. Without it, natural selection cannot operate and organisms cannot adapt to novel environments. Mutation is also generally a source of defect: many mutations are not neutral but cause fitness decreases in the organisms where they arise. In bacteria, another important source of variation is horizontal gene transfer. This source of variation can also cause beneficial or deleterious effects. Determining the distribution of fitness effects of mutations in different environments and genetic backgrounds is an active research field. In bacteria, knowledge of these distributions is key for understanding important traits. For example, for determining the dynamics of microorganisms with a high genomic mutation rate (mutators), and for understanding the evolution of antibiotic resistance, and the emergence of pathogenic traits. All of these characteristics are extremely relevant for human health both at the individual and population levels. Experimental evolution has been a valuable tool to address these questions. Here, we review some of the important findings of mutation effects in bacteria revealed through laboratory experiments.

The genotype-phenotype maps of systems biology and quantitative genetics: distinct and complementary

The processes by which genetic variation in complex traits is generated and maintained in populations has for a long time been treated in abstract and statistical terms. As a consequence, quantitative genetics has provided limited insights into our understanding of the molecular bases of quantitative trait variation. With the developing technological and conceptual tools of systems biology, cellular and molecular processes are being described in greater detail. While we have a good description of how signaling and other molecular networks are organized in the cell, we still do not know how genetic variation affects these pathways, because systems and molecular biology usually ignore the type and extent of genetic variation found in natural populations. Here we discuss the quantitative genetics and systems biology approaches for the study of complex trait architecture and discuss why these two disciplines would synergize with each other to answer questions that neither of the two could answer alone.

High mutation rates in the mitochondrial genomes of Daphnia pulex

Despite the great utility of mitochondrial DNA (mtDNA) sequence data in population genetics and phylogenetics, key parameters describing the process of mitochondrial mutation (e.g., the rate and spectrum of mutational change) are based on few direct estimates. Furthermore, the variation in the mtDNA mutation process within species or between lineages with contrasting reproductive strategies remains poorly understood. In this study, we directly estimate the mtDNA mutation rate and spectrum using Daphnia pulex mutation-accumulation (MA) lines derived from sexual (cyclically parthenogenetic) and asexual (obligately parthenogenetic) lineages. The nearly complete mitochondrial genome sequences of 82 sexual and 47 asexual MA lines reveal high mtDNA mutation rate of 1.37 × 10(-7) and 1.73 × 10(-7) per nucleotide per generation, respectively. The Daphnia mtDNA mutation rate is among the highest in eukaryotes, and its spectrum is dominated by insertions and deletions (70%), largely due to the presence of mutational hotspots at homopolymeric nucleotide stretches. Maximum likelihood estimates of the Daphnia mitochondrial effective population size reveal that between five and ten copies of mitochondrial genomes are transmitted per female per generation. Comparison between sexual and asexual lineages reveals no statistically different mutation rates and highly similar mutation spectra.

Variation in the mutation rate across mammalian genomes

It has been known for many years that the mutation rate varies across the genome. However, only with the advent of large genomic data sets is the full extent of this variation becoming apparent. The mutation rate varies over many different scales, from adjacent sites to whole chromosomes, with the strongest variation seen at the smallest scales. Some of these patterns have clear mechanistic bases, but much of the rate variation remains unexplained, and some of it is deeply perplexing. Variation in the mutation rate has important implications in evolutionary biology and underexplored implications for our understanding of hereditary disease and cancer.

Find the weakest link. A comparison between demographic, genetic and demo-genetic metapopulation extinction times

Background

While the ultimate causes of most species extinctions are environmental, environmental constraints have various secondary consequences on evolutionary and ecological processes. The roles of demographic, genetic mechanisms and their interactions in limiting the viabilities of species or populations have stirred much debate and remain difficult to evaluate in the absence of demography-genetics conceptual and technical framework. Here, I computed projected times to metapopulation extinction using (1) a model focusing on the effects of species properties, habitat quality, quantity and temporal variability on the time to demographic extinction; (2) a genetic model focusing on the dynamics of the drift and inbreeding loads under the same species and habitat constraints; (3) a demo-genetic model accounting for demographic-genetic processes and feedbacks.

Results

Results indicate that a given population may have a high demographic, but low genetic viability or vice versa; and whether genetic or demographic aspects will be the most limiting to overall viability depends on the constraints faced by the species (e.g., reduction of habitat quantity or quality). As a consequence, depending on metapopulation or species characteristics, incorporating genetic considerations to demographically-based viability assessments may either moderately or severely reduce the persistence time. On the other hand, purely genetically-based estimates of species viability may either underestimate (by neglecting demo-genetic interactions) or overestimate (by neglecting the demographic resilience) true viability.

Conclusion

Unbiased assessments of the viabilities of species may only be obtained by identifying and considering the most limiting processes (i.e., demography or genetics), or, preferentially, by integrating them.

Chimeric genes as a source of rapid evolution in Drosophila melanogaster

Chimeric genes form through the combination of portions of existing coding sequences to create a new open reading frame. These new genes can create novel protein structures that are likely to serve as a strong source of novelty upon which selection can act. We have identified 14 chimeric genes that formed through DNA-level mutations in Drosophila melanogaster, and we investigate expression profiles, domain structures, and population genetics for each of these genes to examine their potential to effect adaptive evolution. We find that chimeric gene formation commonly produces mid-domain breaks and unites portions of wholly unrelated peptides, creating novel protein structures that are entirely distinct from other constructs in the genome. These new genes are often involved in selective sweeps. We further find a disparity between chimeric genes that have recently formed and swept to fixation versus chimeric genes that have been preserved over long periods of time, suggesting that preservation and adaptation are distinct processes. Finally, we demonstrate that chimeric gene formation can produce qualitative expression changes that are difficult to mimic through duplicate gene formation, and that extremely young chimeric genes (d(S) < 0.03) are more likely to be associated with selective sweeps than duplicate genes of the same age. Hence, chimeric genes can serve as an exceptional source of genetic novelty that can have a profound influence on adaptive evolution in D. melanogaster.

Kin selection-mutation balance: a model for the origin, maintenance, and consequences of social cheating

Social conflict, in the form of intraspecific selfish "cheating," has been observed in a number of natural systems. However, a formal, evolutionary genetic theory of social cheating that provides an explanatory, predictive framework for these observations is lacking. Here we derive the kin selection-mutation balance, which provides an evolutionary null hypothesis for the statics and dynamics of cheating. When social interactions have linear fitness effects and Hamilton's rule is satisfied, selection is never strong enough to eliminate recurrent cheater mutants from a population, but cheater lineages are transient and do not invade. Instead, cheating lineages are eliminated by kin selection but are constantly reintroduced by mutation, maintaining a stable equilibrium frequency of cheaters. The presence of cheaters at equilibrium creates a "cheater load" that selects for mechanisms of cheater control, such as policing. We find that increasing relatedness reduces the cheater load more efficiently than does policing the costs and benefits of cooperation. Our results provide new insight into the effects of genetic systems, mating systems, ecology, and patterns of sex-limited expression on social evolution. We offer an explanation for the widespread cheater/altruist polymorphism found in nature and suggest that the common fear of conflict-induced social collapse is unwarranted.

Time-dependent rates of molecular evolution

For over half a century, it has been known that the rate of morphological evolution appears to vary with the time frame of measurement. Rates of microevolutionary change, measured between successive generations, were found to be far higher than rates of macroevolutionary change inferred from the fossil record. More recently, it has been suggested that rates of molecular evolution are also time dependent, with the estimated rate depending on the timescale of measurement. This followed surprising observations that estimates of mutation rates, obtained in studies of pedigrees and laboratory mutation-accumulation lines, exceeded long-term substitution rates by an order of magnitude or more. Although a range of studies have provided evidence for such a pattern, the hypothesis remains relatively contentious. Furthermore, there is ongoing discussion about the factors that can cause molecular rate estimates to be dependent on time. Here we present an overview of our current understanding of time-dependent rates. We provide a summary of the evidence for time-dependent rates in animals, bacteria and viruses. We review the various biological and methodological factors that can cause rates to be time dependent, including the effects of natural selection, calibration errors, model misspecification and other artefacts. We also describe the challenges in calibrating estimates of molecular rates, particularly on the intermediate timescales that are critical for an accurate characterization of time-dependent rates. This has important consequences for the use of molecular-clock methods to estimate timescales of recent evolutionary events.

Experimental mutation-accumulation on the X chromosome of Drosophila melanogaster reveals stronger selection on males than females

Background

Sex differences in the magnitude or direction of mutational effect may be important to a variety of population processes, shaping the mutation load and affecting the cost of sex itself. These differences are expected to be greatest after sexual maturity. Mutation-accumulation (MA) experiments provide the most direct way to examine the consequences of new mutations, but most studies have focused on juvenile viability without regard to sex, and on autosomes rather than sex chromosomes; both adult fitness and X-linkage have been little studied. We therefore investigated the effects of 50 generations of X-chromosome mutation accumulation on the fitness of males and females derived from an outbred population of Drosophila melanogaster.

Results

Fitness declined rapidly in both sexes as a result of MA, but adult males showed markedly greater fitness loss relative to their controls compared to females expressing identical genotypes, even when females were made homozygous for the X. We estimate that these mutations are partially additive (h ~ 0.3) in females. In addition, the majority of new mutations appear to harm both males and females.

Conclusions

Our data helps fill a gap in our understanding of the consequences of sexual selection for genetic load, and suggests that stronger selection on males may indeed purge deleterious mutations affecting female fitness.

Use of whole genome sequencing to estimate the mutation rate of Mycobacterium tuberculosis during latent infection

Tuberculosis poses a global health emergency, which has been compounded by the emergence of drug-resistant Mycobacterium tuberculosis (Mtb) strains. We used whole-genome sequencing to compare the accumulation of mutations in Mtb isolated from cynomolgus macaques with active, latent or reactivated disease. We sequenced 33 Mtb isolates from nine macaques with an average genome coverage of 93% and an average read depth of 117×. Based on the distribution of SNPs observed, we calculated the mutation rates for these disease states. We found a similar mutation rate during latency as during active disease or in a logarithmically growing culture over the same period of time. The pattern of polymorphisms suggests that the mutational burden in vivo is because of oxidative DNA damage. We show that Mtb continues to acquire mutations during disease latency, which may explain why isoniazid monotherapy for latent tuberculosis is a risk factor for the emergence of isoniazid resistance.

Molecular population genomics: a short history

Population genomics is the study of the amount and causes of genome-wide variability in natural populations, a topic that has been under discussion since Darwin. This paper first briefly reviews the early development of molecular approaches to the subject: the pioneering unbiased surveys of genetic variability at multiple loci by means of gel electrophoresis and restriction enzyme mapping. The results of surveys of levels of genome-wide variability using DNA resequencing studies are then discussed. Studies of the extent to which variability for different classes of variants (non-synonymous, synonymous and non-coding) are affected by natural selection, or other directional forces such as biased gene conversion, are also described. Finally, the effects of deleterious mutations on population fitness and the possible role of Hill–Robertson interference in shaping patterns of sequence variability are discussed.

The maintenance of genetic variation for oviposition rate in two-spotted spider mites: inferences from artificial selection

Despite the directional selection acting on life-history traits, substantial amounts of standing variation for these traits have frequently been found. This variation may result from balancing selection (e.g., through genetic trade-offs) or from mutation-selection balance. These mechanisms affect allele frequencies in different ways: Under balancing selection alleles are maintained at intermediate frequencies, whereas under mutation-selection balance variation is generated by deleterious mutations and removed by directional selection, which leads to asymmetry in the distribution of allele frequencies. To investigate the importance of these two mechanisms in maintaining heritable variation in oviposition rate of the two-spotted spider mite, we analyzed the response to artificial selection. In three replicate experiments, we selected for higher and lower oviposition rate, compared to control lines. A response to selection only occurred in the downward direction. Selection for lower oviposition rate did not lead to an increase in any other component of fitness, but led to a decline in female juvenile survival. The results suggest standing variation for oviposition rate in this population consists largely of deleterious alleles, as in a mutation-selection balance. Consequently, the standing variation for this trait does not appear to be indicative of its adaptive potential.