Combinatorial Genetics Reveals a Scaling Law for the Effects of Mutations on Splicing.

Despite a wealth of molecular knowledge, quantitative laws for accurate prediction of biological phenomena remain rare. Alternative pre-mRNA splicing is an important regulated step in gene expression frequently perturbed in human disease. To understand the combined effects of mutations during evolution, we quantified the effects of all possible combinations of exonic mutations accumulated during the emergence of an alternatively spliced human exon. This revealed that mutation effects scale non-monotonically with the inclusion level of an exon, with each mutation having maximum effect at a predictable intermediate inclusion level. This scaling is observed genome-wide for cis and trans perturbations of splicing, including for natural and disease-associated variants. Mathematical modeling suggests that competition between alternative splice sites is sufficient to cause this non-linearity in the genotype-phenotype map. Combining the global scaling law with specific pairwise interactions between neighboring mutations allows accurate prediction of the effects of complex genotype changes involving >10 mutations.

Mapping the Genetic Landscape of Human Cells.

Seminal yeast studies have established the value of comprehensively mapping genetic interactions (GIs) for inferring gene function. Efforts in human cells using focused gene sets underscore the utility of this approach, but the feasibility of generating large-scale, diverse human GI maps remains unresolved. We developed a CRISPR interference platform for large-scale quantitative mapping of human GIs. We systematically perturbed 222,784 gene pairs in two cancer cell lines. The resultant maps cluster functionally related genes, assigning function to poorly characterized genes, including TMEM261, a new electron transport chain component. Individual GIs pinpoint unexpected relationships between pathways, exemplified by a specific cholesterol biosynthesis intermediate whose accumulation induces deoxynucleotide depletion, causing replicative DNA damage and a synthetic-lethal interaction with the ATR/9-1-1 DNA repair pathway. Our map provides a broad resource, establishes GI maps as a high-resolution tool for dissecting gene function, and serves as a blueprint for mapping the genetic landscape of human cells.

Perturb-Seq: Dissecting Molecular Circuits with Scalable Single-Cell RNA Profiling of Pooled Genetic Screens.

Genetic screens help infer gene function in mammalian cells, but it has remained difficult to assay complex phenotypes-such as transcriptional profiles-at scale. Here, we develop Perturb-seq, combining single-cell RNA sequencing (RNA-seq) and clustered regularly interspaced short palindromic repeats (CRISPR)-based perturbations to perform many such assays in a pool. We demonstrate Perturb-seq by analyzing 200,000 cells in immune cells and cell lines, focusing on transcription factors regulating the response of dendritic cells to lipopolysaccharide (LPS). Perturb-seq accurately identifies individual gene targets, gene signatures, and cell states affected by individual perturbations and their genetic interactions. We posit new functions for regulators of differentiation, the anti-viral response, and mitochondrial function during immune activation. By decomposing many high content measurements into the effects of perturbations, their interactions, and diverse cell metadata, Perturb-seq dramatically increases the scope of pooled genomic assays.

Combinatorial CRISPR-Cas9 Metabolic Screens Reveal Critical Redox Control Points Dependent on the KEAP1-NRF2 Regulatory Axis.

The metabolic pathways fueling tumor growth have been well characterized, but the specific impact of transforming events on network topology and enzyme essentiality remains poorly understood. To this end, we performed combinatorial CRISPR-Cas9 screens on a set of 51 carbohydrate metabolism genes that represent glycolysis and the pentose phosphate pathway (PPP). This high-throughput methodology enabled systems-level interrogation of metabolic gene dispensability, interactions, and compensation across multiple cell types. The metabolic impact of specific combinatorial knockouts was validated using 13C and 2H isotope tracing, and these assays together revealed key nodes controlling redox homeostasis along the KEAP-NRF2 signaling axis. Specifically, targeting KEAP1 in combination with oxidative PPP genes mitigated the deleterious effects of these knockouts on growth rates. These results demonstrate how our integrated framework, combining genetic, transcriptomic, and flux measurements, can improve elucidation of metabolic network alterations and guide precision targeting of metabolic vulnerabilities based on tumor genetics.

AD-Syn-Net: systematic identification of Alzheimer's disease-associated mutation and co-mutation vulnerabilities via deep learning.

