Traversing the conceptual divide between biological and statistical epistasis: systems biology and a more modern synthesis

The synergism between metabolic and target-site resistance enhances the intensity of resistance to pyrethroids in Spodoptera exigua

The widespread application of insecticides imposes intense selective pressure on pest populations, driving the evolution of high-level resistance and leading to frequent control failures of pest. Insecticide resistance is primarily mediated through two primary mechanisms: target-site insensitivity and enhanced metabolic detoxification. However, the potential interactions and synergistic effects between these mechanisms remain largely unexplored. In this study, we demonstrate a striking cooperative interaction between these two major resistance mechanisms in a field-derived strain of Spodoptera exigua exhibiting extreme resistance (631-fold) to the pyrethroid insecticide lambda-cyhalothrin. Through genetic mapping and linkage analysis, we identified that this resistance phenotype is conferred by the combined effects of overexpression of the P450 CYP9A9 (two copies: CYP9A9a and CYP9A9b) and a target-site mutation (L1014F, kdr) in the voltage-gated sodium channel. Using an introgression approach, we generated two near-isogenic strains: WH-kdr, carrying only the target-site resistance allele (6.2-fold resistance), and WH-CYP9A, harboring only the metabolic resistance genes (79-fold resistance), both compared to the susceptible WH-S strain. CRISPR/Cas9-mediated knockout of both CYP9A9 copies in the QP19 strain dramatically reduced resistance from 631-fold to 19-fold, while transgenic expression of the CYP9A9a variant (containing three amino acid substitutions) from QP19 strain in Helicoverpa armigera conferred 39-fold resistance to lambda-cyhalothrin. These findings provide compelling evidence that target-site resistance can significantly potentiate metabolic resistance, resulting in substantially higher resistance levels than either mechanism alone in S. exigua. These findings enhance the understanding of higher level resistance mechanisms mediated by interactions between resistance genes and provide theoretical basis for devising management strategies of insecticide resistance.

Leveraging hierarchical structures for genetic block interaction studies using the hierarchical transformer

Initially introduced in 1909 by William Bateson, classic epistasis (genetic variant interaction) refers to the phenomenon that one variant prevents another variant from a different locus from manifesting its effects. The potential effects of genetic variant interactions on complex diseases have been recognized for the past decades. Moreover, It has been studied and demonstrated that leveraging the combined SNP effects within the genetic block can significantly increase calculation power, reducing background noise, ultimately leading to novel epistasis discovery that the single SNP statistical epistasis study might overlook. However, it is still an open question how we can best combine gene structure representation modelling and interaction learning into an end-to-end model for gene interaction searching. Here, in the current study, we developed a neural genetic block interaction searching model that can effectively process large SNP chip inputs and output the potential genetic block interaction heatmap. Our model augments a previously published hierarchical transformer architecture (Liu and Lapata, 2019) with the ability to model genetic blocks. The cross-block relationship mapping was achieved via a hierarchical attention mechanism which allows the sharing of information regarding specific phenotypes, as opposed to simple unsupervised dimensionality reduction methods e.g. PCA. Results on both simulation and UK Biobank studies show our model brings substantial improvements compared to traditional exhaustive searching and neural network methods.

Distinct network patterns emerge from Cartesian and XOR epistasis models: a comparative network science analysis

BACKGROUND

Epistasis, the phenomenon where the effect of one gene (or variant) is masked or modified by one or more other genes, significantly contributes to the phenotypic variance of complex traits. Traditionally, epistasis has been modeled using the Cartesian epistatic model, a multiplicative approach based on standard statistical regression. However, a recent study investigating epistasis in obesity-related traits has identified potential limitations of the Cartesian epistatic model, revealing that it likely only detects a fraction of the genetic interactions occurring in natural systems. In contrast, the exclusive-or (XOR) epistatic model has shown promise in detecting a broader range of epistatic interactions and revealing more biologically relevant functions associated with interacting variants. To investigate whether the XOR epistatic model also forms distinct network structures compared to the Cartesian model, we applied network science to examine genetic interactions underlying body mass index (BMI) in rats (Rattus norvegicus).

RESULTS

Our comparative analysis of XOR and Cartesian epistatic models in rats reveals distinct topological characteristics. The XOR model exhibits enhanced sensitivity to epistatic interactions between the network communities found in the Cartesian epistatic network, facilitating the identification of novel trait-related biological functions via community-based enrichment analysis. Additionally, the XOR network features triangle network motifs, indicative of higher-order epistatic interactions. This research also evaluates the impact of linkage disequilibrium (LD)-based edge pruning on network-based epistasis analysis, finding that LD-based edge pruning may lead to increased network fragmentation, which may hinder the effectiveness of network analysis for the investigation of epistasis. We confirmed through network permutation analysis that most XOR and Cartesian epistatic networks derived from the data display distinct structural properties compared to randomly shuffled networks.

CONCLUSIONS

Collectively, these findings highlight the XOR model's ability to uncover meaningful biological associations and higher-order epistasis derived from lower-order network topologies. The introduction of community-based enrichment analysis and motif-based epistatic discovery emphasize network science as a critical approach for advancing epistasis research and understanding complex genetic architectures.

Extreme positive epistasis for fitness in monosomic yeast strains

The loss of a single chromosome in a diploid organism halves the dosage of many genes and is usually accompanied by a substantial decrease in fitness. We asked whether this decrease simply reflects the joint damage caused by individual gene dosage deficiencies. We measured the fitness effects of single heterozygous gene deletions in yeast and combined them for each chromosome. This predicted a negative growth rate, that is, lethality, for multiple monosomies. However, monosomic strains remained alive and grew as if much (often most) of the damage caused by single mutations had disappeared, revealing an exceptionally large and positive epistatic component of fitness. We looked for functional explanations by analyzing the transcriptomes. There was no evidence of increased (compensatory) gene expression on the monosomic chromosomes. Nor were there signs of the cellular stress response that would be expected if monosomy led to protein destabilization and thus cytotoxicity. Instead, all monosomic strains showed extensive upregulation of genes encoding ribosomal proteins, but in an indiscriminate manner that did not correspond to their altered dosage. This response did not restore the stoichiometry required for efficient biosynthesis, which probably became growth limiting, making all other mutation-induced metabolic defects much less important. In general, the modular structure of the cell leads to an effective fragmentation of the total mutational load. Defects outside the module(s) currently defining fitness lose at least some of their relevance, producing the epiphenomenon of positive interactions between individually negative effects.

Assessing the limitations of relief-based algorithms in detecting higher-order interactions

BACKGROUND

Epistasis, the interaction between genetic loci where the effect of one locus is influenced by one or more other loci, plays a crucial role in the genetic architecture of complex traits. However, as the number of loci considered increases, the investigation of epistasis becomes exponentially more complex, making the selection of key features vital for effective downstream analyses. Relief-Based Algorithms (RBAs) are often employed for this purpose due to their reputation as "interaction-sensitive" algorithms and uniquely non-exhaustive approach. However, the limitations of RBAs in detecting interactions, particularly those involving multiple loci, have not been thoroughly defined. This study seeks to address this gap by evaluating the efficiency of RBAs in detecting higher-order epistatic interactions. Motivated by previous findings that suggest some RBAs may rank predictive features involved in higher-order epistasis negatively, we explore the potential of absolute value ranking of RBA feature weights as an alternative approach for capturing complex interactions. In this study, we assess the performance of ReliefF, MultiSURF, and MultiSURFstar on simulated genetic datasets that model various patterns of genotype-phenotype associations, including 2-way to 5-way genetic interactions, and compare their performance to two control methods: a random shuffle and mutual information.

RESULTS

