Integrative modeling of transcriptional regulation in response to antirheumatic therapy

Transcriptome profiling of peripheral blood immune cell populations in multiple sclerosis patients before and during treatment with a sphingosine-1-phosphate receptor modulator

AIMS

Fingolimod is a sphingosine-1-phosphate (S1P) receptor modulator approved for the treatment of the relapsing form of multiple sclerosis (MS). It prevents the egress of lymphocyte subpopulations from lymphoid tissues into the circulation. Here, we explored the broad effects of fingolimod on gene expression in different immune cell subsets.

METHODS

Utilizing 150 high-resolution microarrays from Affymetrix, we obtained the transcriptome profiles of 5 cell populations, which were separated from the peripheral blood of MS patients prior to and following oral administration of fingolimod.

RESULTS

After 3 months of treatment, significant transcriptome shifts were seen in CD4+ and CD8+ cells, which is mainly attributable to the selective homing of naive T cells and central memory T cells. Although the number of B cells was greatly reduced in the blood of fingolimod-treated MS patients, the analysis of differential expression in CD19+ cells identified only a small set of 42 genes, which indicated a slightly higher frequency of transitional B cells. The transcriptome signatures of CD14+ monocytes and CD56+ natural killer cells were not affected.

CONCLUSION

Our study corroborates changes in the composition of circulating immune cells in response to fingolimod and delineates the respective implications at the RNA level. Our data may be valuable for comparing the effects of novel S1P receptor modulating agents, which may be a therapeutic option for patients with secondary progressive MS as well.

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.

Systems pharmacology and enhanced pharmacodynamic models for understanding antibody-based drug action and toxicity

Pharmacokinetic (PK) and pharmacodynamic (PD) models seek to describe the temporal pattern of drug exposures and their associated pharmacological effects produced at micro- and macro-scales of organization. Antibody-based drugs have been developed for a large variety of diseases, with effects exhibited through a comprehensive range of mechanisms of action. Mechanism-based PK/PD and systems pharmacology models can play a major role in elucidating and integrating complex antibody pharmacological properties, such as nonlinear disposition and dynamical intracellular signaling pathways triggered by ligation to their cognate targets. Such complexities can be addressed through the use of robust computational modeling techniques that have proven powerful tools for pragmatic characterization of experimental data and for theoretical exploration of antibody efficacy and adverse effects. The primary objectives of such multi-scale mathematical models are to generate and test competing hypotheses and to predict clinical outcomes. In this review, relevant systems pharmacology and enhanced PD (ePD) models that are used as predictive tools for antibody-based drug action are reported. Their common conceptual features are highlighted, along with approaches used for modeling preclinical and clinically available data. Key examples illustrate how systems pharmacology and ePD models codify the interplay among complex biology, drug concentrations, and pharmacological effects. New hybrid modeling concepts that bridge cutting-edge systems pharmacology models with established PK/ePD models will be needed to anticipate antibody effects on disease in subpopulations and individual patients.

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.

How to Predict Molecular Interactions between Species?

Organisms constantly interact with other species through physical contact which leads to changes on the molecular level, for example the transcriptome. These changes can be monitored for all genes, with the help of high-throughput experiments such as RNA-seq or microarrays. The adaptation of the gene expression to environmental changes within cells is mediated through complex gene regulatory networks. Often, our knowledge of these networks is incomplete. Network inference predicts gene regulatory interactions based on transcriptome data. An emerging application of high-throughput transcriptome studies are dual transcriptomics experiments. Here, the transcriptome of two or more interacting species is measured simultaneously. Based on a dual RNA-seq data set of murine dendritic cells infected with the fungal pathogen Candida albicans, the software tool NetGenerator was applied to predict an inter-species gene regulatory network. To promote further investigations of molecular inter-species interactions, we recently discussed dual RNA-seq experiments for host-pathogen interactions and extended the applied tool NetGenerator (Schulze et al., 2015). The updated version of NetGenerator makes use of measurement variances in the algorithmic procedure and accepts gene expression time series data with missing values. Additionally, we tested multiple modeling scenarios regarding the stimuli functions of the gene regulatory network. Here, we summarize the work by Schulze et al. (2015) and put it into a broader context. We review various studies making use of the dual transcriptomics approach to investigate the molecular basis of interacting species. Besides the application to host-pathogen interactions, dual transcriptomics data are also utilized to study mutualistic and commensalistic interactions. Furthermore, we give a short introduction into additional approaches for the prediction of gene regulatory networks and discuss their application to dual transcriptomics data. We conclude that the application of network inference on dual-transcriptomics data is a promising approach to predict molecular inter-species interactions.

