bbc 2015
BBC2015_booklet
BBC2015_booklet
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
BeNeLux Bioinformatics Conference – Antwerp, December 7-8 <strong>2015</strong><br />
Abstract ID: P<br />
Poster<br />
10th Benelux Bioinformatics Conference <strong>bbc</strong> <strong>2015</strong><br />
P24. DOSE-TIME NETWORK IDENTIFICATION: A NEW METHOD FOR<br />
GENE REGULATORY NETWORK INFERENCE FROM GENE EXPRESSION<br />
DATA WITH MULTIPLE DOSES AND TIME POINTS<br />
Diana M Hendrickx 1* , Danyel G J Jennen 1 & Jos C S Kleinjans 1 .<br />
Department of Toxicogenomics, Maastricht University, The Netherlands 1 .<br />
*d.hendrickx@maastrichtuniversity.nl<br />
Toxicogenomics, the application of ‘omics’ technologies to toxicology, is a rapidly growing field due to the need for<br />
alternatives to animal experiments for toxicity testing of compounds. Identification of gene regulatory networks affected<br />
by compounds is important to gain more insight into the mode of action of a toxic compound. The response to a toxic<br />
compound is both time and dose dependent. Therefore, toxicogenomics data are often measured across several time<br />
points and doses. However, to our knowledge, there does not exist a method for gene regulatory network inference that<br />
takes into account both time and dose dependencies. Here we present Dose-Time Network Identification (DTNI), a novel<br />
gene regulatory network inference algorithm that takes into account both dose and time dependencies in the data. We<br />
show that DTNI can be used to infer gene regulatory networks affected by a group of compounds with the same mode of<br />
action. This is illustrated with gene expression (microarray) data from COX inhibitors, measured in human hepatocytes.<br />
INTRODUCTION<br />
Identifying and understanding gene regulatory networks<br />
(GRN) influenced by chemical compounds is one of the<br />
main challenges of systems toxicology. A GRN affected<br />
by one or more compounds evolves over time and with<br />
dose. The analysis of gene expression data measured at<br />
multiple time points and for multiple doses can provide<br />
more insight in the effects of compounds. Therefore, there<br />
is a need for mathematical approaches for GRN<br />
identification from this type of data.<br />
METHODS<br />
One of the mathematical approaches currently used for<br />
GRN inference is based on ordinary differential equations<br />
(ODE), where changes in gene expression over time are<br />
related to each other and to the external perturbation (i.e.<br />
the dose of the compound). Because gene expression data<br />
usually have less data points than variables (genes), ODE<br />
approaches are often combined with interpolation and/or<br />
dimension reduction techniques (PCA). A current method<br />
that combines ODE with both interpolation and dimension<br />
reduction techniques is Time Series Network<br />
Identification (TSNI) (Bansal et al., 2006).<br />
Here, we present Dose-Time Network Identification<br />
(DTNI), a method that extends TSNI by including ODE<br />
that describe changes in gene expression over dose in<br />
relation to each other and to time. We also adapted the<br />
original method so that it can include data from multiple<br />
perturbations (compounds).<br />
RESULTS & DISCUSSION<br />
By exploiting simulated data, we show that including<br />
ODE for expression changes over dose leads to improved<br />
GRN identification compared with including only ODE<br />
that describe changes over time. Furthermore, we show<br />
that DTNI performs better when including data from<br />
multiple perturbations (compounds) than when applying<br />
DTNI to data from a single perturbation. This suggests<br />
that the method is suitable to infer a GRN affected by<br />
compounds with the same mode of action. As an example,<br />
we infer the network affected by COX inhibitors from<br />
public microarray data of 6 COX inhibitors, measured in<br />
human hepatocytes, available from Open TG-Gates<br />
(http://toxico.nibio.go.jp/english/index.html) (Noriyuki et<br />
al., 2012). The interactions in the inferred network were<br />
compared to interactions from ConsensusPathDB, a<br />
database including interactions from 32 different sources<br />
(Kamburov et al., 2013). The inferred network was<br />
validated by leave-one out cross-validation (LOOCV). Six<br />
datasets were created from the original data by leaving out<br />
the data of one compound. The network constructed from<br />
the whole data set showed large overlap with the networks<br />
constructed from each of the LOOCV datasets. Edges in<br />
the network constructed from the whole data set, but not in<br />
the networks constructed from the LOOCV datasets were<br />
removed from the network. The remaining novel<br />
interactions, i.e. those that are not in ConsensusPathDB,<br />
have to be validated experimentally, e.g. by geneknockdown<br />
experiments.<br />
FIGURE 1. Workflow for identifying a gene regulatory network affected<br />
by a group of compounds with the same mode of action.<br />
REFERENCES<br />
Bansal M et al. Bioinformatics 22, 815-822 (2006).<br />
Noriyuki N et al. J Toxicol Sci 37,791-801 (2012).<br />
Kamburov A et al. Nucl Acids Res 41, D793-D800 (2013).<br />
68