On the evolution of epistasis I: diploids under selection

One interpretation of recent literature on the evolution of phenotypic modularity is that evolution should act to decrease the degree of interaction between genes that contribute to different phenotypes. This issue is addressed directly here using a fitness scheme determined by two genetic loci and a third locus which modifies a measure of statistical interaction between the fitnesses due to the first two. The equilibrium structure of such an epistasis-modifying locus is studied. It is shown that under well-specified conditions a modifying allele that increases epistasis succeeds. In other words, genetic interactions tend to become stronger. It is speculated that this occurs because the mean fitness in such models is locally increasing as a function of the degree of epistasis.     

On the evolution of epistasis III: the haploid case with mutation

Whether interaction between genes is better represented by synergistic or antagonistic epistasis has been a focus of experimental research in bacterial population genetics. Our previous research on evolution of modifiers of epistasis in diploid systems has indicated that the strength of positive or negative epistasis should increase provided linkage disequilibrium is maintained. Here we study a modifier of epistasis in fitness between two loci in a haploid system. Epistasis is modified in the neighborhood of a mutation-selection balance. We show that when linkage in the three-locus system is tight, an increase in the frequency of a modifier allele that induces either more negative or more positive epistasis is possible. Epistasis here can be measured on either an additive or multiplicative scale.     

The effect of epistasis on sexually antagonistic genetic variation

There is increasing evidence of segregating sexually antagonistic (SA) genetic variation for fitness in laboratory and wild populations, yet the conditions for the maintenance of such variation can be restrictive. Epistatic interactions between genes can contribute to the maintenance of genetic variance in fitness and we suggest that epistasis between SA genes should be pervasive. Here, we explore its effect on SA genetic variation in fitness using a two locus model with negative epistasis. Our results demonstrate that epistasis often increases the parameter space showing polymorphism for SA loci. This is because selection in one locus is affected by allele frequencies at the other, which can act to balance net selection in males and females. Increased linkage between SA loci had more marginal effects. We also show that under some conditions, large portions of the parameter space evolve to a state where male benefit alleles are fixed at one locus and female benefit alleles at the other. This novel effect of epistasis on SA loci, which we term the ‘equity effect’, may have important effects on population differentiation and may contribute to speciation. More generally, these results support the suggestion that epistasis contributes to population divergence.     

Marker-Based Inferences about Epistasis for Genes Influencing Inbreeding Depression

We describe a multilocus, marker-based regression method for inferring interactions between genes controlling inbreeding depression in self-fertile organisms. It is based upon selfing a parent heterozygous for several unlinked codominant markers, then analyzing the fitness of progeny marker genotypes. If loci causing inbreeding depression are linked to marker loci, then viability selection is manifested by distorted segregation of markers, and fecundity selection by dependence of the fecundity character upon the marker genotype. To characterize this selection, fitness is regressed on the proportion of loci homozygous for markers linked to deleterious alleles, and epistasis is detected by nonlinearity of the regression. Alternatively, fitness can be regressed on the proportion of heterozygous loci. Other modes of selection can be incorporated with a bivariate regression involving both homozygote and heterozygote marker genotypes. The advantage of this marker-based approach is that ``purging' is minimized and specific chromosomal segments are identified; its disadvantage lies in low statistical power when linkage is not strong and/or the linkage phase between marker and selected loci is uncertain. Using this method, in the wildflower Mimulus guttatus, we found predominant multiplicative gene interaction determining fecundity and some negative synergistic (nonmultiplicative) interaction for viability.     

Red Queen Dynamics with Non-Standard Fitness Interactions