High-Density Peptide Microarray Analysis of IgG Autoantibody Reactivities in Serum and Cerebrospinal Fluid of Multiple Sclerosis Patients

Intrathecal immunoglobulin G (IgG) synthesis and oligoclonal IgG bands in cerebrospinal fluid (CSF) are hallmarks of multiple sclerosis (MS), but the antigen specificities remain enigmatic. Our study is the first investigating the autoantibody repertoire in paired serum and CSF samples from patients with relapsing-remitting MS (RRMS), primary progressive MS (PPMS), and other neurological diseases by the use of high-density peptide microarrays. Protein sequences of 45 presumed MS autoantigens (e.g.MOG, MBP, and MAG) were represented on the microarrays by overlapping 15mer peptides. IgG reactivities were screened against a total of 3991 peptides, including also selected viral epitopes. The measured antibody reactivities were highly individual but correlated for matched serum and CSF samples. We found 54 peptides to be recognized significantly more often by serum or CSF antibodies from MS patients compared with controls (pvalues <0.05). The results for RRMS and PPMS clearly overlapped. However, PPMS patients presented a broader peptide-antibody signature. The highest signals were detected for a peptide mapping to a region of the Epstein-Barr virus protein EBNA1 (amino acids 392-411), which is homologous to the N-terminal part of human crystallin alpha-B. Our data confirmed several known MS-associated antigens and epitopes, and they delivered additional potential linear epitopes, which await further validation. The peripheral and intrathecal humoral immune response in MS is polyspecific and includes antibodies that are also found in serum of patients with other diseases. Further studies are required to assess the pathogenic relevance of autoreactive and anti-EBNA1 antibodies as well as their combinatorial value as biomarkers for MS.

NERI: network-medicine based integrative approach for disease gene prioritization by relative importance

Background

Complex diseases are characterized as being polygenic and multifactorial, so this poses a challenge regarding the search for genes related to them. With the advent of high-throughput technologies for genome sequencing, gene expression measurements (transcriptome), and protein-protein interactions, complex diseases have been sistematically investigated. Particularly, Protein-Protein Interaction (PPI) networks have been used to prioritize genes related to complex diseases according to its topological features. However, PPI networks are affected by ascertainment bias, in which more studied proteins tend to have more connections, degrading the results quality. Additionally, methods using only PPI networks can provide only static and non-specific results, since the topologies of these networks are not specific of a given disease.

Results

The goal of this work is to develop a methodology that integrates PPI networks with disease specific data sources, such as GWAS and gene expression, to find genes more specific of a given complex disease. After the integration of PPI networks and gene expression data, the resulting network is used to connect genes related to the disease through the shortest paths that have the greatest concordance between their gene expressions. Both case and control expression data are used separately and, at the end, the most altered genes between the two conditions are selected. To evaluate the method, schizophrenia was adopted as case study.

Conclusion

Results show that the proposed method successfully retrieves differentially coexpressed genes in two conditions, while avoiding the bias from literature. Moreover we were able to achieve a greater concordance in the selection of important genes from different microarray studies of the same disease and to produce a more specific gene set related to the studied disease.

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.

Master regulators of infiltrate recruitment in autoimmune disease identified through network-based molecular deconvolution

Network-based molecular modeling of physiological behaviors has proven invaluable in the study of complex diseases such as cancer, but these approaches remain largely untested in contexts involving interacting tissues such as autoimmunity. Here, using Alopecia Areata (AA) as a model, we have adapted regulatory network analysis to specifically isolate physiological behaviors in the skin that contribute to the recruitment of immune cells in autoimmune disease. We use context-specific regulatory networks to deconvolve and identify skin-specific regulatory modules with IKZF1 and DLX4 as master regulators (MRs). These MRs are sufficient to induce AA-like molecular states in vitro in three cultured cell lines, resulting in induced NKG2D-dependent cytotoxicity. This work demonstrates the feasibility of a network-based approach for compartmentalizing and targeting molecular behaviors contributing to interactions between tissues in autoimmune disease.

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.

