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.     

The effect of antagonistic pleiotropy on the estimation of the average coefficient of dominance of deleterious mutations

We investigate the impact of antagonistic pleiotropy on the most widely used methods of estimation of the average coefficient of dominance of deleterious mutations from segregating populations. A proportion of the deleterious mutations affecting a given studied fitness component are assumed to have an advantageous effect on another one, generating overdominance on global fitness. Using diffusion approximations and transition matrix methods, we obtain the distribution of gene frequencies for nonpleiotropic and pleiotropic mutations in populations at the mutation-selection-drift balance. From these distributions we build homozygous and heterozygous chromosomes and assess the behavior of the estimators of dominance. A very small number of deleterious mutations with antagonistic pleiotropy produces substantial increases on the estimate of the average degree of dominance of mutations affecting the fitness component under study. For example, estimates are increased three- to fivefold when 2% of segregating loci are over-dominant for fitness. In contrast, strengthening pleiotropy, where pleiotropic effects are assumed to be also deleterious, has little effect on the estimates of the average degree of dominance, supporting previous results. The antagonistic pleiotropy model considered, applied under mutational parameters described in the literature, produces patterns for the distribution of chromosomal viabilities, levels of genetic variance, and homozygous mutation load generally consistent with those observed empirically for viability in Drosophila melanogaster.     

Effects of Single P-Element Insertions on Bristle Number and Viability in Drosophila Melanogaster

Single P-element mutagenesis was used to construct 1094 lines with P[lArB] inserts on all three major chromosomes in an isogenic background previously free of P elements. The effects of insertions on bristle number and on viability were assessed by comparison to 392 control lines. The variance and effects of P-element inserts on bristle number and viability were larger than those inferred from spontaneous mutations. The distributions of effects on bristle number were symmetrical and highly leptokurtic, such that a few inserts with large effects caused most of the increase in variance. The distribution of effects on viability were negatively skewed and platykurtic. On average, the effects of P-element insertions on bristle number were partly recessive and on viability were completely recessive. P-element inserts with large effects on bristle number tended to have reduced viability, but the correlation between the absolute value of the effects on bristle number and on viability was not strong. Fifty P-element inserts tagging quantitative trait loci (QTLs) with large effects on bristle number were mapped cytogenetically. Two P-element-induced scabrous alleles and five extramacrochaetae alleles were generated. Single P-element mutagenesis is a powerful method for identifying QTLs at the level of genetic locus.     

The Impact of Population Demography and Selection on the Genetic Architecture of Complex Traits

Population genetic studies have found evidence for dramatic population growth in recent human history. It is unclear how this recent population growth, combined with the effects of negative natural selection, has affected patterns of deleterious variation, as well as the number, frequency, and effect sizes of mutations that contribute risk to complex traits. Because researchers are performing exome sequencing studies aimed at uncovering the role of low-frequency variants in the risk of complex traits, this topic is of critical importance. Here I use simulations under population genetic models where a proportion of the heritability of the trait is accounted for by mutations in a subset of the exome. I show that recent population growth increases the proportion of nonsynonymous variants segregating in the population, but does not affect the genetic load relative to a population that did not expand. Under a model where a mutation''s effect on a trait is correlated with its effect on fitness, rare variants explain a greater portion of the additive genetic variance of the trait in a population that has recently expanded than in a population that did not recently expand. Further, when using a single-marker test, for a given false-positive rate and sample size, recent population growth decreases the expected number of significant associations with the trait relative to the number detected in a population that did not expand. However, in a model where there is no correlation between a mutation''s effect on fitness and the effect on the trait, common variants account for much of the additive genetic variance, regardless of demography. Moreover, here demography does not affect the number of significant associations detected. These findings suggest recent population history may be an important factor influencing the power of association tests and in accounting for the missing heritability of certain complex traits.     

Effects of P Element Insertions on Quantitative Traits in Drosophila Melanogaster