Our findings indicate that while RBAs effectively identify lower-order (2 to 3-way) interactions, their capability to detect higher-order interactions is significantly limited, primarily by large feature count but also by signal noise. Specifically, we observe that RBAs are successful in detecting fully penetrant 4-way XOR interactions using an absolute value ranking approach, but this is restricted to datasets with only 20 total features.

CONCLUSIONS

These results highlight the inherent limitations of current RBAs and underscore the need for the development of Relief-based approaches with enhanced detection capabilities for the investigation of epistasis, particularly in datasets with large feature counts and complex higher-order interactions.

Synergistic Effect between the APOE ε4 Allele with Genetic Variants of GSK3B and MAPT: Differential Profile between Refractory Epilepsy and Alzheimer Disease

Temporal Lobe Epilepsy (TLE) is a chronic neurological disorder characterized by recurrent focal seizures originating in the temporal lobe. Despite the variety of antiseizure drugs currently available to treat TLE, about 30% of cases continue to have seizures. The etiology of TLE is complex and multifactorial. Increasing evidence indicates that Alzheimer's disease (AD) and drug-resistant TLE present common pathological features that may induce hyperexcitability, especially aberrant hyperphosphorylation of tau protein. Genetic polymorphic variants located in genes of the microtubule-associated protein tau (MAPT) and glycogen synthase kinase-3β (GSK3B) have been associated with the risk of developing AD. The APOE ε4 allele is a major genetic risk factor for AD. Likewise, a gene-dose-dependent effect of ε4 seems to influence TLE. The present study aimed to investigate whether the APOE ɛ4 allele and genetic variants located in the MAPT and GSK3B genes are associated with the risk of developing AD and drug-resistant TLE in a cohort of the Mexican population. A significant association with the APOE ε4 allele was observed in patients with AD and TLE. Additional genetic interactions were identified between this allele and variants of the MAPT and GSK3B genes.

Identification of interactions between genetic risk scores and dietary patterns for personalized prevention of kidney dysfunction in a population-based cohort

BACKGROUND & AIM

Chronic kidney disease (CKD) is a heterogeneous disorder that affects the kidney structure and function. This study investigated the effect of the interaction between genetic factors and dietary pattern on kidney dysfunction in Korean adults.

METHODS

Baseline data were obtained from the Ansan and Ansung Study of the Korean Genome and Epidemiology Study involving 8230 participants aged 40-69 years. Kidney dysfunction was defined as an estimated glomerular filtration rate < 90 mL/minute/1.73 m2. Genomic DNAs genotyped on the Affymetrix® Genome-Wide Human SNP array 5.0 were isolated from peripheral blood. A genome-wide association study using a generalized linear model was performed on 1,590,162 single-nucleotide polymorphisms (SNPs). To select significant SNPs, the threshold criterion was set at P-value < 5 × 10-8. Linkage disequilibrium clumping was performed based on the R2 value, and 94 SNPs had a significant effect. Participants were divided into two groups based on their generic risk score (GRS): the low-GR group had GRS > 0, while the high-GR group had GRS ≤ 0.

RESULTS

Three distinct dietary patterns were extracted, namely, the "prudent pattern," "flour-based and animal food pattern," and "white rice pattern," to analyze the effect of dietary pattern on kidney function. In the "flour-based and animal food pattern," higher pattern scores were associated with a higher prevalence of kidney dysfunction in both the low and high GR groups (P for trend < 0.0001 in the low-, high-GR groups of model 1; 0.0050 and 0.0065 in the low-, high-GR groups of model 2, respectively).

CONCLUSIONS

The results highlight a significant association between the 'flour-based and animal food pattern' and higher kidney dysfunction prevalence in individuals with both low and high GR. These findings suggest that personalized nutritional interventions based on GR profiles may become the basis for presenting GR-based individual dietary patterns for kidney dysfunction.

Interaction models matter: an efficient, flexible computational framework for model-specific investigation of epistasis

PURPOSE

Epistasis, the interaction between two or more genes, is integral to the study of genetics and is present throughout nature. Yet, it is seldom fully explored as most approaches primarily focus on single-locus effects, partly because analyzing all pairwise and higher-order interactions requires significant computational resources. Furthermore, existing methods for epistasis detection only consider a Cartesian (multiplicative) model for interaction terms. This is likely limiting as epistatic interactions can evolve to produce varied relationships between genetic loci, some complex and not linearly separable.

METHODS

We present new algorithms for the interaction coefficients for standard regression models for epistasis that permit many varied models for the interaction terms for loci and efficient memory usage. The algorithms are given for two-way and three-way epistasis and may be generalized to higher order epistasis. Statistical tests for the interaction coefficients are also provided. We also present an efficient matrix based algorithm for permutation testing for two-way epistasis. We offer a proof and experimental evidence that methods that look for epistasis only at loci that have main effects may not be justified. Given the computational efficiency of the algorithm, we applied the method to a rat data set and mouse data set, with at least 10,000 loci and 1,000 samples each, using the standard Cartesian model and the XOR model to explore body mass index.

RESULTS

This study reveals that although many of the loci found to exhibit significant statistical epistasis overlap between models in rats, the pairs are mostly distinct. Further, the XOR model found greater evidence for statistical epistasis in many more pairs of loci in both data sets with almost all significant epistasis in mice identified using XOR. In the rat data set, loci involved in epistasis under the XOR model are enriched for biologically relevant pathways.

CONCLUSION

Our results in both species show that many biologically relevant epistatic relationships would have been undetected if only one interaction model was applied, providing evidence that varied interaction models should be implemented to explore epistatic interactions that occur in living systems.

Uncovering Evolutionary Adaptations in Common Warthogs through Genomic Analyses

In the Suidae family, warthogs show significant survival adaptability and trait specificity. This study offers a comparative genomic analysis between the warthog and other Suidae species, including the Luchuan pig, Duroc pig, and Red River hog. By integrating the four genomes with sequences from the other four species, we identified 8868 single-copy orthologous genes. Based on 8868 orthologous protein sequences, phylogenetic assessments highlighted divergence timelines and unique evolutionary branches within suid species. Warthogs exist on different evolutionary branches compared to DRCs and LCs, with a divergence time preceding that of DRC and LC. Contraction and expansion analyses of warthog gene families have been conducted to elucidate the mechanisms of their evolutionary adaptations. Using GO, KEGG, and MGI databases, warthogs showed a preference for expansion in sensory genes and contraction in metabolic genes, underscoring phenotypic diversity and adaptive evolution direction. Associating genes with the QTLdb-pigSS11 database revealed links between gene families and immunity traits. The overlap of olfactory genes in immune-related QTL regions highlighted their importance in evolutionary adaptations. This work highlights the unique evolutionary strategies and adaptive mechanisms of warthogs, guiding future research into the distinct adaptability and disease resistance in pigs, particularly focusing on traits such as resistance to African Swine Fever Virus.

Gene Interactions in Human Disease Studies-Evidence Is Mounting

Despite monumental advances in molecular technology to generate genome sequence data at scale, there is still a considerable proportion of heritability in most complex diseases that remains unexplained. Because many of the discoveries have been single-nucleotide variants with small to moderate effects on disease, the functional implication of many of the variants is still unknown and, thus, we have limited new drug targets and therapeutics. We, and many others, posit that one primary factor that has limited our ability to identify novel drug targets from genome-wide association studies may be due to gene interactions (epistasis), gene-environment interactions, network/pathway effects, or multiomic relationships. We propose that many of these complex models explain much of the underlying genetic architecture of complex disease. In this review, we discuss the evidence from multiple research avenues, ranging from pairs of alleles to multiomic integration studies and pharmacogenomics, that supports the need for further investigation of gene interactions (or epistasis) in genetic and genomic studies of human disease. Our goal is to catalog the mounting evidence for epistasis in genetic studies and the connections between genetic interactions and human health and disease that could enable precision medicine of the future.

Edge and modular significance assessment in individual-specific networks

