Dissecting genetic networks underlying complex phenotypes: the theoretical framework

Great progress has been made in genetic dissection of quantitative trait variation during the past two decades, but many studies still reveal only a small fraction of quantitative trait loci (QTLs), and epistasis remains elusive. We integrate contemporary knowledge of signal transduction pathways with principles of quantitative and population genetics to characterize genetic networks underlying complex traits, using a model founded upon one-way functional dependency of downstream genes on upstream regulators (the principle of hierarchy) and mutual functional dependency among related genes (functional genetic units, FGU). Both simulated and real data suggest that complementary epistasis contributes greatly to quantitative trait variation, and obscures the phenotypic effects of many 'downstream' loci in pathways. The mathematical relationships between the main effects and epistatic effects of genes acting at different levels of signaling pathways were established using the quantitative and population genetic parameters. Both loss of function and "co-adapted" gene complexes formed by multiple alleles with differentiated functions (effects) are predicted to be frequent types of allelic diversity at loci that contribute to the genetic variation of complex traits in populations. Downstream FGUs appear to be more vulnerable to loss of function than their upstream regulators, but this vulnerability is apparently compensated by different FGUs of similar functions. Other predictions from the model may account for puzzling results regarding responses to selection, genotype by environment interaction, and the genetic basis of heterosis.     

Evolution of pleiotropy: epistatic interaction pattern supports a mechanistic model underlying variation in genotype-phenotype map

The genotype-phenotype (GP) map consists of developmental and physiological mechanisms mapping genetic onto phenotypic variation. It determines the distribution of heritable phenotypic variance on which selection can act. Comparative studies of morphology as well as of gene regulatory networks show that the GP map itself evolves, yet little is known about the actual evolutionary mechanisms involved. The study of such mechanisms requires exploring the variation in GP maps at the population level, which presently is easier to quantify by statistical genetic methods rather than by regulatory network structures. We focus on the evolution of pleiotropy, a major structural aspect of the GP map. Pleiotropic genes affect multiple traits and underlie genetic covariance between traits, often causing evolutionary constraints. Previous quantitative genetic studies have demonstrated population-level variation in pleiotropy in the form of loci, at which genotypes differ in the genetic covariation between traits. This variation can potentially fuel evolution of the GP map under selection and/or drift. Here, we propose a developmental mechanism underlying population genetic variation in covariance and test its predictions. Specifically, the mechanism predicts that the loci identified as responsible for genetic variation in pleiotropy are involved in trait-specific epistatic interactions. We test this prediction for loci affecting allometric relationships between traits in an advanced intercross between inbred mouse strains. The results consistently support the prediction. We further find a high degree of sign epistasis in these interactions, which we interpret as an indication of adaptive gene complexes within the diverged parental lines.     

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

Epistasis plays an important role in the genetic architecture of common human diseases and can be viewed from two perspectives, biological and statistical, each derived from and leading to different assumptions and research strategies. Biological epistasis is the result of physical interactions among biomolecules within gene regulatory networks and biochemical pathways in an individual such that the effect of a gene on a phenotype is dependent on one or more other genes. In contrast, statistical epistasis is defined as deviation from additivity in a mathematical model summarizing the relationship between multilocus genotypes and phenotypic variation in a population. The goal of this essay is to review definitions and examples of biological and statistical epistasis and to explore the relationship between the two. Specifically, we present and discuss the following two questions in the context of human health and disease. First, when does statistical evidence of epistasis in human populations imply underlying biomolecular interactions in the etiology of disease? Second, when do biomolecular interactions produce patterns of statistical epistasis in human populations?Answers to these two reciprocal questions will provide an important framework for using genetic information to improve our ability to diagnose, prevent and treat common human diseases. We propose that systems biology will provide the necessary information for addressing these questions and that model systems such as bacteria, yeast and digital organisms will be a useful place to start. BioEssays 27:637–646, 2005. © 2005 Wiley Periodicals, Inc.     

Epistasis in a model of molecular signal transduction