P element mutagenesis was used to construct 94 third chromosome lines of Drosophila melanogaster which contained on average 3.1 stable P element inserts, in an inbred host strain background previously free of P elements. The homozygous and heterozygous effects of the inserts on viability and abdominal and sternopleural bristle number were ascertained by comparing the chromosome lines with inserts to insert-free control lines of the inbred host strain. P elements reduced average homozygous viability by 12.2% per insert and average heterozygous viability by 5.5% per insert, and induced recessive lethal mutations at a rate of 3.8% per insert. Mutational variation for the bristle traits averaged over both sexes was 0.03Ve per homozygous P insert and 0.003Ve per heterozygous P insert, where Ve is the environmental variance. Mutational variation was greater for the sexes considered separately because inserts had large pleiotropic effects on sex dimorphism of bristle characters. The distributions of homozygous effects of inserts on the bristle traits were asymmetrical, with the largest effects in the direction of reducing bristle number; and highly leptokurtic, with most of the increase in variance contributed by a few lines with large effects. The inserts had partially recessive effects on the bristle traits. Insert lines with extreme bristle effects had on average greatly reduced viability.     

An experimental method for evaluating the contribution of deleterious mutations to quantitative trait variation.

Unconditionally deleterious mutations could be an important source of variation in quantitative traits. Deleterious mutations should be rare (segregating at low frequency in the population) and at least partially recessive. In this paper, I suggest that the contribution of rare, partially recessive alleles to quantitative trait variation can be assessed by comparing the relative magnitudes of two genetic variance components: the covariance of additive and homozygous dominance effects (Cad) and the additive genetic variance (Va). If genetic variation is due to rare recessives, then the ratio of Cad to Va should be equal to or greater than 1. In contrast, Cad/Va should be close to zero or even negative if variation is caused by alleles at intermediate frequencies. The ratio of Cad to Va can be estimated from phenotypic comparisons between inbred and outbred relatives, but such estimates are likely to be highly imprecise. Selection experiments provide an alternative estimator for Cad/Va, one with favourable statistical properties. When combined with other biometrical analyses, the ratio test can provide an incisive test of the deleterious mutation model.     

Pleiotropic model of maintenance of quantitative genetic variation at mutation-selection balance

A pleiotropic model of maintenance of quantitative genetic variation at mutation-selection balance is investigated. Mutations have effects on a metric trait and deleterious effects on fitness, for which a bivariate gamma distribution is assumed. Equations for calculating the strength of apparent stabilizing selection (V(s)) and the genetic variance maintained in segregating populations (V(G)) were derived. A large population can hold a high genetic variance but the apparent stabilizing selection may or may not be relatively strong, depending on other properties such as the distribution of mutation effects. If the distribution of mutation effects on fitness is continuous such that there are few nearly neutral mutants, or a minimum fitness effect is assumed if most mutations are nearly neutral, V(G) increases to an asymptote as the population size increases. Both V(G) and V(s) are strongly affected by the shape of the distribution of mutation effects. Compared with mutants of equal effect, allowing their effects on fitness to vary across loci can produce a much higher V(G) but also a high V(s) (V(s) in phenotypic standard deviation units, which is always larger than the ratio V(P)/V(m)), implying weak apparent stabilizing selection. If the mutational variance V(m) is approximately 10(-3)V(e) (V(e), environmental variance), the model can explain typical values of heritability and also apparent stabilizing selection, provided the latter is quite weak as suggested by a recent review.     

The nature of quantittative genetic variation revisited: Lessons from Drosophila bristles

Most characters that distinguish one individual from another, like height or weight, vary continuously in populations. Continuous variation of these ‘quantitative’ traits is due to the simultaneous segregation of multiple quantitative trait loci (QTLs) as well as environmental influences. A major challenge in human medicine, animal and plant breeding and evolutionary genetics is to identify QTLs and determine their genetic properties. Studies of the classic quantitative traits, abdominal and sternopleural bristle numbers of Drosophila, have shown that: (1) many loci have small effects on bristle number, but a few have large effects and cause most of the genetic variation; (2) ‘candidate’ loci involved in bristle development often have large quantitative effects on bristle number; and (3) alleles at QTLs affecting bristle number have variable degrees of dominance, interact with each other, and affect other quantitative traits, including fitness. Lessons learned from this model system will be applicable to studies of the genetic basis of quantitative variation in other species.     

The Nature of Genetic Variation for Complex Traits Revealed by GWAS and Regional Heritability Mapping Analyses