Individual-specific networks, defined as networks of nodes and connecting edges that are specific to an individual, are promising tools for precision medicine. When such networks are biological, interpretation of functional modules at an individual level becomes possible. An under-investigated problem is relevance or "significance" assessment of each individual-specific network. This paper proposes novel edge and module significance assessment procedures for weighted and unweighted individual-specific networks. Specifically, we propose a modular Cook's distance using a method that involves iterative modeling of one edge versus all the others within a module. Two procedures assessing changes between using all individuals and using all individuals but leaving one individual out (LOO) are proposed as well (LOO-ISN, MultiLOO-ISN), relying on empirically derived edges. We compare our proposals to competitors, including adaptions of OPTICS, kNN, and Spoutlier methods, by an extensive simulation study, templated on real-life scenarios for gene co-expression and microbial interaction networks. Results show the advantages of performing modular versus edge-wise significance assessments for individual-specific networks. Furthermore, modular Cook's distance is among the top performers across all considered simulation settings. Finally, the identification of outlying individuals regarding their individual-specific networks, is meaningful for precision medicine purposes, as confirmed by network analysis of microbiome abundance profiles.

Missing Causality and Heritability of Autoimmune Hepatitis

BACKGROUND

Autoimmune hepatitis has an unknown cause and genetic associations that are not disease-specific or always present. Clarification of its missing causality and heritability could improve prevention and management strategies.

AIMS

Describe the key epigenetic and genetic mechanisms that could account for missing causality and heritability in autoimmune hepatitis; indicate the prospects of these mechanisms as pivotal factors; and encourage investigations of their pathogenic role and therapeutic potential.

METHODS

English abstracts were identified in PubMed using multiple key search phases. Several hundred abstracts and 210 full-length articles were reviewed.

RESULTS

Environmental induction of epigenetic changes is the prime candidate for explaining the missing causality of autoimmune hepatitis. Environmental factors (diet, toxic exposures) can alter chromatin structure and the production of micro-ribonucleic acids that affect gene expression. Epistatic interaction between unsuspected genes is the prime candidate for explaining the missing heritability. The non-additive, interactive effects of multiple genes could enhance their impact on the propensity and phenotype of autoimmune hepatitis. Transgenerational inheritance of acquired epigenetic marks constitutes another mechanism of transmitting parental adaptations that could affect susceptibility. Management strategies could range from lifestyle adjustments and nutritional supplements to precision editing of the epigenetic landscape.

CONCLUSIONS

Autoimmune hepatitis has a missing causality that might be explained by epigenetic changes induced by environmental factors and a missing heritability that might reflect epistatic gene interactions or transgenerational transmission of acquired epigenetic marks. These unassessed or under-evaluated areas warrant investigation.

Learning high-order interactions for polygenic risk prediction

Within the framework of precision medicine, the stratification of individual genetic susceptibility based on inherited DNA variation has paramount relevance. However, one of the most relevant pitfalls of traditional Polygenic Risk Scores (PRS) approaches is their inability to model complex high-order non-linear SNP-SNP interactions and their effect on the phenotype (e.g. epistasis). Indeed, they incur in a computational challenge as the number of possible interactions grows exponentially with the number of SNPs considered, affecting the statistical reliability of the model parameters as well. In this work, we address this issue by proposing a novel PRS approach, called High-order Interactions-aware Polygenic Risk Score (hiPRS), that incorporates high-order interactions in modeling polygenic risk. The latter combines an interaction search routine based on frequent itemsets mining and a novel interaction selection algorithm based on Mutual Information, to construct a simple and interpretable weighted model of user-specified dimensionality that can predict a given binary phenotype. Compared to traditional PRSs methods, hiPRS does not rely on GWAS summary statistics nor any external information. Moreover, hiPRS differs from Machine Learning-based approaches that can include complex interactions in that it provides a readable and interpretable model and it is able to control overfitting, even on small samples. In the present work we demonstrate through a comprehensive simulation study the superior performance of hiPRS w.r.t. state of the art methods, both in terms of scoring performance and interpretability of the resulting model. We also test hiPRS against small sample size, class imbalance and the presence of noise, showcasing its robustness to extreme experimental settings. Finally, we apply hiPRS to a case study on real data from DACHS cohort, defining an interaction-aware scoring model to predict mortality of stage II-III Colon-Rectal Cancer patients treated with oxaliplatin.

Epi-MEIF: detecting higher order epistatic interactions for complex traits using mixed effect conditional inference forests

Understanding the relationship between genetic variations and variations in complex and quantitative phenotypes remains an ongoing challenge. While Genome-wide association studies (GWAS) have become a vital tool for identifying single-locus associations, we lack methods for identifying epistatic interactions. In this article, we propose a novel method for higher-order epistasis detection using mixed effect conditional inference forest (epiMEIF). The proposed method is fitted on a group of single nucleotide polymorphisms (SNPs) potentially associated with the phenotype and the tree structure in the forest facilitates the identification of n-way interactions between the SNPs. Additional testing strategies further improve the robustness of the method. We demonstrate its ability to detect true n-way interactions via extensive simulations in both cross-sectional and longitudinal synthetic datasets. This is further illustrated in an application to reveal epistatic interactions from natural variations of cardiac traits in flies (Drosophila). Overall, the method provides a generalized way to identify higher-order interactions from any GWAS data, thereby greatly improving the detection of the genetic architecture underlying complex phenotypes.

Open problems in human trait genetics

Genetic studies of human traits have revolutionized our understanding of the variation between individuals, and yet, the genetics of most traits is still poorly understood. In this review, we highlight the major open problems that need to be solved, and by discussing these challenges provide a primer to the field. We cover general issues such as population structure, epistasis and gene-environment interactions, data-related issues such as ancestry diversity and rare genetic variants, and specific challenges related to heritability estimates, genetic association studies, and polygenic risk scores. We emphasize the interconnectedness of these problems and suggest promising avenues to address them.

Concurrent outcomes from multiple approaches of epistasis analysis for human body mass index associated loci provide insights into obesity biology

Genome wide association studies (GWAS) have focused on elucidating the genetic architecture of complex traits by assessing single variant effects in additive genetic models, albeit explaining a fraction of the trait heritability. Epistasis has recently emerged as one of the intrinsic mechanisms that could explain part of this missing heritability. We conducted epistasis analysis for genome-wide body mass index (BMI) associated SNPs in Alzheimer's Disease Neuroimaging Initiative (ADNI) and followed up top significant interacting SNPs for replication in the UK Biobank imputed genotype dataset. We report two pairwise epistatic interactions, between rs2177596 (RHBDD1) and rs17759796 (MAPK1), rs1121980 (FTO) and rs6567160 (MC4R), obtained from a consensus of nine different epistatic approaches. Gene interaction maps and tissue expression profiles constructed for these interacting loci highlights co-expression, co-localisation, physical interaction, genetic interaction, and shared pathways emphasising the neuronal influence in obesity and implicating concerted expression of associated genes in liver, pancreas, and adipose tissues insinuating to metabolic abnormalities characterized by obesity. Detecting epistasis could thus be a promising approach to understand the effect of simultaneously interacting multiple genetic loci in disease aetiology, beyond single locus effects.

A systematic analysis of gene-gene interaction in multiple sclerosis

BACKGROUND

For the most part, genome-wide association studies (GWAS) have only partially explained the heritability of complex diseases. One of their limitations is to assume independent contributions of individual variants to the phenotype. Many tools have therefore been developed to investigate the interactions between distant loci, or epistasis. Among them, the recently proposed EpiGWAS models the interactions between a target variant and the rest of the genome. However, applying this approach to studying interactions along all genes of a disease map is not straightforward. Here, we propose a pipeline to that effect, which we illustrate by investigating a multiple sclerosis GWAS dataset from the Wellcome Trust Case Control Consortium 2 through 19 disease maps from the MetaCore pathway database.

