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.     

Modeling quantitative trait Loci and interpretation of models

A quantitative genetic model relates the genotypic value of an individual to the alleles at the loci that contribute to the variation in a population in terms of additive, dominance, and epistatic effects. This partition of genetic effects is related to the partition of genetic variance. A number of models have been proposed to describe this relationship: some are based on the orthogonal partition of genetic variance in an equilibrium population. We compare a few representative models and discuss their utility and potential problems for analyzing quantitative trait loci (QTL) in a segregating population. An orthogonal model implies that estimates of the genetic effects are consistent in a full or reduced model in an equilibrium population and are directly related to the partition of the genetic variance in the population. Linkage disequilibrium does not affect the estimation of genetic effects in a full model, but would in a reduced model. Certainly linkage disequilibrium would complicate the detection of QTL and epistasis. Using different models does not influence the detection of QTL and epistasis. However, it does influence the estimation and interpretation of genetic effects.     

Dissecting genetic networks underlying complex phenotypes: the theoretical framework

Great progress has been made in genetic dissection of quantitative trait variation during the past two decades, but many studies still reveal only a small fraction of quantitative trait loci (QTLs), and epistasis remains elusive. We integrate contemporary knowledge of signal transduction pathways with principles of quantitative and population genetics to characterize genetic networks underlying complex traits, using a model founded upon one-way functional dependency of downstream genes on upstream regulators (the principle of hierarchy) and mutual functional dependency among related genes (functional genetic units, FGU). Both simulated and real data suggest that complementary epistasis contributes greatly to quantitative trait variation, and obscures the phenotypic effects of many 'downstream' loci in pathways. The mathematical relationships between the main effects and epistatic effects of genes acting at different levels of signaling pathways were established using the quantitative and population genetic parameters. Both loss of function and "co-adapted" gene complexes formed by multiple alleles with differentiated functions (effects) are predicted to be frequent types of allelic diversity at loci that contribute to the genetic variation of complex traits in populations. Downstream FGUs appear to be more vulnerable to loss of function than their upstream regulators, but this vulnerability is apparently compensated by different FGUs of similar functions. Other predictions from the model may account for puzzling results regarding responses to selection, genotype by environment interaction, and the genetic basis of heterosis.     

Epistasis,inbreeding depression,and the evolution of self-fertilization

Inbreeding depression resulting from partially recessive deleterious alleles is thought to be the main genetic factor preventing self-fertilizing mutants from spreading in outcrossing hermaphroditic populations. However, deleterious alleles may also generate an advantage to selfers in terms of more efficient purging, while the effects of epistasis among those alleles on inbreeding depression and mating system evolution remain little explored. In this article, we use a general model of selection to disentangle the effects of different forms of epistasis (additive-by-additive, additive-by-dominance, and dominance-by-dominance) on inbreeding depression and on the strength of selection for selfing. Models with fixed epistasis across loci, and models of stabilizing selection acting on quantitative traits (generating distributions of epistasis) are considered as special cases. Besides its effects on inbreeding depression, epistasis may increase the purging advantage associated with selfing (when it is negative on average), while the variance in epistasis favors selfing through the generation of linkage disequilibria that increase mean fitness. Approximations for the strengths of these effects are derived, and compared with individual-based simulation results.     

A simple model of co-evolutionary dynamics caused by epistatic selection

Epistasis is the dependency of the effect of a mutation on the genetic background in which it occurs. Epistasis has been widely documented and implicated in the evolution of species barriers and the evolution of genetic architecture. Here we propose a simple model to formalize the idea that epistasis can also lead to co-evolutionary patterns in molecular evolution of interacting genes. This model epistasis is represented by the influence of one gene substitution on the fitness rank of the resident allele at another locus. We assume that increasing or decreasing fitness rank occur equally likely. In simulations we show that this form of epistasis leads to co-evolution in the sense that the length of an adaptive walk between interacting loci is highly correlated. This effect is caused by episodes of elevated rate of evolution in both loci simultaneously. We find that the influence of epistasis on these measures of co-evolutionary dynamics is relatively robust to the details of the model. The main factor influencing the correlation in evolutionary rates is the probability that a substitution will have an epistatic effect, but the strength of epistasis or the asymmetry of the initial fitness ranks of the alleles have only a minor effect. We suggest that covariance in rates of evolution among loci could be used to detect epistasis among loci.     

Simulations provide support for the common disease-common variant hypothesis