A bayesian framework that integrates heterogeneous data for inferring gene regulatory networks

Reconstruction of gene regulatory networks (GRNs) from experimental data is a fundamental challenge in systems biology. A number of computational approaches have been developed to infer GRNs from mRNA expression profiles. However, expression profiles alone are proving to be insufficient for inferring GRN topologies with reasonable accuracy. Recently, it has been shown that integration of external data sources (such as gene and protein sequence information, gene ontology data, protein-protein interactions) with mRNA expression profiles may increase the reliability of the inference process. Here, I propose a new approach that incorporates transcription factor binding sites (TFBS) and physical protein interactions (PPI) among transcription factors (TFs) in a Bayesian variable selection (BVS) algorithm which can infer GRNs from mRNA expression profiles subjected to genetic perturbations. Using real experimental data, I show that the integration of TFBS and PPI data with mRNA expression profiles leads to significantly more accurate networks than those inferred from expression profiles alone. Additionally, the performance of the proposed algorithm is compared with a series of least absolute shrinkage and selection operator (LASSO) regression-based network inference methods that can also incorporate prior knowledge in the inference framework. The results of this comparison suggest that BVS can outperform LASSO regression-based method in some circumstances.

Inference of the Xenopus tropicalis embryonic regulatory network and spatial gene expression patterns

BACKGROUND

During embryogenesis, signaling molecules produced by one cell population direct gene regulatory changes in neighboring cells and influence their developmental fates and spatial organization. One of the earliest events in the development of the vertebrate embryo is the establishment of three germ layers, consisting of the ectoderm, mesoderm and endoderm. Attempts to measure gene expression in vivo in different germ layers and cell types are typically complicated by the heterogeneity of cell types within biological samples (i.e., embryos), as the responses of individual cell types are intermingled into an aggregate observation of heterogeneous cell types. Here, we propose a novel method to elucidate gene regulatory circuits from these aggregate measurements in embryos of the frog Xenopus tropicalis using gene network inference algorithms and then test the ability of the inferred networks to predict spatial gene expression patterns.

RESULTS

We use two inference models with different underlying assumptions that incorporate existing network information, an ODE model for steady-state data and a Markov model for time series data, and contrast the performance of the two models. We apply our method to both control and knockdown embryos at multiple time points to reconstruct the core mesoderm and endoderm regulatory circuits. Those inferred networks are then used in combination with known dorsal-ventral spatial expression patterns of a subset of genes to predict spatial expression patterns for other genes. Both models are able to predict spatial expression patterns for some of the core mesoderm and endoderm genes, but interestingly of different gene subsets, suggesting that neither model is sufficient to recapitulate all of the spatial patterns, yet they are complementary for the patterns that they do capture.

CONCLUSION

The presented methodology of gene network inference combined with spatial pattern prediction provides an additional layer of validation to elucidate the regulatory circuits controlling the spatial-temporal dynamics in embryonic development.

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.

Glatiramer acetate treatment effects on gene expression in monocytes of multiple sclerosis patients

Background

Glatiramer acetate (GA) is a mixture of synthetic peptides used in the treatment of patients with relapsing-remitting multiple sclerosis (RRMS). The aim of this study was to investigate the effects of GA therapy on the gene expression of monocytes.

Methods

Monocytes were isolated from the peripheral blood of eight RRMS patients. The blood was obtained longitudinally before the start of GA therapy as well as after one day, one week, one month and two months. Gene expression was measured at the mRNA level by microarrays.

Results

More than 400 genes were identified as up-regulated or down-regulated in the course of therapy, and we analyzed their biological functions and regulatory interactions. Many of those genes are known to regulate lymphocyte activation and proliferation, but only a subset of genes was repeatedly differentially expressed at different time points during treatment.

