A network model for the correlation between epistasis and genomic complexity

The study of genetic interactions (epistasis) is central to the understanding of genome organization and evolution. A general correlation between epistasis and genomic complexity has been recently shown, such that in simpler genomes epistasis is antagonistic on average (mutational effects tend to cancel each other out), whereas a transition towards synergistic epistasis occurs in more complex genomes (mutational effects strengthen each other). Here, we use a simple network model to identify basic features explaining this correlation. We show that, in small networks with multifunctional nodes, lack of redundancy, and absence of alternative pathways, epistasis is antagonistic on average. In contrast, lack of multi-functionality, high connectivity, and redundancy favor synergistic epistasis. Moreover, we confirm the previous finding that epistasis is a covariate of mutational robustness: in less robust networks it tends to be antagonistic whereas in more robust networks it tends to be synergistic. We argue that network features associated with antagonistic epistasis are typically found in simple genomes, such as those of viruses and bacteria, whereas the features associated with synergistic epistasis are more extensively exploited by higher eukaryotes.     

The causes of epistasis

Since Bateson's discovery that genes can suppress the phenotypic effects of other genes, gene interactions-called epistasis-have been the topic of a vast research effort. Systems and developmental biologists study epistasis to understand the genotype-phenotype map, whereas evolutionary biologists recognize the fundamental importance of epistasis for evolution. Depending on its form, epistasis may lead to divergence and speciation, provide evolutionary benefits to sex and affect the robustness and evolvability of organisms. That epistasis can itself be shaped by evolution has only recently been realized. Here, we review the empirical pattern of epistasis, and some of the factors that may affect the form and extent of epistasis. Based on their divergent consequences, we distinguish between interactions with or without mean effect, and those affecting the magnitude of fitness effects or their sign. Empirical work has begun to quantify epistasis in multiple dimensions in the context of metabolic and fitness landscape models. We discuss possible proximate causes (such as protein function and metabolic networks) and ultimate factors (including mutation, recombination, and the importance of natural selection and genetic drift). We conclude that, in general, pleiotropy is an important prerequisite for epistasis, and that epistasis may evolve as an adaptive or intrinsic consequence of changes in genetic robustness and evolvability.     

Influence of Gene Interaction on Complex Trait Variation with Multilocus Models

Although research effort is being expended into determining the importance of epistasis and epistatic variance for complex traits, there is considerable controversy about their importance. Here we undertake an analysis for quantitative traits utilizing a range of multilocus quantitative genetic models and gene frequency distributions, focusing on the potential magnitude of the epistatic variance. All the epistatic terms involving a particular locus appear in its average effect, with the number of two-locus interaction terms increasing in proportion to the square of the number of loci and that of third order as the cube and so on. Hence multilocus epistasis makes substantial contributions to the additive variance and does not, per se, lead to large increases in the nonadditive part of the genotypic variance. Even though this proportion can be high where epistasis is antagonistic to direct effects, it reduces with multiple loci. As the magnitude of the epistatic variance depends critically on the heterozygosity, for models where frequencies are widely dispersed, such as for selectively neutral mutations, contributions of epistatic variance are always small. Epistasis may be important in understanding the genetic architecture, for example, of function or human disease, but that does not imply that loci exhibiting it will contribute much genetic variance. Overall we conclude that theoretical predictions and experimental observations of low amounts of epistatic variance in outbred populations are concordant. It is not a likely source of missing heritability, for example, or major influence on predictions of rates of evolution.     

The causes of epistasis in genetic networks

Epistasis refers to the nonadditive interactions between genes in determining phenotypes. Considerable efforts have shown that, even for a given organism, epistasis may vary both in intensity and sign. Recent comparative studies supported that the overall sign of epistasis switches from positive to negative as the complexity of an organism increases, and it has been hypothesized that this change shall be a consequence of the underlying gene network properties. Why should this be the case? What characteristics of genetic networks determine the sign of epistasis? Here we show, by evolving genetic networks that differ in their complexity and robustness against perturbations but that perform the same tasks, that robustness increased with complexity and that epistasis was positive for small nonrobust networks but negative for large robust ones. Our results indicate that robustness and negative epistasis emerge as a consequence of the existence of redundant elements in regulatory structures of genetic networks and that the correlation between complexity and epistasis is a byproduct of such redundancy, allowing for the decoupling of epistasis from the underlying network complexity.     

