The extended TILAR approach: a novel tool for dynamic modeling of the transcription factor network regulating the adaption to in vitro cultivation of murine hepatocytes

Creation of Three-Dimensional Liver Tissue Models from Experimental Images for Systems Medicine

In this chapter, we illustrate how three-dimensional liver tissue models can be created from experimental image modalities by utilizing a well-established processing chain of experiments, microscopic imaging, image processing, image analysis and model construction. We describe how key features of liver tissue architecture are quantified and translated into model parameterizations, and show how a systematic iteration of experiments and model simulations often leads to a better understanding of biological phenomena in systems biology and systems medicine.

LASSIM-A network inference toolbox for genome-wide mechanistic modeling

Recent technological advancements have made time-resolved, quantitative, multi-omics data available for many model systems, which could be integrated for systems pharmacokinetic use. Here, we present large-scale simulation modeling (LASSIM), which is a novel mathematical tool for performing large-scale inference using mechanistically defined ordinary differential equations (ODE) for gene regulatory networks (GRNs). LASSIM integrates structural knowledge about regulatory interactions and non-linear equations with multiple steady state and dynamic response expression datasets. The rationale behind LASSIM is that biological GRNs can be simplified using a limited subset of core genes that are assumed to regulate all other gene transcription events in the network. The LASSIM method is implemented as a general-purpose toolbox using the PyGMO Python package to make the most of multicore computers and high performance clusters, and is available at https://gitlab.com/Gustafsson-lab/lassim. As a method, LASSIM works in two steps, where it first infers a non-linear ODE system of the pre-specified core gene expression. Second, LASSIM in parallel optimizes the parameters that model the regulation of peripheral genes by core system genes. We showed the usefulness of this method by applying LASSIM to infer a large-scale non-linear model of naïve Th2 cell differentiation, made possible by integrating Th2 specific bindings, time-series together with six public and six novel siRNA-mediated knock-down experiments. ChIP-seq showed significant overlap for all tested transcription factors. Next, we performed novel time-series measurements of total T-cells during differentiation towards Th2 and verified that our LASSIM model could monitor those data significantly better than comparable models that used the same Th2 bindings. In summary, the LASSIM toolbox opens the door to a new type of model-based data analysis that combines the strengths of reliable mechanistic models with truly systems-level data. We demonstrate the power of this approach by inferring a mechanistically motivated, genome-wide model of the Th2 transcription regulatory system, which plays an important role in several immune related diseases.

There are excellent methods to mathematically model time-resolved biological data on a small scale using accurate mechanistic models. Despite the rapidly increasing availability of such data, mechanistic models have not been applied on a genome-wide level due to excessive runtimes and the non-identifiability of model parameters. However, genome-wide, mechanistic models could potentially answer key clinical questions, such as finding the best drug combinations to induce an expression change from a disease to a healthy state. We present LASSIM, which is a toolbox built to infer parameters within mechanistic models on a genomic scale. This is made possible due to a property shared across biological systems, namely the existence of a subset of master regulators, here denoted the core system. The introduction of a core system of genes simplifies the network inference into small solvable sub-problems, and implies that all main regulatory actions on peripheral genes come from a small set of regulator genes. This separation allows substantial parts of computations to be solved in parallel, i.e. permitting the use of a computer cluster, which substantially reduces computation time.

Dual-species transcriptional profiling during systemic candidiasis reveals organ-specific host-pathogen interactions

Candida albicans is a common cause of life-threatening fungal bloodstream infections. In the murine model of systemic candidiasis, the kidney is the primary target organ while the fungal load declines over time in liver and spleen. To better understand these organ-specific differences in host-pathogen interaction, we performed gene expression profiling of murine kidney, liver and spleen and determined the fungal transcriptome in liver and kidney. We observed a delayed transcriptional immune response accompanied by late induction of fungal stress response genes in the kidneys. In contrast, early upregulation of the proinflammatory response in the liver was associated with a fungal transcriptome resembling response to phagocytosis, suggesting that phagocytes contribute significantly to fungal control in the liver. Notably, C. albicans hypha-associated genes were upregulated in the absence of visible filamentation in the liver, indicating an uncoupling of gene expression and morphology and a morphology-independent effect by hypha-associated genes in this organ. Consistently, integration of host and pathogen transcriptional data in an inter-species gene regulatory network indicated connections of C. albicans cell wall remodelling and metabolism to the organ-specific immune responses.

Use of systems biology to decipher host-pathogen interaction networks and predict biomarkers

