Bioinformatics Advance Access published online on February 13, 2009
Bioinformatics, doi:10.1093/bioinformatics/btp028
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Gene network reconstruction from transcriptional dynamics un-der kinetic model uncertainty: A case for the second derivative
1Ottawa Institute of Systems Biology; 2 Department of Biochemistry, Microbiology, and Immunology; University of Ottawa; 451 Smyth Road; Ottawa, Ontario K1H 8M5
3 Graduate Institute of Systems Biology and Bioinformatics; National Central University; No. 300, Jhongda Road, Jhongli City, Taoyuan County 32001, Taiwan (R.O.C.)
4 Pioneer Hi-Bred International, Inc.; 7300 NW 62nd Avenue; PO Box 1004; Johnston, Iowa 50131-1004
*To whom correspondence should be addressed. Dr. David R. Bickel, E-mail: dbickel{at}uottawa.ca
| Abstract |
|---|
Motivation: Measurements of gene expression over time enable the reconstruction of transcriptional networks. However, Bayesian net-works and many other current reconstruction methods rely on as-sumptions that conflict with the differential equations that describe transcriptional kinetics. Practical approximations of kinetic models would enable inferring causal relationships between genes from expression data of microarray, tag-based, and conventional plat-forms, but conclusions are sensitive to the assumptions made.
Results: The representation of a sufficiently large portion of genome enables computation of an upper bound on how much confidence one may place in influences between genes on the basis of expression data. Information about which genes encode transcription factors is not necessary but may be incorporated if available. The methodology is generalized to cover cases in which expression measurements are missing for many of the genes that might control the transcription of the genes of interest. The assumption that the gene expression level is roughly proportional to the rate of translation led to better empirical performance than did either the assumption that the gene expression level is roughly proportional to the protein level or the Bayesian model average of both assumptions.
Contact: dbickel{at}uottawa.ca
Availability: {{http://www.oisb.ca}} points to R code implementing the methods (R Development Core Team 2004).
Supplementary Information: {{http://www.davidbickel.com}}
Associate Editor: Prof. Thomas Lengauer
Received on May 9, 2008; revised on December 4, 2008; accepted on January 12, 2008