Travelling the world of gene-gene interactions

Over the last few years, main effect genetic association analysis has proven to be a successful tool to unravel genetic risk components to a variety of complex diseases. In the quest for disease susceptibility factors and the search for the 'missing heritability', supplementary and complementary efforts have been undertaken. These include the inclusion of several genetic inheritance assumptions in model development, the consideration of different sources of information, and the acknowledgement of disease underlying pathways of networks. The search for epistasis or gene-gene interaction effects on traits of interest is marked by an exponential growth, not only in terms of methodological development, but also in terms of practical applications, translation of statistical epistasis to biological epistasis and integration of omics information sources. The current popularity of the field, as well as its attraction to interdisciplinary teams, each making valuable contributions with sometimes rather unique viewpoints, renders it impossible to give an exhaustive review of to-date available approaches for epistasis screening. The purpose of this work is to give a perspective view on a selection of currently active analysis strategies and concerns in the context of epistasis detection, and to provide an eye to the future of gene-gene interaction analysis.     

A Developmental Systems Perspective on Epistasis: Computational Exploration of Mutational Interactions in Model Developmental Regulatory Networks

The way in which the information contained in genotypes is translated into complex phenotypic traits (i.e. embryonic expression patterns) depends on its decoding by a multilayered hierarchy of biomolecular systems (regulatory networks). Each layer of this hierarchy displays its own regulatory schemes (i.e. operational rules such as +/− feedback) and associated control parameters, resulting in characteristic variational constraints. This process can be conceptualized as a mapping issue, and in the context of highly-dimensional genotype-phenotype mappings (GPMs) epistatic events have been shown to be ubiquitous, manifested in non-linear correspondences between changes in the genotype and their phenotypic effects. In this study I concentrate on epistatic phenomena pervading levels of biological organization above the genetic material, more specifically the realm of molecular networks. At this level, systems approaches to studying GPMs are specially suitable to shed light on the mechanistic basis of epistatic phenomena. To this aim, I constructed and analyzed ensembles of highly-modular (fully interconnected) networks with distinctive topologies, each displaying dynamic behaviors that were categorized as either arbitrary or functional according to early patterning processes in the Drosophila embryo. Spatio-temporal expression trajectories in virtual syncytial embryos were simulated via reaction-diffusion models. My in silico mutational experiments show that: 1) the average fitness decay tendency to successively accumulated mutations in ensembles of functional networks indicates the prevalence of positive epistasis, whereas in ensembles of arbitrary networks negative epistasis is the dominant tendency; and 2) the evaluation of epistatic coefficients of diverse interaction orders indicates that, both positive and negative epistasis are more prevalent in functional networks than in arbitrary ones. Overall, I conclude that the phenotypic and fitness effects of multiple perturbations are strongly conditioned by both the regulatory architecture (i.e. pattern of coupled feedback structures) and the dynamic nature of the spatio-temporal expression trajectories displayed by the simulated networks.     

The genetic architecture of disease resistance in plants and the maintenance of recombination by parasites

Parasites represent strong selection on host populations because they are ubiquitous and can drastically reduce host fitness. It has been hypothesized that parasite selection could explain the widespread occurrence of recombination because it is a coevolving force that favours new genetic combinations in the host. A review of deterministic models for the maintenance of recombination reveals that for recombination to be favoured, multiple genes that interact with each other must be under selection. To evaluate whether parasite selection can explain the maintenance of recombination, we review 85 studies that investigated the genetic architecture of plant disease resistance and discuss whether they conform to the requirements that emerge from theoretical models. General characteristics of disease resistance in plants and problems in evaluating resistance experimentally are also discussed. We found strong evidence that disease resistance in plants is determined by multiple loci. Furthermore, in most cases where loci were tested for interactions, epistasis between loci that affect resistance was found. However, we found weak support for the idea that specific allelic combinations determine resistance to different host genotypes and there was little data on whether epistasis between resistance genes is negative or positive. Thus, the current data indicate that it is possible that parasite selection can favour recombination, but more studies in natural populations that specifically address the nature of the interactions between resistance genes are necessary. The data summarized here suggest that disease resistance is a complex trait and that environmental effects and fitness trade-offs should be considered in future models of the coevolutionary dynamics of host and parasites.     

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.     