Alzheimer's disease (AD) is one of the most challenging neurodegenerative diseases because of its complicated and progressive mechanisms, and multiple risk factors. Increasing research evidence demonstrates that genetics may be a key factor responsible for the occurrence of the disease. Although previous reports identified quite a few AD-associated genes, they were mostly limited owing to patient sample size and selection bias. There is a lack of comprehensive research aimed to identify AD-associated risk mutations systematically. To address this challenge, we hereby construct a large-scale AD mutation and co-mutation framework ('AD-Syn-Net'), and propose deep learning models named Deep-SMCI and Deep-CMCI configured with fully connected layers that are capable of predicting cognitive impairment of subjects effectively based on genetic mutation and co-mutation profiles. Next, we apply the customized frameworks to data sets to evaluate the importance scores of the mutations and identified mutation effectors and co-mutation combination vulnerabilities contributing to cognitive impairment. Furthermore, we evaluate the influence of mutation pairs on the network architecture to dissect the genetic organization of AD and identify novel co-mutations that could be responsible for dementia, laying a solid foundation for proposing future targeted therapy for AD precision medicine. Our deep learning model codes are available open access here: https://github.com/Pan-Bio/AD-mutation-effectors.

Chemical-genetic interaction mapping links carbon metabolism and cell wall structure to tuberculosis drug efficacy.

Current chemotherapy against Mycobacterium tuberculosis (Mtb), an important human pathogen, requires a multidrug regimen lasting several months. While efforts have been made to optimize therapy by exploiting drug­drug synergies, testing new drug combinations in relevant host environments remains arduous. In particular, host environments profoundly affect the bacterial metabolic state and drug efficacy, limiting the accuracy of predictions based on in vitro assays alone. In this study, we utilized conditional Mtb knockdown mutants of essential genes as an experimentally tractable surrogate for drug treatment and probe the relationship between Mtb carbon metabolism and chemical­genetic interactions (CGIs). We examined the antitubercular drugs isoniazid, rifampicin, and moxifloxacin and found that CGIs are differentially responsive to the metabolic state, defining both environment-independent and -dependent interactions. Specifically, growth on the in vivo­relevant carbon source, cholesterol, reduced rifampicin efficacy by altering mycobacterial cell surface lipid composition. We report that a variety of perturbations in cell wall synthesis pathways restore rifampicin efficacy during growth on cholesterol, and that both environment-independent and cholesterol-dependent in vitro CGIs could be leveraged to enhance bacterial clearance in the mouse infection model. Our findings present an atlas of chemical­genetic­environmental interactions that can be used to optimize drug­drug interactions, as well as provide a framework for understanding in vitro correlates of in vivo efficacy.

Biomarkers of mitochondrial dysfunction in autism spectrum disorder: A systematic review and meta-analysis.

