Learning Gene Networks under SNP Perturbations Using eQTL Datasets

The standard approach for identifying gene networks is based on experimental perturbations of gene regulatory systems such as gene knock-out experiments, followed by a genome-wide profiling of differential gene expressions. However, this approach is significantly limited in that it is not possible to perturb more than one or two genes simultaneously to discover complex gene interactions or to distinguish between direct and indirect downstream regulations of the differentially-expressed genes. As an alternative, genetical genomics study has been proposed to treat naturally-occurring genetic variants as potential perturbants of gene regulatory system and to recover gene networks via analysis of population gene-expression and genotype data. Despite many advantages of genetical genomics data analysis, the computational challenge that the effects of multifactorial genetic perturbations should be decoded simultaneously from data has prevented a widespread application of genetical genomics analysis. In this article, we propose a statistical framework for learning gene networks that overcomes the limitations of experimental perturbation methods and addresses the challenges of genetical genomics analysis. We introduce a new statistical model, called a sparse conditional Gaussian graphical model, and describe an efficient learning algorithm that simultaneously decodes the perturbations of gene regulatory system by a large number of SNPs to identify a gene network along with expression quantitative trait loci (eQTLs) that perturb this network. While our statistical model captures direct genetic perturbations of gene network, by performing inference on the probabilistic graphical model, we obtain detailed characterizations of how the direct SNP perturbation effects propagate through the gene network to perturb other genes indirectly. We demonstrate our statistical method using HapMap-simulated and yeast eQTL datasets. In particular, the yeast gene network identified computationally by our method under SNP perturbations is well supported by the results from experimental perturbation studies related to DNA replication stress response.     

Systems biology of virus-host signaling network interactions

Viruses have evolved to manipulate the host cell machinery for virus propagation, in part by interfering with the host cellular signaling network. Molecular studies of individual pathways have uncovered many viral host-protein targets; however, it is difficult to predict how viral perturbations will affect the signaling network as a whole. Systems biology approaches rely on multivariate, context-dependent measurements and computational analysis to elucidate how viral infection alters host cell signaling at a network level. Here we describe recent advances in systems analyses of signaling networks in both viral and non-viral biological contexts. These approaches have the potential to uncover virus- mediated changes to host signaling networks, suggest new therapeutic strategies, and assess how cell-to-cell variability affects host responses to infection. We argue that systems approaches will both improve understanding of how individual virus-host protein interactions fit into the progression of viral pathogenesis and help to identify novel therapeutic targets.     

Allele-Specific Behavior of Molecular Networks: Understanding Small-Molecule Drug Response in Yeast

The study of systems genetics is changing the way the genetic and molecular basis of phenotypic variation, such as disease susceptibility and drug response, is being analyzed. Moreover, systems genetics aids in the translation of insights from systems biology into genetics. The use of systems genetics enables greater attention to be focused on the potential impact of genetic perturbations on the molecular states of networks that in turn affects complex traits. In this study, we developed models to detect allele-specific perturbations on interactions, in which a genetic locus with alternative alleles exerted a differing influence on an interaction. We utilized the models to investigate the dynamic behavior of an integrated molecular network undergoing genetic perturbations in yeast. Our results revealed the complexity of regulatory relationships between genetic loci and networks, in which different genetic loci perturb specific network modules. In addition, significant within-module functional coherence was found. We then used the network perturbation model to elucidate the underlying molecular mechanisms of individual differences in response to 100 diverse small molecule drugs. As a result, we identified sub-networks in the integrated network that responded to variations in DNA associated with response to diverse compounds and were significantly enriched for known drug targets. Literature mining results provided strong independent evidence for the effectiveness of these genetic perturbing networks in the elucidation of small-molecule responses in yeast.     

Disentangling the Complexity of HGF Signaling by Combining Qualitative and Quantitative Modeling