The success of mapping genes involved in complex diseases, using association or linkage disequilibrium methods, depends heavily on the number and frequency of susceptibility alleles of these genes. These methods will be economically and statistically feasible if common diseases are usually influenced by one or a few susceptibility alleles at each locus (common disease-common variant, CDCV, hypothesis), but not so if there is a high degree of allelic heterogeneity. Here, we use forward-time population simulations to investigate the impact of various genetic and demographic factors on the allelic spectra of human diseases, on the basis of two models proposed by Reich and Lander and by Pritchard. Factors considered are more complex demographies, a finite-allele mutation model, population structure and migration, and interaction between disease susceptibility loci. The conclusion is that the CDCV hypothesis holds and that the phenomenon is caused by transient effects of demography (population expansion). As a result, we devise a multilocus generalization of the Reich and Lander model and demonstrate how interaction between loci with respect to their response to selection may lead to complex effects. We discuss the implications for mapping of complex diseases.     

Epistasis in a model of molecular signal transduction

Biological functions typically involve complex interacting molecular networks, with numerous feedback and regulation loops. How the properties of the system are affected when one, or several of its parts are modified is a question of fundamental interest, with numerous implications for the way we study and understand biological processes and treat diseases. This question can be rephrased in terms of relating genotypes to phenotypes: to what extent does the effect of a genetic variation at one locus depend on genetic variation at all other loci? Systematic quantitative measurements of epistasis – the deviation from additivity in the effect of alleles at different loci – on a given quantitative trait remain a major challenge. Here, we take a complementary approach of studying theoretically the effect of varying multiple parameters in a validated model of molecular signal transduction. To connect with the genotype/phenotype mapping we interpret parameters of the model as different loci with discrete choices of these parameters as alleles, which allows us to systematically examine the dependence of the signaling output – a quantitative trait – on the set of possible allelic combinations. We show quite generally that quantitative traits behave approximately additively (weak epistasis) when alleles correspond to small changes of parameters; epistasis appears as a result of large differences between alleles. When epistasis is relatively strong, it is concentrated in a sparse subset of loci and in low order (e.g. pair-wise) interactions. We find that focusing on interaction between loci that exhibit strong additive effects is an efficient way of identifying most of the epistasis. Our model study defines a theoretical framework for interpretation of experimental data and provides statistical predictions for the structure of genetic interaction expected for moderately complex biological circuits.     

Epistatic pleiotropy and the genetic architecture of covariation within early and late-developing skull trait complexes in mice

The role of epistasis as a source of trait variation is well established, but its role as a source of covariation among traits (i.e., as a source of "epistatic pleiotropy") is rarely considered. In this study we examine the relative importance of epistatic pleiotropy in producing covariation within early and late-developing skull trait complexes in a population of mice derived from an intercross of the Large and Small inbred strains. Significant epistasis was found for several pairwise combinations of the 21 quantitative trait loci (QTL) affecting early developing traits and among the 20 QTL affecting late-developing traits. The majority of the epistatic effects were restricted to single traits but epistatic pleiotropy still contributed significantly to covariances. Because of their proportionally larger effects on variances than on covariances, epistatic effects tended to reduce within-group correlations of traits and reduce their overall degree of integration. The expected contributions of single-locus and two-locus epistatic pleiotropic QTL effects to the genetic covariance between traits were analyzed using a two-locus population genetic model. The model demonstrates that, for single-locus or epistatic pleiotropy to contribute to trait covariances in the study population, both traits must show the same pattern of single-locus or epistatic effects. As a result, a large number of the cases where loci show pleiotropic effects do not contribute to the covariance between traits in this population because the loci show a different pattern of effect on the different traits. In general, covariance patterns produced by single-locus and epistatic pleiotropy predicted by the model agreed well with actual values calculated from the QTL analysis. Nearly all single-locus and epistatic pleiotropic effects contributed positive components to covariances between traits, suggesting that genetic integration in the skull is achieved by a complex combination of pleiotropic effects.     

Complex genetic architecture of population differences in adult lifespan of a beetle: nonadditive inheritance, gender differences, body size and a large maternal effect

