03.12.2015 Views

bbc 2015

BBC2015_booklet

BBC2015_booklet

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!