Biological functions typically involve complex interacting molecular networks, with numerous feedback and regulation loops. How the properties of the system are affected when one, or several of its parts are modified is a question of fundamental interest, with numerous implications for the way we study and understand biological processes and treat diseases. This question can be rephrased in terms of relating genotypes to phenotypes: to what extent does the effect of a genetic variation at one locus depend on genetic variation at all other loci? Systematic quantitative measurements of epistasis – the deviation from additivity in the effect of alleles at different loci – on a given quantitative trait remain a major challenge. Here, we take a complementary approach of studying theoretically the effect of varying multiple parameters in a validated model of molecular signal transduction. To connect with the genotype/phenotype mapping we interpret parameters of the model as different loci with discrete choices of these parameters as alleles, which allows us to systematically examine the dependence of the signaling output – a quantitative trait – on the set of possible allelic combinations. We show quite generally that quantitative traits behave approximately additively (weak epistasis) when alleles correspond to small changes of parameters; epistasis appears as a result of large differences between alleles. When epistasis is relatively strong, it is concentrated in a sparse subset of loci and in low order (e.g. pair-wise) interactions. We find that focusing on interaction between loci that exhibit strong additive effects is an efficient way of identifying most of the epistasis. Our model study defines a theoretical framework for interpretation of experimental data and provides statistical predictions for the structure of genetic interaction expected for moderately complex biological circuits.     

A stochastic expectation and maximization algorithm for detecting quantitative trait-associated genes

MOTIVATION: Most biological traits may be correlated with the underlying gene expression patterns that are partially determined by DNA sequence variation. The correlations between gene expressions and quantitative traits are essential for understanding the functions of genes and dissecting gene regulatory networks. RESULTS: In the present study, we adopted a novel statistical method, called the stochastic expectation and maximization (SEM) algorithm, to analyze the associations between gene expression levels and quantitative trait values and identify genetic loci controlling the gene expression variations. In the first step, gene expression levels measured from microarray experiments were assigned to two different clusters based on the strengths of their association with the phenotypes of a quantitative trait under investigation. In the second step, genes associated with the trait were mapped to genetic loci of the genome. Because gene expressions are quantitative, the genetic loci controlling the expression traits are called expression quantitative trait loci. We applied the same SEM algorithm to a real dataset collected from a barley genetic experiment with both quantitative traits and gene expression traits. For the first time, we identified genes associated with eight agronomy traits of barley. These genes were then mapped to seven chromosomes of the barley genome. The SEM algorithm and the result of the barley data analysis are useful to scientists in the areas of bioinformatics and plant breeding. Availability and implementation: The R program for the SEM algorithm can be downloaded from our website: http://www.statgen.ucr.edu.     

A Developmental Systems Perspective on Epistasis: Computational Exploration of Mutational Interactions in Model Developmental Regulatory Networks

The way in which the information contained in genotypes is translated into complex phenotypic traits (i.e. embryonic expression patterns) depends on its decoding by a multilayered hierarchy of biomolecular systems (regulatory networks). Each layer of this hierarchy displays its own regulatory schemes (i.e. operational rules such as +/− feedback) and associated control parameters, resulting in characteristic variational constraints. This process can be conceptualized as a mapping issue, and in the context of highly-dimensional genotype-phenotype mappings (GPMs) epistatic events have been shown to be ubiquitous, manifested in non-linear correspondences between changes in the genotype and their phenotypic effects. In this study I concentrate on epistatic phenomena pervading levels of biological organization above the genetic material, more specifically the realm of molecular networks. At this level, systems approaches to studying GPMs are specially suitable to shed light on the mechanistic basis of epistatic phenomena. To this aim, I constructed and analyzed ensembles of highly-modular (fully interconnected) networks with distinctive topologies, each displaying dynamic behaviors that were categorized as either arbitrary or functional according to early patterning processes in the Drosophila embryo. Spatio-temporal expression trajectories in virtual syncytial embryos were simulated via reaction-diffusion models. My in silico mutational experiments show that: 1) the average fitness decay tendency to successively accumulated mutations in ensembles of functional networks indicates the prevalence of positive epistasis, whereas in ensembles of arbitrary networks negative epistasis is the dominant tendency; and 2) the evaluation of epistatic coefficients of diverse interaction orders indicates that, both positive and negative epistasis are more prevalent in functional networks than in arbitrary ones. Overall, I conclude that the phenotypic and fitness effects of multiple perturbations are strongly conditioned by both the regulatory architecture (i.e. pattern of coupled feedback structures) and the dynamic nature of the spatio-temporal expression trajectories displayed by the simulated networks.     