Evolutionary responses to selection can be complicated when there is substantial nonadditivity, which limits our ability to extrapolate from simple models of selection to population differentiation and speciation. Studies of Drosophila melanogaster indicate that lifespan and the rate of senescence are influenced by many genes that have environment- and sex-specific effects. These studies also demonstrate that interactions among alleles (dominance) and loci (epistasis) are common, with the degree of interaction differing between the sexes and among environments. However, little is known about the genetic architecture of lifespan or mortality rates for organisms other than D. melanogaster. We studied genetic architecture of differences in lifespan and shapes of mortality curves between two populations of the seed beetle, Callosobruchus maculatus (South India and Burkina Faso populations). These two populations differ in various traits (such as body size and adult lifespan) that have likely evolved via host-specific selection. We found that the genetic architecture of lifespan differences between populations differs substantially between males and females; there was a large maternal effect on male lifespan (but not on female lifespan), and substantial dominance of long-life alleles in females (but not males). The large maternal effect in males was genetically based (there was no significant cytoplasmic effect) likely due to population differences in maternal effects genes that influence lifespan of progeny. Rearing host did not affect the genetic architecture of lifespan, and there was no evidence that genes on the Y-chromosome influence the population differences in lifespan. Epistatic interactions among loci were detectable for the mortality rate of both males and females, but were detectable for lifespan only after controlling for body size variation among lines. The detection of epistasis, dominance, and sex-specific genetic effects on C. maculatus lifespan is consistent with results from line cross and quantitative trait locus studies of D. melanogaster.     

Weak Epistasis Generally Stabilizes Phenotypes in a Mouse Intercross

The extent and strength of epistasis is commonly unresolved in genetic studies, and observed epistasis is often difficult to interpret in terms of biological consequences or overall genetic architecture. We investigated the prevalence and consequences of epistasis by analyzing four body composition phenotypes—body weight, body fat percentage, femoral density, and femoral circumference—in a large F2 intercross of B6-lit/lit and C3.B6-lit/lit mice. We used Combined Analysis of Pleiotropy and Epistasis (CAPE) to examine interactions for the four phenotypes simultaneously, which revealed an extensive directed network of genetic loci interacting with each other, circulating IGF1, and sex to influence these phenotypes. The majority of epistatic interactions had small effects relative to additive effects of individual loci, and tended to stabilize phenotypes towards the mean of the population rather than extremes. Interactive effects of two alleles inherited from one parental strain commonly resulted in phenotypes closer to the population mean than the additive effects from the two loci, and often much closer to the mean than either single-locus model. Alternatively, combinations of alleles inherited from different parent strains contribute to more extreme phenotypes not observed in either parental strain. This class of phenotype-stabilizing interactions has effects that are close to additive and are thus difficult to detect except in very large intercrosses. Nevertheless, we found these interactions to be useful in generating hypotheses for functional relationships between genetic loci. Our findings suggest that while epistasis is often weak and unlikely to account for a large proportion of heritable variance, even small-effect genetic interactions can facilitate hypotheses of underlying biology in well-powered studies.     

A gene's eye view of epistasis,selection and speciation

In this mini‐review, I discuss the effects of gene interaction or epistasis from a `gene's eye view.' By a `gene's eye view' of epistasis, I mean that I will consider a single, bi‐allelic locus, A , whose effects on fitness result only from its interactions with alleles of another, unknown locus, X . I will show how changes in the frequencies of alleles at the background locus affect the relationship of alleles at the A ‐locus to fitness. Changing the genetic background changes the fundamental characteristics of the A ‐locus, such as the magnitude and sign of allelic effects on fitness, and, consequently, it changes the strength and pattern of selection. I consider each of the four kinds of pair–wise interactions between two loci and show that some kinds of epistasis are more sensitive than others to population genetic subdivision. Lastly, I show that some kinds of epistasis are more likely than others to affect the process of speciation and contribute to or be responsible for general genetic features of interspecific hybrids, such as Haldane's rule.     

A framework for detecting and characterizing genetic background-dependent imprinting effects

Genomic imprinting, where the effects of alleles depend on their parent-of-origin, can be an important component of the genetic architecture of complex traits. Although there has been a rapidly increasing number of studies of genetic architecture that have examined imprinting effects, none have examined whether imprinting effects depend on genetic background. Such effects are critical for the evolution of genomic imprinting because they allow the imprinting state of a locus to evolve as a function of genetic background. Here we develop a two-locus model of epistasis that includes epistatic interactions involving imprinting effects and apply this model to scan the mouse genome for loci that modulate the imprinting effects of quantitative trait loci (QTL). The inclusion of imprinting leads to nine orthogonal forms of epistasis, five of which do not appear in the usual two-locus decomposition of epistasis. Each form represents a change in the imprinting status of one locus across different classes of genotypes at the other locus. Our genome scan identified two different locus pairs that show complex patterns of epistasis, where the imprinting effect at one locus changes across genetic backgrounds at the other locus. Thus, our model provides a framework for the detection of genetic background-dependent imprinting effects that should provide insights into the background dependence and evolution of genomic imprinting. Our application of the model to a genome scan supports this assertion by identifying pairs of loci that show reciprocal changes in their imprinting status as the background provided by the other locus changes.     

