Epistasis and the mutation load: a measurement-theoretical approach

An approximate solution for the mean fitness in mutation-selection balance with arbitrary order of epistatic interaction is derived. The solution is based on the assumptions of coupling equilibrium and that the interaction effects are multilinear. We find that the effect of m-order epistatic interactions (i.e., interactions among groups of m loci) on the load is dependent on the total genomic mutation rate, U, to the mth power. Thus, higher-order gene interactions are potentially important if U is large and the interaction density among loci is not too low. The solution suggests that synergistic epistasis will decrease the mutation load and that variation in epistatic effects will elevate the load. Both of these results, however, are strictly true only if they refer to epistatic interaction strengths measured in the optimal genotype. If gene interactions are measured at mutation-selection equilibrium, only synergistic interactions among even numbers of genes will reduce the load. Odd-ordered synergistic interactions will then elevate the load. There is no systematic relationship between variation in epistasis and load at equilibrium. We argue that empirical estimates of gene interaction must pay attention to the genetic background in which the effects are measured and that it may be advantageous to refer to average interaction intensities as measured in mutation-selection equilibrium. We derive a simple criterion for the strength of epistasis that is necessary to overcome the twofold disadvantage of sex.     

Replication and explorations of high-order epistasis using a large advanced intercross line pedigree

Dissection of the genetic architecture of complex traits persists as a major challenge in biology; despite considerable efforts, much remains unclear including the role and importance of genetic interactions. This study provides empirical evidence for a strong and persistent contribution of both second- and third-order epistatic interactions to long-term selection response for body weight in two divergently selected chicken lines. We earlier reported a network of interacting loci with large effects on body weight in an F(2) intercross between these high- and low-body weight lines. Here, most pair-wise interactions in the network are replicated in an independent eight-generation advanced intercross line (AIL). The original report showed an important contribution of capacitating epistasis to growth, meaning that the genotype at a hub in the network releases the effects of one or several peripheral loci. After fine-mapping of the loci in the AIL, we show that these interactions were persistent over time. The replication of five of six originally reported epistatic loci, as well as the capacitating epistasis, provides strong empirical evidence that the originally observed epistasis is of biological importance and is a contributor in the genetic architecture of this population. The stability of genetic interaction mechanisms over time indicates a non-transient role of epistasis on phenotypic change. Third-order epistasis was for the first time examined in this study and was shown to make an important contribution to growth, which suggests that the genetic architecture of growth is more complex than can be explained by two-locus interactions only. Our results illustrate the importance of designing studies that facilitate exploration of epistasis in populations for obtaining a comprehensive understanding of the genetics underlying a complex trait.     

水稻株高上位性效应和ＱＥ互作效应的ＱＴＬ遗传研究

利用基因混合模型的ＱＴＬ定位方法研究了由籼稻品种ＩＲ６４和粳稻品种Ａzucena杂交衍生的ＤＨ群体在４个环境中的ＱＴＬ上位性效应和环境互作效应,结果表明,上位性是数量性状的重要遗传基础,并揭示了上位性的几个重要特点,所有的ＱＴＬ都参与了上位性效应的形成,６４％的ＱＴＬ还具有本身的加性效应,因此传统方法对ＱＴＬ加性效应的估算会由于上位性的影响而有偏,其他３６％的ＱＴＬ没有本身的加性效应,却参与了４８％的上位性互作用,这些位点可能通过诱发和修饰其他位点而起作用,上位性的特点还包括,经常发现了一个ＱＴＬ与多个ＱＴＬ发生互作;大效应的ＱＴＬ也参与上位性互作;上位性互作受环境影响,ＱＴＬ与环境的互效应比ＱＴＬ的主效应更多地被检测到,表明数量性状基因的表达易受环境影响。     

A General Model for Multilocus Epistatic Interactions in Case-Control Studies

Background

Epistasis, i.e., the interaction of alleles at different loci, is thought to play a central role in the formation and progression of complex diseases. The complexity of disease expression should arise from a complex network of epistatic interactions involving multiple genes.

Methodology

We develop a general model for testing high-order epistatic interactions for a complex disease in a case-control study. We incorporate the quantitative genetic theory of high-order epistasis into the setting of cases and controls sampled from a natural population. The new model allows the identification and testing of epistasis and its various genetic components.

