Time Series Gene Expression Research Articles

Transcriptome dynamics in early zebrafish embryogenesis determined by high-resolution time course analysis of 180 successive, individual zebrafish embryos

BackgroundRecently, much progress has been made in the field of gene-expression in early embryogenesis. However, the dynamic behaviour of transcriptomes in individual embryos has hardly been studied yet and the time points at which pools of embryos are collected are usually still quite far apart. Here, we present a high-resolution gene-expression time series with 180 individual zebrafish embryos, obtained from nine different spawns, developmentally ordered and profiled from late blastula to mid-gastrula stage. On average one embryo per minute was analysed. The focus was on identification and description of the transcriptome dynamics of the expressed genes in this embryonic stage, rather than to biologically interpret profiles in cellular processes and pathways.ResultsIn the late blastula to mid-gastrula stage, we found 6,734 genes being expressed with low variability and rather gradual changes. Ten types of dynamic behaviour were defined, such as genes with continuously increasing or decreasing expression, and all expressed genes were grouped into these types. Also, the exact expression starting and stopping points of several hundred genes during this developmental period could be pinpointed. Although the resolution of the experiment was so high, that we were able to clearly identify four known oscillating genes, no genes were observed with a peaking expression. Additionally, several genes showed expression at two or three distinct levels that strongly related to the spawn an embryo originated from.ConclusionOur unique experimental set-up of whole-transcriptome analysis of 180 individual embryos, provided an unparalleled in-depth insight into the dynamics of early zebrafish embryogenesis. The existence of a tightly regulated embryonic transcriptome program, even between individuals from different spawns is shown. We have made the expression profile of all genes available for domain experts. The fact that we were able to separate the different spawns by their gene-expression variance over all expressed genes, underlines the importance of spawn specificity, as well as the unexpectedly tight gene-expression regulation in early zebrafish embryogenesis.

Systematic identification of an integrative network module during senescence from time-series gene expression

BackgroundCellular senescence irreversibly arrests growth of human diploid cells. In addition, recent studies have indicated that senescence is a multi-step evolving process related to important complex biological processes. Most studies analyzed only the genes and their functions representing each senescence phase without considering gene-level interactions and continuously perturbed genes. It is necessary to reveal the genotypic mechanism inferred by affected genes and their interaction underlying the senescence process.ResultsWe suggested a novel computational approach to identify an integrative network which profiles an underlying genotypic signature from time-series gene expression data. The relatively perturbed genes were selected for each time point based on the proposed scoring measure denominated as perturbation scores. Then, the selected genes were integrated with protein-protein interactions to construct time point specific network. From these constructed networks, the conserved edges across time point were extracted for the common network and statistical test was performed to demonstrate that the network could explain the phenotypic alteration. As a result, it was confirmed that the difference of average perturbation scores of common networks at both two time points could explain the phenotypic alteration. We also performed functional enrichment on the common network and identified high association with phenotypic alteration. Remarkably, we observed that the identified cell cycle specific common network played an important role in replicative senescence as a key regulator.ConclusionsHeretofore, the network analysis from time series gene expression data has been focused on what topological structure was changed over time point. Conversely, we focused on the conserved structure but its context was changed in course of time and showed it was available to explain the phenotypic changes. We expect that the proposed method will help to elucidate the biological mechanism unrevealed by the existing approaches.

The Ubiquitin-like with PHD and Ring Finger Domains 1 (UHRF1)/DNA Methyltransferase 1 (DNMT1) Axis Is a Primary Regulator of Cell Senescence