Antagonistic coevolution between hosts and parasites can involve rapid fluctuations of genotype frequencies that are known as Red Queen dynamics. Under such dynamics, recombination in the hosts may be advantageous because genetic shuffling can quickly produce disproportionately fit offspring (the Red Queen hypothesis). Previous models investigating these dynamics have assumed rather simple models of genetic interactions between hosts and parasites. Here, we assess the robustness of earlier theoretical predictions about the Red Queen with respect to the underlying host-parasite interactions. To this end, we created large numbers of random interaction matrices, analysed the resulting dynamics through simulation, and ascertained whether recombination was favoured or disfavoured. We observed Red Queen dynamics in many of our simulations provided the interaction matrices exhibited sufficient ‘antagonicity’. In agreement with previous studies, strong selection on either hosts or parasites favours selection for increased recombination. However, fast changes in the sign of linkage disequilibrium or epistasis were only infrequently observed and do not appear to be a necessary condition for the Red Queen hypothesis to work. Indeed, recombination was often favoured even though the linkage disequilibrium remained of constant sign throughout the simulations. We conclude that Red Queen-type dynamics involving persistent fluctuations in host and parasite genotype frequencies appear to not be an artefact of specific assumptions about host-parasite fitness interactions, but emerge readily with the general interactions studied here. Our results also indicate that although recombination is often favoured, some of the factors previously thought to be important in this process such as linkage disequilibrium fluctuations need to be reassessed when fitness interactions between hosts and parasites are complex.     

Identifying quantitative trait locus by genetic background interactions in association studies

Association studies are designed to identify main effects of alleles across a potentially wide range of genetic backgrounds. To control for spurious associations, effects of the genetic background itself are often incorporated into the linear model, either in the form of subpopulation effects in the case of structure or in the form of genetic relationship matrices in the case of complex pedigrees. In this context epistatic interactions between loci can be captured as an interaction effect between the associated locus and the genetic background. In this study I developed genetic and statistical models to tie the locus by genetic background interaction idea back to more standard concepts of epistasis when genetic background is modeled using an additive relationship matrix. I also simulated epistatic interactions in four-generation randomly mating pedigrees and evaluated the ability of the statistical models to identify when a biallelic associated locus was epistatic to other loci. Under additive-by-additive epistasis, when interaction effects of the associated locus were quite large (explaining 20% of the phenotypic variance), epistasis was detected in 79% of pedigrees containing 320 individuals. The epistatic model also predicted the genotypic value of progeny better than a standard additive model in 78% of simulations. When interaction effects were smaller (although still fairly large, explaining 5% of the phenotypic variance), epistasis was detected in only 9% of pedigrees containing 320 individuals and the epistatic and additive models were equally effective at predicting the genotypic values of progeny. Epistasis was detected with the same power whether the overall epistatic effect was the result of a single pairwise interaction or the sum of nine pairwise interactions, each generating one ninth of the epistatic variance. The power to detect epistasis was highest (94%) at low QTL minor allele frequency, fell to a minimum (60%) at minor allele frequency of about 0.2, and then plateaued at about 80% as alleles reached intermediate frequencies. The power to detect epistasis declined when the linkage disequilibrium between the DNA marker and the functional polymorphism was not complete.     

A semi-symmetric two-locus model

The two-locus symmetric viability model characterized by its invariance with respect to the exchange of alleles at each locus, is a well-studied model of classical two-locus theory. The symmetric model introduced by Lewontin and Kojima is among the few multi-locus models with epistatic interactions between loci for which a polymorphism with linkage equilibrium can be stable and this happens when recombination is sufficiently large. We show that an analogous property holds true for a different model, in which symmetry need exist at only one locus. The properties of this new semi-symmetric model are compared with those of the classical symmetric model. For tight linkage, two classes of polymorphisms are possible, depending on the magnitude of additive epistasis. The recombination rate above which linkage equilibrium becomes stable is derived analytically. As in the symmetric model, intervals of recombination in which no polymorphism is stable are possible, and stable polymorphisms can coexist with stable fixations.     

Social heterosis and the maintenance of genetic diversity at the genome level

Social heterosis is when individuals in groups or neighbourhoods receive a mutualistic benefit from across‐individual genetic diversity. Although it can be a viable evolutionary mechanism to maintain allelic diversity at a given locus, its efficacy at maintaining genome‐wide diversity is in question when multiple loci are being simultaneously selected. Therefore, we modelled social heterosis in a population of haploid genomes of two‐ or three‐linked loci. With such linkages, social heterosis decreases gametic diversity, but maintains allelic diversity. Genomes tend to survive as complimentary pairs, with alternate alleles at each locus (e.g. the pair AbC and aBc). The outcomes of selection appear similar to fitness epistasis but are novel in the sense that phenotypic interactions occur across rather than within individuals. The model’s results strongly suggest that strong linkage across gene loci actually increases the probability that social heterosis maintains significant genetic diversity at the level of the genome.     

