Epistasis is a phenomenon in genetics in which the effect of a gene mutation is dependent on the presence or absence of mutations in one or more other genes, respectively termed modifier genes. In other words, the effect of the mutation is dependent on the genetic background in which it appears.[1] Epistatic mutations therefore have different effects on their own than when they occur together. Originally, the term epistasis specifically meant that the effect of a gene variant is masked by that of different gene.
The concept of epistasis originated in genetics in 1907[2] but is now used in biochemistry, computational biology and evolutionary biology. The phenomenon arises due to interactions, either between genes (such as mutations also being needed in regulators of gene expression) or within them (multiple mutations being needed before the gene loses function), leading to non-linear effects. Epistasis has a great influence on the shape of evolutionary landscapes, which leads to profound consequences for evolution and for the evolvability of phenotypic traits.
Understanding of epistasis has changed considerably through the history of genetics and so too has the use of the term. The term was first used by William Bateson and his collaborators Florence Durham and Muriel Wheldale Onslow. In early models of natural selection devised in the early 20th century, each gene was considered to make its own characteristic contribution to fitness, against an average background of other genes. Some introductory courses still teach population genetics this way. Because of the way that the science of population genetics was developed, evolutionary geneticists have tended to think of epistasis as the exception. However, in general, the expression of any one allele depends in a complicated way on many other alleles.
In classical genetics, if genes A and B are mutated, and each mutation by itself produces a unique phenotype but the two mutations together show the same phenotype as the gene A mutation, then gene A is epistatic and gene B is hypostatic. For example, the gene for total baldness is epistatic to the gene for brown hair. In this sense, epistasis can be contrasted with genetic dominance, which is an interaction between alleles at the same gene locus. As the study of genetics developed, and with the advent of molecular biology, epistasis started to be studied in relation to quantitative trait loci (QTL) and polygenic inheritance.
The effects of genes are now commonly quantifiable by assaying the magnitude of a phenotype (e.g. height, pigmentation or growth rate) or by biochemically assaying protein activity (e.g. binding or catalysis). Increasingly sophisticated computational and evolutionary biology models aim to describe the effects of epistasis on a genome-wide scale and the consequences of this for evolution.[3] [4] [5] Since identification of epistatic pairs is challenging both computationally and statistically, some studies try to prioritize epistatic pairs.[6] [7]
Terminology about epistasis can vary between scientific fields. Geneticists often refer to wild type and mutant alleles where the mutation is implicitly deleterious and may talk in terms of genetic enhancement, synthetic lethality and genetic suppressors. Conversely, a biochemist may more frequently focus on beneficial mutations and so explicitly state the effect of a mutation and use terms such as reciprocal sign epistasis and compensatory mutation.[8] Additionally, there are differences when looking at epistasis within a single gene (biochemistry) and epistasis within a haploid or diploid genome (genetics). In general, epistasis is used to denote the departure from 'independence' of the effects of different genetic loci. Confusion often arises due to the varied interpretation of 'independence' among different branches of biology.[9] The classifications below attempt to cover the various terms and how they relate to one another.
Two mutations are considered to be purely additive if the effect of the double mutation is the sum of the effects of the single mutations. This occurs when genes do not interact with each other, for example by acting through different metabolic pathways. Simply, additive traits were studied early on in the history of genetics, however they are relatively rare, with most genes exhibiting at least some level of epistatic interaction.[10] [11]
When the double mutation has a fitter phenotype than expected from the effects of the two single mutations, it is referred to as positive epistasis. Positive epistasis between beneficial mutations generates greater improvements in function than expected. Positive epistasis between deleterious mutations protects against the negative effects to cause a less severe fitness drop.
Conversely, when two mutations together lead to a less fit phenotype than expected from their effects when alone, it is called negative epistasis.[12] [13] Negative epistasis between beneficial mutations causes smaller than expected fitness improvements, whereas negative epistasis between deleterious mutations causes greater-than-additive fitness drops.
Independently, when the effect on fitness of two mutations is more radical than expected from their effects when alone, it is referred to as synergistic epistasis. The opposite situation, when the fitness difference of the double mutant from the wild type is smaller than expected from the effects of the two single mutations, it is called antagonistic epistasis. Therefore, for deleterious mutations, negative epistasis is also synergistic, while positive epistasis is antagonistic; conversely, for advantageous mutations, positive epistasis is synergistic, while negative epistasis is antagonistic.
The term genetic enhancement is sometimes used when a double (deleterious) mutant has a more severe phenotype than the additive effects of the single mutants. Strong positive epistasis is sometimes referred to by creationists as irreducible complexity (although most examples are misidentified).
Sign epistasis[14] occurs when one mutation has the opposite effect when in the presence of another mutation. This occurs when a mutation that is deleterious on its own can enhance the effect of a particular beneficial mutation.[9] For example, a large and complex brain is a waste of energy without a range of sense organs, but sense organs are made more useful by a large and complex brain that can better process the information. If a fitness landscape has no sign epistasis then it is called smooth.
At its most extreme, reciprocal sign epistasis[15] occurs when two deleterious genes are beneficial when together. For example, producing a toxin alone can kill a bacterium, and producing a toxin exporter alone can waste energy, but producing both can improve fitness by killing competing organisms. If a fitness landscape has sign epistasis but no reciprocal sign epistasis then it is called semismooth.[16]
Reciprocal sign epistasis also leads to genetic suppression whereby two deleterious mutations are less harmful together than either one on its own, i.e. one compensates for the other. A clear example of genetic suppression was the demonstration that in the assembly of bacteriophage T4 two deleterious mutations, each causing a deficiency in the level of a different morphogenetic protein, could interact positively.[17] If a mutation causes a reduction in a particular structural component, this can bring about an imbalance in morphogenesis and loss of viable virus progeny, but production of viable progeny can be restored by a second (suppressor) mutation in another morphogenetic component that restores the balance of protein components.
The term genetic suppression can also apply to sign epistasis where the double mutant has a phenotype intermediate between those of the single mutants, in which case the more severe single mutant phenotype is suppressed by the other mutation or genetic condition. For example, in a diploid organism, a hypomorphic (or partial loss-of-function) mutant phenotype can be suppressed by knocking out one copy of a gene that acts oppositely in the same pathway. In this case, the second gene is described as a "dominant suppressor" of the hypomorphic mutant; "dominant" because the effect is seen when one wild-type copy of the suppressor gene is present (i.e. even in a heterozygote). For most genes, the phenotype of the heterozygous suppressor mutation by itself would be wild type (because most genes are not haplo-insufficient), so that the double mutant (suppressed) phenotype is intermediate between those of the single mutants.
In non reciprocal sign epistasis, fitness of the mutant lies in the middle of that of the extreme effects seen in reciprocal sign epistasis.
When two mutations are viable alone but lethal in combination, it is called Synthetic lethality or unlinked non-complementation.[18]
In a haploid organism with genotypes (at two loci) ab, Ab, aB or AB, we can think of different forms of epistasis as affecting the magnitude of a phenotype upon mutation individually (Ab and aB) or in combination (AB).
| Interaction type | ab | Ab | aB | AB | ||
| No epistasis (additive) | 0 | 1 | 1 | 2 | AB = Ab + aB + ab | |
| Positive (synergistic) epistasis | 0 | 1 | 1 | 3 | AB > Ab + aB + ab | |
| Negative (antagonistic) epistasis | 0 | 1 | 1 | 1 | AB < Ab + aB + ab | |
| Sign epistasis | 0 | 1 | -1 | 2 | AB has opposite sign to Ab or aB | |
| Reciprocal sign epistasis | 0 | -1 | -1 | 2 | AB has opposite sign to Ab and aB |
Epistasis in diploid organisms is further complicated by the presence of two copies of each gene. Epistasis can occur between loci, but additionally, interactions can occur between the two copies of each locus in heterozygotes. For a two locus, two allele system, there are eight independent types of gene interaction.[19]
| Additive A locus | Additive B locus | Dominance A locus | Dominance B locus | ||||||||||||||||
| aa | aA | AA | aa | aA | AA | aa | aA | AA | aa | aA | AA | ||||||||
| bb | 1 | 0 | –1 | bb | 1 | 1 | 1 | bb | –1 | 1 | –1 | bb | –1 | –1 | –1 | ||||
| bB | 1 | 0 | –1 | bB | 0 | 0 | 0 | bB | –1 | 1 | –1 | bB | 1 | 1 | 1 | ||||
| BB | 1 | 0 | –1 | BB | –1 | –1 | –1 | BB | –1 | 1 | –1 | BB | –1 | –1 | –1 | ||||
| Additive by Additive Epistasis | Additive by Dominance Epistasis | Dominance by Additive Epistasis | Dominance by Dominance Epistasis | ||||||||||||||||
| aa | aA | AA | aa | aA | AA | aa | aA | AA | aa | aA | AA | ||||||||
| bb | 1 | 0 | –1 | bb | 1 | 0 | –1 | bb | 1 | –1 | 1 | bb | –1 | 1 | –1 | ||||
| bB | 0 | 0 | 0 | bB | –1 | 0 | 1 | bB | 0 | 0 | 0 | bB | 1 | –1 | 1 | ||||
| BB | –1 | 0 | 1 | BB | 1 | 0 | –1 | BB | –1 | 1 | –1 | BB | –1 | 1 | –1 | ||||
This can be the case when multiple genes act in parallel to achieve the same effect. For example, when an organism is in need of phosphorus, multiple enzymes that break down different phosphorylated components from the environment may act additively to increase the amount of phosphorus available to the organism. However, there inevitably comes a point where phosphorus is no longer the limiting factor for growth and reproduction and so further improvements in phosphorus metabolism have smaller or no effect (negative epistasis). Some sets of mutations within genes have also been specifically found to be additive.[20] It is now considered that strict additivity is the exception, rather than the rule, since most genes interact with hundreds or thousands of other genes.[10] [11]
Epistasis within the genomes of organisms occurs due to interactions between the genes within the genome. This interaction may be direct if the genes encode proteins that, for example, are separate components of a multi-component protein (such as the ribosome), inhibit each other's activity, or if the protein encoded by one gene modifies the other (such as by phosphorylation). Alternatively the interaction may be indirect, where the genes encode components of a metabolic pathway or network, developmental pathway, signalling pathway or transcription factor network. For example, the gene encoding the enzyme that synthesizes penicillin is of no use to a fungus without the enzymes that synthesize the necessary precursors in the metabolic pathway.
Just as mutations in two separate genes can be non-additive if those genes interact, mutations in two codons within a gene can be non-additive. In genetics this is sometimes called intragenic suppression when one deleterious mutation can be compensated for by a second mutation within that gene. Analysis of bacteriophage T4 mutants that were altered in the rIIB cistron (gene) revealed that certain pairwise combinations of mutations could mutually suppress each other; that is the double mutants had a more nearly wild-type phenotype than either mutant alone.[21] The linear map order of the mutants was established using genetic recombination data, From these sources of information, the triplet nature of the genetic code was logically deduced for the first time in 1961, and other key features of the code were also inferred.[21]
Also intragenic suppression can occur when the amino acids within a protein interact. Due to the complexity of protein folding and activity, additive mutations are rare.
Proteins are held in their tertiary structure by a distributed, internal network of cooperative interactions (hydrophobic, polar and covalent).[22] Epistatic interactions occur whenever one mutation alters the local environment of another residue (either by directly contacting it, or by inducing changes in the protein structure).[23] For example, in a disulphide bridge, a single cysteine has no effect on protein stability until a second is present at the correct location at which point the two cysteines form a chemical bond which enhances the stability of the protein.[24] This would be observed as positive epistasis where the double-cysteine variant had a much higher stability than either of the single-cysteine variants. Conversely, when deleterious mutations are introduced, proteins often exhibit mutational robustness whereby as stabilising interactions are destroyed the protein still functions until it reaches some stability threshold at which point further destabilising mutations have large, detrimental effects as the protein can no longer fold. This leads to negative epistasis whereby mutations that have little effect alone have a large, deleterious effect together.[25]
In enzymes, the protein structure orients a few, key amino acids into precise geometries to form an active site to perform chemistry.[26] Since these active site networks frequently require the cooperation of multiple components, mutating any one of these components massively compromises activity, and so mutating a second component has a relatively minor effect on the already inactivated enzyme. For example, removing any member of the catalytic triad of many enzymes will reduce activity to levels low enough that the organism is no longer viable.[27] [28] [29]
Diploid organisms contain two copies of each gene. If these are different (heterozygous / heteroallelic), the two different copies of the allele may interact with each other to cause epistasis. This is sometimes called allelic complementation, or interallelic complementation. It may be caused by several mechanisms, for example transvection, where an enhancer from one allele acts in trans to activate transcription from the promoter of the second allele. Alternately, trans-splicing of two non-functional RNA molecules may produce a single, functional RNA.
Similarly, at the protein level, proteins that function as dimers may form a heterodimer composed of one protein from each alternate gene and may display different properties to the homodimer of one or both variants. Two bacteriophage T4 mutants defective at different locations in the same gene can undergo allelic complementation during a mixed infection. That is, each mutant alone upon infection cannot produce viable progeny, but upon mixed infection with two complementing mutants, viable phage are formed. Intragenic complementation was demonstrated for several genes that encode structural proteins of the bacteriophage[30] indicating that such proteins function as dimers or even higher order multimers.[31]
See also: fitness landscape and evolvability. In evolutionary genetics, the sign of epistasis is usually more significant than the magnitude of epistasis. This is because magnitude epistasis (positive and negative) simply affects how beneficial mutations are together, however sign epistasis affects whether mutation combinations are beneficial or deleterious.
A fitness landscape is a representation of the fitness where all genotypes are arranged in 2D space and the fitness of each genotype is represented by height on a surface. It is frequently used as a visual metaphor for understanding evolution as the process of moving uphill from one genotype to the next, nearby, fitter genotype.[10]
If all mutations are additive, they can be acquired in any order and still give a continuous uphill trajectory. The landscape is perfectly smooth, with only one peak (global maximum) and all sequences can evolve uphill to it by the accumulation of beneficial mutations in any order. Conversely, if mutations interact with one another by epistasis, the fitness landscape becomes rugged as the effect of a mutation depends on the genetic background of other mutations.[32] At its most extreme, interactions are so complex that the fitness is 'uncorrelated' with gene sequence and the topology of the landscape is random. This is referred to as a rugged fitness landscape and has profound implications for the evolutionary optimisation of organisms. If mutations are deleterious in one combination but beneficial in another, the fittest genotypes can only be accessed by accumulating mutations in one specific order. This makes it more likely that organisms will get stuck at local maxima in the fitness landscape having acquired mutations in the 'wrong' order.[33] [34] For example, a variant of TEM1 β-lactamase with 5 mutations is able to cleave cefotaxime (a third generation antibiotic).[35] However, of the 120 possible pathways to this 5-mutant variant, only 7% are accessible to evolution as the remainder passed through fitness valleys where the combination of mutations reduces activity. In contrast, changes in environment (and therefore the shape of the fitness landscape) have been shown to provide escape from local maxima. In this example, selection in changing antibiotic environments resulted in a "gateway mutation" which epistatically interacted in a positive manner with other mutations along an evolutionary pathway, effectively crossing a fitness valley. This gateway mutation alleviated the negative epistatic interactions of other individually beneficial mutations, allowing them to better function in concert. Complex environments or selections may therefore bypass local maxima found in models assuming simple positive selection.
High epistasis is usually considered a constraining factor on evolution, and improvements in a highly epistatic trait are considered to have lower evolvability. This is because, in any given genetic background, very few mutations will be beneficial, even though many mutations may need to occur to eventually improve the trait. The lack of a smooth landscape makes it harder for evolution to access fitness peaks. In highly rugged landscapes, fitness valleys block access to some genes, and even if ridges exist that allow access, these may be rare or prohibitively long.[36] Moreover, adaptation can move proteins into more precarious or rugged regions of the fitness landscape.[37] These shifting "fitness territories" may act to decelerate evolution and could represent tradeoffs for adaptive traits.
The frustration of adaptive evolution by rugged fitness landscapes was recognized as a potential force for the evolution of evolvability. Michael Conrad in 1972 was the first to propose a mechanism for the evolution of evolvability by noting that a mutation which smoothed the fitness landscape at other loci could facilitate the production of advantageous mutations and hitchhike along with them.[38] [39] Rupert Riedl in 1975 proposed that new genes which produced the same phenotypic effects with a single mutation as other loci with reciprocal sign epistasis would be a new means to attain a phenotype otherwise too unlikely to occur by mutation.[40] [41]
Rugged, epistatic fitness landscapes also affect the trajectories of evolution. When a mutation has a large number of epistatic effects, each accumulated mutation drastically changes the set of available beneficial mutations. Therefore, the evolutionary trajectory followed depends highly on which early mutations were accepted. Thus, repeats of evolution from the same starting point tend to diverge to different local maxima rather than converge on a single global maximum as they would in a smooth, additive landscape.[42] [43]
Negative epistasis and sex are thought to be intimately correlated. Experimentally, this idea has been tested in using digital simulations of asexual and sexual populations. Over time, sexual populations move towards more negative epistasis, or the lowering of fitness by two interacting alleles. It is thought that negative epistasis allows individuals carrying the interacting deleterious mutations to be removed from the populations efficiently. This removes those alleles from the population, resulting in an overall more fit population. This hypothesis was proposed by Alexey Kondrashov, and is sometimes known as the deterministic mutation hypothesis[44] and has also been tested using artificial gene networks.[12]
However, the evidence for this hypothesis has not always been straightforward and the model proposed by Kondrashov has been criticized for assuming mutation parameters far from real world observations.[45] In addition, in those tests which used artificial gene networks, negative epistasis is only found in more densely connected networks,[12] whereas empirical evidence indicates that natural gene networks are sparsely connected,[46] and theory shows that selection for robustness will favor more sparsely connected and minimally complex networks.[46]
Quantitative genetics focuses on genetic variance due to genetic interactions. Any two locus interactions at a particular gene frequency can be decomposed into eight independent genetic effects using a weighted regression. In this regression, the observed two locus genetic effects are treated as dependent variables and the "pure" genetic effects are used as the independent variables. Because the regression is weighted, the partitioning among the variance components will change as a function of gene frequency. By analogy it is possible to expand this system to three or more loci, or to cytonuclear interactions[47]
When assaying epistasis within a gene, site-directed mutagenesis can be used to generate the different genes, and their protein products can be assayed (e.g. for stability or catalytic activity). This is sometimes called a double mutant cycle and involves producing and assaying the wild type protein, the two single mutants and the double mutant. Epistasis is measured as the difference between the effects of the mutations together versus the sum of their individual effects.[48] This can be expressed as a free energy of interaction.The same methodology can be used to investigate the interactions between larger sets of mutations but all combinations have to be produced and assayed. For example, there are 120 different combinations of 5 mutations, some or all of which may show epistasis...
Numerous computational methods have been developed for the detection and characterization of epistasis. Many of these rely on machine learning to detect non-additive effects that might be missed by statistical approaches such as linear regression.[49] For example, multifactor dimensionality reduction (MDR) was designed specifically for nonparametric and model-free detection of combinations of genetic variants that are predictive of a phenotype such as disease status in human populations.[50] [51] Several of these approaches have been broadly reviewed in the literature.[52] Even more recently, methods that utilize insights from theoretical computer science (the Hadamard transform[53] and compressed sensing[54] [55]) or maximum-likelihood inference[56] were shown to distinguish epistatic effects from overall non-linearity in genotype–phenotype map structure,[57] while others used patient survival analysis to identify non-linearity.[58]