We use computer simulations to investigate the amount of genetic variation for complex traits that can be revealed by single-SNP genome-wide association studies (GWAS) or regional heritability mapping (RHM) analyses based on full genome sequence data or SNP chips. We model a large population subject to mutation, recombination, selection, and drift, assuming a pleiotropic model of mutations sampled from a bivariate distribution of effects of mutations on a quantitative trait and fitness. The pleiotropic model investigated, in contrast to previous models, implies that common mutations of large effect are responsible for most of the genetic variation for quantitative traits, except when the trait is fitness itself. We show that GWAS applied to the full sequence increases the number of QTL detected by as much as 50% compared to the number found with SNP chips but only modestly increases the amount of additive genetic variance explained. Even with full sequence data, the total amount of additive variance explained is generally below 50%. Using RHM on the full sequence data, a slightly larger number of QTL are detected than by GWAS if the same probability threshold is assumed, but these QTL explain a slightly smaller amount of genetic variance. Our results also suggest that most of the missing heritability is due to the inability to detect variants of moderate effect (∼0.03–0.3 phenotypic SDs) segregating at substantial frequencies. Very rare variants, which are more difficult to detect by GWAS, are expected to contribute little genetic variation, so their eventual detection is less relevant for resolving the missing heritability problem.     

Polygenic Mutation in Drosophila Melanogaster: Genetic Analysis of Selection Lines

We have conducted genetic analyses of 12 long-term selection lines of Drosophila melanogaster derived from a highly inbred base population, containing new mutations affecting abdominal and sternopleural bristle number. Biometric analysis of the number of effective factors differentiating the selected lines from the base inbred indicated that with the exception of the three lines selected for increased number of abdominal bristles, three or more mutations contributed to the responses of the selection lines. Analysis of the chromosomal distribution of effects revealed that mutations affecting abdominal bristle number occurred on all three major chromosomes. In addition, Y-linked mutations with effects ranging from one to three bristles occurred in all three lines selected for decreased number of abdominal bristles, as well as in one line selected for increased abdominal bristle number. Mutations affecting sternopleural bristle number were mainly on the X and third chromosomes. One abdominal and one sternopleural selection line showed evidence of a segregating lethal with large effects on bristle number. As an indirect test for allelism of mutations occurring in different selection lines, the three lines selected in the same direction for the same trait were crossed in all possible combinations, and selection continued from the F(2) hybrids. Responses of the hybrid lines usually did not exceed those of the most extreme parental lines, indicating that the responses of the parental lines may have been partly due to mutations at the same loci, although other interpretations are possible.     

Dominance Genetic Variance for Traits Under Directional Selection in Drosophila serrata

In contrast to our growing understanding of patterns of additive genetic variance in single- and multi-trait combinations, the relative contribution of nonadditive genetic variance, particularly dominance variance, to multivariate phenotypes is largely unknown. While mechanisms for the evolution of dominance genetic variance have been, and to some degree remain, subject to debate, the pervasiveness of dominance is widely recognized and may play a key role in several evolutionary processes. Theoretical and empirical evidence suggests that the contribution of dominance variance to phenotypic variance may increase with the correlation between a trait and fitness; however, direct tests of this hypothesis are few. Using a multigenerational breeding design in an unmanipulated population of Drosophila serrata, we estimated additive and dominance genetic covariance matrices for multivariate wing-shape phenotypes, together with a comprehensive measure of fitness, to determine whether there is an association between directional selection and dominance variance. Fitness, a trait unequivocally under directional selection, had no detectable additive genetic variance, but significant dominance genetic variance contributing 32% of the phenotypic variance. For single and multivariate morphological traits, however, no relationship was observed between trait–fitness correlations and dominance variance. A similar proportion of additive and dominance variance was found to contribute to phenotypic variance for single traits, and double the amount of additive compared to dominance variance was found for the multivariate trait combination under directional selection. These data suggest that for many fitness components a positive association between directional selection and dominance genetic variance may not be expected.     

Genetic Architecture of Two Fitness-related Traits in Drosophila melanogaster: Ovariole Number and Thorax Length