As senescence develops, cells sequentially acquire diverse senescent phenotypes along with simultaneous multistage gene reprogramming. It remains unclear what acts as the key regulator of the collective changes in gene expression at initiation of senescent reprogramming. Here we analyzed time series gene expression profiles obtained in two different senescence models in human diploid fibroblasts: replicative senescence and H2O2-induced senescence. Our results demonstrate that suppression of DNA methyltransferase 1 (DNMT1)-mediated DNA methylation activity was an initial event prior to the display of senescent phenotypes. We identified seven DNMT1-interacting proteins, ubiquitin-like with PHD and ring finger domains 1 (UHRF1), EZH2, CHEK1, SUV39H1, CBX5, PARP1, and HELLS (also known as LSH (lymphoid-specific helicase) 1), as being commonly down-regulated at the same time point as DNMT1 in both senescence models. Knockdown experiments revealed that, among the DNMT1-interacting proteins, only UHRF1 knockdown suppressed DNMT1 transcription. However, UHRF1 overexpression alone did not induce DNMT1 expression, indicating that UHRF1 was essential but not sufficient for DNMT1 transcription. Although UHRF1 knockdown effectively induced senescence, this was significantly attenuated by DNMT1 overexpression, clearly implicating the UHRF1/DNMT1 axis in senescence. Bioinformatics analysis further identified WNT5A as a downstream effector of UHRF1/DNMT1-mediated senescence. Senescence-associated hypomethylation was found at base pairs -1569 to -1363 from the transcription start site of the WNT5A gene in senescent human diploid fibroblasts. As expected, WNT5A overexpression induced senescent phenotypes. Overall, our results indicate that decreased UHRF1 expression is a key initial event in the suppression of DNMT1-mediated DNA methylation and in the consequent induction of senescence via increasing WNT5A expression.

A novel mutual information-based Boolean network inference method from time-series gene expression data.

BackgroundInferring a gene regulatory network from time-series gene expression data in systems biology is a challenging problem. Many methods have been suggested, most of which have a scalability limitation due to the combinatorial cost of searching a regulatory set of genes. In addition, they have focused on the accurate inference of a network structure only. Therefore, there is a pressing need to develop a network inference method to search regulatory genes efficiently and to predict the network dynamics accurately.ResultsIn this study, we employed a Boolean network model with a restricted update rule scheme to capture coarse-grained dynamics, and propose a novel mutual information-based Boolean network inference (MIBNI) method. Given time-series gene expression data as an input, the method first identifies a set of initial regulatory genes using mutual information-based feature selection, and then improves the dynamics prediction accuracy by iteratively swapping a pair of genes between sets of the selected regulatory genes and the other genes. Through extensive simulations with artificial datasets, MIBNI showed consistently better performance than six well-known existing methods, REVEAL, Best-Fit, RelNet, CST, CLR, and BIBN in terms of both structural and dynamics prediction accuracy. We further tested the proposed method with two real gene expression datasets for an Escherichia coli gene regulatory network and a fission yeast cell cycle network, and also observed better results using MIBNI compared to the six other methods.ConclusionsTaken together, MIBNI is a promising tool for predicting both the structure and the dynamics of a gene regulatory network.

Expectation propagation for large scale Bayesian inference of non-linear molecular networks from perturbation data.

Inferring the structure of molecular networks from time series protein or gene expression data provides valuable information about the complex biological processes of the cell. Causal network structure inference has been approached using different methods in the past. Most causal network inference techniques, such as Dynamic Bayesian Networks and ordinary differential equations, are limited by their computational complexity and thus make large scale inference infeasible. This is specifically true if a Bayesian framework is applied in order to deal with the unavoidable uncertainty about the correct model. We devise a novel Bayesian network reverse engineering approach using ordinary differential equations with the ability to include non-linearity. Besides modeling arbitrary, possibly combinatorial and time dependent perturbations with unknown targets, one of our main contributions is the use of Expectation Propagation, an algorithm for approximate Bayesian inference over large scale network structures in short computation time. We further explore the possibility of integrating prior knowledge into network inference. We evaluate the proposed model on DREAM4 and DREAM8 data and find it competitive against several state-of-the-art existing network inference methods.

Module Anchored Network Inference: A Sequential Module-Based Approach to Novel Gene Network Construction from Genomic Expression Data on Human Disease Mechanism.