The effects of multilocus balancing selection on neutral variability

We studied the effect of multilocus balancing selection on neutral nucleotide variability at linked sites by simulating a model where diallelic polymorphisms are maintained at an arbitrary number of selected loci by means of symmetric overdominance. Different combinations of alleles define different genetic backgrounds that subdivide the population and strongly affect variability. Several multilocus fitness regimes with different degrees of epistasis and gametic disequilibrium are allowed. Analytical results based on a multilocus extension of the structured coalescent predict that the expected linked neutral diversity increases exponentially with the number of selected loci and can become extremely large. Our simulation results show that although variability increases with the number of genetic backgrounds that are maintained in the population, it is reduced by random fluctuations in the frequencies of those backgrounds and does not reach high levels even in very large populations. We also show that previous results on balancing selection in single-locus systems do not extend to the multilocus scenario in a straightforward way. Different patterns of linkage disequilibrium and of the frequency spectrum of neutral mutations are expected under different degrees of epistasis. Interestingly, the power to detect balancing selection using deviations from a neutral distribution of allele frequencies seems to be diminished under the fitness regime that leads to the largest increase of variability over the neutral case. This and other results are discussed in the light of data from the Mhc.     

The effect of epistasis on the structure of hybrid zones

Within hybrid zones that are maintained by a balance between selection and dispersal, linkage disequilibrium is generated by the mixing of divergent populations. This linkage disequilibrium causes selection on each locus to act on all other loci, thereby steepening clines, and generating a barrier to gene flow. Diffusion models predict simple relations between the strength of linkage disequilibrium and the dispersal rate, sigma, and between the barrier to gene flow, B, and the reduction in mean fitness, W. The aim of this paper is to test the accuracy of these predictions by comparison with an exact deterministic model of unlinked loci (r = 0.5). Disruptive selection acts on the proportion of alleles from the parental populations (p,q): W = exp[-S(4pq)beta], such that the least fit genotype has fitness e-s. Where beta < 1, fitness is reduced for a wide range of intermediate genotypes; where beta > 1, fitness is only reduced for those genotypes close to p = 0.5. Even with strong epistasis, linkage disequilibria are close to sigma 2p'ip'j/rij, where p'i, p'j are the gradients in allele frequency at loci i, j. The barrier to gene flow, which is reflected in the steepening of neutral clines, is given by [formula: see text] where r, the harmonic mean recombination rate between the neural and selected loci, is here 0.5. This is a close approximation for weak selection, but underestimates B for strong selection. The barrier is stronger for small beta, because hybrid fitness is then reduced over a wider range of p. The widths of the selected clines are harder to predict: though simple approximations are accurate for beta = 1, they become inaccurate for extreme beta because, then, fitness changes sharply with p. Estimates of gene number, made from neutral clines on the assumption that selection acts against heterozygotes, are accurate for weak selection when beta = 1; however, for strong selection, gene number is overestimated. For beta > 1, gene number is systematically overestimated and, conversely, when beta < 1, it is underestimated.     

Multilocus Behavior in Random Environments I. Random Levene Models

In this paper the consequences of natural selection acting on several loci simultaneously in a spatially fluctuating environment are described. The fitnesses of the genotypes are assumed to be additive both within and between loci. The environment is assumed to be made up of a very large (effectively infinite) number of patches in which fitnesses are assigned at random. The resulting deterministic model is called a Random Levene Model and its properties are approximate by a system of differential equations. The main equilibrium properites are that (1) the linkage disequilibrium is zero and (2) the correlations in fitness between alleles at different loci are the principle determinants of the dynamic inter-locus interactions. Although there is no epistasis as conventionally defined, the equilibrium state at the two loci are highly interdependent, the governing principle being that two alleles at different loci whose fitness are negatively correlated across environments have a higher overall fitness due to the reduction in their variance in fitness through the negative correlation. When a large number of loci are considered, they naturally fall into correlation groupings which lead to an enhanced likelihood for polymorphism over that predicted by single-locus theory.     