Signaling pathways are characterized by crosstalk, feedback and feedforward mechanisms giving rise to highly complex and cell-context specific signaling networks. Dissecting the underlying relations is crucial to predict the impact of targeted perturbations. However, a major challenge in identifying cell-context specific signaling networks is the enormous number of potentially possible interactions. Here, we report a novel hybrid mathematical modeling strategy to systematically unravel hepatocyte growth factor (HGF) stimulated phosphoinositide-3-kinase (PI3K) and mitogen activated protein kinase (MAPK) signaling, which critically contribute to liver regeneration. By combining time-resolved quantitative experimental data generated in primary mouse hepatocytes with interaction graph and ordinary differential equation modeling, we identify and experimentally validate a network structure that represents the experimental data best and indicates specific crosstalk mechanisms. Whereas the identified network is robust against single perturbations, combinatorial inhibition strategies are predicted that result in strong reduction of Akt and ERK activation. Thus, by capitalizing on the advantages of the two modeling approaches, we reduce the high combinatorial complexity and identify cell-context specific signaling networks.     

Robust signal processing in living cells

Cellular signaling networks have evolved an astonishing ability to function reliably and with high fidelity in uncertain environments. A crucial prerequisite for the high precision exhibited by many signaling circuits is their ability to keep the concentrations of active signaling compounds within tightly defined bounds, despite strong stochastic fluctuations in copy numbers and other detrimental influences. Based on a simple mathematical formalism, we identify topological organizing principles that facilitate such robust control of intracellular concentrations in the face of multifarious perturbations. Our framework allows us to judge whether a multiple-input-multiple-output reaction network is robust against large perturbations of network parameters and enables the predictive design of perfectly robust synthetic network architectures. Utilizing the Escherichia coli chemotaxis pathway as a hallmark example, we provide experimental evidence that our framework indeed allows us to unravel the topological organization of robust signaling. We demonstrate that the specific organization of the pathway allows the system to maintain global concentration robustness of the diffusible response regulator CheY with respect to several dominant perturbations. Our framework provides a counterpoint to the hypothesis that cellular function relies on an extensive machinery to fine-tune or control intracellular parameters. Rather, we suggest that for a large class of perturbations, there exists an appropriate topology that renders the network output invariant to the respective perturbations.     

Inferring dynamic architecture of cellular networks using time series of gene expression, protein and metabolite data

MOTIVATION: High-throughput technologies have facilitated the acquisition of large genomics and proteomics datasets. However, these data provide snapshots of cellular behavior, rather than help us reveal causal relations. Here, we propose how these technologies can be utilized to infer the topology and strengths of connections among genes, proteins and metabolites by monitoring time-dependent responses of cellular networks to experimental interventions. RESULTS: We demonstrate that all connections leading to a given network node, e.g. to a particular gene, can be deduced from responses to perturbations none of which directly influences that node, e.g. using strains with knock-outs to other genes. To infer all interactions from stationary data, each node should be perturbed separately or in combination with other nodes. Monitoring time series provides richer information and does not require perturbations to all nodes. Overall, the methods we propose are capable of deducing and quantifying functional interactions within and across cellular gene, signaling and metabolic networks. SUPPLEMENTARY INFORMATION: Supplementary material is available at http://www.dbi.tju.edu/bioinformatics2004.pdf     

Coherent coupling of feedback loops: a design principle of cell signaling networks

MOTIVATION: It is widely accepted that cell signaling networks have been evolved to be robust against perturbations. To investigate the topological characteristics resulting in such robustness, we have examined large-scale signaling networks and found that a number of feedback loops are present mostly in coupled structures. In particular, the coupling was made in a coherent way implying that same types of feedback loops are interlinked together. RESULTS: We have investigated the role of such coherently coupled feedback loops through extensive Boolean network simulations and found that a high proportion of coherent couplings can enhance the robustness of a network against its state perturbations. Moreover, we found that the robustness achieved by coherently coupled feedback loops can be kept evolutionarily stable. All these results imply that the coherent coupling of feedback loops might be a design principle of cell signaling networks devised to achieve the robustness.     

Steady-State Differential Dose Response in Biological Systems