In systems biology, researchers aim to understand complex biological systems as a whole, which is often achieved by mathematical modelling and the analyses of high-throughput data. In this review, we give an overview of medical applications of systems biology approaches with special focus on host-pathogen interactions. After introducing general ideas of systems biology, we focus on (1) the detection of putative biomarkers for improved diagnosis and support of therapeutic decisions, (2) network modelling for the identification of regulatory interactions between cellular molecules to reveal putative drug targets and (3) module discovery for the detection of phenotype-specific modules in molecular interaction networks. Biomarker detection applies supervised machine learning methods utilizing high-throughput data (e.g. single nucleotide polymorphism (SNP) detection, RNA-seq, proteomics) and clinical data. We demonstrate structural analysis of molecular networks, especially by identification of disease modules as a novel strategy, and discuss possible applications to host-pathogen interactions. Pioneering work was done to predict molecular host-pathogen interactions networks based on dual RNA-seq data. However, currently this network modelling is restricted to a small number of genes. With increasing number and quality of databases and data repositories, the prediction of large-scale networks will also be feasible that can used for multidimensional diagnosis and decision support for prevention and therapy of diseases. Finally, we outline further perspective issues such as support of personalized medicine with high-throughput data and generation of multiscale host-pathogen interaction models.

CyNetworkBMA: a Cytoscape app for inferring gene regulatory networks

Background

Inference of gene networks from expression data is an important problem in computational biology. Many algorithms have been proposed for solving the problem efficiently. However, many of the available implementations are programming libraries that require users to write code, which limits their accessibility.

Results

We have developed a tool called CyNetworkBMA for inferring gene networks from expression data that integrates with Cytoscape. Our application offers a graphical user interface for networkBMA, an efficient implementation of Bayesian Model Averaging methods for network construction. The client-server architecture of CyNetworkBMA makes it possible to distribute or centralize computation depending on user needs.

Conclusions

CyNetworkBMA is an easy-to-use tool that makes network inference accessible to non-programmers through seamless integration with Cytoscape. CyNetworkBMA is available on the Cytoscape App Store at http://apps.cytoscape.org/apps/cynetworkbma.

Electronic supplementary material

The online version of this article (doi:10.1186/s13029-015-0043-5) contains supplementary material, which is available to authorized users.

Network Modeling Reveals Cross Talk of MAP Kinases during Adaptation to Caspofungin Stress in Aspergillus fumigatus

Mitogen activated protein kinases (MAPKs) are highly conserved in eukaryotic organisms. In pathogenic fungi, their activities were assigned to different physiological functions including drug adaptation and resistance. Aspergillus fumigatus is a human pathogenic fungus, which causes life-threatening invasive infections. Therapeutic options against invasive mycoses are still limited. One of the clinically used drugs is caspofungin, which specifically targets the fungal cell wall biosynthesis. A systems biology approach, based on comprehensive transcriptome data sets and mathematical modeling, was employed to infer a regulatory network and identify key interactions during adaptation to caspofungin stress in A. fumigatus. Mathematical modeling and experimental validations confirmed an intimate cross talk occurring between the cell wall-integrity and the high osmolarity-glycerol signaling pathways. Specifically, increased concentrations of caspofungin promoted activation of these signalings. Moreover, caspofungin affected the intracellular transport, which caused an additional osmotic stress that is independent of glucan inhibition. High concentrations of caspofungin reduced this osmotic stress, and thus decreased its toxic activity. Our results demonstrated that MAPK signaling pathways play a key role during caspofungin adaptation and are contributing to the paradoxical effect exerted by this drug.

Data- and knowledge-based modeling of gene regulatory networks: an update

Gene regulatory network inference is a systems biology approach which predicts interactions between genes with the help of high-throughput data. In this review, we present current and updated network inference methods focusing on novel techniques for data acquisition, network inference assessment, network inference for interacting species and the integration of prior knowledge. After the advance of Next-Generation-Sequencing of cDNAs derived from RNA samples (RNA-Seq) we discuss in detail its application to network inference. Furthermore, we present progress for large-scale or even full-genomic network inference as well as for small-scale condensed network inference and review advances in the evaluation of network inference methods by crowdsourcing. Finally, we reflect the current availability of data and prior knowledge sources and give an outlook for the inference of gene regulatory networks that reflect interacting species, in particular pathogen-host interactions.

Computational prediction of molecular pathogen-host interactions based on dual transcriptome data