Effect of varying epistasis on the evolution of recombination

Whether recombination decelerates or accelerates a population's response to selection depends, at least in part, on how fitness-determining loci interact. Realistically, all genomes likely contain fitness interactions both with positive and with negative epistasis. Therefore, it is crucial to determine the conditions under which the potential beneficial effects of recombination with negative epistasis prevail over the detrimental effects of recombination with positive epistasis. Here, we examine the simultaneous effects of diverse epistatic interactions with different strengths and signs in a simplified model system with independent pairs of interacting loci and selection acting only on the haploid phase. We find that the average form of epistasis does not predict the average amount of linkage disequilibrium generated or the impact on a recombination modifier when compared to results using the entire distribution of epistatic effects and associated single-mutant effects. Moreover, we show that epistatic interactions of a given strength can produce very different effects, having the greatest impact when selection is weak. In summary, we observe that the evolution of recombination at mutation-selection balance might be driven by a small number of interactions with weak selection rather than by the average epistasis of all interactions. We illustrate this effect with an analysis of published data of Saccharomyces cerevisiae. Thus to draw conclusions on the evolution of recombination from experimental data, it is necessary to consider the distribution of epistatic interactions together with the associated selection coefficients.     

Lewontin and Kojima Meet Fisher: Linkage in a Symmetric Model of Sex Determination

The effect of linkage and epistasis on the evolution of the sex-ratio is studied in a symmetric two-locus model of autosomal sex determination closely related to the symmetric viability model of R. C. Lewontin and K. Kojima. R. A. Fisher's expectation of an even sex ratio for autosomal sex determination by a single gene governs the dynamics when the loci are tightly linked. However, recombination may preclude optimization of the sex ratio just as occurs in viability selection models. Many of the evolutionary phenomena known for the symmetric viability model also occur here. In addition, we exhibit a series of new phenomena related to the presence of surfaces of even sex ratio.     

Impact of epistasis and pleiotropy on evolutionary adaptation

Evolutionary adaptation is often likened to climbing a hill or peak. While this process is simple for fitness landscapes where mutations are independent, the interaction between mutations (epistasis) as well as mutations at loci that affect more than one trait (pleiotropy) are crucial in complex and realistic fitness landscapes. We investigate the impact of epistasis and pleiotropy on adaptive evolution by studying the evolution of a population of asexual haploid organisms (haplotypes) in a model of N interacting loci, where each locus interacts with K other loci. We use a quantitative measure of the magnitude of epistatic interactions between substitutions, and find that it is an increasing function of K. When haplotypes adapt at high mutation rates, more epistatic pairs of substitutions are observed on the line of descent than expected. The highest fitness is attained in landscapes with an intermediate amount of ruggedness that balance the higher fitness potential of interacting genes with their concomitant decreased evolvability. Our findings imply that the synergism between loci that interact epistatically is crucial for evolving genetic modules with high fitness, while too much ruggedness stalls the adaptive process.     

The role of pleiotropy vs signaller-receiver gene epistasis in life history trade-offs: dissecting the genomic architecture of organismal design in social systems

Traditional life history theory ignores trade-offs due to social interactions, yet social systems expand the set of possible trade-offs affecting a species evolution--by introducing asymmetric interactions between the sexes, age classes and invasion of alternative strategies. We outline principles for understanding gene epistasis due to signaller-receiver dynamics, gene interactions between individuals, and impacts on life history trade-offs. Signaller-receiver epistases create trade-offs among multiple correlated traits that affect fitness, and generate multiple fitness optima conditional on frequency of alternative strategies. In such cases, fitness epistasis generated by selection can maintain linkage disequilibrium, even among physically unlinked loci. In reviewing genetic methods for studying life history trade-offs, we conclude that current artificial selection or gene manipulation experiments focus on pleiotropy. Multi-trait selection experiments, multi-gene engineering methods or multiple endocrine manipulations can test for epistasis and circumvent these limitations. In nature, gene mapping in field pedigrees is required to study social gene epistases and associated trade-offs. Moreover, analyses of correlational selection and frequency-dependent selection are necessary to study epistatic social system trade-offs, which can be achieved with group-structured versions of Price's (1970) equation.     