Conclusions

Overall, the observed gene regulatory effects of GA on monocytes were modest and not stable over time. However, our study revealed several genes that are worthy of investigation in future studies on the molecular mechanisms of GA therapy.

MicroRNA expression changes during interferon-beta treatment in the peripheral blood of multiple sclerosis patients

MicroRNAs (miRNAs) are small non-coding RNA molecules acting as post-transcriptional regulators of gene expression. They are involved in many biological processes, and their dysregulation is implicated in various diseases, including multiple sclerosis (MS). Interferon-beta (IFN-beta) is widely used as a first-line immunomodulatory treatment of MS patients. Here, we present the first longitudinal study on the miRNA expression changes in response to IFN-beta therapy. Peripheral blood mononuclear cells (PBMC) were obtained before treatment initiation as well as after two days, four days, and one month, from patients with clinically isolated syndrome (CIS) and patients with relapsing-remitting MS (RRMS). We measured the expression of 651 mature miRNAs and about 19,000 mRNAs in parallel using real-time PCR arrays and Affymetrix microarrays. We observed that the up-regulation of IFN-beta-responsive genes is accompanied by a down-regulation of several miRNAs, including members of the mir-29 family. These differentially expressed miRNAs were found to be associated with apoptotic processes and IFN feedback loops. A network of miRNA-mRNA target interactions was constructed by integrating the information from different databases. Our results suggest that miRNA-mediated regulation plays an important role in the mechanisms of action of IFN-beta, not only in the treatment of MS but also in normal immune responses. miRNA expression levels in the blood may serve as a biomarker of the biological effects of IFN-beta therapy that may predict individual disease activity and progression.

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.

Interferon-beta therapy in multiple sclerosis: the short-term and long-term effects on the patients' individual gene expression in peripheral blood

Therapy with interferon-beta (IFN-beta) is a mainstay in the management of relapsing-remitting multiple sclerosis (MS), with proven long-term effectiveness and safety. Much has been learned about the molecular mechanisms of action of IFN-beta in the past years. Previous studies described more than a hundred genes to be modulated in expression in blood cells in response to the therapy. However, for many of these genes, the precise temporal expression pattern and the therapeutic relevance are unclear. We used Affymetrix microarrays to investigate in more detail the gene expression changes in peripheral blood mononuclear cells from MS patients receiving subcutaneous IFN-beta-1a. The blood samples were obtained longitudinally at five different time points up to 2 years after the start of therapy, and the patients were clinically followed up for 5 years. We examined the functions of the genes that were upregulated or downregulated at the transcript level after short-term or long-term treatment. Moreover, we analyzed their mutual interactions and their regulation by transcription factors. Compared to pretreatment levels, 96 genes were identified as highly differentially expressed, many of them already after the first IFN-beta injection. The interactions between these genes form a large network with multiple feedback loops, indicating the complex crosstalk between innate and adaptive immune responses during therapy. We discuss the genes and biological processes that might be important to reduce disease activity by attenuating the proliferation of autoreactive immune cells and their migration into the central nervous system. In summary, we present novel insights that extend the current knowledge on the early and late pharmacodynamic effects of IFN-beta therapy and describe gene expression differences between the individual patients that reflect clinical heterogeneity.

Global transcriptome sequencing identifies chlamydospore specific markers in Candida albicans and Candida dubliniensis

Candida albicans and Candida dubliniensis are pathogenic fungi that are highly related but differ in virulence and in some phenotypic traits. During in vitro growth on certain nutrient-poor media, C. albicans and C. dubliniensis are the only yeast species which are able to produce chlamydospores, large thick-walled cells of unknown function. Interestingly, only C. dubliniensis forms pseudohyphae with abundant chlamydospores when grown on Staib medium, while C. albicans grows exclusively as a budding yeast. In order to further our understanding of chlamydospore development and assembly, we compared the global transcriptional profile of both species during growth in liquid Staib medium by RNA sequencing. We also included a C. albicans mutant in our study which lacks the morphogenetic transcriptional repressor Nrg1. This strain, which is characterized by its constitutive pseudohyphal growth, specifically produces masses of chlamydospores in Staib medium, similar to C. dubliniensis. This comparative approach identified a set of putatively chlamydospore-related genes. Two of the homologous C. albicans and C. dubliniensis genes (CSP1 and CSP2) which were most strongly upregulated during chlamydospore development were analysed in more detail. By use of the green fluorescent protein as a reporter, the encoded putative cell wall related proteins were found to exclusively localize to C. albicans and C. dubliniensis chlamydospores. Our findings uncover the first chlamydospore specific markers in Candida species and provide novel insights in the complex morphogenetic development of these important fungal pathogens.