In Drosophila melanogaster, ovariole number and thorax length are morphological characters thought to be associated with fitness. Maximum daily egg production in females is positively correlated with ovariole number, while thorax length is correlated with male reproductive success and female fecundity. Though both traits are related to fitness, ovariole number is likely to be under stabilizing selection, while thorax length appears to be under directional selection. Current research has focused on examining the sources of variation for ovariole number in relation to fitness, with a view towards elucidating how segregating variation is maintained in natural populations. Here, we utilize a diallel design to explore the genetic architecture of ovariole number and thorax length in nine isogenic lines derived from a natural population. The full diallel design allows the estimation of general combining ability (GCA), specific combining ability (SCA), and also describes variation due to reciprocal effects (RGCA and RSCA). Ovariole number and thorax length differed with respect to their genetic architecture, reflective of the independent selective forces acting on the traits. For ovariole number, GCA accounted for the majority (67.3%) of variation segregating between the lines, with no evidence of reciprocal effects or inbreeding depression; SCA accounted for a small percentage (3.9%) of the variance, suggesting dominance variation; no reciprocal effects were observed. In contrast, for thorax length, the majority of the non-error variance was accounted for by SCA (17.9%), with only one third as much variance (6.2%) due to GCA. Interestingly, RSCA (nuclear–extranuclear interactions) accounted for slightly more variation (7.5%) than GCA in these data. Thus, genetic variation for thorax length is largely in accord with predictions for a fitness trait under directional selection: little additive genetic variation and substantial dominance variation (including a suggestion of inbreeding depression); while the mechanisms underlying the maintenance of variation for ovariole number are more complex.     

The Response to Artificial Selection from New Mutations in Drosophila Melanogaster

Twenty generations of divergent selection for abdominal bristle number were carried out starting from a completely homozygous population of Drosophila melanogaster. All lines were selected with the same proportion (20%) but at two different numbers of selected parents of each sex (5 or 25). A significant response to selection was detected in eight lines (out of 40) and, in most cases, it could be wholly attributed to a single mutation of relatively large effect (0.5-2 phenotypic standard deviations). The ratio of new mutational variance to environmental variance was estimated to be (0.33 +/- 0.11) X 10(-3). The distribution of mutant effects was asymmetrical, both with respect to bristle number (85% of it was negative) and to fitness (most detected bristle mutations were lethal or semilethal). Moreover, this distribution was leptokurtic, due to the presence of major genes. Gene action on bristles ranged from additive to completely recessive, no epistatic interactions being found. In agreement with theory, larger responses in each direction were achieved by those lines selected at greater effective population sizes. Furthermore, the observed divergence between lines selected in opposite directions was proportional to their effective size, as predicted for mutations of large effect.     

Redistribution of gene frequency and changes of genetic variation following a bottleneck in population size

Although the distribution of frequencies of genes influencing quantitative traits is important to our understanding of their genetic basis and their evolution, direct information from laboratory experiments is very limited. In theory, different models of selection and mutation generate different predictions of frequency distributions. When a large population at mutation-selection balance passes through a rapid bottleneck in size, the frequency distribution of genes is dramatically altered, causing changes in observable quantities such as the mean and variance of quantitative traits. We investigate the gene frequency distribution of a population at mutation-selection balance under a joint-effect model of real stabilizing and pleiotropic selection and its redistribution and thus changes of the genetic properties of metric and fitness traits after the population passes a rapid bottleneck and expands in size. If all genes that affect the trait are neutral with respect to fitness, the additive genetic variance (VA) is always reduced by a bottleneck in population size, regardless of their degree of dominance. For genes that have been under selection, VA increases following a bottleneck if they are (partially) recessive, while the dominance variance increases substantially for any degree of dominance. With typical estimates of mutation parameters, the joint-effect model can explain data from laboratory experiments on the effect of bottlenecking on fitness and morphological traits, providing further support for it as a plausible mechanism for maintenance of quantitative genetic variation.     

Heterosis and outbreeding depression in crosses between natural populations of Arabidopsis thaliana