Autism spectrum disorder (ASD) is a neurodevelopmental disorder affecting 1 in 36 children and is associated with physiological abnormalities, most notably mitochondrial dysfunction, at least in a subset of individuals. This systematic review and meta-analysis discovered 204 relevant articles which evaluated biomarkers of mitochondrial dysfunction in ASD individuals. Significant elevations (all p < 0.01) in the prevalence of lactate (17%), pyruvate (41%), alanine (15%) and creatine kinase (9%) were found in ASD. Individuals with ASD had significant differences (all p < 0.01) with moderate to large effect sizes (Cohen's d' ≥ 0.6) compared to controls in mean pyruvate, lactate-to-pyruvate ratio, ATP, and creatine kinase. Some studies found abnormal TCA cycle metabolites associated with ASD. Thirteen controlled studies reported mitochondrial DNA (mtDNA) deletions or variations in the ASD group in blood, peripheral blood mononuclear cells, lymphocytes, leucocytes, granulocytes, and brain. Meta-analyses discovered significant differences (p < 0.01) in copy number of mtDNA overall and in ND1, ND4 and CytB genes. Four studies linked specific mtDNA haplogroups to ASD. A series of studies found a subgroup of ASD with elevated mitochondrial respiration which was associated with increased sensitivity of the mitochondria to physiological stressors and neurodevelopmental regression. Lactate, pyruvate, lactate-to-pyruvate ratio, carnitine, and acyl-carnitines were associated with clinical features such as delays in language, social interaction, cognition, motor skills, and with repetitive behaviors and gastrointestinal symptoms, although not all studies found an association. Lactate, carnitine, acyl-carnitines, ATP, CoQ10, as well as mtDNA variants, heteroplasmy, haplogroups and copy number were associated with ASD severity. Variability was found across biomarker studies primarily due to differences in collection and processing techniques as well as the intrinsic heterogeneity of the ASD population. Several studies reported alterations in mitochondrial metabolism in mothers of children with ASD and in neonates who develop ASD. Treatments targeting mitochondria, particularly carnitine and ubiquinol, appear beneficial in ASD. The link between mitochondrial dysfunction in ASD and common physiological abnormalities in individuals with ASD including gastrointestinal disorders, oxidative stress, and immune dysfunction is outlined. Several subtypes of mitochondrial dysfunction in ASD are discussed, including one related to neurodevelopmental regression, another related to alterations in microbiome metabolites, and another related to elevations in acyl-carnitines. Mechanisms linking abnormal mitochondrial function with alterations in prenatal brain development and postnatal brain function are outlined. Given the multisystem complexity of some individuals with ASD, this review presents evidence for the mitochondria being central to ASD by contributing to abnormalities in brain development, cognition, and comorbidities such as immune and gastrointestinal dysfunction as well as neurodevelopmental regression. A diagnostic approach to identify mitochondrial dysfunction in ASD is outlined. From this evidence, it is clear that many individuals with ASD have alterations in mitochondrial function which may need to be addressed in order to achieve optimal clinical outcomes. The fact that alterations in mitochondrial metabolism may be found during pregnancy and early in the life of individuals who eventually develop ASD provides promise for early life predictive biomarkers of ASD. Further studies may improve the understanding of the role of the mitochondria in ASD by better defining subgroups and understanding the molecular mechanisms driving some of the unique changes found in mitochondrial function in those with ASD.

A knowledge graph approach to predict and interpret disease-causing gene interactions.

BACKGROUND: Understanding the impact of gene interactions on disease phenotypes is increasingly recognised as a crucial aspect of genetic disease research. This trend is reflected by the growing amount of clinical research on oligogenic diseases, where disease manifestations are influenced by combinations of variants on a few specific genes. Although statistical machine-learning methods have been developed to identify relevant genetic variant or gene combinations associated with oligogenic diseases, they rely on abstract features and black-box models, posing challenges to interpretability for medical experts and impeding their ability to comprehend and validate predictions. In this work, we present a novel, interpretable predictive approach based on a knowledge graph that not only provides accurate predictions of disease-causing gene interactions but also offers explanations for these results. RESULTS: We introduce BOCK, a knowledge graph constructed to explore disease-causing genetic interactions, integrating curated information on oligogenic diseases from clinical cases with relevant biomedical networks and ontologies. Using this graph, we developed a novel predictive framework based on heterogenous paths connecting gene pairs. This method trains an interpretable decision set model that not only accurately predicts pathogenic gene interactions, but also unveils the patterns associated with these diseases. A unique aspect of our approach is its ability to offer, along with each positive prediction, explanations in the form of subgraphs, revealing the specific entities and relationships that led to each pathogenic prediction. CONCLUSION: Our method, built with interpretability in mind, leverages heterogenous path information in knowledge graphs to predict pathogenic gene interactions and generate meaningful explanations. This not only broadens our understanding of the molecular mechanisms underlying oligogenic diseases, but also presents a novel application of knowledge graphs in creating more transparent and insightful predictors for genetic research.

High Lethality of Mycobacterium tuberculosis Infection in Mice Lacking the Phagocyte Oxidase and Caspase1/11.

Immune networks that control antimicrobial and inflammatory mechanisms have overlapping regulation and functions to ensure effective host responses. Genetic interaction studies of immune pathways that compare host responses in single and combined knockout backgrounds are a useful tool to identify new mechanisms of immune control during infection. For disease caused by pulmonary Mycobacterium tuberculosis (Mtb) infections, which currently lacks an effective vaccine, understanding the genetic interactions between protective immune pathways may identify new therapeutic targets or disease-associated genes. Previous studies have suggested a direct link between the activation of NLRP3-Caspase1 inflammasome and the NADPH-dependent phagocyte oxidase complex during Mtb infection. Loss of the phagocyte oxidase complex alone resulted in increased activation of Caspase1 and IL-1ß production during Mtb infection, resulting in failed disease tolerance during the chronic stages of disease. To better understand this interaction, we generated mice lacking both Cybb, a key subunit of the phagocyte oxidase, and Caspase1/11. We found that ex vivo Mtb infection of Cybb-/-Caspase1/11-/- macrophages resulted in the expected loss of IL-1ß secretion but an unexpected change in other inflammatory cytokines and bacterial control. Mtb infected Cybb-/-Caspase1/11-/- mice rapidly progressed to severe TB, succumbing within 4 weeks to disease characterized by high bacterial burden, increased inflammatory cytokines, and the recruitment of granulocytes that associated with Mtb in the lungs. These results uncover a key genetic interaction between the phagocyte oxidase complex and Caspase1/11 that controls protection against TB and highlight the need for a better understanding of the regulation of fundamental immune networks during Mtb infection.