A Population Genetic Signal of Polygenic Adaptation

Adaptation in response to selection on polygenic phenotypes may occur via subtle allele frequencies shifts at many loci. Current population genomic techniques are not well posed to identify such signals. In the past decade, detailed knowledge about the specific loci underlying polygenic traits has begun to emerge from genome-wide association studies (GWAS). Here we combine this knowledge from GWAS with robust population genetic modeling to identify traits that may have been influenced by local adaptation. We exploit the fact that GWAS provide an estimate of the additive effect size of many loci to estimate the mean additive genetic value for a given phenotype across many populations as simple weighted sums of allele frequencies. We use a general model of neutral genetic value drift for an arbitrary number of populations with an arbitrary relatedness structure. Based on this model, we develop methods for detecting unusually strong correlations between genetic values and specific environmental variables, as well as a generalization of comparisons to test for over-dispersion of genetic values among populations. Finally we lay out a framework to identify the individual populations or groups of populations that contribute to the signal of overdispersion. These tests have considerably greater power than their single locus equivalents due to the fact that they look for positive covariance between like effect alleles, and also significantly outperform methods that do not account for population structure. We apply our tests to the Human Genome Diversity Panel (HGDP) dataset using GWAS data for height, skin pigmentation, type 2 diabetes, body mass index, and two inflammatory bowel disease datasets. This analysis uncovers a number of putative signals of local adaptation, and we discuss the biological interpretation and caveats of these results.     

Gene interactions from maternal effects

Theoretical analyses have demonstrated a potential role for epistasis in many of the most important processes in evolution. These analyses generally assume that an individual's genes map directly to its phenotype and epistasis results from interactions among loci that contribute to the same biochemical or developmental pathways (termed physiological, or within-genotype, epistasis). For many characters, particularly those expressed early in life, an individual's phenotype may also be affected by genes expressed by its parents. The presence of these parental effects allows for interactions between the genes present in the parental and offspring genomes. When the phenotypic effect of a locus in the offspring depends on the alleles possessed by its parents, genotype-by-genotype, or among-genotype, epistasis occurs. The among-genotype epistasis resulting from parental effects may contribute to ruggedness of adaptive landscapes because early mortality often accounts for much of the variance in fitness in populations. To demonstrate how parent-offspring interactions can result in among-genotype epistasis, I use a two-locus model, with one maternal effect locus and one direct effect locus, each with two alleles. Dynamical equations are presented for the two-locus model and are directly contrasted with the dynamical equations derived for a model for physiological epistasis. The relationship between the evolutionary dynamics resulting from these two forms of epistasis is discussed. Three scenarios are presented to illustrate systems in which maternal-offspring, genotype-by-genotype epistasis may occur. The implications of maternal-offspring epistasis for quantitative-trait-loci studies are also discussed.     

Perspective: Sign epistasis and genetic constraint on evolutionary trajectories

Epistasis for fitness means that the selective effect of a mutation is conditional on the genetic background in which it appears. Although epistasis is widely observed in nature, our understanding of its consequences for evolution by natural selection remains incomplete. In particular, much attention focuses only on its influence on the instantaneous rate of changes in frequency of selected alleles via epistatic contribution to the additive genetic variance for fitness. Thus, in this framework epistasis only has evolutionary importance if the interacting loci are simultaneously segregating in the population. However, the selective accessibility of mutational trajectories to high fitness genotypes may depend on the genetic background in which novel mutations appear, and this effect is independent of population polymorphism at other loci. Here we explore this second influence of epistasis on evolution by natural selection. We show that it is the consequence of a particular form of epistasis, which we designate sign epistasis. Sign epistasis means that the sign of the fitness effect of a mutation is under epistatic control; thus, such a mutation is beneficial on some genetic backgrounds and deleterious on others. Recent experimental innovations in microbial systems now permit assessment of the fitness effects of individual mutations on multiple genetic backgrounds. We review this literature and identify many examples of sign epistasis, and we suggest that the implications of these results may generalize to other organisms. These theoretical and empirical considerations imply that strong genetic constraint on the selective accessibility of trajectories to high fitness genotypes may exist and suggest specific areas of investigation for future research.     