Data and theory point to mainly additive genetic variance for complex traits

The relative proportion of additive and non-additive variation for complex traits is important in evolutionary biology, medicine, and agriculture. We address a long-standing controversy and paradox about the contribution of non-additive genetic variation, namely that knowledge about biological pathways and gene networks imply that epistasis is important. Yet empirical data across a range of traits and species imply that most genetic variance is additive. We evaluate the evidence from empirical studies of genetic variance components and find that additive variance typically accounts for over half, and often close to 100%, of the total genetic variance. We present new theoretical results, based upon the distribution of allele frequencies under neutral and other population genetic models, that show why this is the case even if there are non-additive effects at the level of gene action. We conclude that interactions at the level of genes are not likely to generate much interaction at the level of variance.     

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.     

Predicting epistasis: an experimental test of metabolic control theory with bacterial transcription and translation

Epistatic interactions between mutations are thought to play a crucial role in a number of evolutionary processes, including adaptation and sex. Evidence for epistasis is abundant, but tests of general theoretical models that can predict epistasis are lacking. In this study, I test the ability of metabolic control theory to predict epistasis using a novel experimental approach that combines phenotypic and genetic perturbations of enzymes involved in gene expression and protein synthesis in the bacterium Pseudomonas aeruginosa. These experiments provide experimental support for two key predictions of metabolic control theory: (i) epistasis between genes involved in the same pathway is antagonistic; (ii) epistasis becomes increasingly antagonistic as mutational severity increases. Metabolic control theory is a general theory that applies to any set of genes that are involved in the same linear processing chain, not just metabolic pathways, and I argue that this theory is likely to have important implications for predicting epistasis between functionally coupled genes, such as those involved in antibiotic resistance. Finally, this study highlights the fact that phenotypic manipulations of gene activity provide a powerful method for studying epistasis that complements existing genetic methods.     

MULTIDIMENSIONAL ADAPTIVE EVOLUTION OF A FEED‐FORWARD NETWORK AND THE ILLUSION OF COMPENSATION

When multiple substitutions affect a trait in opposing ways, they are often assumed to be compensatory, not only with respect to the trait, but also with respect to fitness. This type of compensatory evolution has been suggested to underlie the evolution of protein structures and interactions, RNA secondary structures, and gene regulatory modules and networks. The possibility for compensatory evolution results from epistasis. Yet if epistasis is widespread, then it is also possible that the opposing substitutions are individually adaptive. I term this possibility an adaptive reversal. Although possible for arbitrary phenotype‐fitness mappings, it has not yet been investigated whether such epistasis is prevalent in a biologically realistic setting. I investigate a particular regulatory circuit, the type I coherent feed‐forward loop, which is ubiquitous in natural systems and is accurately described by a simple mathematical model. I show that such reversals are common during adaptive evolution, can result solely from the topology of the fitness landscape, and can occur even when adaptation follows a modest environmental change and the network was well adapted to the original environment. The possibility of adaptive reversals warrants a systems perspective when interpreting substitution patterns in gene regulatory networks.     

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.     

The role of epistatic gene interactions in the response to selection and the evolution of evolvability

It has been argued that the architecture of the genotype-phenotype map determines evolvability, but few studies have attempted to quantify these effects. In this article we use the multilinear epistatic model to study the effects of different forms of epistasis on the response to directional selection. We derive an analytical prediction for the change in the additive genetic variance, and use individual-based simulations to understand the dynamics of evolvability and the evolution of genetic architecture. This shows that the major determinant for the evolution of the additive variance, and thus the evolvability, is directional epistasis. Positive directional epistasis leads to an acceleration of evolvability, while negative directional epistasis leads to canalization. In contrast, pure non-directional epistasis has little effect on the response to selection. One consequence of this is that the classical epistatic variance components, which do not distinguish directional and non-directional effects, are useless as predictors of evolutionary dynamics. The build-up of linkage disequilibrium also has negligible effects. We argue that directional epistasis is likely to have major effects on evolutionary dynamics and should be the focus of empirical studies of epistasis.     