Genome-wide study of early and severe childhood asthma identifies interaction between CDHR3 and GSDMB.

BACKGROUND: Asthma with severe exacerbation is one of the most common causes of hospitalization among young children. Exacerbations are typically triggered by respiratory infections, but the host factors causing recurrent infections and exacerbations in some children are poorly understood. As a result, current treatment options and preventive measures are inadequate. OBJECTIVE: We sought to identify genetic interaction associated with the development of childhood asthma. METHODS: We performed an exhaustive search for pairwise interaction between genetic single nucleotide polymorphisms using 1204 cases of a specific phenotype of early childhood asthma with severe exacerbations in patients aged 2 to 6 years combined with 5328 nonasthmatic controls. Replication was attempted in 3 independent populations, and potential underlying immune mechanisms were investigated in the COPSAC2010 and COPSAC2000 birth cohorts. RESULTS: We found evidence of interaction, including replication in independent populations, between the known childhood asthma loci CDHR3 and GSDMB. The effect of CDHR3 was dependent on the GSDMB genotype, and this interaction was more pronounced for severe and early onset of disease. Blood immune analyses suggested a mechanism related to increased IL-17A production after viral stimulation. CONCLUSIONS: We found evidence of interaction between CDHR3 and GSDMB in development of early childhood asthma, possibly related to increased IL-17A response to viral infections. This study demonstrates the importance of focusing on specific disease subtypes for understanding the genetic mechanisms of asthma.

Exploring a Local Genetic Interaction Network Using Evolutionary Replay Experiments.

Understanding how genes interact is a central challenge in biology. Experimental evolution provides a useful, but underutilized, tool for identifying genetic interactions, particularly those that involve non-loss-of-function mutations or mutations in essential genes. We previously identified a strong positive genetic interaction between specific mutations in KEL1 (P344T) and HSL7 (A695fs) that arose in an experimentally evolved Saccharomyces cerevisiae population. Because this genetic interaction is not phenocopied by gene deletion, it was previously unknown. Using "evolutionary replay" experiments, we identified additional mutations that have positive genetic interactions with the kel1-P344T mutation. We replayed the evolution of this population 672 times from six timepoints. We identified 30 populations where the kel1-P344T mutation reached high frequency. We performed whole-genome sequencing on these populations to identify genes in which mutations arose specifically in the kel1-P344T background. We reconstructed mutations in the ancestral and kel1-P344T backgrounds to validate positive genetic interactions. We identify several genetic interactors with KEL1, we validate these interactions by reconstruction experiments, and we show these interactions are not recapitulated by loss-of-function mutations. Our results demonstrate the power of experimental evolution to identify genetic interactions that are positive, allele specific, and not readily detected by other methods, shedding light on an underexplored region of the yeast genetic interaction network.

The Causes and Consequences of Genetic Interactions (Epistasis).

The same mutation can have different effects in different individuals. One important reason for this is that the outcome of a mutation can depend on the genetic context in which it occurs. This dependency is known as epistasis. In recent years, there has been a concerted effort to quantify the extent of pairwise and higher-order genetic interactions between mutations through deep mutagenesis of proteins and RNAs. This research has revealed two major components of epistasis: nonspecific genetic interactions caused by nonlinearities in genotype-to-phenotype maps, and specific interactions between particular mutations. Here, we provide an overview of our current understanding of the mechanisms causing epistasis at the molecular level, the consequences of genetic interactions for evolution and genetic prediction, and the applications of epistasis for understanding biology and determining macromolecular structures.

Synthetic negative genome screen of the GPN-loop GTPase NPA3 in Saccharomyces cerevisiae.

