it has long been appreciated that the regulation of biological processes involves a complex orchestration of genes in networks .
however , typical analytic frameworks for genome - wide association studies ( gwas ) have assumed a simplified genetic architecture , largely considering independent additive contributions to genetic risk . while recent studies have supported a large contribution of additive genetic variance to complex traits ( yang et al . , 2010
; lee et al . , 2011a ) ,  missing heritability  remains , which has led to a range of potential explanations ( eichler et al . , 2010 ) .
one such explanation is the existence of epistatic ( gene by gene ) effects , which are , in fact , not precluded even in the presence of sizeable additive genetic variance ( hill et al . , 2008 ) .
in this opinion and review , we propose that there is room for additional mining of datasets generated in the gwas era using a paradigm of networks that aggregate gene  gene interactions and main effects .
we argue that network approaches may have utility not only for discovery , but also for further characterizing relationships with genetic effects that may have been identified through standard analytic means .
the recognition that numerous genetic variants may act in concert to modulate disease susceptibility has led to the development of gene set enrichment and pathway approaches ( torkamani et al . , 2008 ; shi et al . , 2009 ; liu et al . , 2010 ) .
the small but consistent effect of many snps in gene sets suggests evidence at the population level of the coordination of genes known to interact through particular biological pathways .
gene set enrichment approaches are able to identify these coordination effects despite typically relying on small effect single - marker association evidence .
gene interaction effects also show small effect sizes , but evidence for coordination from main effects suggests that the aggregation of gene  gene interactions and main effects in epistasis networks may lead to even more consistent pathway enrichment .
epistasis networks may also be used to test the hypothesis of a network paradigm of evolution and disease susceptibility .
a recent study comparing the connectivity changes of arabidopsis networks following gene duplication events suggests a possible model of evolution acting at the level of the interactome network ( dreze et al . , 2011 ) .
these studies suggest improved understanding of disease susceptibility may be achieved by conceptualizing the genotype to phenotype mapping as a network of coupled gene  gene interactions and main effects .
we refer to these as genetic association interaction networks ( gain ) or epistasis networks ; however , the networks we describe model main effects as well as epistasis . in a recent gain analysis of a study of the immune response to smallpox vaccine
, we identified a new association in the rxra gene reflecting a large number of gene  gene interactions ( davis et al . , 2010 ) .
variants in figure 1 after the ability to connect other actors to him in a few jumps in a network constructed by shared movie credit .
such variants may be important to a phenotype , not because of their individual effect , but because of their overall influence in modulating the effect of other variants .
another analogy from popular culture is lady gaga , who has over one million twitter followers .
even a small perturbation from such a node may have a large downstream effect due to a cascade through the subnetworks of followers .
the brando node would be easily found by a standard single - locus statistic , but the kevin bacon node would only be revealed by an epistasis network approach due to its many small gene  gene interactions .
the epistasis network may also be useful for identifying new mechanisms for known effects , such as the connection of the brando node to the pathway represented by the subnetwork below it . returning to acting as an analogy ,  marlon brando
 variants ( the usual target of the gwas approach ) may or may not have a large network centrality , but are deemed important because of their individual ( main ) effect .
