Unsupervised logic-based mechanism inference for network-driven biological processes

Modern analytical techniques enable researchers to collect data about cellular states, before and after perturbations. These states can be characterized using analytical techniques, but the inference of regulatory interactions that explain and predict changes in these states remains a challenge. Here we present a generalizable, unsupervised approach to generate parameter-free, logic-based models of cellular processes, described by multiple discrete states. Our algorithm employs a Hamming-distance based approach to formulate, test, and identify optimized logic rules that link two states. Our approach comprises two steps. First, a model with no prior knowledge except for the mapping between initial and attractor states is built. We then employ biological constraints to improve model fidelity. Our algorithm automatically recovers the relevant dynamics for the explored models and recapitulates key aspects of the biochemical species concentration dynamics in the original model. We present the advantages and limitations of our work and discuss how our approach could be used to infer logic-based mechanisms of signaling, gene-regulatory, or other input-output processes describable by the Boolean formalism.     

Superstability of the yeast cell-cycle dynamics: ensuring causality in the presence of biochemical stochasticity

Gene regulatory dynamics are governed by molecular processes and therefore exhibits an inherent stochasticity. However, for the survival of an organism it is a strict necessity that this intrinsic noise does not prevent robust functioning of the system. It is still an open question how dynamical stability is achieved in biological systems despite the omnipresent fluctuations. In this paper we investigate the cell cycle of the budding yeast Saccharomyces cerevisiae as an example of a well-studied organism. We study a genetic network model of 11 genes that coordinate the cell-cycle dynamics using a modeling framework which generalizes the concept of discrete threshold dynamics. By allowing for fluctuations in the process times, we introduce noise into the model, accounting for the effects of biochemical stochasticity. We study the dynamical attractor of the cell cycle and find a remarkable robustness against fluctuations of this kind. We identify mechanisms that ensure reliability in spite of fluctuations: 'Catcher states' and persistence of activity levels contribute significantly to the stability of the yeast cell cycle despite the inherent stochasticity.     

A unifying framework for metapopulation dynamics

Many biologically important processes, such as genetic differentiation, the spread of disease, and population stability, are affected by the (natural or enforced) subdivision of populations into networks of smaller, partly isolated, subunits. Such "metapopulations" can have extremely complex dynamics. We present a new general model that uses only two functions to capture, at the metapopulation scale, the main behavior of metapopulations. We show how complex, structured metapopulation models can be translated into our generalized framework. The metapopulation dynamics arising from some important biological processes are illustrated: the rescue effect, the Allee effect, and what we term the "antirescue effect." The antirescue effect captures instances where high migration rates are deleterious to population persistence, a phenomenon that has been largely ignored in metapopulation conservation theory. Management regimes that ignore a significant antirescue effect will be inadequate and may actually increase extinction risk. Further, consequences of territoriality and conspecific attraction on metapopulation-level dynamics are investigated. The new, simplified framework can incorporate knowledge from epidemiology, genetics, and population biology in a phenomenological way. It opens up new possibilities to identify and analyze the factors that are important for the evolution and persistence of the many spatially subdivided species.     

Joining forces: feedback and integration in plant development

Feedback and integration of information are of paramount importance for the robust functioning and dynamics of biological systems. In plant developmental biology, experimentation is increasingly combined with computational modeling to obtain a better understanding of how such regulatory interactions shape the systems' behavior. Here we highlight experimental and modeling studies on feedback loops and integration mechanisms involved in plant development. These studies have substantially expanded our understanding of previously characterized gene regulatory networks (GRNs). In addition, they illustrate the pervasiveness of regulatory interactions between seemingly unrelated processes and levels of organization. Modelers in plant development will increasingly face the challenges of what level of detail, which processes and how many levels of organization to incorporate when trying to understand a particular process.     

A general modeling framework for describing spatially structured population dynamics

Variation in movement across time and space fundamentally shapes the abundance and distribution of populations. Although a variety of approaches model structured population dynamics, they are limited to specific types of spatially structured populations and lack a unifying framework. Here, we propose a unified network‐based framework sufficiently novel in its flexibility to capture a wide variety of spatiotemporal processes including metapopulations and a range of migratory patterns. It can accommodate different kinds of age structures, forms of population growth, dispersal, nomadism and migration, and alternative life‐history strategies. Our objective was to link three general elements common to all spatially structured populations (space, time and movement) under a single mathematical framework. To do this, we adopt a network modeling approach. The spatial structure of a population is represented by a weighted and directed network. Each node and each edge has a set of attributes which vary through time. The dynamics of our network‐based population is modeled with discrete time steps. Using both theoretical and real‐world examples, we show how common elements recur across species with disparate movement strategies and how they can be combined under a unified mathematical framework. We illustrate how metapopulations, various migratory patterns, and nomadism can be represented with this modeling approach. We also apply our network‐based framework to four organisms spanning a wide range of life histories, movement patterns, and carrying capacities. General computer code to implement our framework is provided, which can be applied to almost any spatially structured population. This framework contributes to our theoretical understanding of population dynamics and has practical management applications, including understanding the impact of perturbations on population size, distribution, and movement patterns. By working within a common framework, there is less chance that comparative analyses are colored by model details rather than general principles.     