Inference of dynamical gene-regulatory networks based on time-resolved multi-stimuli multi-experiment data applying NetGenerator V2.0

BACKGROUND

Inference of gene-regulatory networks (GRNs) is important for understanding behaviour and potential treatment of biological systems. Knowledge about GRNs gained from transcriptome analysis can be increased by multiple experiments and/or multiple stimuli. Since GRNs are complex and dynamical, appropriate methods and algorithms are needed for constructing models describing these dynamics. Algorithms based on heuristic approaches reduce the effort in parameter identification and computation time.

RESULTS

The NetGenerator V2.0 algorithm, a heuristic for network inference, is proposed and described. It automatically generates a system of differential equations modelling structure and dynamics of the network based on time-resolved gene expression data. In contrast to a previous version, the inference considers multi-stimuli multi-experiment data and contains different methods for integrating prior knowledge. The resulting significant changes in the algorithmic procedures are explained in detail. NetGenerator is applied to relevant benchmark examples evaluating the inference for data from experiments with different stimuli. Also, the underlying GRN of chondrogenic differentiation, a real-world multi-stimulus problem, is inferred and analysed.

CONCLUSIONS

NetGenerator is able to determine the structure and parameters of GRNs and their dynamics. The new features of the algorithm extend the range of possible experimental set-ups, results and biological interpretations. Based upon benchmarks, the algorithm provides good results in terms of specificity, sensitivity, efficiency and model fit.

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

BACKGROUND

Network inference is an important tool to reveal the underlying interactions of biological systems. In the liver, a complex system of transcription factors is active to distribute signals and induce the cellular response following extracellular stimuli. Plenty of information is available about single transcription factors important for the different functions of the liver, but little is known about their causal relations to each other.

RESULTS

Given a DNA microarray time series dataset of collagen monolayers cultured murine hepatocytes, we identified 22 differentially expressed genes for which the corresponding protein is known to exhibit transcription factor activity. We developed the Extended TILAR (ExTILAR) network inference algorithm based on the modeling concept of the previously published TILAR algorithm. Using ExTILAR, we inferred a transcription factor network based on gene expression data which puts these important genes into a functional context. This way, we identified a previously unknown relationship between Tgif1 and Atf3 which we validated experimentally. Beside its known role in metabolic processes, this extends the knowledge about Tgif1 in hepatocytes towards a possible influence of processes such as proliferation and cell cycle. Moreover, two positive (i.e. double negative) regulatory loops were predicted that could give rise to bistable behavior. We further evaluated the performance of ExTILAR by systematic inference of an in silico network.

CONCLUSIONS

We present the ExTILAR algorithm, which combines the advantages of the regression based inference algorithm TILAR, like large network sizes processable and low computational costs, with the advantages of dynamic network models based on ordinary differential equation (i.e. in silico knock-down simulations). Like TILAR, ExTILAR makes use of various prior-knowledge types such as transcription factor binding site information and gene interaction knowledge to infer biologically meaningful gene regulatory networks. Therefore, ExTILAR is especially useful when a large number of genes is modeled using a small number of experimental data points.

Integrating external biological knowledge in the construction of regulatory networks from time-series expression data

BACKGROUND

Inference about regulatory networks from high-throughput genomics data is of great interest in systems biology. We present a Bayesian approach to infer gene regulatory networks from time series expression data by integrating various types of biological knowledge.

RESULTS

We formulate network construction as a series of variable selection problems and use linear regression to model the data. Our method summarizes additional data sources with an informative prior probability distribution over candidate regression models. We extend the Bayesian model averaging (BMA) variable selection method to select regulators in the regression framework. We summarize the external biological knowledge by an informative prior probability distribution over the candidate regression models.

