An official website of the United States government
Cellular fates are determined by genes interacting across large, complex biological networks. A critical question is how to identify causal relationships spanning distinct signaling pathways and underlying organismal phenotypes. Here, we address this question by constructing a Boolean model of a well-studied developmental network and analyzing information flows through the system. Depending on environmental signals Caenorhabditis elegans develop normally to sexual maturity or enter a reproductively delayed, developmentally quiescent ‘dauer’ state, progressing to maturity when the environment changes. The developmental network that starts with environmental signal and ends in the dauer/no dauer fate involves genes across 4 signaling pathways including cyclic GMP, Insulin/IGF-1, TGF- β and steroid hormone synthesis. We identified three stable motifs leading to normal development, each composed of genes interacting across the Insulin/IGF-1, TGF- β and steroid hormone synthesis pathways. Three genes known to influence dauer fate, daf-2 , daf-7 and hsf-1 , acted as driver nodes in the system. Using causal logic analysis, we identified a five gene cyclic subgraph integrating the information flow from environmental signal to dauer fate. Perturbation analysis showed that a multifactorial insulin profile determined the stable motifs the system entered and interacted with daf-12 as the switchpoint driving the dauer/no dauer fate. Our results show that complex organismal systems can be distilled into abstract representations that permit full characterization of the causal relationships driving developmental fates. Analyzing organismal systems from this perspective of logic and function has important implications for studies examining the evolution and conservation of signaling pathways.
more »
« less
Inference of a Boolean Network From Causal Logic Implications
Biological systems contain a large number of molecules that have diverse interactions. A fruitful path to understanding these systems is to represent them with interaction networks, and then describe flow processes in the network with a dynamic model. Boolean modeling, the simplest discrete dynamic modeling framework for biological networks, has proven its value in recapitulating experimental results and making predictions. A first step and major roadblock to the widespread use of Boolean networks in biology is the laborious network inference and construction process. Here we present a streamlined network inference method that combines the discovery of a parsimonious network structure and the identification of Boolean functions that determine the dynamics of the system. This inference method is based on a causal logic analysis method that associates a logic type (sufficient or necessary) to node-pair relationships (whether promoting or inhibitory). We use the causal logic framework to assimilate indirect information obtained from perturbation experiments and infer relationships that have not yet been documented experimentally. We apply this inference method to a well-studied process of hormone signaling in plants, the signaling underlying abscisic acid (ABA)—induced stomatal closure. Applying the causal logic inference method significantly reduces the manual work typically required for network and Boolean model construction. The inferred model agrees with the manually curated model. We also test this method by re-inferring a network representing epithelial to mesenchymal transition based on a subset of the information that was initially used to construct the model. We find that the inference method performs well for various likely scenarios of inference input information. We conclude that our method is an effective approach toward inference of biological networks and can become an efficient step in the iterative process between experiments and computations.
more »
« less
- Award ID(s):
- 1715826
- PAR ID:
- 10333721
- Date Published:
- Journal Name:
- Frontiers in Genetics
- Volume:
- 13
- ISSN:
- 1664-8021
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Signaling networks are at the heart of almost all biological processes. Most of these networks contain large number of components, and often either the connections between these components are not known or the rate equations that govern the dynamics of soluble signaling components are not quantified. This uncertainty in network topology and parameters can make it challenging to formulate detailed mathematical models. Boolean networks, in which all components are either on or off, have emerged as viable alternatives to detailed mathematical models that contain rate constants and other parameters. Therefore, open-source platforms of Boolean models for community use are desirable. Here, we present Boolink, a freely available graphical user interface that allows users to easily construct and analyze existing Boolean networks. Boolink can be applied to any Boolean network. We demonstrate its application using a previously published network for abscisic acid (ABA)-driven stomatal closure in Arabidopsis spp. (Arabidopsis thaliana). We also show how Boolink can be used to generate testable predictions by extending the network to include CO2 regulation of stomatal movements. Predictions of the model were experimentally tested, and the model was iteratively modified based on experiments showing that ABA effectively closes Arabidopsis stomata at near-zero CO2 concentrations (1.5-ppm CO2). Thus, Boolink enables public generation and the use of existing Boolean models, including the prior developed ABA signaling model with added CO2 signaling components.more » « less
-
We present new applications of parity inversion and time reversal to the emergence of complex behavior from simple dynamical rules in stochastic discrete models. Our parity-based encoding of causal relationships and time-reversal construction efficiently reveal discrete analogs of stable and unstable manifolds. We demonstrate their predictive power by studying decision-making in systems biology and statistical physics models. These applications underpin a novel attractor identification algorithm implemented for Boolean networks under stochastic dynamics. Its speed enables resolving a long-standing open question of how attractor count in critical random Boolean networks scales with network size and whether the scaling matches biological observations. Via 80-fold improvement in probed network size ( N = 16,384), we find the unexpectedly low scaling exponent of 0.12 ± 0.05, approximately one-tenth the analytical upper bound. We demonstrate a general principle: A system’s relationship to its time reversal and state-space inversion constrains its repertoire of emergent behaviors.more » « less
-
Coordinated responses to environmental stimuli are critical for multicellular organisms. To overcome the obstacles of cell-to-cell heterogeneity and noisy signaling dynamics within individual cells, cells must effectively exchange information with peers. However, the dynamics and mechanisms of collective information transfer driven by external signals are poorly understood. Here we investigate the calcium dynamics of neuronal cells that form confluent monolayers and respond to cyclic ATP stimuli in microfluidic devices. Using Granger inference to reconstruct the underlying causal relations between the cells, we find that the cells self-organize into spatially decentralized and temporally stationary networks to support information transfer via gap junction channels. The connectivity of the causal networks depends on the temporal profile of the external stimuli, where short periods, or long periods with small duty fractions, lead to reduced connectivity and fractured network topology. We build a theoretical model based on communicating excitable units that reproduces our observations. The model further predicts that connectivity of the causal network is maximal at an optimal communication strength, which is confirmed by the experiments. Together, our results show that information transfer between neuronal cells is externally regulated by the temporal profile of the stimuli and internally regulated by cell–cell communication.more » « less
-
Abstract Despite significant research, discovering causal relationships from fMRI remains a challenge. Popular methods such as Granger causality and dynamic causal modeling fall short in handling contemporaneous effects and latent common causes. Methods from causal structure learning literature can address these limitations but often scale poorly with network size and need acyclicity. In this study, we first provide a taxonomy of existing methods and compare their accuracy and efficiency on simulated fMRI from simple topologies. This analysis demonstrates a pressing need for more accurate and scalable methods, motivating the design of Causal discovery for Large-scale Low-resolution Time-series with Feedback (CaLLTiF). CaLLTiF is a constraint-based method that uses conditional independence between contemporaneous and lagged variables to extract causal relationships. On simulated fMRI from the macaque connectome, CaLLTiF achieves significantly higher accuracy and scalability than all tested alternatives. From resting-state human fMRI, CaLLTiF learns causal connectomes that are highly consistent across individuals, show clear top-down flow of causal effect from attention and default mode to sensorimotor networks, exhibit Euclidean distance dependence in causal interactions, and are highly dominated by contemporaneous effects. Overall, this work takes a major step in enhancing causal discovery from whole-brain fMRI and defines a new standard for future investigations.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3389/fgene.2022.836856