variants could include the cyp2c9 and vkor1c polymorphisms involved in warfarin metabolism ( limdi et al .
, 2010 ) , the wealth of human leukocyte antigen ( hla ) alleles associated with immunologic phenotypes ( lechler and warrens , 2000 ; blackwell et al . ,
2009 ) , or apol1 variants associated with kidney disease in african - americans ( genovese et al . , 2010 ) .
while the majority of variants identified to date do not exhibit such individually strong effects , most are only likely in linkage disequilibrium ( ld ) with causal alleles ; given that the gwas approach relies upon tagsnps to efficiently span the entire genome .
thus , there is an expectation that additional brando snps ( presumably functional variants tied to more penetrant phenotypic effects ) will be found by fine - mapping experiments ( such as the immunochip for autoimmune disorders ( trynka et al . , 2011 ) , trans - racial mapping ( rosenberg et al . , 2010 ; teo et al . ,
recent studies incorporating next - generation sequencing in age - related macular degeneration ( raychaudhuri et al . , 2011 ) and inflammatory bowel disease ( rivas et al . ,
, when considering the issue of the  missing heritability ,  it is important to remember that causal variants are expected to occur at lower frequencies compared to the tagsnps that have been identified thus far ( visscher et al . , 2011 ) .
therefore , the existence of stronger effects will be offset to an unknown extent by a lower prevalence of causal alleles , potentially restricting the total variability that will be explained at the population level .
it seems unlikely that identifying less frequent , causal alleles will immediately fill the missing heritability void , suggesting that complementary approaches , such as understanding variation in a network context , could have important implications in biological mechanistic research and translational medicine .
for example , even for well - known genetic factors such as hla , there is a growing recognition of epistatic considerations , such as the interplay between hla class i polymorphisms and the killer cell immunoglobulin receptor ( kir ) gene complex within the context of infectious disease , autoimmune disorders , cancer , and bone marrow transplantation outcomes ( kulkarni et al . , 2008 ;
thus , network analysis is not restricted to discovery . even in cases where there are known genetic effects , network analysis provides an opportunity to more fully characterize the genetic etiology by identifying the full network of interactions with known susceptibility variants .
the focus of this review is on empirical network approaches , where connections are based on gene  gene interactions ( epistasis ) estimated from gwas data . however , we pause to mention prior knowledge networks ( pkns ) , such as ingenuity pathway analysis ( www.ingenuity.com ) , which represent the most popular network analysis approach applied to gwas to date ( burgner et al . , 2009 ) .
similar to gene set enrichment approaches , pkns typically utilize statistical measures of marginal association to select genes to include in a putative network , with connections then constructed from various knowledge bases including protein 
protein interactions ( ppis ) , co - occurrence in literature mining , or co - expression in microarray experiments .
, 2011b ) uses evidence from pkns based on selected candidate pathways combined with odds ratios from gwas in order to boost the importance of neglected genes .
prior knowledge based strategies have also been used to focus searches for epistatic effects in multiple sclerosis ( bush et al . , 2011 ) and bipolar disorder ( moskvina et al . ,
2011 ) found that a prioritization strategy did not lead to an enrichment of replicable statistical interactions , suggesting a large risk of false positive associations .
these contrasting results could reflect the issue that connections are typically defined independently of the environmental / disease context ; leading to differences in performance as a function of the particular trait under study .
while the knowledge bases for pkns are constantly improving in quality and scope , there is likely a need for empirical network approaches that are able to account for the conditional dependence of gene connections on the specific disease under study .
for example , in studies of infection or vaccination , it is reasonable to hypothesize that genetic effects in networks may be dependent on the particular antigenic perturbation .
there has been a growth in the use of pathway - based approaches to gwas to shed light on biological processes and identify new candidate disease genes .
these approaches identify enriched pathways from a broader set of significant genes rather than focus only on a few of the most significant snps .
for example , pathway - based approaches and pathway clustering have identified important processes in breast cancer ( torkamani et al .
, 2008 ; menashe et al . , 2010 ) and in the diseases from the wtccc gwas ( torkamani et al . , 2008 )
we propose that such pathway - based approaches to gwas may benefit from the prioritization of genes based on the aggregation of gene  gene interactions as well as single - locus effects . a data - driven epistasis network approach
may increase the discovery of enriched disease - specific pathways and new candidate gene targets in gwas over single - locus prioritization alone .
gene interaction networks in this way may require some modifications to pathway analysis because enrichment scores typically rely on permutation to determine the statistical significance of enrichment , and data - driven networks are much more computationally intensive than single - locus calculations .
previous data - driven network approaches for gwas have used shannon information theory for epistasis calculations and network construction ( jakulin and bratko , 2004 ; moore et al . , 2006 ;
mckinney et al . , 2009 ; davis et al . , 2010 ) . however , casting the network in the widely used statistical framework of a general linear model ( glm ) has some advantages over information theory .
for example , use of a glm framework provides the flexibility to handle environmental covariates , longitudinal data , missing data , censoring , and cluster structure ( e.g. , family studies ) through the inclusion of appropriate random effects . as an example
, we present a likelihood ratio test of association between disease and a genetic locus , allowing for the possibility that the genetic effect may be modified by another genetic factor .
control data , but it is straightforward to develop similar tests for non - dichotomous phenotypes and other designs ( e.g. , family studies ) .
the coefficient bb gives the baseline risk of disease and coefficients b1 and b2 correct for main effects in the interaction regression model .
gene edge weights b12 in the gain matrix , we are interested in the b12 regression coefficients that are statistically different from zero .
the statistical framework also allows false discovery rate ( fdr ) procedures to be applied to correct for multiple gene 
fdr would be a more computationally efficient procedure than the permutation approached used in mckinney et al .
the diagonal element bii of the gain is simply the main effect regression coefficient without interactions .
additional terms may be added to the regression equation to define interactions between other factors , such as environment or gene expression , to create a heterogeneous network or to correct for these factors .
heterogeneous interaction networks represent an important frontier in genetics because they may improve our understanding of the interplay between genetic and environmental modifiers of susceptibility .
once an epistasis graph like figure 1 is calculated , one may identify the most influential or central snps or other factors in the network . such snp hubs in the disease - specific network may be potential targets for therapy or diagnostics , and the rankings may be used for a more sensitive pathway enrichment analysis .
various measures of node centrality have been proposed , with one of the most powerful and computationally tractable being eigenvector centrality  the basis of google s pagerank ( page et al . , 1999 ) .
eigenvector centrality calculates the steady state eigenvector of a markov transition matrix whose elements represent the probability to move from one node to another .
we recently designed a transition matrix method based on the main effects and interactions in gain and an eigenvector centrality algorithm we call snprank ( davis et al . , 2010 ) .
in a disease risk interaction network , the network is specific to the disease ( context sensitive ) , and each node and edge contains information about disease risk .
the diagonal of b captures main ( brando ) effects , while the off diagonals represent gene 
the factors dj ( the node degree ) and tr(b ) ( the trace of b ) normalize t to a stochastic matrix .
the elements of b could be calculated using approaches such as information theory or from coefficients of statistical models .
simple recursion is used to compute t s steady - state eigenvector , whose elements are the rankings of each node .
2 ) also includes self - interactions or main effects when i  =  j. in contrast , google s pagerank does not permit self - interactions because this would represent a website including links to itself and would artificially inflate the site s importance .
2 ) includes a parameter  , typically chosen close to 1 so that the main effect is not overwhelmed by the potentially large number of gene  gene interactions a snp may have .
we find  in the range [ 0.8 , 0.9 ] gives the highest internal consistency of the snprank scores , which we estimate by splitting data randomly into halves and calculating the kendall tau ( kendall , 1938 ) rank correlation of the snpranks .
we typically set  to 0.85 ; however ,   =  1 is also a reasonable simplification to ( eq .
the parameter  can also be used to blend prior knowledge from canonical pathways or protein  protein interactions . a helpful way to conceptualize the way the algorithm scores the importance of genes is to think of the snprank centrality algorithm as simulating ants that follow paths through the network that have the most disease risk information , either due to gene  gene interactions or main effects encoded in the edges and nodes of the network .
more pheromones are deposited at more frequently visited genetic or environmental nodes , and the amount of pheromone is proportional to the rank score of the node .
the steady state eigenvector , v ( eigenvalue   =  1 ) , of the snprank transition matrix is given by tv  =  v , where the elements of v are the snprank score of each snp .
the elements of v represent a probability field , and , thus , it is possible to use quantile plots and similar approaches to identify snps with snpranks that deviate from a uniform probability distribution .
one can also see from the eigenvalue equation that the snprank importance of the ith snp , vi , is influenced by the entire network : in other words , the snprank of the ith snp is a function of all interaction coefficients bij ( j    i ) , the main effect ( self interaction ) , bii , of the ith snp , and the recursively estimated importance of each snp vj ( j    i ) .
this approach does not include pure higher order interactions , like three - way etc .
but when ranking an individual snp , this approach models the interaction with every other snp in the network and the snp s main effect .
an important point for future research is the extent to which network approaches offer improved statistical power to detect complex genetic effects . in the case of pkns ,
power critically depends on the extent to which a particular dataset fits with the current state of biological knowledge .
similar to a correctly constructed informative prior distribution in bayesian modeling , pkns will have an advantage in situations where the data is congruent with prior expectations .
much of this is due to limitations in simulation methodology , as it is currently difficult to generate datasets with a complex and  realistic  latent structure . while some inroads have been made ( himmelstein et al . , 2011 ) , because approaches to characterize complex effects are in their infancy , it is also not clear what constitutes a realistic level of complexity . therefore , it is imperative that advancements in network methodologies be pursued in tandem with research into appropriate simulation frameworks .
we do note , however , that approaches like snprank operate on network representation typically defined in a pair - wise fashion .
therefore , they implicitly inherit some of the features of the approach used to construct this pair - wise representation , whether based on information theory , statistical modeling , or other approaches ( cordell , 2009 ) .
the power of network approaches thus could be improved by finding the best approach or combination of approaches for detecting pair - wise interactions ( ritchie et al . , 2001 ; fan et al . , 2011 ) .
in addition to this inheritance , snprank may also boost power to detect associations by aggregating numerous interactions and main effects ( selinger - leneman et al .
strategies for replication of a single snp are well accepted ( kraft et al . , 2009 ) , but the replication strategy for a network is less obvious because of the complex entanglement of all snps in the network .
one level of replication would be to test for replication of a gene set in an independent cohort based on gene set enrichment from the network prioritization .
ld is another area that requires further investigation because the common strategy of ld pruning could run the risk of excluding interactions .
finally , both of these issues simultaneously become pertinent within the context of comparisons across racial / ethnic populations ( teo et al . ,
2010 ) as differential patterns of ld could complicate the interpretation of edges connecting specific haplotype tagging variants .
previous data - driven network approaches for gwas have used shannon information theory for epistasis calculations and network construction ( jakulin and bratko , 2004 ; moore et al . , 2006 ;
however , casting the network in the widely used statistical framework of a general linear model ( glm ) has some advantages over information theory .
for example , use of a glm framework provides the flexibility to handle environmental covariates , longitudinal data , missing data , censoring , and cluster structure ( e.g. , family studies ) through the inclusion of appropriate random effects . as an example
, we present a likelihood ratio test of association between disease and a genetic locus , allowing for the possibility that the genetic effect may be modified by another genetic factor .
control data , but it is straightforward to develop similar tests for non - dichotomous phenotypes and other designs ( e.g. , family studies ) .
the coefficient bb gives the baseline risk of disease and coefficients b1 and b2 correct for main effects in the interaction regression model .
gene edge weights b12 in the gain matrix , we are interested in the b12 regression coefficients that are statistically different from zero .
the statistical framework also allows false discovery rate ( fdr ) procedures to be applied to correct for multiple gene  gene hypotheses .
fdr would be a more computationally efficient procedure than the permutation approached used in mckinney et al .
the diagonal element bii of the gain is simply the main effect regression coefficient without interactions .
additional terms may be added to the regression equation to define interactions between other factors , such as environment or gene expression , to create a heterogeneous network or to correct for these factors .
heterogeneous interaction networks represent an important frontier in genetics because they may improve our understanding of the interplay between genetic and environmental modifiers of susceptibility .
once an epistasis graph like figure 1 is calculated , one may identify the most influential or central snps or other factors in the network . such snp hubs in the disease - specific network may be potential targets for therapy or diagnostics , and the rankings may be used for a more sensitive pathway enrichment analysis .
various measures of node centrality have been proposed , with one of the most powerful and computationally tractable being eigenvector centrality  the basis of google s pagerank ( page et al . , 1999 ) .
eigenvector centrality calculates the steady state eigenvector of a markov transition matrix whose elements represent the probability to move from one node to another .
we recently designed a transition matrix method based on the main effects and interactions in gain and an eigenvector centrality algorithm we call snprank ( davis et al . , 2010 ) .
in a disease risk interaction network , the network is specific to the disease ( context sensitive ) , and each node and edge contains information about disease risk .
the diagonal of b captures main ( brando ) effects , while the off diagonals represent gene 
the factors dj ( the node degree ) and tr(b ) ( the trace of b ) normalize t to a stochastic matrix . the elements of b could be calculated using approaches such as information theory or from coefficients of statistical models .
simple recursion is used to compute t s steady - state eigenvector , whose elements are the rankings of each node .
2 ) also includes self - interactions or main effects when i  =  j. in contrast , google s pagerank does not permit self - interactions because this would represent a website including links to itself and would artificially inflate the site s importance .
2 ) includes a parameter  , typically chosen close to 1 so that the main effect is not overwhelmed by the potentially large number of gene  gene interactions a snp may have .
we find  in the range [ 0.8 , 0.9 ] gives the highest internal consistency of the snprank scores , which we estimate by splitting data randomly into halves and calculating the kendall tau ( kendall , 1938 ) rank correlation of the snpranks .
we typically set  to 0.85 ; however ,   =  1 is also a reasonable simplification to ( eq .
the parameter  can also be used to blend prior knowledge from canonical pathways or protein  protein interactions . a helpful way to conceptualize the way the algorithm scores the importance of genes is to think of the snprank centrality algorithm as simulating ants that follow paths through the network that have the most disease risk information , either due to gene  gene interactions or main effects encoded in the edges and nodes of the network .
more pheromones are deposited at more frequently visited genetic or environmental nodes , and the amount of pheromone is proportional to the rank score of the node .
the steady state eigenvector , v ( eigenvalue   =  1 ) , of the snprank transition matrix is given by tv  =  v , where the elements of v are the snprank score of each snp .
the elements of v represent a probability field , and , thus , it is possible to use quantile plots and similar approaches to identify snps with snpranks that deviate from a uniform probability distribution .
one can also see from the eigenvalue equation that the snprank importance of the ith snp , vi , is influenced by the entire network : in other words , the snprank of the ith snp is a function of all interaction coefficients bij ( j    i ) , the main effect ( self interaction ) , bii , of the ith snp , and the recursively estimated importance of each snp vj ( j    i ) .
this approach does not include pure higher order interactions , like three - way etc .
but when ranking an individual snp , this approach models the interaction with every other snp in the network and the snp s main effect .
an important point for future research is the extent to which network approaches offer improved statistical power to detect complex genetic effects . in the case of pkns ,
power critically depends on the extent to which a particular dataset fits with the current state of biological knowledge .
similar to a correctly constructed informative prior distribution in bayesian modeling , pkns will have an advantage in situations where the data is congruent with prior expectations .
much of this is due to limitations in simulation methodology , as it is currently difficult to generate datasets with a complex and  realistic  latent structure . while some inroads have been made ( himmelstein et al . , 2011 ) , because approaches to characterize complex effects are in their infancy , it is also not clear what constitutes a realistic level of complexity . therefore , it is imperative that advancements in network methodologies be pursued in tandem with research into appropriate simulation frameworks .
we do note , however , that approaches like snprank operate on network representation typically defined in a pair - wise fashion .
therefore , they implicitly inherit some of the features of the approach used to construct this pair - wise representation , whether based on information theory , statistical modeling , or other approaches ( cordell , 2009 ) .
the power of network approaches thus could be improved by finding the best approach or combination of approaches for detecting pair - wise interactions ( ritchie et al . , 2001 ; fan et al . ,
in addition to this inheritance , snprank may also boost power to detect associations by aggregating numerous interactions and main effects ( selinger - leneman et al . , 2003 ) .
strategies for replication of a single snp are well accepted ( kraft et al . , 2009 ) , but the replication strategy for a network is less obvious because of the complex entanglement of all snps in the network .
one level of replication would be to test for replication of a gene set in an independent cohort based on gene set enrichment from the network prioritization .
ld is another area that requires further investigation because the common strategy of ld pruning could run the risk of excluding interactions .
finally , both of these issues simultaneously become pertinent within the context of comparisons across racial / ethnic populations ( teo et al . ,
2010 ) as differential patterns of ld could complicate the interpretation of edges connecting specific haplotype tagging variants .
pathway and gene set enrichment approaches have demonstrated the utility of aggregating information from many moderate - sized single - locus effects .
while such approaches assume implicitly that the targeted genetic architecture reflects a complex interacting system , the prioritization of genes for enrichment typically relies on single - locus effects .
we propose network models of gwas data , built up from many single - locus and gene  gene interactions .
we anticipate systems level network approaches to gwas will reveal new mechanisms and improve our understanding of the complex relationship between genotypes and phenotypes .
however , there are numerous statistical and bioinformatics challenges that remain to be addressed to realize this systems level understanding .
the two network approaches for gwas discussed in this brief review  prior knowledge and gene 
pkns are able to leverage information from different scales , such as ppi or gene co - expression , though these data sources may lack specificity to the disease under study .
specificity could be achieved by calculating the interaction between genes conditional on the phenotype at hand .
thus , the best approach may be to integrate prior knowledge with epistasis networks , perhaps through a bayesian formalism . as the biological connections of pkns and canonical pathways continue to improve in specificity
, the integration of these networks with new interactions discovered in empirical networks will likely be a powerful combination .
previous data - driven epistasis networks have modeled interactions with shannon information theory measures because of the computational efficiency and power to detect interactions .
an advantage of regression - based empirical networks is the ability to incorporate variation from other data types ( gene expression , methylation , copy number ) , covariates , and environmental factors .
another open question is how much of the missing heritability is explained by network / epistatic effects .
answering such questions will require the development of statistical models that provide a coherent predictive mechanism based on a highly interacting network .
there is also a methodological need for tools to construct and understand features of empirical networks , such as subnetwork motifs , hubs , and node degree distributions .
one well - studied degree distribution is the power law , which corresponds to the notion of a small world network ( bassett and bullmore , 2006 ) .
it remains to be seen whether the degree distribution of epistasis networks exhibits a power law and how ld would affect the degree distribution .
a power law edge distribution of small world networks implies that most nodes have a small number of edges , but there are a few nodes with a large number of connections ( hubs ) .
the small world property allows one to traverse from one node to any other in relatively few steps and the network is robust to random attack or mutation .
the need for the application and development of algorithms for characterizing networks is not unique to genomics .
the importance of network concepts and algorithms is recognized throughout biology , notably in neuroscience where there is current interest in resting - state functional connectivity networks from fmri data ( braun et al . , 2011 ) .
thus , there is an opportunity to adapt methods and integrate data from other domains for genomic data .
both data from brain connectivity and gene network analysis for other genomic data seem to exhibit characteristics of small world networks ( barabasi and oltvai , 2004 ) , and it will be important to determine whether this is also reflected in the structure of gwas data .
we are in the early stages of understanding how the network of epistatic and main effects synthesizes with biological networks and pathways .
knowledge of hubs , centralities and other properties of disease risk epistasis networks provide a new path toward identifying critical nodes in the network that may act as therapeutic targets or disease risk predictors . and a data - driven interaction network paradigm of gwas and deeper sequencing may lead to new insights into the mechanisms of evolution and complexity .
the authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest .
 there is growing evidence that much more of the genome than previously thought is required to explain the heritability of complex phenotypes . 
 recent studies have demonstrated that numerous common variants from across the genome explain portions of genetic variability , spawning various avenues of research directed at explaining the remaining heritability . 
 this polygenic structure is also the motivation for the growing application of pathway and gene set enrichment techniques , which have yielded promising results . 
 these findings suggest that the coordination of genes in pathways that are known to occur at the gene regulatory level also can be detected at the population level . 
 although genes in these networks interact in complex ways , most population studies have focused on the additive contribution of common variants and the potential of rare variants to explain additional variation . in this brief review , 
 we discuss the potential to explain additional genetic variation through the agglomeration of multiple gene  
 gene interactions as well as main effects of common variants in terms of a network paradigm . just as is the case for single - locus contributions 
 , we expect each gene  gene interaction edge in the network to have a small effect , but these effects may be reinforced through hubs and other connectivity structures in the network . 
 we discuss some of the opportunities and challenges of network methods for analyzing genome - wide association studies ( gwas ) such as the study of hubs and motifs , and integrating other types of variation and environmental interactions . 
 such network approaches may unveil hidden variation in gwas , improve understanding of mechanisms of disease , and possibly fit into a network paradigm of evolutionary genetics . 