RESULTS

For each disease map, we build an epistatic network by connecting the genes that are deemed to interact. These networks tend to be connected, complementary to the disease maps and contain hubs. In addition, we report 4 epistatic gene pairs involving missense variants, and 25 gene pairs with a deleterious epistatic effect mediated by eQTLs. Among these, we highlight the interaction of GLI-1 and SUFU, and of IP10 and NF-[Formula: see text]B, as they both match known biological interactions. The latter pair is particularly promising for therapeutic development, as both genes have known inhibitors.

CONCLUSIONS

Our study showcases the ability of EpiGWAS to uncover biologically interpretable epistatic interactions that are potentially actionable for the development of combination therapy.

The interactions between interleukin-1 family genes: IL1A, IL1B, IL1RN, and obesity parameters

Background

Obesity has been recognized as a worldwide growing problem, producing many pathologies including the promotion of “proinflammatory state.” The etiology of human obesity is still only partially understood; however, the genetic background has been proved. Its nature is complex, and currently, it appears that the combined effects of the interactions among multiple genes should receive more attention. Due to the fact that obesity promotes proinflammatory conditions, in this study, we investigated the genetic polymorphism of IL-1 family genes in healthy people with normal and elevated body mass index (BMI) and fat %.

Results

The single-nucleotide polymorphisms (SNPs) within the IL1A -889C > T (rs1800587), IL1B + 3954 T > C (rs1143634), and IL1RN -87G > A (rs2234677) genes alone were associated neither with BMI nor fat % values in tested group. The associations between SNP–SNP interaction and BMI for the IL1B × IL1RN interactions were significant for dominant model (p = 0.02) and codominant model (p = 0.03). The same SNP-SNP interaction (IL1B × IL1RN) was associated also with fat % for codominant (p = 0.01) and recessive (p = 0.002) models.

Conclusions

This study further confirmed that IL-1 family genes are involved in genetic background of obesity. It has been shown that interaction IL1B × IL1RN was associated with both BMI and fat % with rare T allele protecting form higher values. Thus, even if certain polymorphisms in single genes of IL-1 family cannot be defined as related to obesity in examined population, the genetic interrelationships should be analyzed.

Interpretable network-guided epistasis detection

BACKGROUND

Detecting epistatic interactions at the gene level is essential to understanding the biological mechanisms of complex diseases. Unfortunately, genome-wide interaction association studies involve many statistical challenges that make such detection hard. We propose a multi-step protocol for epistasis detection along the edges of a gene-gene co-function network. Such an approach reduces the number of tests performed and provides interpretable interactions while keeping type I error controlled. Yet, mapping gene interactions into testable single-nucleotide polymorphism (SNP)-interaction hypotheses, as well as computing gene pair association scores from SNP pair ones, is not trivial.

RESULTS

Here we compare 3 SNP-gene mappings (positional overlap, expression quantitative trait loci, and proximity in 3D structure) and use the adaptive truncated product method to compute gene pair scores. This method is non-parametric, does not require a known null distribution, and is fast to compute. We apply multiple variants of this protocol to a genome-wide association study dataset on inflammatory bowel disease. Different configurations produced different results, highlighting that various mechanisms are implicated in inflammatory bowel disease, while at the same time, results overlapped with known disease characteristics. Importantly, the proposed pipeline also differs from a conventional approach where no network is used, showing the potential for additional discoveries when prior biological knowledge is incorporated into epistasis detection.

Detecting gene-gene interactions from GWAS using diffusion kernel principal components

Genes and gene products do not function in isolation but as components of complex networks of macromolecules through physical or biochemical interactions. Dependencies of gene mutations on genetic background (i.e., epistasis) are believed to play a role in understanding molecular underpinnings of complex diseases such as inflammatory bowel disease (IBD). However, the process of identifying such interactions is complex due to for instance the curse of high dimensionality, dependencies in the data and non-linearity. Here, we propose a novel approach for robust and computationally efficient epistasis detection. We do so by first reducing dimensionality, per gene via diffusion kernel principal components (kpc). Subsequently, kpc gene summaries are used for downstream analysis including the construction of a gene-based epistasis network. We show that our approach is not only able to recover known IBD associated genes but also additional genes of interest linked to this difficult gastrointestinal disease.

Evolution, Chance, and Aging

Aging has provided fruitful challenges for evolutionary theory, and evolutionary theory has deepened our understanding of aging. A great deal of genetic and molecular data now exists concerning mortality regulation and there is a growing body of knowledge concerning the life histories of diverse species. Assimilating all relevant data into a framework for the evolution of aging promises to significantly advance the field. We propose extensions of some key concepts to provide greater precision when applying these concepts to age-structured contexts. Secondary or byproduct effects of mutations are proposed as an important factor affecting survival patterns, including effects that may operate in small populations subject to genetic drift, widening the possibilities for mutation accumulation and pleiotropy. Molecular and genetic studies have indicated a diverse array of mechanisms that can modify aging and mortality rates, while transcriptome data indicate a high level of tissue and species specificity for genes affected by aging. The diversity of mechanisms and gene effects that can contribute to the pattern of aging in different organisms may mirror the complex evolutionary processes behind aging.

PRICKLE1 × FOCAD Interaction Revealed by Genome-Wide vQTL Analysis of Human Facial Traits

The human face is a highly complex and variable structure resulting from the intricate coordination of numerous genetic and non-genetic factors. Hundreds of genomic loci impacting quantitative facial features have been identified. While these associations have been shown to influence morphology by altering the mean size and shape of facial measures, their effect on trait variance remains unclear. We conducted a genome-wide association analysis for the variance of 20 quantitative facial measurements in 2,447 European individuals and identified several suggestive variance quantitative trait loci (vQTLs). These vQTLs guided us to conduct an efficient search for gene-by-gene (G × G) interactions, which uncovered an interaction between PRICKLE1 and FOCAD affecting cranial base width. We replicated this G × G interaction signal at the locus level in an additional 5,128 Korean individuals. We used the hypomorphic Prickle1 Beetlejuice (Prickle1 Bj ) mouse line to directly test the function of Prickle1 on the cranial base and observed wider cranial bases in Prickle1 Bj/Bj . Importantly, we observed that the Prickle1 and Focadhesin proteins co-localize in murine cranial base chondrocytes, and this co-localization is abnormal in the Prickle1 Bj/Bj mutants. Taken together, our findings uncovered a novel G × G interaction effect in humans with strong support from both epidemiological and molecular studies. These results highlight the potential of studying measures of phenotypic variability in gene mapping studies of facial morphology.

Multiomic Big Data Analysis Challenges: Increasing Confidence in the Interpretation of Artificial Intelligence Assessments

The need for holistic molecular measurements to better understand disease initiation, development, diagnosis, and therapy has led to an increasing number of multiomic analyses. The wealth of information available from multiomic assessments, however, requires both the evaluation and interpretation of extremely large data sets, limiting analysis throughput and ease of adoption. Computational methods utilizing artificial intelligence (AI) provide the most promising way to address these challenges, yet despite the conceptual benefits of AI and its successful application in singular omic studies, the widespread use of AI in multiomic studies remains limited. Here, we discuss present and future capabilities of AI techniques in multiomic studies while introducing analytical checks and balances to validate the computational conclusions.

Genetic Context Significantly Influences the Maintenance and Evolution of Degenerate Pathways