Inference of inter-species gene regulatory networks based on gene expression data is an important computational method to predict pathogen-host interactions (PHIs). Both the experimental setup and the nature of PHIs exhibit certain characteristics. First, besides an environmental change, the battle between pathogen and host leads to a constantly changing environment and thus complex gene expression patterns. Second, there might be a delay until one of the organisms reacts. Third, toward later time points only one organism may survive leading to missing gene expression data of the other organism. Here, we account for PHI characteristics by extending NetGenerator, a network inference tool that predicts gene regulatory networks from gene expression time series data. We tested multiple modeling scenarios regarding the stimuli functions of the interaction network based on a benchmark example. We show that modeling perturbation of a PHI network by multiple stimuli better represents the underlying biological phenomena. Furthermore, we utilized the benchmark example to test the influence of missing data points on the inference performance. Our results suggest that PHI network inference with missing data is possible, but we recommend to provide complete time series data. Finally, we extended the NetGenerator tool to incorporate gene- and time point specific variances, because complex PHIs may lead to high variance in expression data. Sample variances are directly considered in the objective function of NetGenerator and indirectly by testing the robustness of interactions based on variance dependent disturbance of gene expression values. We evaluated the method of variance incorporation on dual RNA sequencing (RNA-Seq) data of Mus musculus dendritic cells incubated with Candida albicans and proofed our method by predicting previously verified PHIs as robust interactions.

Dynamic modelling of microRNA regulation during mesenchymal stem cell differentiation

Background

Network inference from gene expression data is a typical approach to reconstruct gene regulatory networks. During chondrogenic differentiation of human mesenchymal stem cells (hMSCs), a complex transcriptional network is active and regulates the temporal differentiation progress. As modulators of transcriptional regulation, microRNAs (miRNAs) play a critical role in stem cell differentiation. Integrated network inference aimes at determining interrelations between miRNAs and mRNAs on the basis of expression data as well as miRNA target predictions. We applied the NetGenerator tool in order to infer an integrated gene regulatory network.

Results

Time series experiments were performed to measure mRNA and miRNA abundances of TGF-beta1+BMP2 stimulated hMSCs. Network nodes were identified by analysing temporal expression changes, miRNA target gene predictions, time series correlation and literature knowledge. Network inference was performed using NetGenerator to reconstruct a dynamical regulatory model based on the measured data and prior knowledge. The resulting model is robust against noise and shows an optimal trade-off between fitting precision and inclusion of prior knowledge. It predicts the influence of miRNAs on the expression of chondrogenic marker genes and therefore proposes novel regulatory relations in differentiation control. By analysing the inferred network, we identified a previously unknown regulatory effect of miR-524-5p on the expression of the transcription factor SOX9 and the chondrogenic marker genes COL2A1, ACAN and COL10A1.

Conclusions

Genome-wide exploration of miRNA-mRNA regulatory relationships is a reasonable approach to identify miRNAs which have so far not been associated with the investigated differentiation process. The NetGenerator tool is able to identify valid gene regulatory networks on the basis of miRNA and mRNA time series data.

The use of transcriptomics to unveil the role of nutrients in Mammalian liver

Liver is the organ primarily responding to diet, and it is crucial in determining plasma carbohydrate, protein, and lipid levels. In addition, it is mainly responsible for transformation of xenobiotics. For these reasons, it has been a target of transcriptomic analyses. In this review, we have covered the works dealing with the response of mammalian liver to different nutritional stimuli such as fasting/feeding, caloric restriction, dietary carbohydrate, cholesterol, fat, protein, bile acid, salt, vitamin, and oligoelement contents. Quality of fats or proteins has been equally addressed, and has the influence of minor dietary components. Other compounds, not purely nutritional as those represented by alcohol and food additives, have been included due to their relevance in processed food. The influence has been studied not only on mRNA but also on miRNA. The wide scope of the technology clearly reflects that any simple intervention has profound changes in many metabolic parameters and that there is a synergy in response when more compounds are included in the intervention. Standardized arrays to systematically test the same genes in all studies and analyzing data to establish patterns of response are required, particularly for RNA sequencing. Moreover, RNA is a valuable, easy-screening ally but always requires further confirmation.

GRN2SBML: automated encoding and annotation of inferred gene regulatory networks complying with SBML

GRN2SBML automatically encodes gene regulatory networks derived from several inference tools in systems biology markup language. Providing a graphical user interface, the networks can be annotated via the simple object access protocol (SOAP)-based application programming interface of BioMart Central Portal and minimum information required in the annotation of models registry. Additionally, we provide an R-package, which processes the output of supported inference algorithms and automatically passes all required parameters to GRN2SBML. Therefore, GRN2SBML closes a gap in the processing pipeline between the inference of gene regulatory networks and their subsequent analysis, visualization and storage.

AVAILABILITY

GRN2SBML is freely available under the GNU Public License version 3 and can be downloaded from http://www.hki-jena.de/index.php/0/2/490.

SUPPLEMENTARY INFORMATION

General information on GRN2SBML, examples and tutorials are available at the tool's web page.