Higher-order genetic interactions and their contribution to complex traits

Assessing cell-specific effects of genetic variations using tRNA microarrays

Background

By definition, effect of synonymous single-nucleotide variants (SNVs) on protein folding and function are neutral, as they alter the codon and not the encoded amino acid. Recent examples indicate tissue-specific and transfer RNA (tRNA)-dependent effects of some genetic variations arguing against neutrality of synonymous SNVs for protein biogenesis.

Results

We performed systematic analysis of tRNA abunandance across in various models used in cystic fibrosis (CF) research and drug development, including Fischer rat thyroid (FRT) cells, patient-derived primary human bronchial epithelia (HBE) from lung biopsies, primary human nasal epithelia (HNE) from nasal curettage, intestinal organoids, and airway progenitor-directed differentiation of human induced pluripotent stem cells (iPSCs). These were compared to an immortalized CF bronchial cell model (CFBE41o⁻) and two widely used laboratory cell lines, HeLa and HEK293. We discovered that specific synonymous SNVs exhibited differential effects which correlated with variable concentrations of cognate tRNAs.

Conclusions

Our results highlight ways in which the presence of synonymous SNVs may alter local kinetics of mRNA translation; and thus, impact protein biogenesis and function. This effect is likely to influence results from mechansistic analysis and/or drug screeining efforts, and establishes importance of cereful model system selection based on genetic variation profile.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5864-1) contains supplementary material, which is available to authorized users.

Oligogenic inheritance of a human heart disease involving a genetic modifier

Complex genetic mechanisms are thought to underlie many human diseases, yet experimental proof of this model has been elusive. Here, we show that a human cardiac anomaly can be caused by a combination of rare, inherited heterozygous mutations. Whole-exome sequencing of a nuclear family revealed that three offspring with childhood-onset cardiomyopathy had inherited three missense single-nucleotide variants in the MKL2, MYH7, and NKX2-5 genes. The MYH7 and MKL2 variants were inherited from the affected, asymptomatic father and the rare NKX2-5 variant (minor allele frequency, 0.0012) from the unaffected mother. We used CRISPR-Cas9 to generate mice encoding the orthologous variants and found that compound heterozygosity for all three variants recapitulated the human disease phenotype. Analysis of murine hearts and human induced pluripotent stem cell-derived cardiomyocytes provided histologic and molecular evidence for the NKX2-5 variant's contribution as a genetic modifier.

Layers of Cryptic Genetic Variation Underlie a Yeast Complex Trait

To better understand cryptic genetic variation, Lee et al. comprehensively map the genetic basis of a trait that is typically suppressed in a yeast cross. By determining how three different genetic perturbations give rise...

Cryptic genetic variation may be an important contributor to heritable traits, but its extent and regulation are not fully understood. Here, we investigate the cryptic genetic variation underlying a Saccharomyces cerevisiae colony phenotype that is typically suppressed in a cross of the laboratory strain BY4716 (BY) and a derivative of the clinical isolate 322134S (3S). To do this, we comprehensively dissect the trait’s genetic basis in the BYx3S cross in the presence of three different genetic perturbations that enable its expression. This allows us to detect and compare the specific loci that interact with each perturbation to produce the trait. In total, we identify 21 loci, all but one of which interact with just a subset of the perturbations. Beyond impacting which loci contribute to the trait, the genetic perturbations also alter the extent of additivity, epistasis, and genotype–environment interaction among the detected loci. Additionally, we show that the single locus interacting with all three perturbations corresponds to the coding region of the cell surface gene FLO11. While nearly all of the other remaining loci influence FLO11 transcription in cis or trans, the perturbations tend to interact with loci in different pathways and subpathways. Our work shows how layers of cryptic genetic variation can influence complex traits. Here, these layers mainly represent different regulatory inputs into the transcription of a single key gene.

Time-resolved mapping of genetic interactions to model rewiring of signaling pathways

Context-dependent changes in genetic interactions are an important feature of cellular pathways and their varying responses under different environmental conditions. However, methodological frameworks to investigate the plasticity of genetic interaction networks over time or in response to external stresses are largely lacking. To analyze the plasticity of genetic interactions, we performed a combinatorial RNAi screen in Drosophila cells at multiple time points and after pharmacological inhibition of Ras signaling activity. Using an image-based morphology assay to capture a broad range of phenotypes, we assessed the effect of 12768 pairwise RNAi perturbations in six different conditions. We found that genetic interactions form in different trajectories and developed an algorithm, termed MODIFI, to analyze how genetic interactions rewire over time. Using this framework, we identified more statistically significant interactions compared to end-point assays and further observed several examples of context-dependent crosstalk between signaling pathways such as an interaction between Ras and Rel which is dependent on MEK activity.

Editorial note

This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

The Genetic Architecture of Intra-Species Hybrid Mito-Nuclear Epistasis

Genetic variants that are neutral within, but deleterious between, populations (Dobzhansky-Muller Incompatibilities) are thought to initiate hybrid dysfunction and then to accumulate and complete the speciation process. To identify the types of genetic differences that might initiate speciation, it is useful to study inter-population (intra-species) hybrids that exhibit reduced fitness. In Caenorhabditis briggsae, a close relative of the nematode C. elegans, such minor genetic incompatibilities have been identified. One incompatibility between the mitochondrial and nuclear genomes reduces the fitness of some hybrids. To understand the nuclear genetic architecture of this epistatic interaction, we constructed two sets of recombinant inbred lines by hybridizing two genetically diverse wild populations. In such lines, selection is able to eliminate deleterious combinations of alleles derived from the two parental populations. The genotypes of surviving hybrid lines thus reveal favorable allele combinations at loci experiencing selection. Our genotype data from the resulting lines are consistent with the interpretation that the X alleles participate in epistatic interactions with autosomes and the mitochondrial genome. We evaluate this possibility given predictions that mitochondria-X epistasis should be more prevalent than between mitochondria and autosomes. Our empirical identification of inter-genomic linkage disequilibrium supports the body of literature indicating that the accumulation of mito-nuclear genetic incompatibilities might initiate the speciation process through the generation of less-fit inter-population hybrids.

Analysis of Epistasis in Natural Traits Using Model Organisms