Population and quantitative genetics of regulatory networks

I evolved boolean regulatory networks in a computer simulation. I varied mutation, recombination, the size of the network, and the number of connections per node. I measured the performance of networks and the heritability and epistasis of genetic effects. Networks of intermediate connectivity performed best. The distinction between metabolic and quantitative genetic additivity explained some of the variation in performance. Metabolic additivity describes the interaction between changes in a single network, whereas quantitative genetic additivity measures the consistency of phenotypic effect caused by gene substitution in randomly chosen members of the population. I analysed metabolic additivity by the distribution of epistatic effects of pairs of mutations in individual networks. I measured quantitative genetic additivity by heritability. Highly connected networks had greater metabolic additivity for perturbations to individual networks, but had lower additivity when measured by the average effect of a gene substitution (heritability). The lower heritability of highly connected nets appeared to reduce the effectiveness of recombination in searching evolutionary space.     

Quantitative trait locus mapping using human pedigrees

In the past decade phenomenal progress has been made in molecular and statistical genetic methods for localizing quantitative trait loci. Because of these advances, we can anticipate a long period of active genetic research in which the genes influencing human quantitative variability will be mapped and their effects accurately evaluated. Here, we review the current state of the science in statistical genetic methods for quantitative trait linkage analysis. In particular, we detail a variance component-based framework for localizing quantitative trait loci and for accurately estimating their relative effect sizes. Attention is paid to the optimal design of human family studies for localizing genes of small to moderate effect. In addition, methods and strategies are described for dealing with the most important complications of quantitative variation, including the assessment of genotype x environment interaction and epistasis.     

Improved power of sib-pair linkage analysis using measures of complex trait dynamics

The influence of epistasis on a quantitative trait can reduce the power of linkage analysis to identify the underlying loci. In the present study, we simulated a complex trait derived from a dynamic one-locus gene expression system with epistasis arising from feedback regulation and tested the power of sib-pair linkage analysis methods for detecting the underlying quantitative trait locus (QTL). Using this simple genetic architecture, we demonstrate that the power of sib-pair linkage analysis can be greatly improved if measures of complex trait dynamics are considered.     

Quantitative Trait Loci Underlying Gene Product Variation: A Novel Perspective for Analyzing Regulation of Genome Expression

A methodology to dissect the genetic architecture of quantitative variation of numerous gene products simultaneously is proposed. For each individual of a segregating progeny, proteins extracted from a given organ are separated using two-dimensional electrophoresis, and their amounts are estimated with a computer-assisted system for spot quantification. Provided a complete genetic map is available, statistical procedures allow determination of the number, effects and chromosomal locations of factors controlling the amounts of individual proteins. This approach was applied to anonymous proteins of etiolated coleoptiles of maize, in an F(2) progeny between two distant lines. The genetic map included both restriction fragment length polymorphism and protein markers. Minimum estimates of one to five unlinked regulatory factors were found for 42 of the 72 proteins analyzed, with a large diversity of effects. Dominance and epistasis interactions were involved in the control of 38% and 14% of the 72 proteins, respectively. Such a methodology might help understanding the architecture of regulatory networks and the possible adaptive or phenotypic significance of the polymorphism of the genes involved.     

Identifying the molecular basis of QTLs: eQTLs add a new dimension

Natural genetic variation within plant species is at the core of plant science ranging from agriculture to evolution. Whereas much progress has been made in mapping quantitative trait loci (QTLs) controlling this natural variation, the elucidation of the underlying molecular mechanisms has remained a bottleneck. Recent systems biology tools have significantly shortened the time required to proceed from a mapped locus to testing of candidate genes. These tools enable research on natural variation to move from simple reductionistic studies focused on individual genes to integrative studies connecting molecular variation at multiple loci with physiological consequences. This review focuses on recent examples that demonstrate how expression QTL data can be used for gene discovery and exploited to untangle complex regulatory networks.     

Gene regulatory networks generating the phenomena of additivity, dominance and epistasis