The GPN-loop GTPase Npa3 is encoded by an essential gene in the yeast Saccharomyces cerevisiae. Npa3 plays a critical role in the assembly and nuclear accumulation of RNA polymerase II (RNAPII), a function that may explain its essentiality. Genetic interactions describe the extent to which a mutation in a particular gene affects a specific phenotype when co-occurring with an alteration in a second gene. Discovering synthetic negative genetic interactions has long been used as a tool to delineate the functional relatedness between pairs of genes participating in common or compensatory biological pathways. Previously, our group showed that nuclear targeting and transcriptional activity of RNAPII were unaffected in cells expressing exclusively a C-terminal truncated mutant version of Npa3 (npa3∆C) lacking the last 106 residues naturally absent from the single GPN protein in Archaea, but universally conserved in all Npa3 orthologs of eukaryotes. To gain insight into novel cellular functions for Npa3, we performed here a genome-wide Synthetic Genetic Array (SGA) study coupled to bulk fluorescence monitoring to identify negative genetic interactions of NPA3 by crossing an npa3∆C strain with a 4,389 nonessential gene-deletion collection. This genetic screen revealed previously unknown synthetic negative interactions between NPA3 and 15 genes. Our results revealed that the Npa3 C-terminal tail extension regulates the participation of this essential GTPase in previously unknown biological processes related to mitochondrial homeostasis and ribosome biogenesis.

Natural variants suppress mutations in hundreds of essential genes.

The consequence of a mutation can be influenced by the context in which it operates. For example, loss of gene function may be tolerated in one genetic background, and lethal in another. The extent to which mutant phenotypes are malleable, the architecture of modifiers and the identities of causal genes remain largely unknown. Here, we measure the fitness effects of ~ 1,100 temperature-sensitive alleles of yeast essential genes in the context of variation from ten different natural genetic backgrounds and map the modifiers for 19 combinations. Altogether, fitness defects for 149 of the 580 tested genes (26%) could be suppressed by genetic variation in at least one yeast strain. Suppression was generally driven by gain-of-function of a single, strong modifier gene, and involved both genes encoding complex or pathway partners suppressing specific temperature-sensitive alleles, as well as general modifiers altering the effect of many alleles. The emerging frequency of suppression and range of possible mechanisms suggest that a substantial fraction of monogenic diseases could be managed by modulating other gene products.

The effects of epistasis and linkage on the invasion of locally beneficial mutations and the evolution of genomic islands.

We study local adaptation of a peripheral population by investigating the fate of new mutations using a haploid two-locus two-allele continent-island migration model. We explore how linkage, epistasis, and maladaptive gene flow affect the invasion probability of weakly beneficial de-novo mutations that arise on the island at an arbitrary physical distance to a locus that already maintains a stable migration-selection polymorphism. By assuming a slightly supercritical branching process, we deduce explicit conditions on the parameters that permit a positive invasion probability and we derive approximations for it. They show how the invasion probability depends on the additive and epistatic effects of the mutant, on its linkage to the polymorphism, and on the migration rate. We use these approximations together with empirically motivated distributions of epistatic effects to analyze the influence of epistasis on the expected invasion probability if mutants are drawn randomly from such a distribution and occur at a random physical distance to the existing polymorphism. We find that the invasion probability generally increases as the epistasis parameter increases or the migration rate decreases, but not necessarily as the recombination rate decreases. Finally, we shed light on the size of emerging genomic islands of divergence by exploring the size of the chromosomal neighborhood of the already established polymorphism in which 50% or 90% of the successfully invading mutations become established. These 'window sizes' always decrease in a reverse sigmoidal way with stronger migration and typically increase with increasing epistatic effect.

Systematic analysis of bypass suppression of essential genes.

Essential genes tend to be highly conserved across eukaryotes, but, in some cases, their critical roles can be bypassed through genetic rewiring. From a systematic analysis of 728 different essential yeast genes, we discovered that 124 (17%) were dispensable essential genes. Through whole-genome sequencing and detailed genetic analysis, we investigated the genetic interactions and genome alterations underlying bypass suppression. Dispensable essential genes often had paralogs, were enriched for genes encoding membrane-associated proteins, and were depleted for members of protein complexes. Functionally related genes frequently drove the bypass suppression interactions. These gene properties were predictive of essential gene dispensability and of specific suppressors among hundreds of genes on aneuploid chromosomes. Our findings identify yeast's core essential gene set and reveal that the properties of dispensable essential genes are conserved from yeast to human cells, correlating with human genes that display cell line-specific essentiality in the Cancer Dependency Map (DepMap) project.

