Establishing large-scale gene regulatory networks using a gene-knowledge-embedded evolutionary computation method

With the increased availability of DNA microarray time-series data, establishing large-scale gene regulatory networks (GRNs) plays an important role in studying systems biology. S-system is well recognized as a suitable model to reconstruct the kinetic model of GRNs. However, the inference of an S-system model of N-gene genetic networks has 2N(N+1) parameters in a set of nonlinear differential equations to be optimized. The existing method iTEA has shown to be effective in reconstructing GRNs. In this work, an improved version of iTEA is proposed (named iTEA2) for efficiently establishing large-scale GRNs by incorporating the domain knowledge of gene regulation into the proposed evolutionary computation method. High performance of iTEA2 arises mainly from the novel method of encoding chromosomes in using the intelligent genetic algorithm. Let I be the maximal number of genes that directly regulated each gene. iTEA2 uses a hybrid encoding method which consists of regulation strength, gene number regulated, and binary control parameters in a chromosome. The value of the strength parameter is the kinetic order and control parameters indicate that the up- and down-regulated kinetic orders are active or not. The total number of parameters encoded in a chromosome for one gene is 5I+2, which is independent upon N. The experimental results show that iTEA2 is significantly better than iTEA in terms of accuracy, convergence speed, and robustness in establishing large-scale GRN.

[1]  Shinn-Ying Ho,et al.  Intelligent evolutionary algorithms for large parameter optimization problems , 2004, IEEE Trans. Evol. Comput..

[2]  Feng-Sheng Wang,et al.  Inference of biochemical network models in S-system using multiobjective optimization approach , 2008, Bioinform..

[3]  Araceli M. Huerta,et al.  From specific gene regulation to genomic networks: a global analysis of transcriptional regulation in Escherichia coli. , 1998, BioEssays : news and reviews in molecular, cellular and developmental biology.

[4]  Shuhei Kimura,et al.  Inference of S-system models of genetic networks using a cooperative coevolutionary algorithm , 2005, Bioinform..

[5]  W. Marsden I and J , 2012 .

[6]  Chih-Hung Hsieh,et al.  An Intelligent Two-Stage Evolutionary Algorithm for Dynamic Pathway Identification From Gene Expression Profiles , 2007, IEEE/ACM Transactions on Computational Biology and Bioinformatics.

[7]  Shinn-Ying Ho,et al.  OSA: orthogonal simulated annealing algorithm and its application to designing mixed H2/H∞ optimal controllers , 2004, IEEE Trans. Syst. Man Cybern. Part A.

[8]  Takanori Ueda,et al.  Inference of Genetic Network Using the Expression Profile Time Course Data of Mouse P19 Cells , 2002 .

[9]  Masaru Tomita,et al.  Dynamic modeling of genetic networks using genetic algorithm and S-system , 2003, Bioinform..

[10]  S. Kimura,et al.  Inference of S-system Models of Genetic Networks from Noisy Time-series Data , 2004 .

[11]  Jonas S. Almeida,et al.  Decoupling dynamical systems for pathway identification from metabolic profiles , 2004, Bioinform..

[12]  Feng-Sheng Wang,et al.  Evolutionary optimization with data collocation for reverse engineering of biological networks , 2005, Bioinform..

[13]  Jonas S. Almeida,et al.  Parameter optimization in S-system models , 2008, BMC Systems Biology.

[14]  P. Brazhnik,et al.  Gene networks: how to put the function in genomics. , 2002, Trends in biotechnology.

[15]  E. Winzeler,et al.  Genomics, gene expression and DNA arrays , 2000, Nature.