CANTATA - prediction of missing links in Boolean networks using genetic programming

Motivation

Biological processes are complex systems with distinct behaviour. Despite the growing amount of available data, knowledge is sparse and often insufficient to investigate the complex regulatory behaviour of these systems. Moreover, different cellular phenotypes are possible under varying conditions. Mathematical models attempt to unravel these mechanisms by investigating the dynamics of regulatory networks. Therefore, a major challenge is to combine regulations and phenotypical information as well as the underlying mechanisms. To predict regulatory links in these models, we established an approach called CANTATA to support the integration of information into regulatory networks and retrieve potential underlying regulations. This is achieved by optimizing both static and dynamic properties of these networks.

Results

Initial results show that the algorithm predicts missing interactions by recapitulating the known phenotypes while preserving the original topology and optimizing the robustness of the model. The resulting models allow for hypothesizing about the biological impact of certain regulatory dependencies.

Availability and implementation

Source code of the application, example files and results are available at https://github.com/sysbio-bioinf/Cantata.

Supplementary information

Supplementary data are available at Bioinformatics online.

Robust network structure reconstruction based on Bayesian compressive sensing

Complex network has proven to be a general model to characterize interactions of practical complex systems. Recently, reconstructing the structure of complex networks with limited and noisy data attracts much research attention and has gradually become a hotspot. However, the collected data are often contaminated by unknown outliers inevitably, which seriously affects the accuracy of network reconstruction. Unfortunately, the existence of outliers is hard to be identified and always ignored in the network structure reconstruction task. To address this issue, here we propose a novel method which involves the outliers from the Bayesian perspective. The accuracy and the robustness of the proposed method have been verified in network reconstruction with payoff data contaminated with some outliers on both artificial networks and empirical networks. Extensive simulation results demonstrate the superiority of the proposed method. Thus, it can be concluded that since the proposed method can identify and get rid of outliers in observation data, it is conducive to improve the performance of network reconstruction.

Integrated network model provides new insights into castration-resistant prostate cancer

Castration-resistant prostate cancer (CRPC) is the main challenge for prostate cancer treatment. Recent studies have indicated that extending the treatments to simultaneously targeting different pathways could provide better approaches. To better understand the regulatory functions of different pathways, a system-wide study of CRPC regulation is necessary. For this purpose, we constructed a comprehensive CRPC regulatory network by integrating multiple pathways such as the MEK/ERK and the PI3K/AKT pathways. We studied the feedback loops of this network and found that AKT was involved in all detected negative feedback loops. We translated the network into a predictive Boolean model and analyzed the stable states and the control effects of genes using novel methods. We found that the stable states naturally divide into two obvious groups characterizing PC3 and DU145 cells respectively. Stable state analysis further revealed that several critical genes, such as PTEN, AKT, RAF, and CDKN2A, had distinct expression behaviors in different clusters. Our model predicted the control effects of many genes. We used several public datasets as well as FHL2 overexpression to verify our finding. The results of this study can help in identifying potential therapeutic targets, especially simultaneous targets of multiple pathways, for CRPC.

Boolean networks with multiexpressions and parameters

To model biological systems using networks, it is desirable to allow more than two levels of expression for the nodes and to allow the introduction of parameters. Various modeling and simulation methods addressing these needs using Boolean models, both synchronous and asynchronous, have been proposed in the literature. However, analytical study of these more general Boolean networks models is lagging. This paper aims to develop a concise theory for these different Boolean logic-based modeling methods. Boolean models for networks where each node can have more than two levels of expression and Boolean models with parameters are defined algebraically with examples provided. Certain classes of random asynchronous Boolean networks and deterministic moduli asynchronous Boolean networks are investigated in detail using the setting introduced in this paper. The derived theorems provide a clear picture for the attractor structures of these asynchronous Boolean networks.

Finding Missing Interactions of the Arabidopsis thaliana Root Stem Cell Niche Gene Regulatory Network

Over the last few decades, the Arabidopsis thaliana root stem cell niche (RSCN) has become a model system for the study of plant development and stem cell niche dynamics. Currently, many of the molecular mechanisms involved in RSCN maintenance and development have been described. A few years ago, we published a gene regulatory network (GRN) model integrating this information. This model suggested that there were missing components or interactions. Upon updating the model, the observed stable gene configurations of the RSCN could not be recovered, indicating that there are additional missing components or interactions in the model. In fact, due to the lack of experimental data, GRNs inferred from published data are usually incomplete. However, predicting the location and nature of the missing data is a not trivial task. Here, we propose a set of procedures for detecting and predicting missing interactions in Boolean networks. We used these procedures to predict putative missing interactions in the A. thaliana RSCN network model. Using our approach, we identified three necessary interactions to recover the reported gene activation configurations that have been experimentally uncovered for the different cell types within the RSCN: (1) a regulation of PHABULOSA to restrict its expression domain to the vascular cells, (2) a self-regulation of WOX5, possibly by an indirect mechanism through the auxin signaling pathway, and (3) a positive regulation of JACKDAW by MAGPIE. The procedures proposed here greatly reduce the number of possible Boolean functions that are biologically meaningful and experimentally testable and that do not contradict previous data. We believe that these procedures can be used on any Boolean network. However, because the procedures were designed for the specific case of the RSCN, formal demonstrations of the procedures should be shown in future efforts.

Optimization-based inference for temporally evolving networks with applications in biology

The problem of identifying dynamics of biological networks is of critical importance in order to understand biological systems. In this article, we propose a data-driven inference scheme to identify temporally evolving network representations of genetic networks. In the formulation of the optimization problem, we use an adjacency map as a priori information and define a cost function that both drives the connectivity of the graph to match biological data as well as generates a sparse and robust network at corresponding time intervals. Through simulation studies of simple examples, it is shown that this optimization scheme can help capture the topological change of a biological signaling pathway, and furthermore, might help to understand the structure and dynamics of biological genetic networks.

A systems biology perspective on plant-microbe interactions: biochemical and structural targets of pathogen effectors

Plants have biochemical defences against stresses from predators, parasites and pathogens. In this review we discuss the interaction of plant defences with microbial pathogens such as bacteria, fungi and oomycetes, and viruses. We examine principles of complex dynamic networks that allow identification of network components that are differentially and predictably sensitive to perturbation, thus making them likely effector targets. We relate these principles to recent developments in our understanding of known effector targets in plant-pathogen systems, and propose a systems-level framework for the interpretation and modelling of host-microbe interactions mediated by effectors. We describe this framework briefly, and conclude by discussing useful experimental approaches for populating this framework.