Different computational approaches have been examined and compared for inferring network relationships from time-series genomic data on human disease mechanisms under the recent Dialogue on Reverse Engineering Assessment and Methods (DREAM) challenge. Many of these approaches infer all possible relationships among all candidate genes, often resulting in extremely crowded candidate network relationships with many more False Positives than True Positives. To overcome this limitation, we introduce a novel approach, Module Anchored Network Inference (MANI), that constructs networks by analyzing sequentially small adjacent building blocks (modules). Using MANI, we inferred a 7-gene adipogenesis network based on time-series gene expression data during adipocyte differentiation. MANI was also applied to infer two 10-gene networks based on time-course perturbation datasets from DREAM3 and DREAM4 challenges. MANI well inferred and distinguished serial, parallel, and time-dependent gene interactions and network cascades in these applications showing a superior performance to other in silico network inference techniques for discovering and reconstructing gene network relationships.

An Advanced Omic Approach to Identify Co-Regulated Clusters and Transcription Regulation Network with AGCT and SHOE Methods.

To obtain the global picture of genetic machinery for massive high-throughput gene expression data, novel data-driven unsupervised learning approaches are becoming essentially important. For this purpose, basic analytic workflow has been established and should include two steps: first, unsupervised clustering to identify genes with similar behavior upon exposure to a signal, and second, identification of transcription factors regulating those genes. In this chapter, we will describe an advanced tool that can be used for analyzing and characterizing large-scale time-series gene expression composed of a two-step approach. For the first step, we developed an original method "A Geometric Clustering Tool" (AGCT) that unveils the complex architecture of large-scale time-series gene expression data in a real-time manner using cutting edge techniques of low dimension manifold learning, data clustering, and visualization. For the second step, we established an original method "Sequence Homology in Eukaryotes" (SHOE) executing comparative genomic analysis on humans, mice, and rats.

Integration of Pathway Knowledge and Dynamic Bayesian Networks for the Prediction of Oral Cancer Recurrence.

Oral squamous cell carcinoma has been characterized as a complex disease which involves dynamic genomic changes at the molecular level. These changes indicate the worth to explore the interactions of the molecules and especially of differentially expressed genes that contribute to cancer progression. Moreover, based on this knowledge the identification of differentially expressed genes and related molecular pathways is of great importance. In the present study, we exploit differentially expressed genes in order to further perform pathway enrichment analysis. According to our results we found significant pathways in which the disease associated genes have been identified as strongly enriched. Furthermore, based on the results of the pathway enrichment analysis we propose a methodology for predicting oral cancer recurrence using dynamic Bayesian networks. The methodology takes into consideration time series gene expression data in order to predict a disease recurrence. Subsequently, we are able to conjecture about the causal interactions between genes in consecutive time intervals. Concerning the performance of the predictive models, the overall accuracy of the algorithm is 81.8% and the area under the ROC curve 89.2% regarding the knowledge from the overrepresented pre-NOTCH Expression and processing pathway.

Reverse engineering gene regulatory networks from measurement with missing values.

BackgroundGene expression time series data are usually in the form of high-dimensional arrays. Unfortunately, the data may sometimes contain missing values: for either the expression values of some genes at some time points or the entire expression values of a single time point or some sets of consecutive time points. This significantly affects the performance of many algorithms for gene expression analysis that take as an input, the complete matrix of gene expression measurement. For instance, previous works have shown that gene regulatory interactions can be estimated from the complete matrix of gene expression measurement. Yet, till date, few algorithms have been proposed for the inference of gene regulatory network from gene expression data with missing values.ResultsWe describe a nonlinear dynamic stochastic model for the evolution of gene expression. The model captures the structural, dynamical, and the nonlinear natures of the underlying biomolecular systems. We present point-based Gaussian approximation (PBGA) filters for joint state and parameter estimation of the system with one-step or two-step missing measurements. The PBGA filters use Gaussian approximation and various quadrature rules, such as the unscented transform (UT), the third-degree cubature rule and the central difference rule for computing the related posteriors. The proposed algorithm is evaluated with satisfying results for synthetic networks, in silico networks released as a part of the DREAM project, and the real biological network, the in vivo reverse engineering and modeling assessment (IRMA) network of yeast Saccharomyces cerevisiae.ConclusionPBGA filters are proposed to elucidate the underlying gene regulatory network (GRN) from time series gene expression data that contain missing values. In our state-space model, we proposed a measurement model that incorporates the effect of the missing data points into the sequential algorithm. This approach produces a better inference of the model parameters and hence, more accurate prediction of the underlying GRN compared to when using the conventional Gaussian approximation (GA) filters ignoring the missing data points.Electronic supplementary materialThe online version of this article (doi:10.1186/s13637-016-0055-8) contains supplementary material, which is available to authorized users.