Overview of QTL detection in plants and tests for synergistic epistatic interactions

Improvements in the usefulness of QTL analysis arise from better statistical methods applied to the problem, ability to analyze more complex mating designs, and the fitting of less simplified genetic models. Here we review the advantages of different plant mating designs in QTL analysis and conclude that diallel designs have several favorable properties. We then turn to the detection of systematic genome-wide synergistic epistasis. This form of epistasis has important implications from evolutionary (maintenance of sexual reproduction and concealment of cryptic genetic variation) and practical perspectives (response to pyramided favorable alleles). We develop two methods for detecting systematic synergistic epistasis, one based on analyzing interactions between locus effects and predicted individual genotypic values and one based on analyzing pairwise locus interactions. Using the first method we detect synergistic epistasis in a barley and a wheat dataset but not in a maize dataset. We fail to detect synergistic epistasis with the second method. We discuss our results in the light of theoretical questions concerning the mechanisms of synergistic epistasis.     

The Effect of an Experimental Bottleneck upon Quantitative Genetic Variation in the Housefly

Effects of a population bottleneck (founder-flush cycle) upon quantitative genetic variation of morphometric traits were examined in replicated experimental lines of the housefly founded with one, four or 16 pairs of flies. Heritability and additive genetic variances for eight morphometric traits generally increased as a result of the bottleneck, but the pattern of increase among bottleneck sizes differed among traits. Principal axes of the additive genetic correlation matrix for the control line yielded two suites of traits, one associated with general body size and another set largely independent of body size. In the former set containing five of the traits, additive genetic variance was greatest in the bottleneck size of four pairs, whereas in the latter set of two traits the largest additive genetic variance occurred in the smallest bottleneck size of one pair. One trait exhibited changes in additive genetic variance intermediate between these two major responses. These results were inconsistent with models of additive effects of alleles within loci or of additive effects among loci. An observed decline in viability measures and body size in the bottleneck lines also indicated that there was nonadditivity of allelic effects for these traits. Several possible nonadditive models were explored that increased additive genetic variance as a result of a bottleneck. These included a model with complete dominance, a model with overdominance and a model incorporating multiplicative epistasis.     

From classical genetics to quantitative genetics to systems biology: modeling epistasis

Gene expression data has been used in lieu of phenotype in both classical and quantitative genetic settings. These two disciplines have separate approaches to measuring and interpreting epistasis, which is the interaction between alleles at different loci. We propose a framework for estimating and interpreting epistasis from a classical experiment that combines the strengths of each approach. A regression analysis step accommodates the quantitative nature of expression measurements by estimating the effect of gene deletions plus any interaction. Effects are selected by significance such that a reduced model describes each expression trait. We show how the resulting models correspond to specific hierarchical relationships between two regulator genes and a target gene. These relationships are the basic units of genetic pathways and genomic system diagrams. Our approach can be extended to analyze data from a variety of experiments, multiple loci, and multiple environments.     

Statistical modeling of interlocus interactions in a complex disease: rejection of the multiplicative model of epistasis in type 1 diabetes

In general, common diseases do not follow a Mendelian inheritance pattern. To identify disease mechanisms and etiology, their genetic dissection may be assisted by evaluation of linkage in mouse models of human disease. Statistical modeling of multiple-locus linkage data from the nonobese diabetic (NOD) mouse model of type 1 diabetes has previously provided evidence for epistasis between alleles of several Idd (insulin-dependent diabetes) loci. The construction of NOD congenic strains containing selected segments of the diabetes-resistant strain genome allows analysis of the joint effects of alleles of different loci in isolation, without the complication of other segregating Idd loci. In this article, we analyze data from congenic strains carrying two chromosome intervals (a double congenic strain) for two pairs of loci: Idd3 and Idd10 and Idd3 and Idd5. The joint action of both pairs is consistent with models of additivity on either the log odds of the penetrance, or the liability scale, rather than with the previously proposed multiplicative model of epistasis. For Idd3 and Idd5 we would also not reject a model of additivity on the penetrance scale, which might indicate a disease model mediated by more than one pathway leading to beta-cell destruction and development of diabetes. However, there has been confusion between different definitions of interaction or epistasis as used in the biological, statistical, epidemiological, and quantitative and human genetics fields. The degree to which statistical analyses can elucidate underlying biologic mechanisms may be limited and may require prior knowledge of the underlying etiology.