We show how the phenomena of genetic dominance, overdominance, additivity, and epistasis are generic features of simple diploid gene regulatory networks. These regulatory network models are together sufficiently complex to catch most of the suggested molecular mechanisms responsible for generating dominant mutations. These include reduced gene dosage, expression or protein activity (haploinsufficiency), increased gene dosage, ectopic or temporarily altered mRNA expression, increased or constitutive protein activity, and dominant negative effects. As classical genetics regards the phenomenon of dominance to be generated by intralocus interactions, we have studied two one-locus models, one with a negative autoregulatory feedback loop, and one with a positive autoregulatory feedback loop. To include the phenomena of epistasis and downstream regulatory effects, a model of a three-locus signal transduction network is also analyzed. It is found that genetic dominance as well as overdominance may be an intra- as well as interlocus interaction phenomenon. In the latter case the dominance phenomenon is intimately connected to either feedback-mediated epistasis or downstream-mediated epistasis. It appears that in the intra- as well as the interlocus case there is considerable room for additive gene action, which may explain to some degree the predictive power of quantitative genetic theory, with its emphasis on this type of gene action. Furthermore, the results illuminate and reconcile the prevailing explanations of heterosis, and they support the old conjecture that the phenomenon of dominance may have an evolutionary explanation related to life history strategy.     

Genome-Wide Variation of Cytosine Modifications Between European and African Populations and the Implications for Complex Traits

Elucidating cytosine modification differences between human populations can enhance our understanding of ethnic specificity in complex traits. In this study, cytosine modification levels in 133 HapMap lymphoblastoid cell lines derived from individuals of European or African ancestry were profiled using the Illumina HumanMethylation450 BeadChip. Approximately 13% of the analyzed CpG sites showed differential modification between the two populations at a false discovery rate of 1%. The CpG sites with greater modification levels in European descent were enriched in the proximal regulatory regions, while those greater in African descent were biased toward gene bodies. More than half of the detected population-specific cytosine modifications could be explained primarily by local genetic variation. In addition, a substantial proportion of local modification quantitative trait loci exhibited population-specific effects, suggesting that genetic epistasis and/or genotype × environment interactions could be common. Distinct correlations were observed between gene expression levels and cytosine modifications in proximal regions and gene bodies, suggesting epigenetic regulation of interindividual expression variation. Furthermore, quantitative trait loci associated with population-specific modifications can be colocalized with expression quantitative trait loci and single nucleotide polymorphisms previously identified for complex traits with known racial disparities. Our findings revealed abundant population-specific cytosine modifications and the underlying genetic basis, as well as the relatively independent contribution of genetic and epigenetic variations to population differences in gene expression.     

The molecular basis of quantitative genetic variation in natural populations

DNA markers allow us to study quantitative trait loci (QTL) - the genes that control adaptation and quantitative variation. Experiments can map the genes responsible for quantitative variation and address the evolutionary and ecological significance of this variation. Recent studies suggest that major genes segregate within and among natural populations. It is now feasible to study the genes that cause morphological variation, life history trade-offs, heterosis and speciation. These methods can determine the role of epistasis and genotype-by-environment interaction in maintaining genetic variation. QTL mapping is an important tool used to address evolutionary and ecological questions of long-standing interest.     

The causes of epistasis in genetic networks

Epistasis refers to the nonadditive interactions between genes in determining phenotypes. Considerable efforts have shown that, even for a given organism, epistasis may vary both in intensity and sign. Recent comparative studies supported that the overall sign of epistasis switches from positive to negative as the complexity of an organism increases, and it has been hypothesized that this change shall be a consequence of the underlying gene network properties. Why should this be the case? What characteristics of genetic networks determine the sign of epistasis? Here we show, by evolving genetic networks that differ in their complexity and robustness against perturbations but that perform the same tasks, that robustness increased with complexity and that epistasis was positive for small nonrobust networks but negative for large robust ones. Our results indicate that robustness and negative epistasis emerge as a consequence of the existence of redundant elements in regulatory structures of genetic networks and that the correlation between complexity and epistasis is a byproduct of such redundancy, allowing for the decoupling of epistasis from the underlying network complexity.