SwitchFinder - a novel method and query facility for discovering dynamic gene expression patterns.

BackgroundBiological systems and processes are highly dynamic. To gain insights into their functioning time-resolved measurements are necessary. Time-resolved gene expression data captures temporal behaviour of the genes genome-wide under various biological conditions: in response to stimuli, during cell cycle, differentiation or developmental programs. Dissecting dynamic gene expression patterns from this data may shed light on the functioning of the gene regulatory system. The present approach facilitates this discovery. The fundamental idea behind it is the following: there are change-points (switches) in the gene behaviour separating intervals of increasing and decreasing activity, whereas the intervals may have different durations. Elucidating the switch-points is important for the identification of biologically meanigfull features and patterns of the gene dynamics.ResultsWe developed a statistical method, called SwitchFinder, for the analysis of time-series data, in particular gene expression data, based on a change-point model. Fitting the model to the gene expression time-courses indicates switch-points between increasing and decreasing activities of each gene. Two types of the model - based on linear and on generalized logistic function - were used to capture the data between the switch-points. Model inference was facilitated with the Bayesian methodology using Markov chain Monte Carlo (MCMC) technique Gibbs sampling. Further on, we introduced features of the switch-points: growth, decay, spike and cleft, which reflect important dynamic aspects. With this, the gene expression profiles are represented in a qualitative manner - as sets of the dynamic features at their onset-times. We developed a Web application of the approach, enabling to put queries to the gene expression time-courses and to deduce groups of genes with common dynamic patterns.SwitchFinder was applied to our original data - the gene expression time-series measured in neuroblastoma cell line upon treatment with all-trans retinoic acid (ATRA). The analysis revealed eight patterns of the gene expression responses to ATRA, indicating the induction of the BMP, WNT, Notch, FGF and NTRK-receptor signaling pathways involved in cell differentiation, as well as the repression of the cell-cycle related genes.ConclusionsSwitchFinder is a novel approach to the analysis of biological time-series data, supporting inference and interactive exploration of its inherent dynamic patterns, hence facilitating biological discovery process. SwitchFinder is freely available at https://newbioinformatics.eu/switchfinder.Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-016-1391-0) contains supplementary material, which is available to authorized users.

Prioritizing biological pathways by recognizing context in time-series gene expression data.

BackgroundThe primary goal of pathway analysis using transcriptome data is to find significantly perturbed pathways. However, pathway analysis is not always successful in identifying pathways that are truly relevant to the context under study. A major reason for this difficulty is that a single gene is involved in multiple pathways. In the KEGG pathway database, there are 146 genes, each of which is involved in more than 20 pathways. Thus activation of even a single gene will result in activation of many pathways. This complex relationship often makes the pathway analysis very difficult. While we need much more powerful pathway analysis methods, a readily available alternative way is to incorporate the literature information.ResultsIn this study, we propose a novel approach for prioritizing pathways by combining results from both pathway analysis tools and literature information. The basic idea is as follows. Whenever there are enough articles that provide evidence on which pathways are relevant to the context, we can be assured that the pathways are indeed related to the context, which is termed as relevance in this paper. However, if there are few or no articles reported, then we should rely on the results from the pathway analysis tools, which is termed as significance in this paper. We realized this concept as an algorithm by introducing Context Score and Impact Score and then combining the two into a single score. Our method ranked truly relevant pathways significantly higher than existing pathway analysis tools in experiments with two data sets.ConclusionsOur novel framework was implemented as ContextTRAP by utilizing two existing tools, TRAP and BEST. ContextTRAP will be a useful tool for the pathway based analysis of gene expression data since the user can specify the context of the biological experiment in a set of keywords. The web version of ContextTRAP is available at http://biohealth.snu.ac.kr/software/contextTRAP.Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-016-1335-8) contains supplementary material, which is available to authorized users.