Understanding the causes and architecture of genetic differentiation between natural populations is of central importance in evolutionary biology. Crosses between natural populations can result in heterosis if recessive or nearly recessive deleterious mutations have become fixed within populations because of genetic drift. Divergence between populations can also result in outbreeding depression because of genetic incompatibilities. The net fitness consequences of between-population crosses will be a balance between heterosis and outbreeding depression. We estimated the magnitude of heterosis and outbreeding depression in the highly selfing model plant Arabidopsis thaliana, by crossing replicate line pairs from two sets of natural populations (C↔R, B↔S) separated by similar geographic distances (Italy↔Sweden). We examined the contribution of different modes of gene action to overall differences in estimates of lifetime fitness and fitness components using joint scaling tests with parental, reciprocal F₁ and F₂, and backcross lines. One of these population pairs (C↔R) was previously demonstrated to be locally adapted, but locally maladaptive quantitative trait loci were also found, suggesting a role for genetic drift in shaping adaptive variation. We found markedly different genetic architectures for fitness and fitness components in the two sets of populations. In one (C↔R), there were consistently positive effects of dominance, indicating the masking of recessive or nearly recessive deleterious mutations that had become fixed by genetic drift. The other set (B↔S) exhibited outbreeding depression because of negative dominance effects. Additional studies are needed to explore the molecular genetic basis of heterosis and outbreeding depression, and how their magnitudes vary across environments.     

Polygenic Mutation in Drosophila Melanogaster: The Causal Relationship of Bristle Number to Fitness

The association between sternopleural and abdominal bristle number and fitness in Drosophila melanogaster was determined for sublines of an initially highly inbred strain that were maintained by divergent artificial selection for 150 generations or by random mating for 180 generations. Replicate selection lines had more extreme bristle numbers than those that were maintained without artificial selection at the same census size for approximately the same number of generations. The average fitness, estimated by a single generation of competition against a compound autosome strain, was 0.17 for lines selected for high and low abdominal bristle numbers and 0.19 for lines selected for high and low sternopleural bristle number. The average fitness of unselected lines, 0.46, was significantly higher than that of the selection lines. The fitnesses and the relationships of bristle number to fitness in progeny of all possible crosses of high X high (H X H), high X low (H X L) and low X low (L X L) selection lines were examined to determine whether the observed intermediate optima were caused by direct stabilizing selection on bristle number or by apparent stabilizing selection mediated through deleterious pleiotropic fitness effects of mutations affecting bristle number. Although bristle number was nearly additive for progeny of H X H, H X L and L X L crosses among sternopleural bristle selection lines, their mean fitnesses were not significantly different from each other, or from the mean fitness of the unselected lines, suggesting partly or completely recessive pleiotropic fitness effects cause apparent stabilizing selection. The average fitness of the progeny of H X H abdominal bristle selection lines was not significantly different from the fitness of unselected lines, but the mean fitness of the progeny of L X L crosses was not significantly different from that of the pure low lines. This is consistent with direct selection against low but not high abdominal bristle number, but the interpretation is confounded by variation in average degree of dominance for fitness (on average recessive in the high abdominal bristle selection lines and additive in the low abdominal bristle selection lines). Neither direct stabilizing selection nor pleiotropy, therefore, can account for all the observations.     

Deleterious Mutations, Apparent Stabilizing Selection and the Maintenance of Quantitative Variation

Apparent stabilizing selection on a quantitative trait that is not causally connected to fitness can result from the pleiotropic effects of unconditionally deleterious mutations, because as N. Barton noted, "...individuals with extreme values of the trait will tend to carry more deleterious alleles...." We use a simple model to investigate the dependence of this apparent selection on the genomic deleterious mutation rate, U; the equilibrium distribution of K, the number of deleterious mutations per genome; and the parameters describing directional selection against deleterious mutations. Unlike previous analyses, we allow for epistatic selection against deleterious alleles. For various selection functions and realistic parameter values, the distribution of K, the distribution of breeding values for a pleiotropically affected trait, and the apparent stabilizing selection function are all nearly Gaussian. The additive genetic variance for the quantitative trait is kQa2, where k is the average number of deleterious mutations per genome, Q is the proportion of deleterious mutations that affect the trait, and a2 is the variance of pleiotropic effects for individual mutations that do affect the trait. In contrast, when the trait is measured in units of its additive standard deviation, the apparent fitness function is essentially independent of Q and a2; and beta, the intensity of selection, measured as the ratio of additive genetic variance to the "variance" of the fitness curve, is very close to s = U/k, the selection coefficient against individual deleterious mutations at equilibrium. Therefore, this model predicts appreciable apparent stabilizing selection if s exceeds about 0.03, which is consistent with various data. However, the model also predicts that beta must equal Vm/VG, the ratio of new additive variance for the trait introduced each generation by mutation to the standing additive variance. Most, although not all, estimates of this ratio imply apparent stabilizing selection weaker than generally observed. A qualitative argument suggests that even when direct selection is responsible for most of the selection observed on a character, it may be essentially irrelevant to the maintenance of variation for the character by mutation-selection balance. Simple experiments can indicate the fraction of observed stabilizing selection attributable to the pleiotropic effects of deleterious mutations.     