Understanding the evolution of novel physiological traits is highly relevant for expanding the characterization and manipulation of biological systems. Acquisition of new traits can be achieved through horizontal gene transfer (HGT). Here, we investigate drivers that promote or deter the maintenance of HGT-driven degeneracy, occurring when processes accomplish identical functions through nonidentical components. Subsequent evolution can optimize newly acquired functions; for example, beneficial alleles identified in an engineered Methylorubrum extorquens strain allowed it to utilize a “Foreign” formaldehyde oxidation pathway substituted for its Native pathway for methylotrophic growth. We examined the fitness consequences of interactions between these alleles when they were combined with the Native pathway or both (Dual) pathways. Unlike the Foreign pathway context where they evolved, these alleles were often neutral or deleterious when moved into these alternative genetic backgrounds. However, there were instances where combinations of multiple alleles resulted in higher fitness outcomes than individual allelic substitutions could provide. Importantly, the genetic context accompanying these allelic substitutions significantly altered the fitness landscape, shifting local fitness peaks and restricting the set of accessible evolutionary trajectories. These findings highlight how genetic context can negatively impact the probability of maintaining native and HGT-introduced functions together, making it difficult for degeneracy to evolve. However, in cases where the cost of maintaining degeneracy was mitigated by adding evolved alleles impacting the function of these pathways, we observed rare opportunities for pathway coevolution to occur. Together, our results highlight the importance of genetic context and resulting epistasis in retaining or losing HGT-acquired degenerate functions.

A model and test for coordinated polygenic epistasis in complex traits

Interactions between genetic variants-epistasis-is pervasive in model systems and can profoundly impact evolutionary adaption, population disease dynamics, genetic mapping, and precision medicine efforts. In this work, we develop a model for structured polygenic epistasis, called coordinated epistasis (CE), and prove that several recent theories of genetic architecture fall under the formal umbrella of CE. Unlike standard epistasis models that assume epistasis and main effects are independent, CE captures systematic correlations between epistasis and main effects that result from pathway-level epistasis, on balance skewing the penetrance of genetic effects. To test for the existence of CE, we propose the even-odd (EO) test and prove it is calibrated in a range of realistic biological models. Applying the EO test in the UK Biobank, we find evidence of CE in 18 of 26 traits spanning disease, anthropometric, and blood categories. Finally, we extend the EO test to tissue-specific enrichment and identify several plausible tissue-trait pairs. Overall, CE is a dimension of genetic architecture that can capture structured, systemic forms of epistasis in complex human traits.

A comparison of methods for interpreting random forest models of genetic association in the presence of non-additive interactions

BACKGROUND

Non-additive interactions among genes are frequently associated with a number of phenotypes, including known complex diseases such as Alzheimer's, diabetes, and cardiovascular disease. Detecting interactions requires careful selection of analytical methods, and some machine learning algorithms are unable or underpowered to detect or model feature interactions that exhibit non-additivity. The Random Forest method is often employed in these efforts due to its ability to detect and model non-additive interactions. In addition, Random Forest has the built-in ability to estimate feature importance scores, a characteristic that allows the model to be interpreted with the order and effect size of the feature association with the outcome. This characteristic is very important for epidemiological and clinical studies where results of predictive modeling could be used to define the future direction of the research efforts. An alternative way to interpret the model is with a permutation feature importance metric which employs a permutation approach to calculate a feature contribution coefficient in units of the decrease in the model's performance and with the Shapely additive explanations which employ cooperative game theory approach. Currently, it is unclear which Random Forest feature importance metric provides a superior estimation of the true informative contribution of features in genetic association analysis.

RESULTS

To address this issue, and to improve interpretability of Random Forest predictions, we compared different methods for feature importance estimation in real and simulated datasets with non-additive interactions. As a result, we detected a discrepancy between the metrics for the real-world datasets and further established that the permutation feature importance metric provides more precise feature importance rank estimation for the simulated datasets with non-additive interactions.

CONCLUSIONS

By analyzing both real and simulated data, we established that the permutation feature importance metric provides more precise feature importance rank estimation in the presence of non-additive interactions.

Personizing the prediction of future susceptibility to a specific disease