CONCLUSIONS

We demonstrate our method on simulated data and a set of time-series microarray experiments measuring the effect of a drug perturbation on gene expression levels, and show that it outperforms leading regression-based methods in the literature.

An Interspecies Regulatory Network Inferred from Simultaneous RNA-seq of Candida albicans Invading Innate Immune Cells

The ability to adapt to diverse micro-environmental challenges encountered within a host is of pivotal importance to the opportunistic fungal pathogen Candida albicans. We have quantified C. albicans and M. musculus gene expression dynamics during phagocytosis by dendritic cells in a genome-wide, time-resolved analysis using simultaneous RNA-seq. A robust network inference map was generated from this dataset using NetGenerator, predicting novel interactions between the host and the pathogen. We experimentally verified predicted interdependent sub-networks comprising Hap3 in C. albicans, and Ptx3 and Mta2 in M. musculus. Remarkably, binding of recombinant Ptx3 to the C. albicans cell wall was found to regulate the expression of fungal Hap3 target genes as predicted by the network inference model. Pre-incubation of C. albicans with recombinant Ptx3 significantly altered the expression of Mta2 target cytokines such as IL-2 and IL-4 in a Hap3-dependent manner, further suggesting a role for Mta2 in host-pathogen interplay as predicted in the network inference model. We propose an integrated model for the functionality of these sub-networks during fungal invasion of immune cells, according to which binding of Ptx3 to the C. albicans cell wall induces remodeling via fungal Hap3 target genes, thereby altering the immune response to the pathogen. We show the applicability of network inference to predict interactions between host-pathogen pairs, demonstrating the usefulness of this systems biology approach to decipher mechanisms of microbial pathogenesis.

Genome-Wide Scale-Free Network Inference for Candida albicans

Discovery of essential genes in pathogenic organisms is an important step in the development of new medication. Despite a growing number of genome data available, little is known about C. albicans, a major fungal pathogen. Most of the human population carries C. albicans as commensal, but it can cause systemic infection that may lead to the death of the host if the immune system has deteriorated. In many organisms central nodes in the interaction network (hubs) play a crucial role for information and energy transport. Knock-outs of such hubs often lead to lethal phenotypes making them interesting drug targets. To identify these central genes via topological analysis, we inferred gene regulatory networks that are sparse and scale-free. We collected information from various sources to complement the limited expression data available. We utilized a linear regression algorithm to infer genome-wide gene regulatory interaction networks. To evaluate the predictive power of our approach, we used an automated text-mining system that scanned full-text research papers for known interactions. With the help of the compendium of known interactions, we also optimize the influence of the prior knowledge and the sparseness of the model to achieve the best results. We compare the results of our approach with those of other state-of-the-art network inference methods and show that we outperform those methods. Finally we identify a number of hubs in the genome of the fungus and investigate their biological relevance.

Regulatory interactions for iron homeostasis in Aspergillus fumigatus inferred by a Systems Biology approach

BACKGROUND

In System Biology, iterations of wet-lab experiments followed by modelling approaches and model-inspired experiments describe a cyclic workflow. This approach is especially useful for the inference of gene regulatory networks based on high-throughput gene expression data. Experiments can verify or falsify the predicted interactions allowing further refinement of the network model. Aspergillus fumigatus is a major human fungal pathogen. One important virulence trait is its ability to gain sufficient amounts of iron during infection process. Even though some regulatory interactions are known, we are still far from a complete understanding of the way iron homeostasis is regulated.

RESULTS