The ability to detect and understand epistasis in natural populations is important for understanding how biological traits are influenced by genetic variation. However, identification and characterization of epistasis in natural populations remains difficult due to statistical issues that arise as a result of multiple comparisons, and the fact that most genetic variants segregate at low allele frequencies. In this review, we discuss how model organisms may be used to manipulate genotypic combinations to power the detection of epistasis as well as test interactions between specific genes. Findings from a number of species indicate that statistical epistasis is pervasive between natural genetic variants. However, the properties of experimental systems that enable analysis of epistasis also constrain extrapolation of these results back into natural populations.

The complex underpinnings of genetic background effects

Genetic interactions between mutations and standing polymorphisms can cause mutations to show distinct phenotypic effects in different individuals. To characterize the genetic architecture of these so-called background effects, we genotype 1411 wild-type and mutant yeast cross progeny and measure their growth in 10 environments. Using these data, we map 1086 interactions between segregating loci and 7 different gene knockouts. Each knockout exhibits between 73 and 543 interactions, with 89% of all interactions involving higher-order epistasis between a knockout and multiple loci. Identified loci interact with as few as one knockout and as many as all seven knockouts. In mutants, loci interacting with fewer and more knockouts tend to show enhanced and reduced phenotypic effects, respectively. Cross–environment analysis reveals that most interactions between the knockouts and segregating loci also involve the environment. These results illustrate the complicated interactions between mutations, standing polymorphisms, and the environment that cause background effects.

Mutations often show distinct phenotypic effects across different genetic backgrounds. Here the authors describe the genetic basis of these so-called background effects using data on genotype and growth in 10 environments from 1411 segregants from a cross of two strains of budding yeast.

Genetic variants in major depressive disorder: From pathophysiology to therapy

In spite of promising preclinical results there is a decreasing number of new registered medications in major depression. The main reason behind this fact is the lack of confirmation in clinical studies for the assumed, and in animals confirmed, therapeutic results. This suggests low predictive value of animal studies for central nervous system disorders. One solution for identifying new possible targets is the application of genetics and genomics, which may pinpoint new targets based on the effect of genetic variants in humans. The present review summarizes such research focusing on depression and its therapy. The inconsistency between most genetic studies in depression suggests, first of all, a significant role of environmental stress. Furthermore, effect of individual genes and polymorphisms is weak, therefore gene x gene interactions or complete biochemical pathways should be analyzed. Even genes encoding target proteins of currently used antidepressants remain non-significant in genome-wide case control investigations suggesting no main effect in depression, but rather an interaction with stress. The few significant genes in GWASs are related to neurogenesis, neuronal synapse, cell contact and DNA transcription and as being nonspecific for depression are difficult to harvest pharmacologically. Most candidate genes in replicable gene x environment interactions, on the other hand, are connected to the regulation of stress and the HPA axis and thus could serve as drug targets for depression subgroups characterized by stress-sensitivity and anxiety while other risk polymorphisms such as those related to prominent cognitive symptoms in depression may help to identify additional subgroups and their distinct treatment. Until these new targets find their way into therapy, the optimization of current medications can be approached by pharmacogenomics, where metabolizing enzyme polymorphisms remain prominent determinants of therapeutic success.

An Ultrahigh-Dimensional Mapping Model of High-order Epistatic Networks for Complex Traits

BACKGROUND

Genetic interactions involving more than two loci have been thought to affect quantitatively inherited traits and diseases more pervasively than previously appreciated. However, the detection of such high-order interactions to chart a complete portrait of genetic architecture has not been well explored.

METHODS

We present an ultrahigh-dimensional model to systematically characterize genetic main effects and interaction effects of various orders among all possible markers in a genetic mapping or association study. The model was built on the extension of a variable selection procedure, called iFORM, derived from forward selection. The model shows its unique power to estimate the magnitudes and signs of high-order epistatic effects, in addition to those of main effects and pairwise epistatic effects.

RESULTS

The statistical properties of the model were tested and validated through simulation studies. By analyzing a real data for shoot growth in a mapping population of woody plant, mei (Prunus mume), we demonstrated the usefulness and utility of the model in practical genetic studies. The model has identified important high-order interactions that contribute to shoot growth for mei.

CONCLUSION

The model provides a tool to precisely construct genotype-phenotype maps for quantitative traits by identifying any possible high-order epistasis which is often ignored in the current genetic literature.

Strategies for integrated analysis in imaging genetics studies

Imaging Genetics (IG) integrates neuroimaging and genomic data from the same individual, deepening our knowledge of the biological mechanisms behind neurodevelopmental domains and neurological disorders. Although the literature on IG has exponentially grown over the past years, the majority of studies have mainly analyzed associations between candidate brain regions and individual genetic variants. However, this strategy is not designed to deal with the complexity of neurobiological mechanisms underlying behavioral and neurodevelopmental domains. Moreover, larger sample sizes and increased multidimensionality of this type of data represents a challenge for standardizing modeling procedures in IG research. This review provides a systematic update of the methods and strategies currently used in IG studies, and serves as an analytical framework for researchers working in this field. To complement the functionalities of the Neuroconductor framework, we also describe existing R packages that implement these methodologies. In addition, we present an overview of how these methodological approaches are applied in integrating neuroimaging and genetic data.

Effect of genetic architecture on the prediction accuracy of quantitative traits in samples of unrelated individuals

Predicting complex phenotypes from genomic data is a fundamental aim of animal and plant breeding, where we wish to predict genetic merits of selection candidates; and of human genetics, where we wish to predict disease risk. While genomic prediction models work well with populations of related individuals and high linkage disequilibrium (LD) (e.g., livestock), comparable models perform poorly for populations of unrelated individuals and low LD (e.g., humans). We hypothesized that low prediction accuracies in the latter situation may occur when the genetics architecture of the trait departs from the infinitesimal and additive architecture assumed by most prediction models. We used simulated data for 10,000 lines based on sequence data from a population of unrelated, inbred Drosophila melanogaster lines to evaluate this hypothesis. We show that, even in very simplified scenarios meant as a stress test of the commonly used Genomic Best Linear Unbiased Predictor (G-BLUP) method, using all common variants yields low prediction accuracy regardless of the trait genetic architecture. However, prediction accuracy increases when predictions are informed by the genetic architecture inferred from mapping the top variants affecting main effects and interactions in the training data, provided there is sufficient power for mapping. When the true genetic architecture is largely or partially due to epistatic interactions, the additive model may not perform well, while models that account explicitly for interactions generally increase prediction accuracy. Our results indicate that accounting for genetic architecture can improve prediction accuracy for quantitative traits.