Quantitative Genetics of Sperm Precedence in Drosophila Melanogaster

To assess the genetic basis of sperm competition under conditions in which it occurs, I estimated additive, dominance, homozygous and environmental variance components, the effects of inbreeding, and the weighted average dominance of segregating alleles for two measures of sperm precedence in a large, outbred laboratory population. Both first and second male precedence show significant decline on inbreeding. Second male precedence demonstrates significant dominance variance and homozygous genetic variance, but the additive variance is low and not significantly different from zero. For first male precedence, the variance among homozygous lines is again significant, and dominance variance is larger than the additive variance, but is not statistically significant. In contrast, male mating success and other fitness components in Drosophila generally exhibit significant additive variance and little or no dominance variance. Other recent experiments have shown significant genotypic variation for sperm precedence and have associated it with allelic variants of accessory-gland proteins. The contrast between sperm precedence and other male fitness traits in the structure of quantitative genetic variation suggests that different mechanisms may be responsible for the maintenance of variation in these traits. The pattern of genetic variation and inbreeding decline shown in this experiment suggests that one or a few genes with major effects on sperm precedence may be segregating in this population.     

Effect of E(sev) and Su(Raf) Hsp83 mutants and trans-heterozygotes on bristle trait means and variation in Drosophila melanogaster

The Hsp90 protein encoded by the Hsp83 gene is required for the development of many traits in Drosophila. Hsp83 is also thought to play a role in the expression of phenotypic and genetic variability for subsequent selection and evolutionary change. Here we examine the impact of different E(sev) and Su(Raf) Hsp83 mutants on means and phenotypic variances of invariant and variable bristle traits. One of the mutants influenced the normally invariant thoracic bristle number, while none affected invariant scutellar bristle number. E(sev) alleles consistently influenced variable bristle traits while there were fewer effects of the Su(Raf) alleles. For the variable traits, none of the Hsp83 alleles had any effect on phenotypic variance, environmental variance, or developmental stability of the bristle traits. When alleles were combined in trans-heterozygotes, there were both cumulative and complementary effects on thoracic and variable bristle trait numbers, depending on the allelic combination. Overall, the results suggest that Hsp83 mutants do not have detectable effects on the phenotypic or environmental variance of bristle traits and that complementation of E(sev) and Su(Raf) Hsp83 mutants can extend to thoracic bristles as well as previously reported effects on viability. Some allelic combinations lead to more severe effects on variable bristle trait means than do single Hsp83 mutations.     

Epistatic interactions among herbicide resistances in Arabidopsis thaliana: the fitness cost of multiresistance

The type of interactions among deleterious mutations is considered to be crucial in numerous areas of evolutionary biology, including the evolution of sex and recombination, the evolution of ploidy, the evolution of selfing, and the conservation of small populations. Because the herbicide resistance genes could be viewed as slightly deleterious mutations in the absence of the pesticide selection pressure, the epistatic interactions among three herbicide resistance genes (acetolactate synthase CSR, cellulose synthase IXR1, and auxin-induced AXR1 target genes) were estimated in both the homozygous and the heterozygous states, giving 27 genotype combinations in the model plant Arabidopsis thaliana. By analyzing eight quantitative traits in a segregating population for the three herbicide resistances in the absence of herbicide, we found that most interactions in both the homozygous and the heterozygous states were best explained by multiplicative effects (each additional resistance gene causes a comparable reduction in fitness) rather than by synergistic effects (each additional resistance gene causes a disproportionate fitness reduction). Dominance coefficients of the herbicide resistance cost ranged from partial dominance to underdominance, with a mean dominance coefficient of 0.07. It was suggested that the csr1-1, ixr1-2, and axr1-3 resistance alleles are nearly fully recessive for the fitness cost. More interestingly, the dominance of a specific resistance gene in the absence of herbicide varied according to, first, the presence of the other resistance genes and, second, the quantitative trait analyzed. These results and their implications for multiresistance evolution are discussed in relation to the maintenance of polymorphism at resistance loci in a heterogeneous environment.