In this study, we make use of a reverse engineering strategy to infer a regulatory network controlling iron homeostasis in A. fumigatus. The inference approach utilizes the temporal change in expression data after a change from iron depleted to iron replete conditions. The modelling strategy is based on a set of linear differential equations and offers the possibility to integrate known regulatory interactions as prior knowledge. Moreover, it makes use of important selection criteria, such as sparseness and robustness. By compiling a list of known regulatory interactions for iron homeostasis in A. fumigatus and softly integrating them during network inference, we are able to predict new interactions between transcription factors and target genes. The proposed activation of the gene expression of hapX by the transcriptional regulator SrbA constitutes a so far unknown way of regulating iron homeostasis based on the amount of metabolically available iron. This interaction has been verified by Northern blots in a recent experimental study. In order to improve the reliability of the predicted network, the results of this experimental study have been added to the set of prior knowledge. The final network includes three SrbA target genes. Based on motif searching within the regulatory regions of these genes, we identify potential DNA-binding sites for SrbA. Our wet-lab experiments demonstrate high-affinity binding capacity of SrbA to the promoters of hapX, hemA and srbA.

CONCLUSIONS

This study presents an application of the typical Systems Biology circle and is based on cooperation between wet-lab experimentalists and in silico modellers. The results underline that using prior knowledge during network inference helps to predict biologically important interactions. Together with the experimental results, we indicate a novel iron homeostasis regulating system sensing the amount of metabolically available iron and identify the binding site of iron-related SrbA target genes. It will be of high interest to study whether these regulatory interactions are also important for close relatives of A. fumigatus and other pathogenic fungi, such as Candida albicans.

Regulatory network modelling of iron acquisition by a fungal pathogen in contact with epithelial cells

BACKGROUND

Reverse engineering of gene regulatory networks can be used to predict regulatory interactions of an organism faced with environmental changes, but can prove problematic, especially when focusing on complicated multi-factorial processes. Candida albicans is a major human fungal pathogen. During the infection process, this fungus is able to adapt to conditions of very low iron availability. Such adaptation is an important virulence attribute of virtually all pathogenic microbes. Understanding the regulation of iron acquisition genes will extend our knowledge of the complex regulatory changes during the infection process and might identify new potential drug targets. Thus, there is a need for efficient modelling approaches predicting key regulatory events of iron acquisition genes during the infection process.

RESULTS

This study deals with the regulation of C. albicans iron uptake genes during adhesion to and invasion into human oral epithelial cells. A reverse engineering strategy is presented, which is able to infer regulatory networks on the basis of gene expression data, making use of relevant selection criteria such as sparseness and robustness. An exhaustive use of available knowledge from different data sources improved the network prediction. The predicted regulatory network proposes a number of new target genes for the transcriptional regulators Rim101, Hap3, Sef1 and Tup1. Furthermore, the molecular mode of action for Tup1 is clarified. Finally, regulatory interactions between the transcription factors themselves are proposed. This study presents a model describing how C. albicans may regulate iron acquisition during contact with and invasion of human oral epithelial cells. There is evidence that some of the proposed regulatory interactions might also occur during oral infection.

CONCLUSIONS

This study focuses on a typical problem in Systems Biology where an interesting biological phenomenon is studied using a small number of available experimental data points. To overcome this limitation, a special modelling strategy was used which identifies sparse and robust networks. The data is augmented by an exhaustive search for additional data sources, helping to make proposals on regulatory interactions and to guide the modelling approach. The proposed modelling strategy is capable of finding known regulatory interactions and predicts a number of yet unknown biologically relevant regulatory interactions.

Network analysis of transcriptional regulation in response to intramuscular interferon-β-1a multiple sclerosis treatment

Interferon-β (IFN-β) is one of the major drugs for multiple sclerosis (MS) treatment. The purpose of this study was to characterize the transcriptional effects induced by intramuscular IFN-β-1a therapy in patients with relapsing-remitting form of MS. By using Affymetrix DNA microarrays, we obtained genome-wide expression profiles of peripheral blood mononuclear cells of 24 MS patients within the first 4 weeks of IFN-β administration. We identified 121 genes that were significantly up- or downregulated compared with baseline, with stronger changed expression at 1 week after start of therapy. Eleven transcription factor-binding sites (TFBS) are overrepresented in the regulatory regions of these genes, including those of IFN regulatory factors and NF-κB. We then applied TFBS-integrating least angle regression, a novel integrative algorithm for deriving gene regulatory networks from gene expression data and TFBS information, to reconstruct the underlying network of molecular interactions. An NF-κB-centered sub-network of genes was highly expressed in patients with IFN-β-related side effects. Expression alterations were confirmed by real-time PCR and literature mining was applied to evaluate network inference accuracy.