Decanalizing thinking on genetic canalization

The concept of genetic canalization has had an abiding influence on views of complex-trait evolution. A genetically canalized system has evolved to become less sensitive to the effects of mutation. When a gene product that supports canalization is compromised, the phenotypic impacts of a mutation should be more pronounced. This expected increase in mutational effects not only has important consequences for evolution, but has also motivated strategies to treat disease. However, recent studies demonstrate that, when putative agents of genetic canalization are impaired, systems do not behave as expected. Here, we review the evidence that is used to infer whether particular gene products are agents of genetic canalization. Then we explain how such inferences often succumb to a converse error. We go on to show that several candidate agents of genetic canalization increase the phenotypic impacts of some mutations while decreasing the phenotypic impacts of others. These observations suggest that whether a gene product acts as a 'buffer' (lessening mutational effects) or a 'potentiator' (increasing mutational effects) is not a fixed property of the gene product but instead differs for the different mutations with which it interacts. To investigate features of genetic interactions that might predispose them toward buffering versus potentiation, we explore simulated gene-regulatory networks. Similarly to putative agents of genetic canalization, the gene products in simulated networks also modify the phenotypic effects of mutations in other genes without a strong overall tendency towards lessening or increasing these effects. In sum, these observations call into question whether complex traits have evolved to become less sensitive (i.e., are canalized) to genetic change, and the degree to which trends exist that predict how one genetic change might alter another's impact. We conclude by discussing approaches to address these and other open questions that are brought into focus by re-thinking genetic canalization.

Haplotype-Based Genome-Wide Prediction Models Exploit Local Epistatic Interactions Among Markers

Genome-wide prediction approaches represent versatile tools for the analysis and prediction of complex traits. Mostly they rely on marker-based information, but scenarios have been reported in which models capitalizing on closely-linked markers that were combined into haplotypes outperformed marker-based models. Detailed comparisons were undertaken to reveal under which circumstances haplotype-based genome-wide prediction models are superior to marker-based models. Specifically, it was of interest to analyze whether and how haplotype-based models may take local epistatic effects between markers into account. Assuming that populations consisted of fully homozygous individuals, a marker-based model in which local epistatic effects inside haplotype blocks were exploited (LEGBLUP) was linearly transformable into a haplotype-based model (HGBLUP). This theoretical derivation formally revealed that haplotype-based genome-wide prediction models capitalize on local epistatic effects among markers. Simulation studies corroborated this finding. Due to its computational efficiency the HGBLUP model promises to be an interesting tool for studies in which ultra-high-density SNP data sets are studied. Applying the HGBLUP model to empirical data sets revealed higher prediction accuracies than for marker-based models for both traits studied using a mouse panel. In contrast, only a small subset of the traits analyzed in crop populations showed such a benefit. Cases in which higher prediction accuracies are observed for HGBLUP than for marker-based models are expected to be of immediate relevance for breeders, due to the tight linkage a beneficial haplotype will be preserved for many generations. In this respect the inheritance of local epistatic effects very much resembles the one of additive effects.

Systematic analysis of complex genetic interactions

To systematically explore complex genetic interactions, we constructed ~200,000 yeast triple mutants and scored negative trigenic interactions. We selected double-mutant query genes across a broad spectrum of biological processes, spanning a range of quantitative features of the global digenic interaction network and tested for a genetic interaction with a third mutation. Trigenic interactions often occurred among functionally related genes, and essential genes were hubs on the trigenic network. Despite their functional enrichment, trigenic interactions tended to link genes in distant bioprocesses and displayed a weaker magnitude than digenic interactions. We estimate that the global trigenic interaction network is ~100 times as large as the global digenic network, highlighting the potential for complex genetic interactions to affect the biology of inheritance, including the genotype-to-phenotype relationship.

Facilitating Anti-Cancer Combinatorial Drug Discovery by Targeting Epistatic Disease Genes

Due to synergistic effects, combinatorial drugs are widely used for treating complex diseases. However, combining drugs and making them synergetic remains a challenge. Genetic disease genes are considered a promising source of drug targets with important implications for navigating the drug space. Most diseases are not caused by a single pathogenic factor, but by multiple disease genes, in particular, interacting disease genes. Thus, it is reasonable to consider that targeting epistatic disease genes may enhance the therapeutic effects of combinatorial drugs. In this study, synthetic lethality gene pairs of tumors, similar to epistatic disease genes, were first targeted by combinatorial drugs, resulting in the enrichment of the combinatorial drugs with cancer treatment, which verified our hypothesis. Then, conventional epistasis detection software was used to identify epistatic disease genes from the genome wide association studies (GWAS) dataset. Furthermore, combinatorial drugs were predicted by targeting these epistatic disease genes, and five combinations were proven to have synergistic anti-cancer effects on MCF-7 cells through cell cytotoxicity assay. Combined with the three-dimensional (3D) genome-based method, the epistatic disease genes were filtered and were more closely related to disease. By targeting the filtered gene pairs, the efficiency of combinatorial drug discovery has been further improved.

Gene-gene and gene-environment interactions in complex traits in yeast

One of the fundamental questions in biology is how the genotype regulates the phenotype. An increasing number of studies indicate that, in most cases, the effect of a genetic locus on the phenotype is context-dependent, i.e. it is influenced by the genetic background and the environment in which the phenotype is measured. Still, the majority of the studies, in both model organisms and humans, that map the genetic regulation of phenotypic variation in complex traits primarily identify additive loci with independent effects. This does not reflect an absence of the contribution of genetic interactions to phenotypic variation, but instead is a consequence of the technical limitations in mapping gene-gene interactions (GGI) and gene-environment interactions (GEI). Yeast, with its detailed molecular understanding, diverse population genomics and ease of genetic manipulation, is a unique and powerful resource to study the contributions of GGI and GEI in the regulation of phenotypic variation. Here we review studies in yeast that have identified GGI and GEI that regulate phenotypic variation, and discuss the contribution of these findings in explaining missing heritability of complex traits, and how observations from these GGI and GEI studies enhance our understanding of the mechanisms underlying genetic robustness and adaptability that shape the architecture of the genotype-phenotype map.