From yeast to hypha: defining transcriptomic signatures of the morphological switch in the dimorphic fungal pathogen Ophiostoma novo-ulmi.

BackgroundYeast-to-hypha transition is a major morphological change in fungi. Molecular regulators and pathways that are involved in this process have been extensively studied in model species, including Saccharomyces cerevisiae. The Mitogen-Actived Protein Kinase (MAPK) cascade, for example, is known to be involved in the yeast-to-pseudohypha switch. Yet the conservation of mechanisms regulating such morphological changes in non-model fungi is still poorly understood. Here, we investigate cell remodeling and transcriptomic modifications that occur during this morphological switch in the highly aggressive ascomycete fungus Ophiostoma novo-ulmi, the causal agent of Dutch elm disease.ResultsUsing a combination of light microscopy, scanning electron microscopy and flow cytometry, we demonstrate that the morphological switch occurs in less than 27 h, with phenotypic cell modifications being detected within the first 4 h. Using RNAseq, we found that over 22% of the genome of O. novo-ulmi is differentially expressed during the transition. By performing clustering analyses of time series gene expression data, we identified several sets of genes that are differentially expressed according to distinct and representative temporal profiles. Further, we found that several genes that are homologous to S. cerevisiae MAPK genes are regulated during the yeast-to-hypha transition in O. novo-ulmi and mostly over-expressed, suggesting convergence in gene expression regulation.ConclusionsOur results are the first report of a time-course experiment monitoring the morphological transition in a non-model Sordariomycota species and reveal many genes of interest for further functional investigations of fungal dimorphism.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3251-8) contains supplementary material, which is available to authorized users.

DREISS: Using State-Space Models to Infer the Dynamics of Gene Expression Driven by External and Internal Regulatory Networks.

Gene expression is controlled by the combinatorial effects of regulatory factors from different biological subsystems such as general transcription factors (TFs), cellular growth factors and microRNAs. A subsystem’s gene expression may be controlled by its internal regulatory factors, exclusively, or by external subsystems, or by both. It is thus useful to distinguish the degree to which a subsystem is regulated internally or externally–e.g., how non-conserved, species-specific TFs affect the expression of conserved, cross-species genes during evolution. We developed a computational method (DREISS, dreiss.gerteinlab.org) for analyzing the Dynamics of gene expression driven by Regulatory networks, both External and Internal based on State Space models. Given a subsystem, the “state” and “control” in the model refer to its own (internal) and another subsystem’s (external) gene expression levels. The state at a given time is determined by the state and control at a previous time. Because typical time-series data do not have enough samples to fully estimate the model’s parameters, DREISS uses dimensionality reduction, and identifies canonical temporal expression trajectories (e.g., degradation, growth and oscillation) representing the regulatory effects emanating from various subsystems. To demonstrate capabilities of DREISS, we study the regulatory effects of evolutionarily conserved vs. divergent TFs across distant species. In particular, we applied DREISS to the time-series gene expression datasets of C. elegans and D. melanogaster during their embryonic development. We analyzed the expression dynamics of the conserved, orthologous genes (orthologs), seeing the degree to which these can be accounted for by orthologous (internal) versus species-specific (external) TFs. We found that between two species, the orthologs have matched, internally driven expression patterns but very different externally driven ones. This is particularly true for genes with evolutionarily ancient functions (e.g. the ribosomal proteins), in contrast to those with more recently evolved functions (e.g., cell-cell communication). This suggests that despite striking morphological differences, some fundamental embryonic-developmental processes are still controlled by ancient regulatory systems.

