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.     

Epistasis Creates Invariant Sites and Modulates the Rate of Molecular Evolution

Invariant sites are a common feature of amino acid sequence evolution. The presence of invariant sites is frequently attributed to the need to preserve function through site-specific conservation of amino acid residues. Amino acid substitution models without a provision for invariant sites often fit the data significantly worse than those that allow for an excess of invariant sites beyond those predicted by models that only incorporate rate variation among sites (e.g., a Gamma distribution). An alternative is epistasis between sites to preserve residue interactions that can create invariant sites. Through computer-simulated sequence evolution, we evaluated the relative effects of site-specific preferences and site-site couplings in the generation of invariant sites and the modulation of the rate of molecular evolution. In an analysis of ten major families of protein domains with diverse sequence and functional properties, we find that the negative selection imposed by epistasis creates many more invariant sites than site-specific residue preferences alone. Further, epistasis plays an increasingly larger role in creating invariant sites over longer evolutionary periods. Epistasis also dictates rates of domain evolution over time by exerting significant additional purifying selection to preserve site couplings. These patterns illuminate the mechanistic role of epistasis in the processes underlying observed site invariance and evolutionary rates.     

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.     

Influence of metabolic network structure and function on enzyme evolution

Background  

Most studies of molecular evolution are focused on individual genes and proteins. However, understanding the design principles and evolutionary properties of molecular networks requires a system-wide perspective. In the present work we connect molecular evolution on the gene level with system properties of a cellular metabolic network. In contrast to protein interaction networks, where several previous studies investigated the molecular evolution of proteins, metabolic networks have a relatively well-defined global function. The ability to consider fluxes in a metabolic network allows us to relate the functional role of each enzyme in a network to its rate of evolution.     

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 intramolecular chaperone-mediated protein folding

Some proteins have evolved to contain a specific sequence as an intramolecular chaperone, which is essential for protein folding but not required for protein function, as it is removed after the protein is folded by autoprocessing or by an exogenous protease. To date, a large number of sequences encoded as N-terminal or C-terminal extensions have been identified to function as intramolecular chaperones. An increasing amount of evidence has revealed that these intramolecular chaperones play an important role in protein folding both in vivo and in vitro. Here, we summarize recent studies on intramolecular chaperone-assisted protein folding and discuss the mechanisms as to how intramolecular chaperones play roles in protein folding.     

The length scale of selection in protein evolution

Central to the study of molecular evolution, and an area of long-standing debate, is the appropriate model for the fitness landscape of proteins. Much of this debate has focused on the strength and frequency of positive and purifying selection, but the form and frequency of selective correlations is also a vital element. The constituent amino acids within a protein generically interact and share selective pressures in predictable ways, which conflicts with the selective independence assumed by common caricatures of the fitness landscape. Here, I discuss a recent study by myself and coauthors that used whole-genome comparisons of orthologous molecular sequences from closely related Drosophilids to explore the form of the selective correlations and selective interactions (epistasis) between the amino acids within a protein. I outline our results and highlight our finding of a selective length scale of ten amino acids within which individual amino acids are substantially and generically more likely to share selective pressures and interact epistatically. I then focus on the evidence presented in our study supporting a substantial role for epistasis in the process of molecular evolution, and discuss further the implications of this widespread epistasis on the overdispersion of the molecular clock and the efficacy of common tests for positive selection.     

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.     

Epistatic Adaptive Evolution of Human Color Vision

Establishing genotype-phenotype relationship is the key to understand the molecular mechanism of phenotypic adaptation. This initial step may be untangled by analyzing appropriate ancestral molecules, but it is a daunting task to recapitulate the evolution of non-additive (epistatic) interactions of amino acids and function of a protein separately. To adapt to the ultraviolet (UV)-free retinal environment, the short wavelength-sensitive (SWS1) visual pigment in human (human S1) switched from detecting UV to absorbing blue light during the last 90 million years. Mutagenesis experiments of the UV-sensitive pigment in the Boreoeutherian ancestor show that the blue-sensitivity was achieved by seven mutations. The experimental and quantum chemical analyses show that 4,008 of all 5,040 possible evolutionary trajectories are terminated prematurely by containing a dehydrated nonfunctional pigment. Phylogenetic analysis further suggests that human ancestors achieved the blue-sensitivity gradually and almost exclusively by epistasis. When the final stage of spectral tuning of human S1 was underway 45–30 million years ago, the middle and long wavelength-sensitive (MWS/LWS) pigments appeared and so-called trichromatic color vision was established by interprotein epistasis. The adaptive evolution of human S1 differs dramatically from orthologous pigments with a major mutational effect used in achieving blue-sensitivity in a fish and several mammalian species and in regaining UV vision in birds. These observations imply that the mechanisms of epistatic interactions must be understood by studying various orthologues in different species that have adapted to various ecological and physiological environments.     

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.     

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.     

How enzymes adapt: lessons from directed evolution

Enzymes that are adapted to widely different temperature niches are being used to investigate the molecular basis of protein stability and enzyme function. However, natural evolution is complex: random noise, historical accidents and ignorance of the selection pressures at work during adaptation all cloud comparative studies. Here, we review how adaptation in the laboratory by directed evolution can complement studies of natural enzymes in the effort to understand stability and function. Laboratory evolution experiments can attempt to mimic natural evolution and identify different adaptive mechanisms. However, laboratory evolution might make its biggest contribution in explorations of nonnatural functions, by allowing us to distinguish the properties nutured by evolution from those dictated by the laws of physical chemistry.     

