Effects of genetic drift on variance components under a general model of epistasis

We analyze the changes in the mean and variance components of a quantitative trait caused by changes in allele frequencies, concentrating on the effects of genetic drift. We use a general representation of epistasis and dominance that allows an arbitrary relation between genotype and phenotype for any number of diallelic loci. We assume initial and final Hardy-Weinberg and linkage equilibrium in our analyses of drift-induced changes. Random drift generates transient linkage disequilibria that cause correlations between allele frequency fluctuations at different loci. However, we show that these have negligible effects, at least for interactions among small numbers of loci. Our analyses are based on diffusion approximations that summarize the effects of drift in terms of F, the inbreeding coefficient, interpreted as the expected proportional decrease in heterozygosity at each locus. For haploids, the variance of the trait mean after a population bottleneck is var(Delta(z) over bar) = Sigma(k)(=1)(n) (FVA(k))-V-k, where n is the number of loci contributing to the trait variance, V-A(1) = V-A is the additive genetic variance, and V-A(k) is the kth-order additive epistatic variance. The expected additive genetic variance after the bottleneck, denoted <V-A(*)> is closely related to var(Delta(z) over bar); <V-A(*)> = (1 - F) Sigma(k)(=1)(n) kF(k-1) V-A(k). Thus, epistasis inflates the expected additive variance above V-A(1 - F), the expectation under additivity. For haploids (and diploids without dominance), the expected value of every variance component is inflated by the existence of higher order interactions (e.g., third-order epistasis inflates <V-AA(*)>). This is not true in general with diploidy, because dominance alone can reduce <VA(*)> below V-A(1 - F) (e.g., when dominant alleles are rare). Without A dominance, diploidy produces simple expressions: var(Delta(z) over bar) = Sigma(k)(=1)(n) (2F)(k) V-k(A(k)) and <V-A(*)> = (1 - F) Sigma(k)(=1)(n) k(2F)V-k-1(A(k).) k With dominance (and even without epistasis), var(Delta(z) over bar) and no longer depend solely on the variance components in the base population. For small F, the expected additive variance simplifies to <V-A(*)> similar or equal to (1 - F)V-A +4FV(AA) + 2FV(D) + 2FC(AD), where C-AD is a sum of two terms describing covariances between additive effects and dominance and additive X dominance interactions. Whether population bottlenecks lead to expected increases in additive variance depends primarily on the ratio of nonadditive to additive genetic variance in the base population, but dominance precludes simple predictions based solely on variance components. We illustrate these results using a model in which genotypic values are drawn at random, allowing extreme and erratic epistatic interactions. Although our analyses clarify the conditions under which drift is expected to increase V-A, we question the evolutionary importance of such increases.