Dominance of mutations affecting viability in Drosophila melanogaster

Why are deleterious mutations maintained in selfing populations? An analysis of the effects of early- and late-acting mutations by a two-locus two-allele model

Frequencies of deleterious mutations are higher than expected in many plants. Here, by developing a two-locus two-allele model, I examine the effects of differential timing of the expression of deleterious mutations (two-stage effects) on the maintenance of mutations. I assume early- and late-acting loci to distinguish whether maintenance of mutations in populations with high selfing rates is explained better by two-stage effects of single mutations, or by separate mutations in both early- and late-acting loci. I found that, when ovules are overproduced, the stable frequency of early-acting mutations is higher if mutations also occur in a late-acting locus than if a late-acting mutation is lacking. The stable frequency of late-acting mutations is higher if mutations also occur in an early-acting locus than if an early-acting mutation is lacking. Selective interference does not account for these results because analyses in which the number of loci subject to mutations is equalized are included. Overproduction of ovules has little effect on maintenance if either early- or late-acting mutations are lacking, whereas when ovules are not overproduced, the two-stage effect does not enhance the maintenance of mutations. Hence, mutations occurring in both loci coupled with overproduction of ovules enhances the maintenance of mutations in populations with high selfing rates. The detailed mechanisms underlying the two-stage effect were also analyzed.

Meiotic drive does not cause condition-dependent reduction of the sexual ornament in stalk-eyed flies

Meiotic drive systems are associated with low-frequency chromosomal inversions. These are expected to accumulate deleterious mutations due to reduced recombination and low effective population size. We test this prediction using the 'sex-ratio' (SR) meiotic drive system of the Malaysian stalk-eyed fly Teleopsis dalmanni. SR is associated with a large inversion (or inversions) on the X chromosome. In particular, we study eyespan in males carrying the SR chromosome, as this trait is a highly exaggerated, sexually dimorphic trait, known to have heightened condition-dependent expression. Larvae were raised in low and high larval food stress environments. SR males showed reduced eyespan under the low and high stress treatments, but there was no evidence of a condition-dependent decrease in eyespan under high stress. Similar but more complex patterns were observed for female eyespan, with evidence of additivity under low stress and heterosis under high stress. These results do not support the hypothesis that reduced sexual ornament size in meiotic drive males is due to a condition-dependent response to the putative increase in mutation load. Instead, reduced eyespan likely reflects compensatory resource allocation to different traits in response to drive-mediated destruction of sperm.

Maintenance of high inbreeding depression in selfing populations: Two-stage effect of early- and late-acting mutations

High estimates of inbreeding depression have been obtained in many plant populations with high selfing rates. However, deleterious mutations might be purged from such populations as a result of selfing. I developed a simulation model assuming the presence of mutations at two sets of loci, namely, early- and late-acting loci, and the selective abortion of embryos coupled with ovule overproduction. In the model, early-acting loci are expressed during embryo initiation, and less vigorous embryos are aborted. Late-acting loci are expressed after selective abortion ends; the surviving embryos (seeds) compete, and some of them form the next generation. If mutations are allowed to occur in both early- and late-acting loci, both types increase in frequency in populations with high selfing rates. However, this phenomenon does not occur if mutations occur only in the early- or only in the late-acting loci. Consistent results are observed even if the total number of loci in which mutations are allowed to occur is the same for simulations with both early- and late-acting loci, only early-acting loci, or only late-acting loci, indicating that the presence of both types of loci is the causal factor. Thus, the two-stage effect, or occurrence of both early- and late-acting mutations, promotes the maintenance of these mutations in populations with high selfing rates.

Obstruction of adaptation in diploids by recessive, strongly deleterious alleles

Recessive deleterious mutations are common, causing many genetic disorders in humans and producing inbreeding depression in the majority of sexually reproducing diploids. The abundance of recessive deleterious mutations in natural populations suggests they are likely to be present on a chromosome when a new adaptive mutation occurs, yet the dynamics of recessive deleterious hitchhikers and their impact on adaptation remains poorly understood. Here we model how a recessive deleterious mutation impacts the fate of a genetically linked dominant beneficial mutation. The frequency trajectory of the adaptive mutation in this case is dramatically altered and results in what we have termed a "staggered sweep." It is named for its three-phased trajectory: (i) Initially, the two linked mutations have a selective advantage while rare and will increase in frequency together, then (ii), at higher frequencies, the recessive hitchhiker is exposed to selection and can cause a balanced state via heterozygote advantage (the staggered phase), and (iii) finally, if recombination unlinks the two mutations, then the beneficial mutation can complete the sweep to fixation. Using both analytics and simulations, we show that strongly deleterious recessive mutations can substantially decrease the probability of fixation for nearby beneficial mutations, thus creating zones in the genome where adaptation is suppressed. These mutations can also significantly prolong the number of generations a beneficial mutation takes to sweep to fixation, and cause the genomic signature of selection to resemble that of soft or partial sweeps. We show that recessive deleterious variation could impact adaptation in humans and Drosophila.

