Bioinformatics Advance Access originally published online on July 4, 2006
Bioinformatics 2006 22(17):2107-2113; doi:10.1093/bioinformatics/btl361
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Probe-level measurement error improves accuracy in detecting differential gene expression
1 School of Computer Science, University of Manchester Oxford Road, Manchester M13 9PL, UK
2 Department of Computer Science, University of Sheffield Regent Court 211 Portobello Street, Sheffield S1 4DP, UK
*To whom correspondence should be addressed.
Motivation: Finding differentially expressed genes is a fundamental objective of a microarray experiment. Numerous methods have been proposed to perform this task. Existing methods are based on point estimates of gene expression level obtained from each microarray experiment. This approach discards potentially useful information about measurement error that can be obtained from an appropriate probe-level analysis. Probabilistic probe-level models can be used to measure gene expression and also provide a level of uncertainty in this measurement. This probe-level measurement error provides useful information which can help in the identification of differentially expressed genes.
Results: We propose a Bayesian method to include probe-level measurement error into the detection of differentially expressed genes from replicated experiments. A variational approximation is used for efficient parameter estimation. We compare this approximation with MAP and MCMC parameter estimation in terms of computational efficiency and accuracy. The method is used to calculate the probability of positive log-ratio (PPLR) of expression levels between conditions. Using the measurements from a recently developed Affymetrix probe-level model, multi-mgMOS, we test PPLR on a spike-in dataset and a mouse time-course dataset. Results show that the inclusion of probe-level measurement error improves accuracy in detecting differential gene expression.
Availability: The MAP approximation and variational inference described in this paper have been implemented in an R package pplr. The MCMC method is implemented in Matlab. Both software are available from http://umber.sbs.man.ac.uk/resources/puma
Contact: magnus.rattray{at}manchester.ac.uk
Supplementary Information: Supplementary data are available at Bioinformatics Online.
Received on March 30, 2006; revised on June 5, 2006; accepted on June 26, 2006
This article has been cited by other articles:
![]() |
D. G. Morris, S. M. Waters, S. D. McCarthy, J. Patton, B. Earley, R. Fitzpatrick, J. J. Murphy, M. G. Diskin, D. A. Kenny, A. Brass, et al. Pleiotropic effects of negative energy balance in the postpartum dairy cow on splenic gene expression: repercussions for innate and adaptive immunity Physiol Genomics, September 1, 2009; 39(1): 28 - 37. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. M. Packham, C. Gray, P. R. Heath, P. G. Hellewell, P. W. Ingham, D. C. Crossman, M. Milo, and T. J. A. Chico Microarray profiling reveals CXCR4a is downregulated by blood flow in vivo and mediates collateral formation in zebrafish embryos Physiol Genomics, August 7, 2009; 38(3): 319 - 327. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Geeleher, D. Morris, J. P. Hinde, and A. Golden BioconductorBuntu: a Linux distribution that implements a web-based DNA microarray analysis server Bioinformatics, June 1, 2009; 25(11): 1438 - 1439. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Jayawardhana, D. B. Kell, and M. Rattray Bayesian inference of the sites of perturbations in metabolic pathways via Markov chain Monte Carlo Bioinformatics, May 1, 2008; 24(9): 1191 - 1197. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. B. Silkworth, E. A. Carlson, C. McCulloch, K. Illouz, S. Goodwin, and T. R. Sutter Toxicogenomic Analysis of Gender, Chemical, and Dose Effects in Livers of TCDD- or Aroclor 1254-Exposed Rats Using a Multifactor Linear Model Toxicol. Sci., April 1, 2008; 102(2): 291 - 309. [Abstract] [Full Text] [PDF] |
||||