Conclusions

Simulation studies were used to examine the power and false positive rates of the model under different sampling strategies. The model was used to detect epistasis in a case-control study of inflammatory bowel disease, in which five SNPs at a candidate gene were typed, leading to the identification of a significant three-locus epistasis.     

Heterosis, sex and maternal interactions in crosses among inbred lines of mice.

All possible crosses and reciprocals were made among four inbred lines (F = 92%) developed from 12 generations of full-sib mating. All lines originated from a common outbred base population of ICR-albino mice. Data were obtained from 356 litter containing 2,734 mice to evaluate heterosis, reciprocal effects, sex effects and their interactions as they affect body weight and weight gain. Heterosis was significant for most of the postweaning traits (42- and 56-day weight and gain from 21 to 42 days). Nonadditive gene action may have included overdominance and epistasis since both reciprocal linecrosses were generally heavier than those of the better inbred lines. Although significant differences in reciprocals and inbred lines were not frequent, there were sufficient differences to indicate that lines varied in the fixation of loci during inbreeding. Sex-heterosis interactions were significant for 12 of 30 possible cases. However, eight of the 12 significant interactions occurred in crosses involving only one of the lines. The interactions were of the divergent type and arose from males exhibiting more heterosis than females. Overdominance in genes on the sex chromosomes modified by other loci (epistasis) was proposed as a possible explanation for these results. Some sex-linkage affecting growth was evident from the interaction of sex with reciprocal effects.     

Mapping genetically compensatory pathways from synthetic lethal interactions in yeast

Background

Synthetic lethal genetic interaction analysis has been successfully applied to predicting the functions of genes and their pathway identities. In the context of synthetic lethal interaction data alone, the global similarity of synthetic lethal interaction patterns between two genes is used to predict gene function. With physical interaction data, such as protein-protein interactions, the enrichment of physical interactions within subsets of genes and the enrichment of synthetic lethal interactions between those subsets of genes are used as an indication of compensatory pathways.

Result

In this paper, we propose a method of mapping genetically compensatory pathways from synthetic lethal interactions. Our method is designed to discover pairs of gene-sets in which synthetic lethal interactions are depleted among the genes in an individual set and where such gene-set pairs are connected by many synthetic lethal interactions. By its nature, our method could select compensatory pathway pairs that buffer the deleterious effect of the failure of either one, without the need of physical interaction data. By focusing on compensatory pathway pairs where genes in each individual pathway have a highly homogenous cellular function, we show that many cellular functions have genetically compensatory properties.

Conclusion

We conclude that synthetic lethal interaction data are a powerful source to map genetically compensatory pathways, especially in systems lacking physical interaction information, and that the cellular function network contains abundant compensatory properties.     

水稻籽粒淀粉分支酶活性的遗传分析

用由247个株系组成的珍汕97B/密阳46重组自交系群体及其含207个分子标记的连锁图谱,在2002年和2003年分别测定亲本和重组自交系群体开花后10 d和20 d籽粒的淀粉分支酶的活性,检测到3个控制开花后10 d Q酶活性的主效应QTL(qnantitative trait loci),联合贡献率为10%,其中qQ10-6与环境发生显著的互作;分别检测到5对和2对染色体区间对开花后10 d、20 d Q酶活性的影响具有加性×加性上位性作用,其中开花后10 d的3对染色体区间具有显著的上位性×环境互作效应.由此可见,水稻籽粒Q酶活性相关基因的表达,受到环境因子的极大影响.     

Epistatic interaction is an important genetic basis of grain yield and its components in maize

A population of 294 recombinant inbred lines (RIL) derived from Yuyu22, an elite maize hybrid extending broadly in China, has been constructed to investigate the genetic basis of grain yield, and associated yield components in maize. The main-effect quantitative trait loci (QTL), digenic epistatic interactions, and their interactions with the environment for grain yield and its three components were identified by using the mixed linear model approach. Thirty-two main-effect QTL and forty-four pairs of digenic epistatic interactions were detected for the four measured traits in four environments. Our results suggest that both additive effects and epistasis (additive × additive) effects are important genetic bases of grain yield and its components in the RIL population. Only 30.4% of main-effect QTL for ear length were involved in epistatic interactions. This implies that many loci in epistatic interactions may not have significant effects for traits alone but may affect trait expression by epistatic interaction with the other loci.     