Epistasis interaction of QTL effects as a genetic parameter influencing estimation of the genetic additive effect

Epistasis, an additive-by-additive interaction between quantitative trait loci, has been defined as a deviation from the sum of independent effects of individual genes. Epistasis between QTLs assayed in populations segregating for an entire genome has been found at a frequency close to that expected by chance alone. Recently, epistatic effects have been considered by many researchers as important for complex traits. In order to understand the genetic control of complex traits, it is necessary to clarify additive-by-additive interactions among genes. Herein we compare estimates of a parameter connected with the additive gene action calculated on the basis of two models: a model excluding epistasis and a model with additive-by-additive interaction effects. In this paper two data sets were analysed: 1) 150 barley doubled haploid lines derived from the Steptoe × Morex cross, and 2) 145 DH lines of barley obtained from the Harrington × TR306 cross. The results showed that in cases when the effect of epistasis was different from zero, the coefficient of determination was larger for the model with epistasis than for the one excluding epistasis. These results indicate that epistatic interaction plays an important role in controlling the expression of complex traits.     

Review. Neurogenetic studies of alcohol addiction

Neurogenetic studies of alcohol dependence have relied substantially on genetic animal models, particularly rodents. Studies of inbred strains, selectively bred lines and mutants bearing genes whose function has been targeted for over or under expression are reviewed. Studies focused on gene expression changes are the most recent contributors to this literature, and some genetic effects may work through epigenetic mechanisms. In a few instances, interesting parallels have been revealed between genetic risk in humans and studies in non-human animal models. Future approaches are likely to be increasingly complex.     

A unified model for functional and statistical epistasis and its application in quantitative trait Loci analysis

Interaction between genes, or epistasis, is found to be common and it is a key concept for understanding adaptation and evolution of natural populations, response to selection in breeding programs, and determination of complex disease. Currently, two independent classes of models are used to study epistasis. Statistical models focus on maintaining desired statistical properties for detection and estimation of genetic effects and for the decomposition of genetic variance using average effects of allele substitutions in populations as parameters. Functional models focus on the evolutionary consequences of the attributes of the genotype-phenotype map using natural effects of allele substitutions as parameters. Here we provide a new, general and unified model framework: the natural and orthogonal interactions (NOIA) model. NOIA implements tools for transforming genetic effects measured in one population to the ones of other populations (e.g., between two experimental designs for QTL) and parameters of statistical and functional epistasis into each other (thus enabling us to obtain functional estimates of QTL), as demonstrated numerically. We develop graphical interpretations of functional and statistical models as regressions of the genotypic values on the gene content, which illustrates the difference between the models--the constraint on the slope of the functional regression--and when the models are equivalent. Furthermore, we use our theoretical foundations to conceptually clarify functional and statistical epistasis, discuss the advantages of NOIA over previous theory, and stress the importance of linking functional and statistical models.     

Gene regulatory networks generating the phenomena of additivity, dominance and epistasis

We show how the phenomena of genetic dominance, overdominance, additivity, and epistasis are generic features of simple diploid gene regulatory networks. These regulatory network models are together sufficiently complex to catch most of the suggested molecular mechanisms responsible for generating dominant mutations. These include reduced gene dosage, expression or protein activity (haploinsufficiency), increased gene dosage, ectopic or temporarily altered mRNA expression, increased or constitutive protein activity, and dominant negative effects. As classical genetics regards the phenomenon of dominance to be generated by intralocus interactions, we have studied two one-locus models, one with a negative autoregulatory feedback loop, and one with a positive autoregulatory feedback loop. To include the phenomena of epistasis and downstream regulatory effects, a model of a three-locus signal transduction network is also analyzed. It is found that genetic dominance as well as overdominance may be an intra- as well as interlocus interaction phenomenon. In the latter case the dominance phenomenon is intimately connected to either feedback-mediated epistasis or downstream-mediated epistasis. It appears that in the intra- as well as the interlocus case there is considerable room for additive gene action, which may explain to some degree the predictive power of quantitative genetic theory, with its emphasis on this type of gene action. Furthermore, the results illuminate and reconcile the prevailing explanations of heterosis, and they support the old conjecture that the phenomenon of dominance may have an evolutionary explanation related to life history strategy.     

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.