On the relationship between epistasis and genetic variance heterogeneity

Epistasis and genetic variance heterogeneity are two non-additive genetic inheritance patterns that are often, but not always, related. Here we use theoretical examples and empirical results from earlier analyses of experimental data to illustrate the connection between the two. This includes an introduction to the relationship between epistatic gene action, statistical epistasis, and genetic variance heterogeneity, and a brief discussion about how genetic processes other than epistasis can also give rise to genetic variance heterogeneity.

How well do you know your mutation? Complex effects of genetic background on expressivity, complementation, and ordering of allelic effects

For a given gene, different mutations influence organismal phenotypes to varying degrees. However, the expressivity of these variants not only depends on the DNA lesion associated with the mutation, but also on factors including the genetic background and rearing environment. The degree to which these factors influence related alleles, genes, or pathways similarly, and whether similar developmental mechanisms underlie variation in the expressivity of a single allele across conditions and among alleles is poorly understood. Besides their fundamental biological significance, these questions have important implications for the interpretation of functional genetic analyses, for example, if these factors alter the ordering of allelic series or patterns of complementation. We examined the impact of genetic background and rearing environment for a series of mutations spanning the range of phenotypic effects for both the scalloped and vestigial genes, which influence wing development in Drosophila melanogaster. Genetic background and rearing environment influenced the phenotypic outcome of mutations, including intra-genic interactions, particularly for mutations of moderate expressivity. We examined whether cellular correlates (such as cell proliferation during development) of these phenotypic effects matched the observed phenotypic outcome. While cell proliferation decreased with mutations of increasingly severe effects, surprisingly it did not co-vary strongly with the degree of background dependence. We discuss these findings and propose a phenomenological model to aid in understanding the biology of genes, and how this influences our interpretation of allelic effects in genetic analysis.

Different mutations in a gene, or in genes with related functions, can have effects of varying severity. Studying sets of mutations and analyzing how they interact are essential components of a geneticist's toolkit. However, the effects caused by a mutation depend not only on the mutation itself, but on additional genetic variation throughout an organism's genome and on the environment that organism has experienced. Therefore, identifying how the genomic and environmental context alter the expression of mutations is critical for making reliable inferences about how genes function. Yet studies on this context dependence have largely been limited to single mutations in single genes. We examined how the genomic and environmental context influence the expression of multiple mutations in two related genes affecting the fruit fly wing. Our results show that the genetic and environmental context generally affect the expression of related mutations in similar ways. However, the interactions between two different mutations in a single gene sometimes depended strongly on context. In addition, cell proliferation in the developing wing and adult wing size were not affected by the genetic and environmental context in similar ways in mutant flies, suggesting that variation in cell growth cannot fully explain how mutations affect wings. Overall, our findings show that context can have a big impact on the interpretation of genetic experiments, including how we draw conclusions about gene function and cause-and-effect relationships.

Towards Precision Medicine for Hypertension: A Review of Genomic, Epigenomic, and Microbiomic Effects on Blood Pressure in Experimental Rat Models and Humans

Compelling evidence for the inherited nature of essential hypertension has led to extensive research in rats and humans. Rats have served as the primary model for research on the genetics of hypertension resulting in identification of genomic regions that are causally associated with hypertension. In more recent times, genome-wide studies in humans have also begun to improve our understanding of the inheritance of polygenic forms of hypertension. Based on the chronological progression of research into the genetics of hypertension as the "structural backbone," this review catalogs and discusses the rat and human genetic elements mapped and implicated in blood pressure regulation. Furthermore, the knowledge gained from these genetic studies that provide evidence to suggest that much of the genetic influence on hypertension residing within noncoding elements of our DNA and operating through pervasive epistasis or gene-gene interactions is highlighted. Lastly, perspectives on current thinking that the more complex "triad" of the genome, epigenome, and the microbiome operating to influence the inheritance of hypertension, is documented. Overall, the collective knowledge gained from rats and humans is disappointing in the sense that major hypertension-causing genes as targets for clinical management of essential hypertension may not be a clinical reality. On the other hand, the realization that the polygenic nature of hypertension prevents any single locus from being a relevant clinical target for all humans directs future studies on the genetics of hypertension towards an individualized genomic approach.

Widespread epistasis regulates glucose homeostasis and gene expression

The relative contributions of additive versus non-additive interactions in the regulation of complex traits remains controversial. This may be in part because large-scale epistasis has traditionally been difficult to detect in complex, multi-cellular organisms. We hypothesized that it would be easier to detect interactions using mouse chromosome substitution strains that simultaneously incorporate allelic variation in many genes on a controlled genetic background. Analyzing metabolic traits and gene expression levels in the offspring of a series of crosses between mouse chromosome substitution strains demonstrated that inter-chromosomal epistasis was a dominant feature of these complex traits. Epistasis typically accounted for a larger proportion of the heritable effects than those due solely to additive effects. These epistatic interactions typically resulted in trait values returning to the levels of the parental CSS host strain. Due to the large epistatic effects, analyses that did not account for interactions consistently underestimated the true effect sizes due to allelic variation or failed to detect the loci controlling trait variation. These studies demonstrate that epistatic interactions are a common feature of complex traits and thus identifying these interactions is key to understanding their genetic regulation.

Most complex traits and diseases are regulated by the combined influence of multiple genetic variants. However, it remains controversial whether these genetic variants independently influence complex traits, and therefore the impact of each variant could be simply added together (additivity), or whether the variants work together to influence trait variation, in which case the combined impact of multiple variants would differ from the summed impact of each individual variant (epistasis). In this study in mice, we discovered that the genetic regulation of blood sugar levels and gene expression in the liver were predominantly controlled by non-additive interactions, whereas body weight was predominantly controlled by additive interactions. Remarkably, the expression level of nearly 25% of all genes in the liver was controlled by non-additive interactions. The non-additive interactions typically acted to return trait values to the levels detected in control mice, thus contributing to a reduction in trait variation. We also demonstrated that not accounting for non-additive interactions significantly underestimated the phenotypic effect of a genetic variant on a particular genetic background, suggesting that many previously identified risk loci may have significantly larger effects on disease susceptibility in a subset of individuals. These studies highlight the importance of understanding interactions between genetic variants to better understand disease risk and personalize clinical care.