The dominance effect of the adaptive transposable element insertion Bari-Jheh depends on the genetic background

Although adaptive mutations are often considered to be dominant, it has been recently shown that a substantial proportion of adaptive mutations should display heterozygote advantage. In this work, we take advantage of a recently characterized transposable element insertion mediating oxidative stress response in Drosophila melanogaster to test the dominance effect of an adaptive mutation. The comparison of the survival curves of heterozygous and the two corresponding homozygous flies indicated that the dominance effect of Bari-Jheh depends on the genetic background. Both in homozygous and in heterozygous flies, Bari-Jheh was associated with upregulation of Jheh1 (Juvenile Hormone Epoxyde Hydrolase 1) and/or Jheh2 genes. Our results add to the limited number of studies in which the dominance effect of adaptive mutations has been empirically estimated and highlights the complexity of their inheritance.

On the retention of gene duplicates prone to dominant deleterious mutations

Recent studies have shown that gene families from different functional categories have been preferentially expanded either by small scale duplication (SSD) or by whole-genome duplication (WGD). In particular, gene families prone to dominant deleterious mutations and implicated in cancers and other genetic diseases in human have been greatly expanded through two rounds of WGD dating back from early vertebrates. Here, we strengthen this intriguing observation, showing that human oncogenes involved in different primary tumors have retained many WGD duplicates compared to other human genes. In order to rationalize this evolutionary outcome, we propose a consistent population genetics model to analyze the retention of SSD and WGD duplicates taking into account their propensity to acquire dominant deleterious mutations. We solve a deterministic haploid model including initial duplicated loci, their retention through sub-functionalization or their neutral loss-of-function or deleterious gain-of-function at one locus. Extensions to diploid genotypes are presented and population size effects are analyzed using stochastic simulations. The only difference between the SSD and WGD scenarios is the initial number of individuals with duplicated loci. While SSD duplicates need to spread through the entire population from a single individual to reach fixation, WGD duplicates are de facto fixed in the small initial post-WGD population arising through the ploidy incompatibility between post-WGD individuals and the rest of the pre-WGD population. WGD duplicates prone to dominant deleterious mutations are then shown to be indirectly selected through purifying selection in post-WGD species, whereas SSD duplicates typically require positive selection. These results highlight the long-term evolution mechanisms behind the surprising accumulation of WGD duplicates prone to dominant deleterious mutations and are shown to be consistent with cancer genome data on the prevalence of human oncogenes with WGD duplicates.

Fitness landscapes: an alternative theory for the dominance of mutation

Deleterious mutations tend to be recessive. Several theories, notably those of Fisher (based on selection) and Wright (based on metabolism), have been put forward to explain this pattern. Despite a long-lasting debate, the matter remains unresolved. This debate has focused on the average dominance of mutations. However, we also know very little about the distribution of dominance coefficients among mutations, and about its variation across environments. In this article we present a new approach to predicting this distribution. Our approach is based on a phenotypic fitness landscape model. First, we show that under a very broad range of conditions (and environments), the average dominance of mutation of small effects should be approximately one-quarter as long as adaptation of organisms to their environment can be well described by stabilizing selection on an arbitrary set of phenotypic traits. Second, the theory allows predicting the whole distribution of dominance coefficients among mutants. Because it provides quantitative rather than qualitative predictions, this theory can be directly compared to data. We found that its prediction on mean dominance (average dominance close to 0.25) agreed well with the data, based on a meta-analysis of dominance data for mildly deleterious mutations. However, a simple landscape model does not account for the dominance of mutations of large effects and we provide possible extension of the theory for this class of mutations. Because dominance is a central parameter for evolutionary theory, and because these predictions are quantitative, they set the stage for a wide range of applications and further empirical tests.

