EINVis: a visualization tool for analyzing and exploring genetic interactions in large-scale association studies

Polygenic risk scores: the future of cancer risk prediction, screening, and precision prevention

Genome-wide association studies (GWASs) have shown that the genetic architecture of cancers are highly polygenic and enabled researchers to identify genetic risk loci for cancers. The genetic variants associated with a cancer can be combined into a polygenic risk score (PRS), which captures part of an individual's genetic susceptibility to cancer. Recently, PRSs have been widely used in cancer risk prediction and are shown to be capable of identifying groups of individuals who could benefit from the knowledge of their probabilistic susceptibility to cancer, which leads to an increased interest in understanding the potential utility of PRSs that might further refine the assessment and management of cancer risk. In this context, we provide an overview of the major discoveries from cancer GWASs. We then review the methodologies used for PRS construction, and describe steps for the development and evaluation of risk prediction models that include PRS and/or conventional risk factors. Potential utility of PRSs in cancer risk prediction, screening, and precision prevention are illustrated. Challenges and practical considerations relevant to the implementation of PRSs in health care settings are discussed.

Genetic associations of the vitamin D and antiviral pathways with natural resistance to HIV-1 infection are influenced by interpopulation variability

Vitamin D (VitD) may modulate anti-HIV-1 responses modifying the risk to acquire the HIV-1-infection. We performed a nested case-control exploratory study involving 413 individuals; HIV-1-exposed seropositives (cases) and seronegatives (HESN) (controls) from three cohorts: sexually-exposed from Colombia and Italy and parenterally-exposed from Spain. The association and interactions of 139 variants in 9 VitD pathway genes, and in 14 antiviral genes with resistance/susceptibility (R/S) to HIV-1 infection was evaluated. Associations between variants and mRNA levels were also analyzed in the Colombian samples. Variants and haplotypes in genes of VitD and antiviral pathways were associated with R/S, but specific associations were not reproduced in all cohorts. Allelic heterogeneity could explain such inconsistency since the associations found in all cohorts were consistently in the same genes: VDR and RXRA of the VitD pathway genes and in TLR2 and RNASE4. Remarkably, the multi-locus genotypes (interacting variants) observed in genes of VitD and antiviral pathways were present in most HESNs of all cohorts. Finally, HESNs carrying resistance-associated variants had higher levels of VitD in plasma, of VDR mRNA in blood cells, and of ELAFIN and defensins mRNA in the oral mucosa. In conclusion, despite allelic heterogeneity, most likely due to differences in the genetic history of the populations, the associations were locus dependent suggesting that genes of the VitD pathway might act in concert with antiviral genes modulating the resistance phenotype of the HESNs. Although these associations were significant after permutation test, only haplotype results remained statistically significant after Bonferroni test, requiring further replications in larger cohorts and functional analyzes to validate these conclusions.

HC-HDSD: A method of hypergraph construction and high-density subgraph detection for inferring high-order epistatic interactions

Detecting epistatic interactions, or nonlinear interactive effects of Single Nucleotide Polymorphisms (SNPs), has gained increasing attention in explaining the "missing heritability" of complex diseases. Though much work has been done in mapping SNPs underlying diseases, most of them constrain to 2-order epistatic interactions. In this paper, a method of hypergraph construction and high-density subgraph detection, named HC-HDSD, is proposed for detecting high-order epistatic interactions. The hypergraph is constructed by low-order epistatic interactions that identified using the normalized co-information measure and the exhaustive search. The hypergraph consists of two types of vertices: real ones representing main effects of SNPs and virtual ones denoting interactive effects of epistatic interactions. Then, both maximal clique centrality algorithm and near-clique mining algorithm are employed to detect high-density subgraphs from the constructed hypergraph. These high-density subgraphs are inferred as high-order epistatic interactions in the HC-HDSD. Experiments are performed on several simulation data sets, results of which show that HC-HDSD is promising in inferring high-order epistatic interactions while substantially reducing the computation cost. In addition, the application of HC-HDSD on a real Age-related Macular Degeneration (AMD) data set provides several new clues for the exploration of causative factors of AMD.

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

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

CINOEDV: a co-information based method for detecting and visualizing n-order epistatic interactions

BACKGROUND