Functional genomics and rheumatoid arthritis: where have we been and where should we go?

Studies in model organisms and humans have begun to reveal the complexity of the transcriptome. In addition to serving as passive templates from which genes are translated, RNA molecules are active, functional elements of the cell whose products can detect, interact with, and modify other transcripts. Gene expression profiling is the method most commonly used thus far to enrich our understanding of the molecular basis of rheumatoid arthritis in adults and juvenile idiopathic arthritis in children. The feasibility of this approach for patient classification (for example, active versus inactive disease, disease subsets) and improving prognosis (for example, response to therapy) has been demonstrated over the past 7 years. Mechanistic understanding of disease-related differences in gene expression must be interpreted in the context of interactions with transcriptional regulatory molecules and epigenetic alterations of the genome. Ongoing work regarding such functional complexities in the human genome will likely bring both insight and surprise to our understanding of rheumatoid arthritis.

Regulatory impact factors: unraveling the transcriptional regulation of complex traits from expression data

MOTIVATION

Although transcription factors (TF) play a central regulatory role, their detection from expression data is limited due to their low, and often sparse, expression. In order to fill this gap, we propose a regulatory impact factor (RIF) metric to identify critical TF from gene expression data.

RESULTS

To substantiate the generality of RIF, we explore a set of experiments spanning a wide range of scenarios including breast cancer survival, fat, gonads and sex differentiation. We show that the strength of RIF lies in its ability to simultaneously integrate three sources of information into a single measure: (i) the change in correlation existing between the TF and the differentially expressed (DE) genes; (ii) the amount of differential expression of DE genes; and (iii) the abundance of DE genes. As a result, RIF analysis assigns an extreme score to those TF that are consistently most differentially co-expressed with the highly abundant and highly DE genes (RIF1), and to those TF with the most altered ability to predict the abundance of DE genes (RIF2). We show that RIF analysis alone recovers well-known experimentally validated TF for the processes studied. The TF identified confirm the importance of PPAR signaling in adipose development and the importance of transduction of estrogen signals in breast cancer survival and sexual differentiation. We argue that RIF has universal applicability, and advocate its use as a promising hypotheses generating tool for the systematic identification of novel TF not yet documented as critical.

Long-term genome-wide blood RNA expression profiles yield novel molecular response candidates for IFN-β-1b treatment in relapsing remitting MS

Aims: In multiple sclerosis patients, treatment with recombinant IFN-β (rIFN-β) is partially efficient in reducing clinical exacerbations. However, its molecular mechanism of action is still under scrutiny. Materials & methods: We used DNA microarrays (Affymetrix, CA, USA) and peripheral mononuclear blood cells from 25 relapsing remitting multiple sclerosis patients to analyze the longitudinal transcriptional profile within 2 years of rIFN-β administration. Sets of differentially expressed genes were attained by applying a combination of independent criteria, thereby providing efficient data curation and gene filtering that accounted for technical and biological noise. Gene ontology term-association analysis and scientific literature text mining were used to explore evidence of gene interaction. Results: Post-therapy initiation, we identified 42 (day 2), 175 (month 1), 103 (month 12) and 108 (month 24) differentially expressed genes. Increased expression of established IFN-β marker genes, as well as differential expression of circulating IFN-β-responsive candidate genes, were observed. MS4A1 (CD20), a known target of B-cell depletion therapy, was significantly downregulated after one month. CMPK2, FCER1A, and FFAR2 appeared as hitherto unrecognized multiple sclerosis treatment-related differentially expressed genes that were consistently modulated over time. Overall, 84 interactions between 54 genes were attained, of which two major gene networks were identified at an earlier stage of therapy: the first (n = 15 genes) consisted of mostly known IFN-β-activated genes, whereas the second (n = 12) mainly contained downregulated genes that to date have not been associated with IFN-β effects in multiple sclerosis array research. Conclusion: We achieved both a broadening of the knowledge of IFN-β mechanism-of-action-related constituents and the identification of time-dependent interactions between IFN-β regulated genes.