Inferences about the distribution of dominance drawn from yeast gene knockout data

Data from several thousand knockout mutations in yeast (Saccharomyces cerevisiae) were used to estimate the distribution of dominance coefficients. We propose a new unbiased likelihood approach to measuring dominance coefficients. On average, deleterious mutations are partially recessive, with a mean dominance coefficient ~0.2. Alleles with large homozygous effects are more likely to be more recessive than are alleles of weaker effect. Our approach allows us to quantify, for the first time, the substantial variance and skew in the distribution of dominance coefficients. This heterogeneity is so great that many population genetic processes analyses based on the mean dominance coefficient alone will be in substantial error. These results are applied to the debate about various mechanisms for the evolution of dominance, and we conclude that they are most consistent with models that depend on indirect selection on homeostatic gene expression or on the ability to perform well under periods of high demand for a protein.

Inferring the strength of selection in Drosophila under complex demographic models

Transposable elements (TEs) constitute a substantial fraction of the genomes of many species, and it is thus important to understand their population dynamics. The strength of natural selection against TEs is a key parameter in understanding these dynamics. In principle, the strength of selection can be inferred from the frequencies of a sample of TEs. However, complicated demographic histories, such as found in Drosophila melanogaster, could lead to a substantial distortion of the TE frequency distribution compared with that expected for a panmictic, constant-sized population. The current methodology for the estimation of selection intensity acting against TEs does not take into account demographic history and might generate erroneous estimates especially for TE families under weak selection. Here, we develop a flexible maximum likelihood methodology that explicitly accounts both for demographic history and for the ascertainment biases of identifying TEs. We apply this method to the newly generated frequency data of the BS family of non-long terminal repeat retrotransposons in D. melanogaster in concert with two recent models of the demographic history of the species to infer the intensity of selection against this family. We find the estimate to differ substantially compared with a prior estimate that was made assuming a model of constant population size. Further, we find there to be relatively little information about selection intensity present in the derived non-African frequency data and that the ancestral African subpopulation is much more informative in this respect. These findings highlight the importance of accounting for demographic history and bear on study design for the inference of selection coefficients generally.

Increase of the spontaneous mutation rate in a long-term experiment with Drosophila melanogaster

In a previous experiment, the effect of 255 generations of mutation accumulation (MA) on the second chromosome viability of Drosophila melanogaster was studied using 200 full-sib MA1 lines and a large C1 control, both derived from a genetically homogeneous base population. At generation 265, one of those MA1 lines was expanded to start 150 new full-sib MA2 lines and a new C2 large control. After 46 generations, the rate of decline in mean viability in MA2 was approximately 2.5 times that estimated in MA1, while the average degree of dominance of mutations was small and nonsignificant by generation 40 and moderate by generation 80. In parallel, the inbreeding depression rate for viability and the amount of additive variance for two bristle traits in C2 were 2-3 times larger than those in C1. The results are consistent with a mutation rate in the line from which MA2 and C2 were derived about 2.5 times larger than that in MA1. The mean viability of C2 remained roughly similar to that of C1, but the rate of MA2 line extinction increased progressively, leading to mutational collapse, which can be ascribed to accelerated mutation and/or synergy after important deleterious accumulation.

Gene action of new mutations in Arabidopsis thaliana

For a newly arising mutation affecting a trait under selection, its degree of dominance relative to the preexisting allele(s) strongly influences its evolutionary impact. We have estimated dominance parameters for spontaneous mutations in a subset of lines derived from a highly inbred founder of Arabidopsis thaliana by at least 17 generations of mutation accumulation (MA). The labor-intensive nature of the crosses and the anticipated subtlety of effects limited the number of MA lines included in this study to 8. Each MA line was selfed and reciprocally crossed to plants representing the founder genotype, and progeny were assayed in the greenhouse. Significant mutational effects on reproductive fitness included a recessive fitness-enhancing effect in one line and fitness-reducing effects, one additive and the other slightly recessive. Mutations conferring earlier phenology or smaller leaves were significantly recessive. For effects increasing leaf number and reducing height at flowering, additive gene action accounted for the expression of the traits. The sole example of a significantly dominant mutational effect delayed phenology. Our findings of recessive action of a fitness-enhancing mutational effect and additive action of a deleterious effect counter a common expectation of (partial) dominance of alleles that increase fitness, but the frequency of occurrence of such mutations is unknown.