The Local Edge Machine: inference of dynamic models of gene regulation.

We present a novel approach, the Local Edge Machine, for the inference of regulatory interactions directly from time-series gene expression data. We demonstrate its performance, robustness, and scalability on in silico datasets with varying behaviors, sizes, and degrees of complexity. Moreover, we demonstrate its ability to incorporate biological prior information and make informative predictions on a well-characterized in vivo system using data from budding yeast that have been synchronized in the cell cycle. Finally, we use an atlas of transcription data in a mammalian circadian system to illustrate how the method can be used for discovery in the context of large complex networks.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-016-1076-z) contains supplementary material, which is available to authorized users.

Application of dynamic topic models to toxicogenomics data.

BackgroundAll biological processes are inherently dynamic. Biological systems evolve transiently or sustainably according to sequential time points after perturbation by environment insults, drugs and chemicals. Investigating the temporal behavior of molecular events has been an important subject to understand the underlying mechanisms governing the biological system in response to, such as, drug treatment. The intrinsic complexity of time series data requires appropriate computational algorithms for data interpretation. In this study, we propose, for the first time, the application of dynamic topic models (DTM) for analyzing time-series gene expression data.ResultsA large time-series toxicogenomics dataset was studied. It contains over 3144 microarrays of gene expression data corresponding to rat livers treated with 131 compounds (most are drugs) at two doses (control and high dose) in a repeated schedule containing four separate time points (4-, 8-, 15- and 29-day). We analyzed, with DTM, the topics (consisting of a set of genes) and their biological interpretations over these four time points. We identified hidden patterns embedded in this time-series gene expression profiles. From the topic distribution for compound-time condition, a number of drugs were successfully clustered by their shared mode-of-action such as PPARɑ agonists and COX inhibitors. The biological meaning underlying each topic was interpreted using diverse sources of information such as functional analysis of the pathways and therapeutic uses of the drugs. Additionally, we found that sample clusters produced by DTM are much more coherent in terms of functional categories when compared to traditional clustering algorithms.ConclusionsWe demonstrated that DTM, a text mining technique, can be a powerful computational approach for clustering time-series gene expression profiles with the probabilistic representation of their dynamic features along sequential time frames. The method offers an alternative way for uncovering hidden patterns embedded in time series gene expression profiles to gain enhanced understanding of dynamic behavior of gene regulation in the biological system.Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-016-1225-0) contains supplementary material, which is available to authorized users.

Parameter estimation in tree graph metabolic networks.

We study the glycosylation processes that convert initially toxic substrates to nutritionally valuable metabolites in the flavonoid biosynthesis pathway of tomato (Solanum lycopersicum) seedlings. To estimate the reaction rates we use ordinary differential equations (ODEs) to model the enzyme kinetics. A popular choice is to use a system of linear ODEs with constant kinetic rates or to use Michaelis–Menten kinetics. In reality, the catalytic rates, which are affected among other factors by kinetic constants and enzyme concentrations, are changing in time and with the approaches just mentioned, this phenomenon cannot be described. Another problem is that, in general these kinetic coefficients are not always identifiable. A third problem is that, it is not precisely known which enzymes are catalyzing the observed glycosylation processes. With several hundred potential gene candidates, experimental validation using purified target proteins is expensive and time consuming. We aim at reducing this task via mathematical modeling to allow for the pre-selection of most potential gene candidates. In this article we discuss a fast and relatively simple approach to estimate time varying kinetic rates, with three favorable properties: firstly, it allows for identifiable estimation of time dependent parameters in networks with a tree-like structure. Secondly, it is relatively fast compared to usually applied methods that estimate the model derivatives together with the network parameters. Thirdly, by combining the metabolite concentration data with a corresponding microarray data, it can help in detecting the genes related to the enzymatic processes. By comparing the estimated time dynamics of the catalytic rates with time series gene expression data we may assess potential candidate genes behind enzymatic reactions. As an example, we show how to apply this method to select prominent glycosyltransferase genes in tomato seedlings.