Exploiting the proteome to improve the genome-wide genetic analysis of epistasis in common human diseases

One of the central goals of human genetics is the identification of loci with alleles or genotypes that confer increased susceptibility. The availability of dense maps of single-nucleotide polymorphisms (SNPs) along with high-throughput genotyping technologies has set the stage for routine genome-wide association studies that are expected to significantly improve our ability to identify susceptibility loci. Before this promise can be realized, there are some significant challenges that need to be addressed. We address here the challenge of detecting epistasis or gene–gene interactions in genome-wide association studies. Discovering epistatic interactions in high dimensional datasets remains a challenge due to the computational complexity resulting from the analysis of all possible combinations of SNPs. One potential way to overcome the computational burden of a genome-wide epistasis analysis would be to devise a logical way to prioritize the many SNPs in a dataset so that the data may be analyzed more efficiently and yet still retain important biological information. One of the strongest demonstrations of the functional relationship between genes is protein-protein interaction. Thus, it is plausible that the expert knowledge extracted from protein interaction databases may allow for a more efficient analysis of genome-wide studies as well as facilitate the biological interpretation of the data. In this review we will discuss the challenges of detecting epistasis in genome-wide genetic studies and the means by which we propose to apply expert knowledge extracted from protein interaction databases to facilitate this process. We explore some of the fundamentals of protein interactions and the databases that are publicly available.     

On the evolution of epistasis II: a generalized Wright-Kimura framework