Shaping a fruit: Developmental pathways that impact growth patterns.

Angiosperms produce seeds as their progeny enclosed in maternally-derived structures called fruits. Evolutionarily, fruits have contributed enormously to the success of the Angiosperms phylum by providing protection and nutrition to the developing seeds, while ensuring the efficient dispersal upon maturity. Fruits vary massively in both size and shape and certain species have been targeted for domestication due to their nutritional value and delicious taste. Among the vast array of 3D fruit shapes that exist in nature, the mechanism by which growth is oriented and coordinated to generate this diversity of forms is unclear. In this review, we discuss the latest results in identifying components that control fruit morphology and their effect on isotropic and anisotropic growth. Moreover, we will compare the current knowledge on the mechanisms that control fruit growth, size and shape between the domesticated Solanaceae species, tomato and members of the large family of Brassicaceae.

Multiple functionalities of molecular chaperones revealed through systematic mapping of their interaction networks.

Chaperones are a highly interactive group of proteins that function globally in many cellular processes involved in maintaining protein homeostasis. Traditional biochemical assays typically do not provide a complete view of the intricate networks through which chaperones collaborate to promote proteostasis. Recent advances in high-throughput systematic analyses of chaperone interactions have uncovered that chaperones display a remarkable cooperativity in their interactions with numerous client proteins. This cooperativity has been found to be a fundamental aspect of a properly functioning cell. Aberrant formation or improper regulation of these interactions can easily lead to disease states. Herein, we provide an overview of the use of large-scale interaction assays, whether physical (protein-protein) or genetic (epistatic), to study chaperone interaction networks. Importantly, we discuss the ongoing need for such studies to determine the mechanisms by which protein homeostasis is controlled in the cell.

Non-additive (dominance) effects of genetic variants associated with refractive error and myopia.

Genome-wide association studies (GWAS) have revealed that the genetic contribution to certain complex diseases is well-described by Fisher's infinitesimal model in which a vast number of polymorphisms each confer a small effect. Under Fisher's model, variants have additive effects both across loci and within loci. However, the latter assumption is at odds with the common observation of dominant or recessive rare alleles responsible for monogenic disorders. Here, we searched for evidence of non-additive (dominant or recessive) effects for GWAS variants known to confer susceptibility to the highly heritable quantitative trait, refractive error. Of 146 GWAS variants examined in a discovery sample of 228,423 individuals whose refractive error phenotype was inferred from their age-of-onset of spectacle wear, only 8 had even nominal evidence (p < 0.05) of non-additive effects. In a replication sample of 73,577 individuals who underwent direct assessment of refractive error, 1 of these 8 variants had robust independent evidence of non-additive effects (rs7829127 within ZMAT4, p = 4.76E-05) while a further 2 had suggestive evidence (rs35337422 in RD3L, p = 7.21E-03 and rs12193446 in LAMA2, p = 2.57E-02). Accounting for non-additive effects had minimal impact on the accuracy of a polygenic risk score for refractive error (R2 = 6.04% vs. 6.01%). Our findings demonstrate that very few GWAS variants for refractive error show evidence of a departure from an additive mode of action and that accounting for non-additive risk variants offers little scope to improve the accuracy of polygenic risk scores for myopia.

Transferring entropy to the realm of GxG interactions.

Genome-wide association studies are moving to genome-wide interaction studies, as the genetic background of many diseases appears to be more complex than previously supposed. Thus, many statistical approaches have been proposed to detect gene-gene (GxG) interactions, among them numerous information theory-based methods, inspired by the concept of entropy. These are suggested as particularly powerful and, because of their nonlinearity, as better able to capture nonlinear relationships between genetic variants and/or variables. However, the introduced entropy-based estimators differ to a surprising extent in their construction and even with respect to the basic definition of interactions. Also, not every entropy-based measure for interaction is accompanied by a proper statistical test. To shed light on this, a systematic review of the literature is presented answering the following questions: (1) How are GxG interactions defined within the framework of information theory? (2) Which entropy-based test statistics are available? (3) Which underlying distribution do the test statistics follow? (4) What are the given strengths and limitations of these test statistics?