Seasonal Invasion Dynamics in a Spatially Heterogeneous River with Fluctuating Flows

A key problem in environmental flow assessment is the explicit linking of the flow regime with ecological dynamics. We present a hybrid modeling approach to couple hydrodynamic and biological processes, focusing on the combined impact of spatial heterogeneity and temporal variability on population dynamics. Studying periodically alternating pool-riffle rivers that are subjected to seasonally varying flows, we obtain an invasion ratchet mechanism. We analyze the ratchet process for a caricature model and a hybrid physical–biological model. The water depth and current are derived from a hydrodynamic equation for variable stream bed water flows and these quantities feed into a reaction-diffusion-advection model that governs population dynamics of a river species. We establish the existence of spreading speeds and the invasion ratchet phenomenon, using a mixture of mathematical approximations and numerical computations. Finally, we illustrate the invasion ratchet phenomenon in a spatially two-dimensional hydraulic simulation model of a meandering river structure. Our hybrid modeling approach strengthens the ecological component of stream hydraulics and allows us to gain a mechanistic understanding as to how flow patterns affect population survival.     

Probabilistic discovery of overlapping cellular processes and their regulation.

In this paper, we explore modeling overlapping biological processes. We discuss a probabilistic model of overlapping biological processes, gene membership in those processes, and an addition to that model that identifies regulatory mechanisms controlling process activation. A key feature of our approach is that we allow genes to participate in multiple processes, thus providing a more biologically plausible model for the process of gene regulation. We present algorithms to learn each model automatically from data, using only genomewide measurements of gene expression as input. We compare our results to those obtained by other approaches and show that significant benefits can be gained by modeling both the organization of genes into overlapping cellular processes and the regulatory programs of these processes. Moreover, our method successfully grouped genes known to function together, recovered many regulatory relationships that are known in the literature, and suggested novel hypotheses regarding the regulatory role of previously uncharacterized proteins.     

Xenopus as a model organism in developmental chemical genetic screens

Chemical genetics is a potentially powerful tool for studying developmental processes in vertebrate systems. We present data showing Xenopus laevis as a model organism in which systematic chemical genetic screens can be carried out. Previous forward chemical genetic screens, including those with developing zebrafish embryos, have demonstrated the nature and value of biological information gained with this approach. We show how amenable Xenopus is to chemical genetics by investigating a series of compounds either with known biochemical effects, or previously identified to give developmental phenotypes, on a range of biological functions, including the development of pigmentation, the heart and the central nervous system in zebrafish. We have found that the compounds give comparable phenotypes when applied to developing Xenopus embryos. We have also studied the penetrance and expressivity of these chemical genetic phenotypes in relation to genetic variation and the developmental window during which the compound is present. Finally, we assess the feasibility and the potential throughput of a screen in this vertebrate species.     

Network-based auto-probit modeling for protein function prediction

Predicting the functional roles of proteins based on various genome-wide data, such as protein-protein association networks, has become a canonical problem in computational biology. Approaching this task as a binary classification problem, we develop a network-based extension of the spatial auto-probit model. In particular, we develop a hierarchical Bayesian probit-based framework for modeling binary network-indexed processes, with a latent multivariate conditional autoregressive Gaussian process. The latter allows for the easy incorporation of protein-protein association network topologies-either binary or weighted-in modeling protein functional similarity. We use this framework to predict protein functions, for functions defined as terms in the Gene Ontology (GO) database, a popular rigorous vocabulary for biological functionality. Furthermore, we show how a natural extension of this framework can be used to model and correct for the high percentage of false negative labels in training data derived from GO, a serious shortcoming endemic to biological databases of this type. Our method performance is evaluated and compared with standard algorithms on weighted yeast protein-protein association networks, extracted from a recently developed integrative database called Search Tool for the Retrieval of INteracting Genes/proteins (STRING). Results show that our basic method is competitive with these other methods, and that the extended method-incorporating the uncertainty in negative labels among the training data-can yield nontrivial improvements in predictive accuracy.     

Anticipatory Functions, Digital-Analog Forms and Biosemiotics: Integrating the Tools to Model Information and Normativity in Autonomous Biological Agents