The evolution of fitness interactions between genes at two major loci is studied where the alleles at a third locus modify the epistatic interaction between the two major loci. The epistasis is defined by a parameter epsilon and a matrix structure that specifies the nature of the interactions. When epsilon=0 the two major loci have additive fitnesses, and when these are symmetric the interaction matrices studied here produce symmetric viabilities of the Wright [1952. The genetics of quantitative variability. In: Reeve, E.C.R., Waddington, C.H. (Eds.), Quantitative Inheritance. Her Majesty's Stationary Office, London]-Kimura [1956. A model of a genetic system which leads to closer linkage by natural selection. Evolution 10, 278-281] form. Two such interaction matrices are studied, for one of which epistasis as measured by |epsilon| always increases, and for the other it increases when the linkage between the major loci is tight enough and there is initial linkage disequilibrium. Increase of epistasis does not necessarily coincide with increase in equilibrium mean fitness.     

Rules of nonallelic noncomplementation at the synapse in Caenorhabditis elegans

Nonallelic noncomplementation occurs when recessive mutations in two different loci fail to complement one another, in other words, the double heterozygote exhibits a phenotype. We observed that mutations in the genes encoding the physically interacting synaptic proteins UNC-13 and syntaxin/UNC-64 failed to complement one another in the nematode Caenorhabditis elegans. Noncomplementation was not observed between null alleles of these genes and thus this genetic interaction does not occur with a simple decrease in dosage at the two loci. However, noncomplementation was observed if at least one gene encoded a partially functional gene product. Thus, this genetic interaction requires a poisonous gene product to sensitize the genetic background. Nonallelic noncomplementation was not limited to interacting proteins: Although the strongest effects were observed between loci encoding gene products that bind to one another, interactions were also observed between proteins that do not directly interact but are members of the same complex. We also observed noncomplementation between genes that function at distant points in the same pathway, implying that physical interactions are not required for nonallelic noncomplementation. Finally, we observed that mutations in genes that function in different processes such as neurotransmitter synthesis or synaptic development complement one another. Thus, this genetic interaction is specific for genes acting in the same pathway, that is, for genes acting in synaptic vesicle trafficking.     

Characterisation of genome-wide association epistasis signals for serum uric acid in human population isolates

Genome-wide association (GWA) studies have identified a number of loci underlying variation in human serum uric acid (SUA) levels with the SLC2A9 gene having the largest effect identified so far. Gene-gene interactions (epistasis) are largely unexplored in these GWA studies. We performed a full pair-wise genome scan in the Italian MICROS population (n = 1201) to characterise epistasis signals in SUA levels. In the resultant epistasis profile, no SNP pairs reached the Bonferroni adjusted threshold for the pair-wise genome-wide significance. However, SLC2A9 was found interacting with multiple loci across the genome, with NFIA - SLC2A9 and SLC2A9 - ESRRAP2 being significant based on a threshold derived for interactions between GWA significant SNPs and the genome and jointly explaining 8.0% of the phenotypic variance in SUA levels (3.4% by interaction components). Epistasis signal replication in a CROATIAN population (n = 1772) was limited at the SNP level but improved dramatically at the gene ontology level. In addition, gene ontology terms enriched by the epistasis signals in each population support links between SUA levels and neurological disorders. We conclude that GWA epistasis analysis is useful despite relatively low power in small isolated populations.     

Regulatory hierarchy of photomorphogenic loci: allele-specific and light-dependent interaction between the HY5 and COP1 loci.

Previous studies suggested that the CONSTITUTIVE PHOTOMORPHOGENIC 1 (COP1) gene product represses photomorphogenic development in darkness and that light signals reverse this action. In this report, we used genetic analysis to investigate the regulatory hierarchical relationship of COP1 and the loci encoding the photoreceptors and other signaling components. Our results showed that cop1 mutations are epistatic to the long hypocotyl mutations hy1, hy2, hy3, and hy4, suggesting that COP1 acts downstream of the phytochromes and a blue light receptor. Although epistasis of a putative null cop1-5 mutation over a hy5 mutation implied that COP1 acts downstream of HY5, the same hy5 mutation can suppress the dark photomorphogenic phenotypes (including hypocotyl elongation and cotyledon cellular differentiation) of the weak cop1-6 mutation. This, and other allele-specific interactions between COP1 and HY5, may suggest direct physical contact of their gene products. In addition, the synthetic lethality of the weak deetiolated1 (det1) and cop1 mutations and the fact that the cop1-6 mutation is epistatic to the det1-1 mutation with respect to light control of seed germination and dark-adaptative gene expression suggested that DET1 and COP1 may act in the same pathway, with COP1 being downstream. These results, together with previous epistasis studies, support models in which light signals, once perceived by different photoreceptors, converge downstream and act through a common cascade(s) of regulatory steps, as defined by DET1, HY5, COP1, and likely others, to derepress photomorphogenic development.     

Detection of epistatic interactions between exotic alleles introgressed from wild barley (H. vulgare ssp. spontaneum)

The expression of a quantitative phenotype can be controlled through genotype, environment and genotype by environment interaction effects. Further, genotype effects can be attributed to major genes, quantitative trait loci (QTL) and gene by gene interactions, which are also termed epistatic interactions. The present study demonstrates that two-way epistatic interactions can play an important role for the expression of domestication-related traits like heading date, plant height and yield. In the BC₂DH population S42, carrying wild barley introgressions in the genetic background of the spring barley cultivar Scarlett, 13, 8 and 12 marker by marker interaction effects could be detected for the traits heading date, plant height and yield, respectively. Significant allelic combinations at interacting loci coincided for heading date, plant height and yield suggesting the presence of pleiotropic effects rather than several linked QTL. The mode of epistasis observed was primarily characterised by either (1) compensatory effects, where allelic combinations from the same genotype buffered the phenotype, or (2) augmented effects, where only the combination of the exotic allele at both interacting loci caused an altered phenotype. The present study shows that estimates of main effects of QTL can be confounded by interactions with background loci, suggesting that the identification of epistatic effects is important for gene cloning and marker-assisted selection. Furthermore, interaction effects between loci and putative candidate genes detected in the present study reveal potential functional relationships, which can be used to further elucidate gene networks in barley.     

Association Mapping for Epistasis and Environmental Interaction of Yield Traits in 323 Cotton Cultivars under 9 Different Environments

Improving yield is a major objective for cotton breeding schemes, and lint yield and its three component traits (boll number, boll weight and lint percentage) are complex traits controlled by multiple genes and various environments. Association mapping was performed to detect markers associated with these four traits using 651 simple sequence repeats (SSRs). A mixed linear model including epistasis and environmental interaction was used to screen the loci associated with these four yield traits by 323 accessions of Gossypium hirsutum L. evaluated in nine different environments. 251 significant loci were detected to be associated with lint yield and its three components, including 69 loci with individual effects and all involved in epistasis interactions. These significant loci explain ∼ 62.05% of the phenotypic variance (ranging from 49.06% ∼ 72.29% for these four traits). It was indicated by high contribution of environmental interaction to the phenotypic variance for lint yield and boll numbers, that genetic effects of SSR loci were susceptible to environment factors. Shared loci were also observed among these four traits, which may be used for simultaneous improvement in cotton breeding for yield traits. Furthermore, consistent and elite loci were screened with −Log₁₀ (P-value) >8.0 based on predicted effects of loci detected in different environments. There was one locus and 6 pairs of epistasis for lint yield, 4 loci and 10 epistasis for boll number, 15 loci and 2 epistasis for boll weight, and 2 loci and 5 epistasis for lint percentage, respectively. These results provided insights into the genetic basis of lint yield and its components and may be useful for marker-assisted breeding to improve cotton production.     

Synthetic Cytotoxicity: Digenic Interactions with TEL1/ATM Mutations Reveal Sensitivity to Low Doses of Camptothecin

Many tumors contain mutations that confer defects in the DNA-damage response and genome stability. DNA-damaging agents are powerful therapeutic tools that can differentially kill cells with an impaired DNA-damage response. The response to DNA damage is complex and composed of a network of coordinated pathways, often with a degree of redundancy. Tumor-specific somatic mutations in DNA-damage response genes could be exploited by inhibiting the function of a second gene product to increase the sensitivity of tumor cells to a sublethal concentration of a DNA-damaging therapeutic agent, resulting in a class of conditional synthetic lethality we call synthetic cytotoxicity. We used the Saccharomyces cerevisiae nonessential gene-deletion collection to screen for synthetic cytotoxic interactions with camptothecin, a topoisomerase I inhibitor, and a null mutation in TEL1, the S. cerevisiae ortholog of the mammalian tumor-suppressor gene, ATM. We found and validated 14 synthetic cytotoxic interactions that define at least five epistasis groups. One class of synthetic cytotoxic interaction was due to telomere defects. We also found that at least one synthetic cytotoxic interaction was conserved in Caenorhabditis elegans. We have demonstrated that synthetic cytotoxicity could be a useful strategy for expanding the sensitivity of certain tumors to DNA-damaging therapeutics.     

dSLAM analysis of genome-wide genetic interactions in Saccharomyces cerevisiae

Analysis of genetic interactions has been extensively exploited to study gene functions and to dissect pathway structures. One such genetic interaction is synthetic lethality, in which the combination of two non-lethal mutations leads to loss of organism viability. We have developed a dSLAM (heterozygote diploid-based synthetic lethality analysis with microarrays) technology that effectively studies synthetic lethality interactions on a genome-wide scale in the budding yeast Saccharomyces cerevisiae. Typically, a query mutation is introduced en masse into a population of approximately 6000 haploid-convertible heterozygote diploid Yeast Knockout (YKO) mutants via integrative transformation. Haploid pools of single and double mutants are freshly generated from the resultant heterozygote diploid double mutant pool after meiosis and haploid selection and studied for potential growth defects of each double mutant combination by microarray analysis of the "molecular barcodes" representing each YKO. This technology has been effectively adapted to study other types of genome-wide genetic interactions including gene-compound synthetic lethality, secondary mutation suppression, dosage-dependent synthetic lethality and suppression.     

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.     

基于单片段代换系的水稻穗长QTL加性及其上位性效应

穗长是影响水稻(Oryza sativa)产量的重要因子之一, 研究水稻穗长QTL间的上位性效应对于发掘水稻产量潜力具有重要意义。该研究以16个单片段代换系(single segment substitution lines, SSSLs)和15个双片段聚合系(double segment pyramiding lines, DSPLs)为材料研究了水稻穗长QTL的加性及上位性效应。以P<0.01为阈值, 共检测到6个穗长QTL和9对基因互作座位。其中2个(Pl-2和Pl-10)是尚未报道的穗长QTL。穗长QTL互作后, 一些互作对的上位性效应与单个QTL的作用方式及效应大小各不相同, 预示着基因聚合后会产生不同的互作效应。该研究结果对于通过分子聚合育种手段改良穗长具有重要意义。