The first steps of transposable elements invasion: parasitic strategy vs. genetic drift

Transposable elements are often considered as selfish DNA sequences able to invade the genome of their host species. Their evolutive dynamics are complex, due to the interaction between their intrinsic amplification capacity, selection at the host level, transposition regulation, and genetic drift. Here, we propose modeling the first steps of TE invasion, i.e., just after a horizontal transfer, when a single copy is present in the genome of one individual. If the element has a constant transposition rate, it will disappear in most cases: the elements with low-transposition rate are frequently lost through genetic drift, while those with high-transposition rate may amplify, leading to the sterility of their host. Elements whose transposition rate is regulated are able to successfully invade the populations, thanks to an initial transposition burst followed by a strong limitation of their activity. Self-regulation or hybrid dysgenesis may thus represent some genome-invasion parasitic strategies.

Estimates of the genomic mutation rate for detrimental alleles in Drosophila melanogaster

The net rate of mutation to deleterious but nonlethal alleles and the sizes of effects of these mutations are of great significance for many evolutionary questions. Here we describe three replicate experiments in which mutations have been accumulated on chromosome 3 of Drosophila melanogaster by means of single-male backcrosses of heterozygotes for a wild-type third chromosome. Egg-to-adult viability was assayed for nonlethal homozygous chromosomes. The rates of decline in mean and increase in variance (DM and DV, respectively) were estimated. Scaled up to the diploid whole genome, the mean DM for homozygous detrimental mutations over the three experiments was between 0.8 and 1.8%. The corresponding DV estimate was approximately 0.11%. Overall, the results suggest a lower bound estimate of at least 12% for the diploid per genome mutation rate for detrimentals. The upper bound estimates for the mean selection coefficient were between 2 and 10%, depending on the method used. Mutations with selection coefficients of at least a few percent must be the major contributors to the effects detected here and are likely to be caused mostly by transposable element insertions or indels.

Analysis of the estimators of the average coefficient of dominance of deleterious mutations

We investigate the sources of bias that affect the most commonly used methods of estimation of the average degree of dominance (h) of deleterious mutations, focusing on estimates from segregating populations. The main emphasis is on the effect of the finite size of the populations, but other sources of bias are also considered. Using diffusion approximations to the distribution of gene frequencies in finite populations as well as stochastic simulations, we assess the behavior of the estimators obtained from populations at mutation-selection-drift balance under different mutational scenarios and compare averages of h for newly arisen and segregating mutations. Because of genetic drift, the inferences concerning newly arisen mutations based on the mutation-selection balance theory can have substantial upward bias depending upon the distribution of h. In addition, estimates usually refer to h weighted by the homozygous deleterious effect in different ways, so that inferences are complicated when these two variables are negatively correlated. Due to both sources of bias, the widely used regression of heterozygous on homozygous means underestimates the arithmetic mean of h for segregating mutations, in contrast to their repeatedly assumed equality in the literature. We conclude that none of the estimators from segregating populations provides, under general conditions, a useful tool to ascertain the properties of the degree of dominance, either for segregating or for newly arisen deleterious mutations. Direct estimates of the average h from mutation-accumulation experiments are shown to suffer some bias caused by purging selection but, because they do not require assumptions on the causes maintaining segregating variation, they appear to give a more reliable average dominance for newly arisen mutations.

On the rate and linearity of viability declines in Drosophila mutation-accumulation experiments: genomic mutation rates and synergistic epistasis revisited

High rates of deleterious mutations could severely reduce the fitness of populations, even endangering their persistence; these effects would be mitigated if mutations synergize each others' effects. An experiment by Mukai in the 1960s gave evidence that in Drosophila melanogaster, viability-depressing mutations occur at the surprisingly high rate of around one per zygote and that the mutations interact synergistically. A later experiment by Ohnishi seemed to support the high mutation rate, but gave no evidence for synergistic epistasis. Both of these studies, however, were flawed by the lack of suitable controls for assessing viability declines of the mutation-accumulation (MA) lines. By comparing homozygous viability of the MA lines to simultaneously estimated heterozygous viability and using estimates of the dominance of mutations in the experiments, I estimate the viability declines relative to an appropriate control. This approach yields two unexpected conclusions. First, in Ohnishi's experiment as well as in Mukai's, MA lines showed faster-than-linear declines in viability, indicative of synergistic epistasis. Second, while Mukai's estimate of the genomic mutation rate is supported, that from Ohnishi's experiment is an order of magnitude lower. The different results of the experiments most likely resulted from differences in the starting genotypes; even within Mukai's experiment, a subset of MA lines, which I argue probably resulted from a contamination event, showed much slower viability declines than did the majority of lines. Because different genotypes may show very different mutational behavior, only studies using many founding genotypes can determine the average rate and distribution of effects of mutations relevant to natural populations.