We argue that living systems process information such that functionality emerges in them on a continuous basis. We then provide a framework that can explain and model the normativity of biological functionality. In addition we offer an explanation of the anticipatory nature of functionality within our overall approach. We adopt a Peircean approach to Biosemiotics, and a dynamical approach to Digital-Analog relations and to the interplay between different levels of functionality in autonomous systems, taking an integrative approach. We then apply the underlying biosemiotic logic to a particular biological system, giving a model of the B-Cell Receptor signaling system, in order to demonstrate how biosemiotic concepts can be used to build an account of biological information and functionality. Next we show how this framework can be used to explain and model more complex aspects of biological normativity, for example, how cross-talk between different signaling pathways can be avoided. Overall, we describe an integrated theoretical framework for the emergence of normative functions and, consequently, for the way information is transduced across several interconnected organizational levels in an autonomous system, and we demonstrate how this can be applied in real biological phenomena. Our aim is to open the way towards realistic tools for the modeling of information and normativity in autonomous biological agents.     

The role of mitochondria in the life of the nematode,Caenorhabditis elegans

Mitochondria are essential organelles involved in energy metabolism via oxidative phosphorylation. They play a vital role in diverse biological processes such as aging and apoptosis. In humans, defects in the mitochondrial respiratory chain (MRC) are responsible for or associated with a bewildering variety of diseases. The nematode Caenorhabditis elegans is a simple animal and a powerful genetic and developmental model system. In this review, we discuss how the nematode model system has contributed to our understanding of mitochondrial dynamics, of the genetics and inheritance of the mitochondrial genome, and of the consequences of nuclear and mitochondrial DNA (mtDNA) mutations. Mitochondrial respiration is vital to energy metabolism but also to other aspects of multicellular life such as aging and development. We anticipate that further significant contributions to our understanding of mitochondrial function in animal biology are forthcoming with the C. elegans model system.     

A dynamic model for functional mapping of biological rhythms

Functional mapping is a statistical method for mapping quantitative trait loci (QTLs) that regulate the dynamic pattern of a biological trait. This method integrates mathematical aspects of biological complexity into a mixture model for genetic mapping and tests the genetic effects of QTLs by comparing genotype-specific curve parameters. As a way of quantitatively specifying the dynamic behavior of a system, differential equations have proven to be powerful for modeling and unraveling the biochemical, molecular, and cellular mechanisms of a biological process, such as biological rhythms. The equipment of functional mapping with biologically meaningful differential equations provides new insights into the genetic control of any dynamic processes. We formulate a new functional mapping framework for a dynamic biological rhythm by incorporating a group of ordinary differential equations (ODE). The Runge-Kutta fourth order algorithm was implemented to estimate the parameters that define the system of ODE. The new model will find its implications for understanding the interplay between gene interactions and developmental pathways in complex biological rhythms.     

Biological applications of the theory of birth-and-death processes

In this review, we discuss applications of the theory of birth-and-death processes to problems in biology, primarily, those of evolutionary genomics. The mathematical principles of the theory of these processes are briefly described. Birth-and-death processes, with some straightforward additions such as innovation, are a simple, natural and formal framework for modeling a vast variety of biological processes such as population dynamics, speciation, genome evolution, including growth of paralogous gene families and horizontal gene transfer and somatic evolution of cancers. We further describe how empirical data, e.g. distributions of paralogous gene family size, can be used to choose the model that best reflects the actual course of evolution among different versions of birth-death-and-innovation models. We conclude that birth-and-death processes, thanks to their mathematical transparency, flexibility and relevance to fundamental biological processes, are going to be an indispensable mathematical tool for the burgeoning field of systems biology.     

Natural selection governs local, but not global, evolutionary gene coexpression networks in Caenorhabditis elegans

Background

Signal duration (e.g. the time over which an active signaling intermediate persists) is a key regulator of biological decisions in myriad contexts such as cell growth, proliferation, and developmental lineage commitments. Accompanying differences in signal duration are numerous downstream biological processes that require multiple steps of biochemical regulation.

Results

Here we present an analysis that investigates how simple biochemical motifs that involve multiple stages of regulation can be constructed to differentially process signals that persist at different time scales. We compute the dynamic, frequency dependent gain within these networks and resulting power spectra to better understand how biochemical networks can integrate signals at different time scales. We identify topological features of these networks that allow for different frequency dependent signal processing properties.

Conclusion

We show that multi-staged cascades are effective in integrating signals of long duration whereas multi-staged cascades that operate in the presence of negative feedback are effective in integrating signals of short duration. Our studies suggest principles for why signal duration in connection with multiple steps of downstream regulation is a ubiquitous motif in biochemical systems.     