xSyn: A Software Tool for Identifying Sophisticated 3-Way Interactions From Cancer Expression Data

Background:

Constructing gene co-expression networks from cancer expression data is important for investigating the genetic mechanisms underlying cancer. However, correlation coefficients or linear regression models are not able to model sophisticated relationships among gene expression profiles. Here, we address the 3-way interaction that 2 genes’ expression levels are clustered in different space locations under the control of a third gene’s expression levels.

Results:

We present xSyn, a software tool for identifying such 3-way interactions from cancer gene expression data based on an optimization procedure involving the usage of UPGMA (Unweighted Pair Group Method with Arithmetic Mean) and synergy. The effectiveness is demonstrated by application to 2 real gene expression data sets.

Conclusions:

xSyn is a useful tool for decoding the complex relationships among gene expression profiles. xSyn is available at http://www.bdxconsult.com/xSyn.html.

Mapping a diversity of genetic interactions in yeast

Genetic interactions occur when the combination of multiple mutations yields an unexpected phenotype, and they may confound our ability to fully understand the genetic mechanisms underlying complex diseases. Genetic interactions are challenging to study because there are millions of possible different variant combinations within a given genome. Consequently, they have primarily been systematically explored in unicellular model organisms, such as yeast, with a focus on pairwise genetic interactions between loss-of-function alleles. However, there are many different types of genetic interactions, such as those occurring between gain-of-function or heterozygous mutations. Here, we review recent advances made in the systematic analysis of such diverse genetic interactions in yeast, and briefly discuss how similar studies could be undertaken in human cells.

Transcriptional Profiling of Biofilm Regulators Identified by an Overexpression Screen in Saccharomyces cerevisiae

Biofilm formation by microorganisms is a major cause of recurring infections and removal of biofilms has proven to be extremely difficult given their inherent drug resistance . Understanding the biological processes that underlie biofilm formation is thus extremely important and could lead to the development of more effective drug therapies, resulting in better infection outcomes. Using the yeast Saccharomyces cerevisiae as a biofilm model, overexpression screens identified DIG1, SFL1, HEK2, TOS8, SAN1, and ROF1/YHR177W as regulators of biofilm formation. Subsequent RNA-seq analysis of biofilm and nonbiofilm-forming strains revealed that all of the overexpression strains, other than DIG1 and TOS8, were adopting a single differential expression profile, although induced to varying degrees. TOS8 adopted a separate profile, while the expression profile of DIG1 reflected the common pattern seen in most of the strains, plus substantial DIG1-specific expression changes. We interpret the existence of the common transcriptional pattern seen across multiple, unrelated overexpression strains as reflecting a transcriptional state, that the yeast cell can access through regulatory signaling mechanisms, allowing an adaptive morphological change between biofilm-forming and nonbiofilm states.

Mechanisms of suppression: The wiring of genetic resilience

Recent analysis of genome sequences has identified individuals that are healthy despite carrying severe disease-associated mutations. A possible explanation is that these individuals carry a second genomic perturbation that can compensate for the detrimental effects of the disease allele, a phenomenon referred to as suppression. In model organisms, suppression interactions are generally divided into two classes: genomic suppressors which are secondary mutations in the genome that bypass a mutant phenotype, and dosage suppression interactions in which overexpression of a suppressor gene rescues a mutant phenotype. Here, we describe the general properties of genomic and dosage suppression, with an emphasis on the budding yeast. We propose that suppression interactions between genetic variants are likely relevant for determining the penetrance of human traits. Consequently, an understanding of suppression mechanisms may guide the discovery of protective variants in healthy individuals that carry disease alleles, which could direct the rational design of new therapeutics.

HiSeeker: Detecting High-Order SNP Interactions Based on Pairwise SNP Combinations

Detecting single nucleotide polymorphisms’ (SNPs) interaction is one of the most popular approaches for explaining the missing heritability of common complex diseases in genome-wide association studies. Many methods have been proposed for SNP interaction detection, but most of them only focus on pairwise interactions and ignore high-order ones, which may also contribute to complex traits. Existing methods for high-order interaction detection can hardly handle genome-wide data and suffer from low detection power, due to the exponential growth of search space. In this paper, we proposed a flexible two-stage approach (called HiSeeker) to detect high-order interactions. In the screening stage, HiSeeker employs the chi-squared test and logistic regression model to efficiently obtain candidate pairwise combinations, which have intermediate or significant associations with the phenotype for interaction detection. In the search stage, two different strategies (exhaustive search and ant colony optimization-based search) are utilized to detect high-order interactions from candidate combinations. The experimental results on simulated datasets demonstrate that HiSeeker can more efficiently and effectively detect high-order interactions than related representative algorithms. On two real case-control datasets, HiSeeker also detects several significant high-order interactions, whose individual SNPs and pairwise interactions have no strong main effects or pairwise interaction effects, and these high-order interactions can hardly be identified by related algorithms.

Genetic suppression: Extending our knowledge from lab experiments to natural populations

Many mutations have deleterious phenotypic effects that can be alleviated by suppressor mutations elsewhere in the genome. High-throughput approaches have facilitated the large-scale identification of these suppressors and have helped shed light on core functional mechanisms that give rise to suppression. Following reports that suppression occurs naturally within species, it is important to determine how our understanding of this phenomenon based on lab experiments extends to genetically diverse natural populations. Although suppression is typically mediated by individual genetic changes in lab experiments, recent studies have shown that suppression in natural populations can involve combinations of genetic variants. This difference in complexity suggests that sets of variants can exhibit similar functional effects to individual suppressors found in lab experiments. In this review, we discuss how characterizing the way in which these variants jointly lead to suppression could provide important insights into the genotype-phenotype map that are relevant to evolution and health.

REVIEW: Epistasis and dominance in the emergence of catalytic function as exemplified by the evolution of plant terpene synthases