Influence of dominance, leptokurtosis and pleiotropy of deleterious mutations on quantitative genetic variation at mutation-selection balance

In models of maintenance of genetic variance (V (G)) it has often been assumed that mutant alleles act additively. However, experimental data show that the dominance coefficient varies among mutant alleles and those of large effect tend to be recessive. On the basis of empirical knowledge of mutations, a joint-effect model of pleiotropic and real stabilizing selection that includes dominance is constructed and analyzed. It is shown that dominance can dramatically alter the prediction of equilibrium V (G). Analysis indicates that for the situations where mutations are more recessive for fitness than for a quantitative trait, as supported by the available data, the joint-effect model predicts a significantly higher V (G) than does an additive model. Importantly, for what seem to be realistic distributions of mutational effects (i.e., many mutants may not affect the quantitative trait substantially but are likely to affect fitness), the observed high levels of genetic variation in the quantitative trait under strong apparent stabilizing selection can be generated. This investigation supports the hypothesis that most V (G) comes from the alleles nearly neutral for fitness in heterozygotes while apparent stabilizing selection is contributed mainly by the alleles of large effect on the quantitative trait. Thus considerations of dominance coefficients of mutations lend further support to our previous conclusion that mutation-selection balance is a plausible mechanism of the maintenance of the genetic variance in natural populations.

On the Rate and Linearity of Viability Declines in Drosophila Mutation-Accumulation Experiments: Genomic Mutation Rates and Synergistic Epistasis Revisited

High rates of deleterious mutations could severely reduce the fitness of populations, even endangering their persistence; these effects would be mitigated if mutations synergize each others’ effects. An experiment by Mukai in the 1960s gave evidence that in Drosophila melanogaster, viability-depressing mutations occur at the surprisingly high rate of around one per zygote and that the mutations interact synergistically. A later experiment by Ohnishi seemed to support the high mutation rate, but gave no evidence for synergistic epistasis. Both of these studies, however, were flawed by the lack of suitable controls for assessing viability declines of the mutation-accumulation (MA) lines. By comparing homozygous viability of the MA lines to simultaneously estimated heterozygous viability and using estimates of the dominance of mutations in the experiments, I estimate the viability declines relative to an appropriate control. This approach yields two unexpected conclusions. First, in Ohnishi’s experiment as well as in Mukai’s, MA lines showed faster-than-linear declines in viability, indicative of synergistic epistasis. Second, while Mukai’s estimate of the genomic mutation rate is supported, that from Ohnishi’s experiment is an order of magnitude lower. The different results of the experiments most likely resulted from differences in the starting genotypes; even within Mukai’s experiment, a subset of MA lines, which I argue probably resulted from a contamination event, showed much slower viability declines than did the majority of lines. Because different genotypes may show very different mutational behavior, only studies using many founding genotypes can determine the average rate and distribution of effects of mutations relevant to natural populations.

Dominance and Overdominance of Mildly Deleterious Induced Mutations for Fitness Traits inCaenorhabditis elegans

We estimated the average dominance coefficient of mildly deleterious mutations (h, the proportion by which mutations in the heterozygous state reduce fitness components relative to those in the homozygous state) in the nematode Caenorhabditis elegans. From 56 worm lines that carry mutations induced by the point mutagen ethyl methanesulfonate (EMS), we selected 19 lines that are relatively high in fitness and estimated the viabilities, productivities, and relative fitnesses of heterozygotes and homozygotes compared to the ancestral wild type. There was very little effect of homozygous or heterozygous mutations on egg-to-adult viability. For productivity and relative fitness, we found that the average dominance coefficient, h, was ∼0.1, suggesting that mildly deleterious mutations are on average partially recessive. These estimates were not significantly different from zero (complete recessivity) but were significantly different from 0.5 (additivity). In addition, there was a significant amount of variation in h among lines, and analysis of average dominance coefficients of individual lines suggested that several lines showed overdominance for fitness. Further investigation of two of these lines partially confirmed this finding.