Defective discoidin domain structure,subunit assembly,and endoplasmic reticulum processing of retinoschisin are primary mechanisms responsible for X-linked retinoschisis

Retinoschisin is a 24-kDa discoidin domain-containing protein that is secreted from photoreceptor and bipolar cells as a large disulfide-linked multisubunit complex. It functions as a cell adhesion protein to maintain the cellular organization and synaptic structure of the retina. Over 125 different mutations in the RS1 gene are associated with X-linked juvenile retinoschisis, the most common form of early onset macular degeneration in males. To identify molecular determinants important for retinoschisin structure and function and elucidate molecular and cellular mechanisms responsible for X-linked juvenile retinoschisis, we have analyzed the expression, protein folding, disulfide-linked subunit assembly, intracellular localization, and secretion of wild-type retinoschisin, 15 Cys-to-Ser variants and 12 disease-linked mutants. Our studies, together with molecular modeling of the discoidin domain, identify Cys residues involved in intramolecular and intermolecular disulfide bonds essential for protein folding and subunit assembly. We show that misfolding of the discoidin domain, defective disulfide-linked subunit assembly, and inability of retinoschisin to insert into the endoplasmic reticulum membrane as part of the protein secretion process are three primary mechanisms responsible for the loss in the function of retinoschisin as a cell adhesion protein and the pathogenesis of X-linked juvenile retinoschisis.     

Social interaction in synthetic and natural microbial communities

Social interaction among cells is essential for multicellular complexity. But how do molecular networks within individual cells confer the ability to interact? And how do those same networks evolve from the evolutionary conflict between individual‐ and population‐level interests? Recent studies have dissected social interaction at the molecular level by analyzing both synthetic and natural microbial populations. These studies shed new light on the role of population structure for the evolution of cooperative interactions and revealed novel molecular mechanisms that stabilize cooperation among cells. New understanding of populations is changing our view of microbial processes, such as pathogenesis and antibiotic resistance, and suggests new ways to fight infection by exploiting social interaction. The study of social interaction is also challenging established paradigms in cancer evolution and immune system dynamics. Finding similar patterns in such diverse systems suggests that the same ‘social interaction motifs’ may be general to many cell populations.     

Epistasis from functional dependence of fitness on underlying traits

Epistasis between mutations in two genes is thought to reflect an interdependence of their functions. While sometimes epistasis is predictable using mechanistic models, its roots seem, in general, hidden in the complex architecture of biological networks. Here, we ask how epistasis can be quantified based on the mathematical dependence of a system-level trait (e.g. fitness) on lower-level traits (e.g. molecular or cellular properties). We first focus on a model in which fitness is the difference between a benefit and a cost trait, both pleiotropically affected by mutations. We show that despite its simplicity, this model can be used to analytically predict certain properties of the ensuing distribution of epistasis, such as a global negative bias, resulting in antagonism between beneficial mutations, and synergism between deleterious ones. We next extend these ideas to derive a general expression for epistasis given an arbitrary functional dependence of fitness on other traits. This expression demonstrates how epistasis relative to fitness can emerge despite the absence of epistasis relative to lower level traits, leading to a formalization of the concept of independence between biological processes. Our results suggest that epistasis may be largely shaped by the pervasiveness of pleiotropic effects and modular organization in biological networks.     

Recreating ancestral proteins

Tracing the history of molecular changes using phylogenetic methods can provide powerful insights into how and why molecules work the way they do. It is now possible to recreate inferred ancestral proteins in the laboratory and study the function of these molecules. This provides a unique opportunity to study the paths and the mechanisms of functional change during molecular evolution. What insights have already emerged from such phylogenetic studies of protein evolution and function, what are the impediments to progress and what are the prospects for the future?     

The proteomics of feather development in pied flycatchers (Ficedula hypoleuca) with different plumage coloration

The genetic theory of morphological evolution postulates that form evolves largely by changing the expression proteins that are functionally conserved. It follows that understanding the function of proteins during different phases of development as well as the mechanisms by which the functions are modified is a prerequisite for understanding evolutionary change. Male pied flycatchers exhibit marked phenotypic variation in their breeding plumage. This variation has repeatedly been shown to have adaptive significance, but the molecular basis of this variation is not known. Here, we characterize the proteome of developing pied flycatcher feathers from differently pigmented males and also introduce a new method for examining the effect sizes of expression differences in protein interaction networks. Approximately 300 proteins were identified in the developing feathers of males. Gene products associated with cellular transport, cell metabolism and protein synthesis formed a large part of the developing feather proteome. Sixty‐five proteins associated with the development of the epidermis and/or pigmentation were detected in the data. The examination of expression level differences of protein–protein interaction networks revealed an immunological signalling–related network to exhibit significantly higher expression in black compared to brown males. Additionally, indications of differences in energy balance and oxidative stress related characteristics were detected. Together, these results provide new insight into the molecular mechanisms and evolutionary significance of plumage colour variation.