Detecting and visualizing nonlinear interaction effects of single nucleotide polymorphisms (SNPs) or epistatic interactions are important topics in bioinformatics since they play an important role in unraveling the mystery of "missing heritability". However, related studies are almost limited to pairwise epistatic interactions due to their methodological and computational challenges.

RESULTS

We develop CINOEDV (Co-Information based N-Order Epistasis Detector and Visualizer) for the detection and visualization of epistatic interactions of their orders from 1 to n (n ≥ 2). CINOEDV is composed of two stages, namely, detecting stage and visualizing stage. In detecting stage, co-information based measures are employed to quantify association effects of n-order SNP combinations to the phenotype, and two types of search strategies are introduced to identify n-order epistatic interactions: an exhaustive search and a particle swarm optimization based search. In visualizing stage, all detected n-order epistatic interactions are used to construct a hypergraph, where a real vertex represents the main effect of a SNP and a virtual vertex denotes the interaction effect of an n-order epistatic interaction. By deeply analyzing the constructed hypergraph, some hidden clues for better understanding the underlying genetic architecture of complex diseases could be revealed.

CONCLUSIONS

Experiments of CINOEDV and its comparison with existing state-of-the-art methods are performed on both simulation data sets and a real data set of age-related macular degeneration. Results demonstrate that CINOEDV is promising in detecting and visualizing n-order epistatic interactions. CINOEDV is implemented in R and is freely available from R CRAN: http://cran.r-project.org and https://sourceforge.net/projects/cinoedv/files/ .

An interaction quantitative trait loci tool implicates epistatic functional variants in an apoptosis pathway in smallpox vaccine eQTL data

Expression quantitative trait loci (eQTL) studies have functionalized nucleic acid variants through the regulation of gene expression. Although most eQTL studies only examine the effects of single variants on transcription, a more complex process of variant-variant interaction (epistasis) may regulate transcription. Herein, we describe a tool called interaction QTL (iQTL) designed to efficiently detect epistatic interactions that regulate gene expression. To maximize biological relevance and minimize the computational and hypothesis testing burden, iQTL restricts interactions such that one variant is within a user-defined proximity of the transcript (cis-regulatory). We apply iQTL to a data set of 183 smallpox vaccine study participants with genome-wide association study and gene expression data from unstimulated samples and samples stimulated by inactivated vaccinia virus. While computing only 0.15% of possible interactions, we identify 11 probe sets whose expression is regulated through a variant-variant interaction. We highlight the functional epistatic interactions among apoptosis-related genes, DIABLO, TRAPPC4 and FADD, in the context of smallpox vaccination. We also use an integrative network approach to characterize these iQTL interactions in a posterior network of known prior functional interactions. iQTL is an efficient, open-source tool to analyze variant interactions in eQTL studies, providing better understanding of the function of epistasis in immune response and other complex phenotypes.

Functional dyadicity and heterophilicity of gene-gene interactions in statistical epistasis networks

BACKGROUND

The interaction effect among multiple genetic factors, i.e. epistasis, plays an important role in explaining susceptibility on common human diseases and phenotypic traits. The uncertainty over the number of genetic attributes involved in interactions poses great challenges in genetic association studies and calls for advanced bioinformatics methodologies. Network science has gained popularity in modeling genetic interactions thanks to its structural characterization of large numbers of entities and their complex relationships. However, little has been done on functionally interpreting statistically inferred epistatic interactions using networks.

RESULTS

In this study, we propose to characterize gene functional properties in the context of interaction network structure. We used Gene Ontology (GO) to functionally annotate genes as vertices in a statistical epistasis network, and quantitatively characterize the correlation between the distribution of gene functional properties and the network structure by measuring dyadicity and heterophilicity of each functional category in the network. These two parameters quantify whether genetic interactions tend to occur more frequently for genes from the same functional category, i.e. dyadic effect, or more frequently for genes from across different functional categories, i.e. heterophilic effect.

CONCLUSIONS

By applying this framework to a population-based bladder cancer dataset, we were able to identify several GO categories that have significant dyadicity or heterophilicity associated with bladder cancer susceptibility. Thus, our informatics framework suggests a new methodology for embedding functional analysis in network modeling of statistical epistasis in genetic association studies.