A traceable biomarker is a member of a disease's molecular pathway. A disease may be associated with several molecular pathways. Each different combination of these molecular pathways, to which detected traceable biomarkers belong, may serve as an indicative of the elicitation of the disease at a different time frame in the future. Based on this notion, we introduce a novel methodology for personalizing an individual's degree of future susceptibility to a specific disease. We implemented the methodology in a working system called Susceptibility Degree to a Disease Predictor (SDDP). For a specific disease d, let S be the set of molecular pathways, to which traceable biomarkers detected from most patients of d belong. For the same disease d, let S' be the set of molecular pathways, to which traceable biomarkers detected from a certain individual belong. SDDP is able to infer the subset S'' ⊆{S-S'} of undetected molecular pathways for the individual. Thus, SDDP can infer undetected molecular pathways of a disease for an individual based on few molecular pathways detected from the individual. SDDP can also help in inferring the combination of molecular pathways in the set {S'+S''}, whose traceable biomarkers collectively is an indicative of the disease. SDDP is composed of the following four components: information extractor, interrelationship between molecular pathways modeler, logic inferencer, and risk indicator. The information extractor takes advantage of the exponential increase of biomedical literature to automatically extract the common traceable biomarkers for a specific disease. The interrelationship between molecular pathways modeler models the hierarchical interrelationships between the molecular pathways of the traceable biomarkers. The logic inferencer transforms the hierarchical interrelationships between the molecular pathways into rule-based specifications. It employs the specification rules and the inference rules for predicate logic to infer as many as possible undetected molecular pathways of a disease for an individual. The risk indicator outputs a risk indicator value that reflects the individual's degree of future susceptibility to the disease. We evaluated SDDP by comparing it experimentally with other methods. Results revealed marked improvement.

A Framework for Efficient N-Way Interaction Testing in Case/Control Studies With Categorical Data

Goal: Most common diseases are influenced by multiple gene interactions and interactions with the environment. Performing an exhaustive search to identify such interactions is computationally expensive and needs to address the multiple testing problem. A four-step framework is proposed for the efficient identification of n-Way interactions. Methods: The framework was applied on a Multiple Sclerosis dataset with 725 subjects and 147 tagging SNPs. The first two steps of the framework are quality control and feature selection. The next step uses clustering and binary encodes the features. The final step performs the n-Way interaction testing. Results: The feature space was reduced to 7 SNPs and using the proposed binary encoding, more 2-SNP and 3-SNP interactions were identified compared to using the initial encoding. Conclusions: The framework selects informative features and with the proposed binary encoding it is able to identify more n-way interactions by increasing the power of the statistical analysis.

Novel Alzheimer risk genes determine the microglia response to amyloid-β but not to TAU pathology

Polygenic risk scores have identified that genetic variants without genome-wide significance still add to the genetic risk of developing Alzheimer's disease (AD). Whether and how subthreshold risk loci translate into relevant disease pathways is unknown. We investigate here the involvement of AD risk variants in the transcriptional responses of two mouse models: APPswe/PS1L166P and Thy-TAU22. A unique gene expression module, highly enriched for AD risk genes, is specifically responsive to Aβ but not TAU pathology. We identify in this module 7 established AD risk genes (APOE, CLU, INPP5D, CD33, PLCG2, SPI1, and FCER1G) and 11 AD GWAS genes below the genome-wide significance threshold (GPC2, TREML2, SYK, GRN, SLC2A5, SAMSN1, PYDC1, HEXB, RRBP1, LYN, and BLNK), that become significantly upregulated when exposed to Aβ. Single microglia sequencing confirms that Aβ, not TAU, pathology induces marked transcriptional changes in microglia, including increased proportions of activated microglia. We conclude that genetic risk of AD functionally translates into different microglia pathway responses to Aβ pathology, placing AD genetic risk downstream of the amyloid pathway but upstream of TAU pathology.

Genetic Modifiers and Rare Mendelian Disease

Despite advances in high-throughput sequencing that have revolutionized the discovery of gene defects in rare Mendelian diseases, there are still gaps in translating individual genome variation to observed phenotypic outcomes. While we continue to improve genomics approaches to identify primary disease-causing variants, it is evident that no genetic variant acts alone. In other words, some other variants in the genome (genetic modifiers) may alleviate (suppress) or exacerbate (enhance) the severity of the disease, resulting in the variability of phenotypic outcomes. Thus, to truly understand the disease, we need to consider how the disease-causing variants interact with the rest of the genome in an individual. Here, we review the current state-of-the-field in the identification of genetic modifiers in rare Mendelian diseases and discuss the potential for future approaches that could bridge the existing gap.

How Computational Experiments Can Improve Our Understanding of the Genetic Architecture of Common Human Diseases

Susceptibility to common human diseases such as cancer is influenced by many genetic and environmental factors that work together in a complex manner. The state of the art is to perform a genome-wide association study (GWAS) that measures millions of single-nucleotide polymorphisms (SNPs) throughout the genome followed by a one-SNP-at-a-time statistical analysis to detect univariate associations. This approach has identified thousands of genetic risk factors for hundreds of diseases. However, the genetic risk factors detected have very small effect sizes and collectively explain very little of the overall heritability of the disease. Nonetheless, it is assumed that the genetic component of risk is due to many independent risk factors that contribute additively. The fact that many genetic risk factors with small effects can be detected is taken as evidence to support this notion. It is our working hypothesis that the genetic architecture of common diseases is partly driven by non-additive interactions. To test this hypothesis, we developed a heuristic simulation-based method for conducting experiments about the complexity of genetic architecture. We show that a genetic architecture driven by complex interactions is highly consistent with the magnitude and distribution of univariate effects seen in real data. We compare our results with measures of univariate and interaction effects from two large-scale GWASs of sporadic breast cancer and find evidence to support our hypothesis that is consistent with the results of our computational experiment.

Epistasis between phenylethanolamine N-methyltransferase and β2-adrenergic receptor influences extracellular epinephrine level and associates with the susceptibility to allergic asthma

BACKGROUND

Reduced extracellular epinephrine level often associates with asthma-related symptoms; however, the correlation between asthma and genetic variants in genes participating in the epinephrine signalling pathway remains unclear.

OBJECTIVE

To characterize the functions of single nucleotide polymorphisms (SNPs) in phenylethanolamine N-methyltransferase (PNMT) and β2-adrenergic receptor (ADRB2), and to study the effects, including both direct and epistatic, of these SNPs on serum epinephrine level and asthma susceptibility.

METHODS

Single nucleotide polymorphisms functions were characterized through in vitro luciferase assay. ADRB2 gene expression level in peripheral blood mononuclear cell (PBMC) was measured by transcriptome sequencing and expression microarray on two separate Asian cohorts (NUS-UTAR, n = 278 and NUS-TA, n = 58). Serum epinephrine level was assessed on a Singapore Chinese cohort (NUS-SH, n = 314) with 155 asthmatic and 159 non-asthmatic subjects. A separate Singapore Chinese cohort (NUS-G, n = 3009) was genotyped to show disease association (direct and epistatic effect) of functional SNPs in PNMT and ADRB2.

RESULTS

Reduced serum epinephrine level was associated with increased asthma risk in Singapore Chinese. The minor allele of rs876493 was shown to increase PNMT promoter activity and reduce asthma risk. Multiple SNPs in ADRB2 forms a haplotype that was associated with the differential promoter activity of this gene. In this haplotype, rs11168070 was associated directly with ADRB2 expression in PBMCs. Both minor alleles from rs876493 and rs11168070 contribute synergistically to reduce asthma risk and increase serum epinephrine level.

CONCLUSION AND CLINICAL RELEVANCE

Epistatic interaction between genetic variants from PNMT (rs876493) and ADRB2 (rs11168070) is associated with serum epinephrine level and the susceptibility of asthma. Our findings improved the current understanding of the genetic basis of this disease, while genotypic states of these SNPs may serve as potential biomarkers to predict susceptibility to the disease.

Selecting Closely-Linked SNPs Based on Local Epistatic Effects for Haplotype Construction Improves Power of Association Mapping

Genome-wide association studies (GWAS) have gained central importance for the identification of candidate loci underlying complex traits. Single nucleotide polymorphism (SNP) markers are mostly used as genetic variants for the analysis of genotype-phenotype associations in populations, but closely linked SNPs that are grouped into haplotypes are also exploited. The benefit of haplotype-based GWAS approaches vs. SNP-based approaches is still under debate because SNPs in high linkage disequilibrium provide redundant information. To overcome some constraints of the commonly-used haplotype-based GWAS in which only consecutive SNPs are considered for haplotype construction, we propose a new method called functional haplotype-based GWAS (FH GWAS). FH GWAS is featured by combining SNPs into haplotypes based on the additive and epistatic effects among SNPs. Such haplotypes were termed functional haplotypes (FH). As shown by simulation studies, the FH GWAS approach clearly outperformed the SNP-based approach unless the minor allele frequency of the SNPs making up the haplotypes is low and the linkage disequilibrium between them is high. Applying FH GWAS for the trait flowering time in a large Arabidopsis thaliana population with whole-genome sequencing data revealed its potential empirically. FH GWAS identified all candidate regions which were detected in SNP-based and two other haplotype-based GWAS approaches. In addition, a novel region on chromosome 4 was solely detected by FH GWAS. Thus both the results of our simulation and empirical studies demonstrate that FH GWAS is a promising method and superior to the SNP-based approach even if almost complete genotype information is available.

Gene-gene interaction: the curse of dimensionality

Identified genetic variants from genome wide association studies frequently show only modest effects on the disease risk, leading to the "missing heritability" problem. An avenue, to account for a part of this "missingness" is to evaluate gene-gene interactions (epistasis) thereby elucidating their effect on complex diseases. This can potentially help with identifying gene functions, pathways, and drug targets. However, the exhaustive evaluation of all possible genetic interactions among millions of single nucleotide polymorphisms (SNPs) raises several issues, otherwise known as the "curse of dimensionality". The dimensionality involved in the epistatic analysis of such exponentially growing SNPs diminishes the usefulness of traditional, parametric statistical methods. With the immense popularity of multifactor dimensionality reduction (MDR), a non-parametric method, proposed in 2001, that classifies multi-dimensional genotypes into one- dimensional binary approaches, led to the emergence of a fast-growing collection of methods that were based on the MDR approach. Moreover, machine-learning (ML) methods such as random forests and neural networks (NNs), deep-learning (DL) approaches, and hybrid approaches have also been applied profusely, in the recent years, to tackle this dimensionality issue associated with whole genome gene-gene interaction studies. However, exhaustive searching in MDR based approaches or variable selection in ML methods, still pose the risk of missing out on relevant SNPs. Furthermore, interpretability issues are a major hindrance for DL methods. To minimize this loss of information, Python based tools such as PySpark can potentially take advantage of distributed computing resources in the cloud, to bring back smaller subsets of data for further local analysis. Parallel computing can be a powerful resource that stands to fight this "curse". PySpark supports all standard Python libraries and C extensions thus making it convenient to write codes to deliver dramatic improvements in processing speed for extraordinarily large sets of data.

Flux balance analysis with or without molecular crowding fails to predict two thirds of experimentally observed epistasis in yeast

Computational predictions of double gene knockout effects by flux balance analysis (FBA) have been used to characterized genome-wide patterns of epistasis in microorganisms. However, it is unclear how in silico predictions are related to in vivo epistasis, as FBA predicted only a minority of experimentally observed genetic interactions between non-essential metabolic genes in yeast. Here, we perform a detailed comparison of yeast experimental epistasis data to predictions generated with different constraint-based metabolic modeling algorithms. The tested methods comprise standard FBA; a variant of MOMA, which was specifically designed to predict fitness effects of non-essential gene knockouts; and two alternative implementations of FBA with macro-molecular crowding, which account approximately for enzyme kinetics. The number of interactions uniquely predicted by one method is typically larger than its overlap with any alternative method. Only 20% of negative and 10% of positive interactions jointly predicted by all methods are confirmed by the experimental data; almost all unique predictions appear to be false. More than two thirds of epistatic interactions are undetectable by any of the tested methods. The low prediction accuracies indicate that the physiology of yeast double metabolic gene knockouts is dominated by processes not captured by current constraint-based analysis methods.

Confounding of linkage disequilibrium patterns in large scale DNA based gene-gene interaction studies

Background

In Genome-Wide Association Studies (GWAS), the concept of linkage disequilibrium is important as it allows identifying genetic markers that tag the actual causal variants. In Genome-Wide Association Interaction Studies (GWAIS), similar principles hold for pairs of causal variants. However, Linkage Disequilibrium (LD) may also interfere with the detection of genuine epistasis signals in that there may be complete confounding between Gametic Phase Disequilibrium (GPD) and interaction. GPD may involve unlinked genetic markers, even residing on different chromosomes. Often GPD is eliminated in GWAIS, via feature selection schemes or so-called pruning algorithms, to obtain unconfounded epistasis results. However, little is known about the optimal degree of GPD/LD-pruning that gives a balance between false positive control and sufficient power of epistasis detection statistics. Here, we focus on Model-Based Multifactor Dimensionality Reduction as one large-scale epistasis detection tool. Its performance has been thoroughly investigated in terms of false positive control and power, under a variety of scenarios involving different trait types and study designs, as well as error-free and noisy data, but never with respect to multicollinear SNPs.

Results

Using real-life human LD patterns from a homogeneous subpopulation of British ancestry, we investigated the impact of LD-pruning on the statistical sensitivity of MB-MDR. We considered three different non-fully penetrant epistasis models with varying effect sizes. There is a clear advantage in pre-analysis pruning using sliding windows at r² of 0.75 or lower, but using a threshold of 0.20 has a detrimental effect on the power to detect a functional interactive SNP pair (power < 25%). Signal sensitivity, directly using LD-block information to determine whether an epistasis signal is present or not, benefits from LD-pruning as well (average power across scenarios: 87%), but is largely hampered by functional loci residing at the boundaries of an LD-block.

Conclusions

Our results confirm that LD patterns and the position of causal variants in LD blocks do have an impact on epistasis detection, and that pruning strategies and LD-blocks definitions combined need careful attention, if we wish to maximize the power of large-scale epistasis screenings.

Assessing the fitness consequences of mitonuclear interactions in natural populations

Metazoans exist only with a continuous and rich supply of chemical energy from oxidative phosphorylation in mitochondria. The oxidative phosphorylation machinery that mediates energy conservation is encoded by both mitochondrial and nuclear genes, and hence the products of these two genomes must interact closely to achieve coordinated function of core respiratory processes. It follows that selection for efficient respiration will lead to selection for compatible combinations of mitochondrial and nuclear genotypes, and this should facilitate coadaptation between mitochondrial and nuclear genomes (mitonuclear coadaptation). Herein, we outline the modes by which mitochondrial and nuclear genomes may coevolve within natural populations, and we discuss the implications of mitonuclear coadaptation for diverse fields of study in the biological sciences. We identify five themes in the study of mitonuclear interactions that provide a roadmap for both ecological and biomedical studies seeking to measure the contribution of intergenomic coadaptation to the evolution of natural populations. We also explore the wider implications of the fitness consequences of mitonuclear interactions, focusing on central debates within the fields of ecology and biomedicine.

How to increase our belief in discovered statistical interactions via large-scale association studies?

The understanding that differences in biological epistasis may impact disease risk, diagnosis, or disease management stands in wide contrast to the unavailability of widely accepted large-scale epistasis analysis protocols. Several choices in the analysis workflow will impact false-positive and false-negative rates. One of these choices relates to the exploitation of particular modelling or testing strategies. The strengths and limitations of these need to be well understood, as well as the contexts in which these hold. This will contribute to determining the potentially complementary value of epistasis detection workflows and is expected to increase replication success with biological relevance. In this contribution, we take a recently introduced regression-based epistasis detection tool as a leading example to review the key elements that need to be considered to fully appreciate the value of analytical epistasis detection performance assessments. We point out unresolved hurdles and give our perspectives towards overcoming these.

Generalized multifactor dimensionality reduction approaches to identification of genetic interactions underlying ordinal traits

The manifestation of complex traits is influenced by gene-gene and gene-environment interactions, and the identification of multifactor interactions is an important but challenging undertaking for genetic studies. Many complex phenotypes such as disease severity are measured on an ordinal scale with more than two categories. A proportional odds model can improve statistical power for these outcomes, when compared to a logit model either collapsing the categories into two mutually exclusive groups or limiting the analysis to pairs of categories. In this study, we propose a proportional odds model-based generalized multifactor dimensionality reduction (GMDR) method for detection of interactions underlying polytomous ordinal phenotypes. Computer simulations demonstrated that this new GMDR method has a higher power and more accurate predictive ability than the GMDR methods based on a logit model and a multinomial logit model. We applied this new method to the genetic analysis of low-density lipoprotein (LDL) cholesterol, a causal risk factor for coronary artery disease, in the Multi-Ethnic Study of Atherosclerosis, and identified a significant joint action of the CELSR2, SERPINA12, HPGD, and APOB genes. This finding provides new information to advance the limited knowledge about genetic regulation and gene interactions in metabolic pathways of LDL cholesterol. In conclusion, the proportional odds model-based GMDR is a useful tool that can boost statistical power and prediction accuracy in studying multifactor interactions underlying ordinal traits.

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.

Benchmarking relief-based feature selection methods for bioinformatics data mining

Modern biomedical data mining requires feature selection methods that can (1) be applied to large scale feature spaces (e.g. 'omics' data), (2) function in noisy problems, (3) detect complex patterns of association (e.g. gene-gene interactions), (4) be flexibly adapted to various problem domains and data types (e.g. genetic variants, gene expression, and clinical data) and (5) are computationally tractable. To that end, this work examines a set of filter-style feature selection algorithms inspired by the 'Relief' algorithm, i.e. Relief-Based algorithms (RBAs). We implement and expand these RBAs in an open source framework called ReBATE (Relief-Based Algorithm Training Environment). We apply a comprehensive genetic simulation study comparing existing RBAs, a proposed RBA called MultiSURF, and other established feature selection methods, over a variety of problems. The results of this study (1) support the assertion that RBAs are particularly flexible, efficient, and powerful feature selection methods that differentiate relevant features having univariate, multivariate, epistatic, or heterogeneous associations, (2) confirm the efficacy of expansions for classification vs. regression, discrete vs. continuous features, missing data, multiple classes, or class imbalance, (3) identify previously unknown limitations of specific RBAs, and (4) suggest that while MultiSURF∗ performs best for explicitly identifying pure 2-way interactions, MultiSURF yields the most reliable feature selection performance across a wide range of problem types.

Leveraging epigenomics and contactomics data to investigate SNP pairs in GWAS

Although Genome Wide Association Studies (GWAS) have led to many valuable insights into the genetic bases of common diseases over the past decade, the issue of missing heritability has surfaced, as the discovered main effect genetic variants found to date do not account for much of a trait's predicted genetic component. We present a workflow, integrating epigenomics and topologically associating domain data, aimed at discovering trait-associated SNP pairs from GWAS where neither SNP achieved independent genome-wide significance. Each analyzed SNP pair consists of one SNP in a putative active enhancer and another SNP in a putative physically interacting gene promoter in a trait-relevant tissue. As a proof-of-principle case study, we used this approach to identify focused collections of SNP pairs that we analyzed in three independent Type 2 diabetes (T2D) GWAS. This approach led us to discover 35 significant SNP pairs, encompassing both novel signals and signals for which we have found orthogonal support from other sources. Nine of these pairs are consistent with eQTL results, two are consistent with our own capture C experiments, and seven involve signals supported by recent T2D literature.

The Mitonuclear Dimension of Neanderthal and Denisovan Ancestry in Modern Human Genomes

Some human populations interbred with Neanderthals and Denisovans, resulting in substantial contributions to modern-human genomes. Therefore, it is now possible to use genomic data to investigate mechanisms that shaped historical gene flow between humans and our closest hominin relatives. More generally, in eukaryotes, mitonuclear interactions have been argued to play a disproportionate role in generating reproductive isolation. There is no evidence of mtDNA introgression into modern human populations, which means that all introgressed nuclear alleles from archaic hominins must function on a modern-human mitochondrial background. Therefore, mitonuclear interactions are also potentially relevant to hominin evolution. We performed a detailed accounting of mtDNA divergence among hominin lineages and used population-genomic data to test the hypothesis that mitonuclear incompatibilities have preferentially restricted the introgression of nuclear genes with mitochondrial functions. We found a small but significant underrepresentation of introgressed Neanderthal alleles at such nuclear loci. Structural analyses of mitochondrial enzyme complexes revealed that these effects are unlikely to be mediated by physically interacting sites in mitochondrial and nuclear gene products. We did not detect any underrepresentation of introgressed Denisovan alleles at mitochondrial-targeted loci, but this may reflect reduced power because locus-specific estimates of Denisovan introgression are more conservative. Overall, we conclude that genes involved in mitochondrial function may have been subject to distinct selection pressures during the history of introgression from archaic hominins but that mitonuclear incompatibilities have had, at most, a small role in shaping genome-wide introgression patterns, perhaps because of limited functional divergence in mtDNA and interacting nuclear genes.

Alzheimer's disease susceptibility genes modify the risk of Parkinson disease and Parkinson's disease-associated cognitive impairment

The pathogenic mechanism underlying Parkinson's disease (PD) and PD- Cognitive impairment (CI) remains elusive. Its potential link to the risk factors in Alzheimer's disease (AD) is unclear. In this study, we analyzed 16 CE-associated single nucleotide polymorphisms (SNPs) in twelve genes in a Chinese cohort of 450 PD cases and 449 controls. Among our 298 cases clinically evaluated for CI, 113 cases did not show CI signs (PD-NC), 86 cases had mildly cognitive impairment (PD-MCI) and 99 cases had dementia (PD-D). We found that the APOE ε4 allele is associated with a higher risk for PD-D. Gene-gene interaction analysis revealed that three significant gene-gene interactions, including BDNF and CLU, APOE and CR1, and DYRK1A and CD2AP increase the risk for PD. Because these SNPs are known genetic risk factors for AD, their contribution to PD and PD-D shown in this study suggests that PD/PD-D and AD may share convergent pathways in their pathogenesis through gene-gene interactions.

Improving machine learning reproducibility in genetic association studies with proportional instance cross validation (PICV)

Background

Machine learning methods and conventions are increasingly employed for the analysis of large, complex biomedical data sets, including genome-wide association studies (GWAS). Reproducibility of machine learning analyses of GWAS can be hampered by biological and statistical factors, particularly so for the investigation of non-additive genetic interactions. Application of traditional cross validation to a GWAS data set may result in poor consistency between the training and testing data set splits due to an imbalance of the interaction genotypes relative to the data as a whole. We propose a new cross validation method, proportional instance cross validation (PICV), that preserves the original distribution of an independent variable when splitting the data set into training and testing partitions.

Results

We apply PICV to simulated GWAS data with epistatic interactions of varying minor allele frequencies and prevalences and compare performance to that of a traditional cross validation procedure in which individuals are randomly allocated to training and testing partitions. Sensitivity and positive predictive value are significantly improved across all tested scenarios for PICV compared to traditional cross validation. We also apply PICV to GWAS data from a study of primary open-angle glaucoma to investigate a previously-reported interaction, which fails to significantly replicate; PICV however improves the consistency of testing and training results.

Conclusions

Application of traditional machine learning procedures to biomedical data may require modifications to better suit intrinsic characteristics of the data, such as the potential for highly imbalanced genotype distributions in the case of epistasis detection. The reproducibility of genetic interaction findings can be improved by considering this variable imbalance in cross validation implementation, such as with PICV. This approach may be extended to problems in other domains in which imbalanced variable distributions are a concern.

The search for gene-gene interactions in genome-wide association studies: challenges in abundance of methods, practical considerations, and biological interpretation

One of the primary goals in this era of precision medicine is to understand the biology of human diseases and their treatment, such that each individual patient receives the best possible treatment for their disease based on their genetic and environmental exposures. One way to work towards achieving this goal is to identify the environmental exposures and genetic variants that are relevant to each disease in question, as well as the complex interplay between genes and environment. Genome-wide association studies (GWAS) have allowed for a greater understanding of the genetic component of many complex traits. However, these genetic effects are largely small and thus, our ability to use these GWAS finding for precision medicine is limited. As more and more GWAS have been performed, rather than focusing only on common single nucleotide polymorphisms (SNPs) and additive genetic models, many researchers have begun to explore alternative heritable components of complex traits including rare variants, structural variants, epigenetics, and genetic interactions. While genetic interactions are a plausible reality that could explain some of the heritabliy that has not yet been identified, especially when one considers the identification of genetic interactions in model organisms as well as our understanding of biological complexity, still there are significant challenges and considerations in identifying these genetic interactions. Broadly, these can be summarized in three categories: abundance of methods, practical considerations, and biological interpretation. In this review, we will discuss these important elements in the search for genetic interactions along with some potential solutions. While genetic interactions are theoretically understood to be important for complex human disease, the body of evidence is still building to support this component of the underlying genetic architecture of complex human traits. Our hope is that more sophisticated modeling approaches and more robust computational techniques will enable the community to identify these important genetic interactions and improve our ability to implement precision medicine in the future.

A Review on Methods for Detecting SNP Interactions in High-Dimensional Genomic Data

In this era of genome-wide association studies (GWAS), the quest for understanding the genetic architecture of complex diseases is rapidly increasing more than ever before. The development of high throughput genotyping and next generation sequencing technologies enables genetic epidemiological analysis of large scale data. These advances have led to the identification of a number of single nucleotide polymorphisms (SNPs) responsible for disease susceptibility. The interactions between SNPs associated with complex diseases are increasingly being explored in the current literature. These interaction studies are mathematically challenging and computationally complex. These challenges have been addressed by a number of data mining and machine learning approaches. This paper reviews the current methods and the related software packages to detect the SNP interactions that contribute to diseases. The issues that need to be considered when developing these models are addressed in this review. The paper also reviews the achievements in data simulation to evaluate the performance of these models. Further, it discusses the future of SNP interaction analysis.

Another Round of "Clue" to Uncover the Mystery of Complex Traits

A plethora of genetic association analyses have identified several genetic risk loci. Technological and statistical advancements have now led to the identification of not only common genetic variants, but also low-frequency variants, structural variants, and environmental factors, as well as multi-omics variations that affect the phenotypic variance of complex traits in a population, thus referred to as complex trait architecture. The concept of heritability, or the proportion of phenotypic variance due to genetic inheritance, has been studied for several decades, but its application is mainly in addressing the narrow sense heritability (or additive genetic component) from Genome-Wide Association Studies (GWAS). In this commentary, we reflect on our perspective on the complexity of understanding heritability for human traits in comparison to model organisms, highlighting another round of clues beyond GWAS and an alternative approach, investigating these clues comprehensively to help in elucidating the genetic architecture of complex traits.