Accounting for epistasis in linkage analysis of general pedigrees

Complex traits are generally believed to be influenced by multiple loci. Identification of loci involved in complex traits is more difficult for interacting than for additive loci. Here we describe an extension of the program lm_twoqtl in the package MORGAN to handle two quantitative trait loci (QTLs) with gene-gene interaction. We investigate whether parametric linkage analysis that accounts for such epistasis improves prospects for linkage detection and accuracy of localization of QTLs. Through use of simulated data we show that analysis that accounts for epistasis provides higher lod scores and better localization than does analysis without epistasis. In addition, we demonstrate that the difference between lod scores in the presence vs. absence of use of an interaction model in analysis is greater in extended than in nuclear pedigrees.     

Evolutionary dynamics, epistatic interactions, and biological information

We investigate a definition of biological information that connects population genetics with the tools of information theory by focusing on the distribution of genotypes found in a population. Previous research has treated loci as non-interacting by making specific approximations in the calculation of information-theoretic quantities. We expand earlier mathematical forms to include epistasis, or interactions between mutations at all pairs of loci. Application of our improved measure of biological information to evolution on two-locus, two-allele fitness landscapes demonstrates that mutual information between loci reflects epistatic interaction of mutations. Finally, we consider four-locus, two-allele fitness landscapes with modular structure. As modular interactions are inherently epistatic, we demonstrate that our refined approximation provides insight into the underlying structure of these non-trivial fitness landscapes.     

Sexual and parental antagonism shape genomic architecture

Populations with two sexes are vulnerable to a pair of genetic conflicts: sexual antagonism that can arise when alleles have opposing fitness effects on females and males; and parental antagonism that arises when alleles have opposing fitness effects when maternally and paternally inherited. This paper extends previous theoretical work that found stable linkage disequilibrium (LD) between sexually antagonistic loci. We find that LD is also generated between parentally antagonistic loci, and between sexually and parentally antagonistic loci, without any requirement of epistasis. We contend that the LD in these models arises from the admixture of gene pools subject to different selective histories. We also find that polymorphism maintained by parental antagonism at one locus expands the opportunity for polymorphism at a linked locus experiencing parental or sexual antagonism. Taken together, our results predict the chromosomal clustering of loci that segregate for sexually and parentally antagonistic alleles. Thus, genetic conflict may play a role in the evolution of genomic architecture.     

Foundations of a mathematical theory of darwinism

This paper pursues the ‘formal darwinism’ project of Grafen, whose aim is to construct formal links between dynamics of gene frequencies and optimization programmes, in very abstract settings with general implications for biologically relevant situations. A major outcome is the definition, within wide assumptions, of the ubiquitous but problematic concept of ‘fitness’. This paper is the first to present the project for mathematicians. Within the framework of overlapping generations in discrete time and no social interactions, the current model shows links between fitness maximization and gene frequency change in a class-structured population, with individual-level uncertainty but no uncertainty in the class projection operator, where individuals are permitted to observe and condition their behaviour on arbitrary parts of the uncertainty. The results hold with arbitrary numbers of loci and alleles, arbitrary dominance and epistasis, and make no assumptions about linkage, linkage disequilibrium or mating system. An explicit derivation is given of Fisher’s Fundamental Theorem of Natural Selection in its full generality.     

Limits to the Relationship among Recombination, Disequilibrium and Epistasis in Two-Locus Models

I determine limits to the equilibrium relationship among epistasis, recombination and disequilibrium in two-locus, two-allele models using linear programming techniques. I show that when allele frequencies are one-half at each locus, the symmetric model is the fitness pattern that generates the most disequilibrium for the smallest level of epistasis. When allele frequencies deviate from one-half much larger levels of epistasis are required to generate similar levels of disequilibrium. I determine the level of epistasis required to generate observed significant levels of disequilibrium in natural populations. The overall implication is that disequilibrium will be large at equilibrium only between strongly interacting, closely linked loci.