Epistasis, the interaction between mutations and the genetic background, is a pervasive force in evolution that is difficult to predict yet derives from a simple principle - biological systems are interconnected. Therefore, one effect may be intimately linked to another, hence interdependent. Untangling epistatic interactions between and within genes is a vibrant area of research. Deriving a mechanistic understanding of epistasis is a major challenge. Particularly, elucidating how epistasis can attenuate the effects of otherwise dominant mutations that control phenotypes. Using the emergence of terpene cyclization in specialized metabolism as an excellent example, this review describes the process of discovery and interpretation of dominance and epistasis in relation to current efforts. Specifically, we outline experimental approaches to isolating epistatic networks of mutations in protein structure, formally quantifying epistatic interactions, then building biochemical models with chemical mechanisms in efforts to achieve an understanding of the physical basis for epistasis. From these models we describe informed conjectures about past evolutionary events that underlie the emergence, divergence and specialization of terpene synthases to illustrate key principles of the constraining forces of epistasis in enzyme function.

Detecting High-Order Epistasis in Nonlinear Genotype-Phenotype Maps

High-order epistasis has been observed in many genotype-phenotype maps. These multi-way interactions between mutations may be useful for dissecting complex traits and could have profound implications for evolution. Alternatively, they could be a statistical artifact. High-order epistasis models assume the effects of mutations should add, when they could in fact multiply or combine in some other nonlinear way. A mismatch in the “scale” of the epistasis model and the scale of the underlying map would lead to spurious epistasis. In this article, we develop an approach to estimate the nonlinear scales of arbitrary genotype-phenotype maps. We can then linearize these maps and extract high-order epistasis. We investigated seven experimental genotype-phenotype maps for which high-order epistasis had been reported previously. We find that five of the seven maps exhibited nonlinear scales. Interestingly, even after accounting for nonlinearity, we found statistically significant high-order epistasis in all seven maps. The contributions of high-order epistasis to the total variation ranged from 2.2 to 31.0%, with an average across maps of 12.7%. Our results provide strong evidence for extensive high-order epistasis, even after nonlinear scale is taken into account. Further, we describe a simple method to estimate and account for nonlinearity in genotype-phenotype maps.

Spiraling Complexity: A Test of the Snowball Effect in a Computational Model of RNA Folding

Genetic incompatibilities can emerge as a byproduct of genetic divergence. According to Dobzhansky and Muller, an allele that fixes in one population may be incompatible with an allele at a different locus in another population when the two alleles are brought together in hybrids. Orr showed that the number of Dobzhansky-Muller incompatibilities (DMIs) should accumulate faster than linearly-i.e., snowball-as two lineages diverge. Several studies have attempted to test the snowball effect using data from natural populations. One limitation of these studies is that they have focused on predictions of the Orr model, but not on its underlying assumptions. Here, we use a computational model of RNA folding to test both predictions and assumptions of the Orr model. Two populations are allowed to evolve in allopatry on a holey fitness landscape. We find that the number of inviable introgressions (an indicator for the number of DMIs) snowballs, but does so more slowly than expected. We show that this pattern is explained, in part, by the fact that DMIs can disappear after they have arisen, contrary to the assumptions of the Orr model. This occurs because DMIs become progressively more complex (i.e., involve alleles at more loci) as a result of later substitutions. We also find that most DMIs involve >2 loci, i.e., they are complex. Reproductive isolation does not snowball because DMIs do not act independently of each other. We conclude that the RNA model supports the central prediction of the Orr model that the number of DMIs snowballs, but challenges other predictions, as well as some of its underlying assumptions.

A global test for gene-gene interactions based on random matrix theory

Statistical interactions between markers of genetic variation, or gene-gene interactions, are believed to play an important role in the etiology of many multifactorial diseases and other complex phenotypes. Unfortunately, detecting gene-gene interactions is extremely challenging due to the large number of potential interactions and ambiguity regarding marker coding and interaction scale. For many data sets, there is insufficient statistical power to evaluate all candidate gene-gene interactions. In these cases, a global test for gene-gene interactions may be the best option. Global tests have much greater power relative to multiple individual interaction tests and can be used on subsets of the markers as an initial filter prior to testing for specific interactions. In this paper, we describe a novel global test for gene-gene interactions, the global epistasis test (GET), that is based on results from random matrix theory. As we show via simulation studies based on previously proposed models for common diseases including rheumatoid arthritis, type 2 diabetes, and breast cancer, our proposed GET method has superior performance characteristics relative to existing global gene-gene interaction tests. A glaucoma GWAS data set is used to demonstrate the practical utility of the GET method.

Modifiers of the Genotype-Phenotype Map: Hsp90 and Beyond

Disruption of certain genes alters the heritable phenotypic variation among individuals. Research on the chaperone Hsp90 has played a central role in determining the genetic basis of this phenomenon, which may be important to evolution and disease. Key studies have shown that Hsp90 perturbation modifies the effects of many genetic variants throughout the genome. These modifications collectively transform the genotype–phenotype map, often resulting in a net increase or decrease in heritable phenotypic variation. Here, we summarize some of the foundational work on Hsp90 that led to these insights, discuss a framework for interpreting this research that is centered upon the standard genetics concept of epistasis, and propose major questions that future studies in this area should address.

Gene-Environment Interactions in Stress Response Contribute Additively to a Genotype-Environment Interaction

How combinations of gene-environment interactions collectively give rise to genotype-environment interactions is not fully understood. To shed light on this problem, we genetically dissected an environment-specific poor growth phenotype in a cross of two budding yeast strains. This phenotype is detectable when certain segregants are grown on ethanol at 37°C (‘E37’), a condition that differs from the standard culturing environment in both its carbon source (ethanol as opposed to glucose) and temperature (37°C as opposed to 30°C). Using recurrent backcrossing with phenotypic selection, we identified 16 contributing loci. To examine how these loci interact with each other and the environment, we focused on a subset of four loci that together can lead to poor growth in E37. We measured the growth of all 16 haploid combinations of alleles at these loci in all four possible combinations of carbon source (ethanol or glucose) and temperature (30 or 37°C) in a nearly isogenic population. This revealed that the four loci act in an almost entirely additive manner in E37. However, we also found that these loci have weaker effects when only carbon source or temperature is altered, suggesting that their effect magnitudes depend on the severity of environmental perturbation. Consistent with such a possibility, cloning of three causal genes identified factors that have unrelated functions in stress response. Thus, our results indicate that polymorphisms in stress response can show effects that are intensified by environmental stress, thereby resulting in major genotype-environment interactions when multiple of these variants co-occur.