In pharmacology and systems biology, it is a fundamental problem to determine how biological systems change their dose-response behavior upon perturbations. In particular, it is unclear how topologies, reactions, and parameters (differentially) affect the dose response. Because parameters are often unknown, systematic approaches should directly relate network structure and function. Here, we outline a procedure to compare general non-monotone dose-response curves and subsequently develop a comprehensive theory for differential dose responses of biochemical networks captured by non-equilibrium steady-state linear framework models. Although these models are amenable to analytical derivations of non-equilibrium steady states in principle, their size frequently increases (super) exponentially with model size. We extract general principles of differential responses based on a model’s graph structure and thereby alleviate the combinatorial explosion. This allows us, for example, to determine reactions that affect differential responses, to identify classes of networks with equivalent differential, and to reject hypothetical models reliably without needing to know parameter values. We exemplify such applications for models of insulin signaling.     

Mutualistic networks: moving closer to a predictive theory

Plant–animal mutualistic networks sustain terrestrial biodiversity and human food security. Global environmental changes threaten these networks, underscoring the urgency for developing a predictive theory on how networks respond to perturbations. Here, I synthesise theoretical advances towards predicting network structure, dynamics, interaction strengths and responses to perturbations. I find that mathematical models incorporating biological mechanisms of mutualistic interactions provide better predictions of network dynamics. Those mechanisms include trait matching, adaptive foraging, and the dynamic consumption and production of both resources and services provided by mutualisms. Models incorporating species traits better predict the potential structure of networks (fundamental niche), while theory based on the dynamics of species abundances, rewards, foraging preferences and reproductive services can predict the extremely dynamic realised structures of networks, and may successfully predict network responses to perturbations. From a theoretician's standpoint, model development must more realistically represent empirical data on interaction strengths, population dynamics and how these vary with perturbations from global change. From an empiricist's standpoint, theory needs to make specific predictions that can be tested by observation or experiments. Developing models using short‐term empirical data allows models to make longer term predictions of community dynamics. As more longer term data become available, rigorous tests of model predictions will improve.     

The aging stress response

Aging is the outcome of a balance between damage and repair. The rate of aging and the appearance of age-related pathology are modulated by stress response and repair pathways that gradually decline, including the proteostasis and DNA damage repair networks and mitochondrial respiratory metabolism. Highly conserved insulin/IGF-1, TOR, and sirtuin signaling pathways in turn control these critical cellular responses. The coordinated action of these signaling pathways maintains cellular and organismal homeostasis in the face of external perturbations, such as changes in nutrient availability, temperature, and oxygen level, as well as internal perturbations, such as protein misfolding and DNA damage. Studies in model organisms suggest that changes in signaling can augment these critical stress response systems, increasing life span and reducing age-related pathology. The systems biology of stress response signaling thus provides a new approach to the understanding and potential treatment of age-related diseases.     

Inference of signaling and gene regulatory networks by steady-state perturbation experiments: structure and accuracy

One of the fundamental problems of cell biology is the understanding of complex regulatory networks. Such networks are ubiquitous in cells and knowledge of their properties is essential for the understanding of cellular behavior. In earlier work (Kholodenko et al. (PNAS 99: 12841), it was shown how the structure of biological networks can be quantified from experimental measurements of steady-state concentrations of key intermediates as a result of perturbations using a simple algorithm called "unravelling". Here, we study the effect of experimental uncertainty on the accuracy of the inferred structure (i.e. whether interactions are excitatory or inhibitory) of the networks determined using the unravelling algorithm. We show that the accuracy of the network structure depends not only on the noise level but on the strength of the interactions within the network. In particular, both very small and very large values of the connection strengths lead to large uncertainty in the inferred network. We describe a powerful geometric tool for the intuitive understanding of the effect of experimental error on the qualitative accuracy of the inferred network. In addition, we show that the use of additional data beyond that needed to minimally constrain the network not only improves the accuracy of the inferred network, but also may allow the detection of situations in which the initial assumptions of unravelling with respect to the network and the perturbations have been violated. Our ideas are illustrated using the mitogen-activated protein kinase (MAPK) signaling network as an example.