Prediction of oral cancer recurrence using dynamic Bayesian networks.

We propose a methodology for predicting oral cancer recurrence using Dynamic Bayesian Networks. The methodology takes into consideration time series gene expression data collected at the follow-up study of patients that had or had not suffered a disease relapse. Based on that knowledge, our aim is to infer the corresponding dynamic Bayesian networks and subsequently conjecture about the causal relationships among genes within the same time-slice and between consecutive time-slices. Moreover, the proposed methodology aims to (i) assess the prognosis of patients regarding oral cancer recurrence and at the same time, (ii) provide important information about the underlying biological processes of the disease.

Inferring gene regulatory networks using a time-delayed mass action model.

This paper demonstrates a new time-delayed mass action model which applies a set of delay differential equations (DDEs) to represent the dynamics of gene regulatory networks (GRNs). The mass action model is a classical model which is often used to describe the kinetics of biochemical processes, so it is fit for GRN modeling. The ability to incorporate time-delayed parameters in this model enables different time delays of interaction between genes. Moreover, an efficient learning method which employs population-based incremental learning (PBIL) algorithm and trigonometric differential evolution (TDE) algorithm TDE is proposed to automatically evolve the structure of the network and infer the optimal parameters from observed time-series gene expression data. Experiments on three well-known motifs of GRN and a real budding yeast cell cycle network show that the proposal can not only successfully infer the network structure and parameters but also has a strong anti-noise ability. Compared with other works, this method also has a great improvement in performances.

Batch Mode TD($\lambda$ ) for Controlling Partially Observable Gene Regulatory Networks.

External control of gene regulatory networks (GRNs) has received much attention in recent years. The aim is to find a series of actions to apply to a gene regulation system making it avoid its diseased states. In this work, we propose a novel method for controlling partially observable GRNs combining batch mode reinforcement learning (Batch RL) and TD() algorithms. Unlike the existing studies inferring a computational model from gene expression data, and obtaining a control policy over the constructed model, our idea is to interpret the time series gene expression data as a sequence of observations that the system produced, and obtain an approximate stochastic policy directly from the gene expression data without estimation of the internal states of the partially observable environment. Thereby, we get rid of the most time consuming phases of the existing studies, inferring a model and running the model for the control. Results show that our method is able to provide control solutions for regulation systems of several thousands of genes only in seconds, whereas existing studies cannot solve control problems of even a few dozens of genes. Results also show that our approximate stochastic policies are almost as good as the policies generated by the existing studies.

S100A6 Regulates Endothelial Cell Cycle Progression by Attenuating Antiproliferative Signal Transducers and Activators of Transcription 1 Signaling.

S100A6, a member of the S100 protein family, has been described as relevant for cell cycle entry and progression in endothelial cells. The molecular mechanism conferring S100A6's proliferative actions, however, remained elusive. Originating from the clinically relevant observation of enhanced S100A6 protein expression in proliferating endothelial cells in remodeling coronary and carotid arteries, our study unveiled S100A6 as a suppressor of antiproliferative signal transducers and activators of transcription 1 signaling. Discovery of the molecular liaison was enabled by combining gene expression time series analysis with bioinformatic pathway modeling in S100A6-silenced human endothelial cells stimulated with vascular endothelial growth factor A. This unbiased approach led to successful identification and experimental validation of interferon-inducible transmembrane protein 1 and protein inhibitors of activated signal transducers and activators of transcription as key components of the link between S100A6 and signal transducers and activators of transcription 1. Given the important role of coordinated endothelial cell cycle activity for integrity and reconstitution of the inner lining of arterial blood vessels in health and disease, signal transducers and activators of transcription 1 suppression by S100A6 may represent a promising therapeutic target to facilitate reendothelialization in damaged vessels.