Determining the genetic and molecular mechanisms that give rise to genotype-environment interaction (‘GxE’) is important for many areas of biology, including agriculture, evolution, and medicine. To help advance knowledge regarding this topic, we dissect the genetic basis of an example of GxE in which certain Saccharomyces cerevisiae cross progeny show extremely poor growth specifically on ethanol at 37°C. This environment differs from the standard condition used for culturing budding yeast in both its carbon source (ethanol as opposed to glucose) and temperature (37°C as opposed to 30°C). We provide evidence that poor growth on ethanol at 37°C is caused by a number of predominantly additive loci that individually exhibit gene-environment interactions with both carbon source and temperature. These loci show their largest effects when carbon source and temperature are simultaneously modified, indicating their effect magnitudes may be influenced by the severity of environmental stress. Consistent with this possibility, we clone three causal genes and find they encode functionally unrelated components of stress response. Our work suggests that polymorphisms in stress response can contribute additively to genotype-environment interactions that vary in intensity across conditions in a stress level-dependent manner.

Diverse genetic architectures lead to the same cryptic phenotype in a yeast cross

Cryptic genetic variants that do not typically influence traits can interact epistatically with each other and mutations to cause unexpected phenotypes. To improve understanding of the genetic architectures and molecular mechanisms that underlie these interactions, we comprehensively dissected the genetic bases of 17 independent instances of the same cryptic colony phenotype in a yeast cross. In eight cases, the phenotype resulted from a genetic interaction between a de novo mutation and one or more cryptic variants. The number and identities of detected cryptic variants depended on the mutated gene. In the nine remaining cases, the phenotype arose without a de novo mutation due to two different classes of higher-order genetic interactions that only involve cryptic variants. Our results may be relevant to other species and disease, as most of the mutations and cryptic variants identified in our study reside in components of a partially conserved and oncogenic signalling pathway.

Cryptic genetic variants may not individually show discernible phenotypic effects, but collectively, these polymorphisms can lead to unexpected, genetically complex traits that might be relevant to evolution and disease. Here, the authors use large yeast populations to comprehensively dissect the genetic bases of 17 independent occurrences of a phenotype that arises due to combinations of epistatically interacting cryptic variants.

Local fitness landscape of the green fluorescent protein

Fitness landscapes1,2, depictions of how genotypes manifest at the phenotypic level, form the basis for our understanding of many areas of biology2–7 yet their properties remain elusive. Studies addressing this issue often consider specific genes and their function as proxy for fitness2,4, experimentally assessing the impact on function of single mutations and their combinations in a specific sequence2,5,8–15 or in different sequences2,3,5,16–18. However, systematic high-throughput studies of the local fitness landscape of an entire protein have not yet been reported. Here, we chart an extensive region of the local fitness landscape of the green fluorescent protein from Aequorea victoria (avGFP) by measuring the native function, fluorescence, of tens of thousands of derivative genotypes of avGFP. We find that its fitness landscape is narrow, with half of genotypes with two mutations showing reduced fluorescence and half of genotypes with five mutations being completely non-fluorescent. The narrowness is enhanced by epistasis, which was detected in up to 30% of genotypes with multiple mutations arising mostly through the cumulative impact of slightly deleterious mutations causing a threshold-like decrease of protein stability and concomitant loss of fluorescence. A model of orthologous sequence divergence spanning hundreds of millions of years predicted the extent of epistasis in our data, indicating congruence between the fitness landscape properties at the local and global scales. The characterization of the local fitness landscape of avGFP has important implications for a number of fields including molecular evolution, population genetics and protein design.

Multi-locus Genotypes Underlying Temperature Sensitivity in a Mutationally Induced Trait

Determining how genetic variation alters the expression of heritable phenotypes across conditions is important for agriculture, evolution, and medicine. Central to this problem is the concept of genotype-by-environment interaction (or 'GxE'), which occurs when segregating genetic variation causes individuals to show different phenotypic responses to the environment. While many studies have sought to identify individual loci that contribute to GxE, obtaining a deeper understanding of this phenomenon may require defining how sets of loci collectively alter the relationship between genotype, environment, and phenotype. Here, we identify combinations of alleles at seven loci that control how a mutationally induced colony phenotype is expressed across a range of temperatures (21, 30, and 37 °C) in a panel of yeast recombinants. We show that five predominant multi-locus genotypes involving the detected loci result in trait expression with varying degrees of temperature sensitivity. By comparing these genotypes and their patterns of trait expression across temperatures, we demonstrate that the involved alleles contribute to temperature sensitivity in different ways. While alleles of the transcription factor MSS11 specify the potential temperatures at which the trait can occur, alleles at the other loci modify temperature sensitivity within the range established by MSS11 in a genetic background- and/or temperature-dependent manner. Our results not only represent one of the first characterizations of GxE at the resolution of multi-locus genotypes, but also provide an example of the different roles that genetic variants can play in altering trait expression across conditions.

Prediction of Synergism from Chemical-Genetic Interactions by Machine Learning

The structure of genetic interaction networks predicts that, analogous to synthetic lethal interactions between non-essential genes, combinations of compounds with latent activities may exhibit potent synergism. To test this hypothesis, we generated a chemical-genetic matrix of 195 diverse yeast deletion strains treated with 4,915 compounds. This approach uncovered 1,221 genotype-specific inhibitors, which we termed cryptagens. Synergism between 8,128 structurally disparate cryptagen pairs was assessed experimentally and used to benchmark predictive algorithms. A model based on the chemical-genetic matrix and the genetic interaction network failed to accurately predict synergism. However, a combined random forest and Naive Bayesian learner that associated chemical structural features with genotype-specific growth inhibition had strong predictive power. This approach identified previously unknown compound combinations that exhibited species-selective toxicity toward human fungal pathogens. This work demonstrates that machine learning methods trained on unbiased chemical-genetic interaction data may be widely applicable for the discovery of synergistic combinations in different species.

A new conceptual framework for investigating complex genetic disease