A computer-based model for the regulation of mitogen activated protein kinase (MAPK) activation

Computer simulations and mathematical modeling of biological processes are becoming increasingly popular, and yet the complexity of the biochemical systems or the differences between experimental setups make it very difficult to establish a standard formula for these modeling projects. Before we can start using computer-based models for predictions or targeted experiment designs, it is very important to establish a reliable model on which those predictions can be based and experimentally tested. Here we attempt to present a computer model for the mitogen-activated protein kinase (MAPK) signaling cascade which is consistent with previously published experimental results. In this study we have focused our attention to a generic MAPK ERK (extracellular signal-regulated kinase) pathway activated by epidermal growth factor (EGF) in an attempt to understand how receptors may achieve different activation kinetics of the MAPK signaling. We successfully show that the level of receptor expression is one key determinant in this regulation, and that the binding affinity of the active receptor to adaptor proteins can have a small but albeit direct effect on the downstream activation.     

A scenario-based approach to modeling development: a prototype model of C. elegans vulval fate specification

Studies of developmental biology are often facilitated by diagram “models” that summarize the current understanding of underlying mechanisms. The increasing complexity of our understanding of development necessitates computational models that can extend these representations to include their dynamic behavior. Here we present a prototype model of Caenorhabditis elegans vulval precursor cell fate specification that represents many processes crucial for this developmental event but that are hard to integrate using other modeling methodologies. We demonstrate the integrative capabilities of our methodology by comprehensively incorporating the contents of three seminal papers, showing that this methodology can lead to comprehensive models of developmental biology. The prototype computational model was built and is run using a language (Live Sequence Charts) and tool (the Play-Engine) that facilitate the same conceptual processes biologists use to construct and probe diagram-type models. We demonstrate that this modeling approach permits rigorous tests of mutual consistency between experimental data and mechanistic hypotheses and can identify specific conflicting results, providing a useful approach to probe developmental systems.     

A process‐based metacommunity framework linking local and regional scale community ecology

The metacommunity concept has the potential to integrate local and regional dynamics within a general community ecology framework. To this end, the concept must move beyond the discrete archetypes that have largely defined it (e.g. neutral vs. species sorting) and better incorporate local scale species interactions and coexistence mechanisms. Here, we present a fundamental reconception of the framework that explicitly links local coexistence theory to the spatial processes inherent to metacommunity theory, allowing for a continuous range of competitive community dynamics. These dynamics emerge from the three underlying processes that shape ecological communities: (1) density‐independent responses to abiotic conditions, (2) density‐dependent biotic interactions and (3) dispersal. Stochasticity is incorporated in the demographic realisation of each of these processes. We formalise this framework using a simulation model that explores a wide range of competitive metacommunity dynamics by varying the strength of the underlying processes. Using this model and framework, we show how existing theories, including the traditional metacommunity archetypes, are linked by this common set of processes. We then use the model to generate new hypotheses about how the three processes combine to interactively shape diversity, functioning and stability within metacommunities.     

Functional mapping of ontogeny in flowering plants

All organisms face the problem of how to perform a sequence of developmental changes and transitions during ontogeny. We revise functional mapping, a statistical model originally derived to map genes that determine developmental dynamics, to take into account the entire process of ontogenetic growth from embryo to adult and from the vegetative to reproductive phase. The revised model provides a framework that reconciles the genetic architecture of development at different stages and elucidates a comprehensive picture of the genetic control mechanisms of growth that change gradually from a simple to a more complex level. We use an annual flowering plant, as an example, to demonstrate our model by which to map genes and their interactions involved in embryo and postembryonic growth. The model provides a useful tool to study the genetic control of ontogenetic growth in flowering plants and any other organisms through proper modifications based on their biological characteristics.     

Evolution in stage-structured populations

For many organisms, stage is a better predictor of demographic rates than age. Yet no general theoretical framework exists for understanding or predicting evolution in stage-structured populations. Here, we provide a general modeling approach that can be used to predict evolution and demography of stage-structured populations. This advances our ability to understand evolution in stage-structured populations to a level previously available only for populations structured by age. We use this framework to provide the first rigorous proof that Lande's theorem, which relates adaptive evolution to population growth, applies to stage-classified populations, assuming only normality and that evolution is slow relative to population dynamics. We extend this theorem to allow for different means or variances among stages. Our next major result is the formulation of Price's theorem, a fundamental law of evolution, for stage-structured populations. In addition, we use data from Trillium grandiflorum to demonstrate how our models can be applied to a real-world population and thereby show their practical potential to generate accurate projections of evolutionary and population dynamics. Finally, we use our framework to compare rates of evolution in age- versus stage-structured populations, which shows how our methods can yield biological insights about evolution in stage-structured populations.