Some common diseases are known to have an inherited component, however, their population- and familial-incidence patterns do not conform to any known monogenic Mendelian pattern of inheritance and instead they are currently much better explained if an underlying polygenic architecture is posited. Studies that have attempted to identify the causative genetic factors have been designed on this polygenic framework, but so far the yield has been largely unsatisfactory. Based on accumulating recent observations concerning the roles of somatic mosaicism in disease, in this article a second framework which posits a single gene-two hit model which can be modulated by a mutator/anti-mutator genetic background is suggested. I discuss whether such a model can be considered a viable alternative based on current knowledge, its advantages over the current polygenic framework, and describe practical routes via which the new framework can be investigated.

The complex genetic and molecular basis of a model quantitative trait

Sixty‐four genomic loci and seven genes that contribute to heritable variation in a model quantitative trait—resistance to oxidative stress—are identified across three yeast strains. The high‐resolution understanding of this phenotype provides new insight into the genetic and molecular basis of quantitative traits.

Quantitative traits are often influenced by many loci with small effects. Identifying most of these loci and resolving them to specific genes or genetic variants is challenging. Yet, achieving such a detailed understanding of quantitative traits is important, as it can improve our knowledge of the genetic and molecular basis of heritable phenotypic variation. In this study, we use a genetic mapping strategy that involves recurrent backcrossing with phenotypic selection to obtain new insights into an ecologically, industrially, and medically relevant quantitative trait—tolerance of oxidative stress, as measured based on resistance to hydrogen peroxide. We examine the genetic basis of hydrogen peroxide resistance in three related yeast crosses and detect 64 distinct genomic loci that likely influence the trait. By precisely resolving or cloning a number of these loci, we demonstrate that a broad spectrum of cellular processes contribute to hydrogen peroxide resistance, including DNA repair, scavenging of reactive oxygen species, stress-induced MAPK signaling, translation, and water transport. Consistent with the complex genetic and molecular basis of hydrogen peroxide resistance, we show two examples where multiple distinct causal genetic variants underlie what appears to be a single locus. Our results improve understanding of the genetic and molecular basis of a highly complex, model quantitative trait.

Transcriptional Derepression Uncovers Cryptic Higher-Order Genetic Interactions

Disruption of certain genes can reveal cryptic genetic variants that do not typically show phenotypic effects. Because this phenomenon, which is referred to as ‘phenotypic capacitance’, is a potential source of trait variation and disease risk, it is important to understand how it arises at the genetic and molecular levels. Here, we use a cryptic colony morphology trait that segregates in a yeast cross to explore the mechanisms underlying phenotypic capacitance. We find that the colony trait is expressed when a mutation in IRA2, a negative regulator of the Ras pathway, co-occurs with specific combinations of cryptic variants in six genes. Four of these genes encode transcription factors that act downstream of the Ras pathway, indicating that the phenotype involves genetically complex changes in the transcriptional regulation of Ras targets. We provide evidence that the IRA2 mutation reveals the phenotypic effects of the cryptic variants by disrupting the transcriptional silencing of one or more genes that contribute to the trait. Supporting this role for the IRA2 mutation, deletion of SFL1, a repressor that acts downstream of the Ras pathway, also reveals the phenotype, largely due to the same cryptic variants that were detected in the IRA2 mutant cross. Our results illustrate how higher-order genetic interactions among mutations and cryptic variants can result in phenotypic capacitance in specific genetic backgrounds, and suggests these interactions might reflect genetically complex changes in gene expression that are usually suppressed by negative regulation.

Some genetic polymorphisms have phenotypic effects that are masked under most conditions, but can be revealed by mutations or environmental change. The genetic and molecular mechanisms that suppress and uncover these cryptic genetic variants are important to understand. Here, we show that a single mutation in a yeast cross causes a major phenotypic change through its genetic interactions with two specific combinations of cryptic variants in six genes. This result suggests that in some cases cryptic variants themselves play roles in revealing their own phenotypic effects through their genetic interactions with each other and the mutations that reveal them. We also demonstrate that most of the genes harboring cryptic variation in our system are transcription factors, a finding that supports an important role for perturbation of gene regulatory networks in the uncovering of cryptic variation. As a final part of our study, we interrogate how a mutation exposes combinations of cryptic variants and obtain evidence that it does so by disrupting the silencing of one or more genes that must be expressed for the cryptic variants to exert their effects. To prove this point, we delete the transcriptional repressor that mediates this silencing and demonstrate that this deletion reveals a similar set of cryptic variants to the ones that were discovered in the initial mutant background. These findings advance our understanding of the genetic and molecular mechanisms that reveal cryptic variation.

A survey about methods dedicated to epistasis detection

During the past decade, findings of genome-wide association studies (GWAS) improved our knowledge and understanding of disease genetics. To date, thousands of SNPs have been associated with diseases and other complex traits. Statistical analysis typically looks for association between a phenotype and a SNP taken individually via single-locus tests. However, geneticists admit this is an oversimplified approach to tackle the complexity of underlying biological mechanisms. Interaction between SNPs, namely epistasis, must be considered. Unfortunately, epistasis detection gives rise to analytic challenges since analyzing every SNP combination is at present impractical at a genome-wide scale. In this review, we will present the main strategies recently proposed to detect epistatic interactions, along with their operating principle. Some of these methods are exhaustive, such as multifactor dimensionality reduction, likelihood ratio-based tests or receiver operating characteristic curve analysis; some are non-exhaustive, such as machine learning techniques (random forests, Bayesian networks) or combinatorial optimization approaches (ant colony optimization, computational evolution system).

Review: High-performance computing to detect epistasis in genome scale data sets

It is becoming clear that most human diseases have a complex etiology that cannot be explained by single nucleotide polymorphisms (SNPs) or simple additive combinations; the general consensus is that they are caused by combinations of multiple genetic variations. The limited success of some genome-wide association studies is partly a result of this focus on single genetic markers. A more promising approach is to take into account epistasis, by considering the association of multiple SNP interactions with disease. However, as genomic data continues to grow in resolution, and genome and exome sequencing become more established, the number of combinations of variants to consider increases rapidly. Two potential solutions should be considered: the use of high-performance computing, which allows us to consider a larger number of variables, and heuristics to make the solution more tractable, essential in the case of genome sequencing. In this review, we look at different computational methods to analyse epistatic interactions within disease-related genetic data sets created by microarray technology. We also review efforts to use epistatic analysis results to produce biomarkers for diagnostic tests and give our views on future directions in this field in light of advances in sequencing technology and variants in non-coding regions.