Bioinformatics Advance Access originally published online on February 10, 2004
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Bioinformatics 20(8) © Oxford University Press 2004; all rights reserved.
Bayesian mixture model based clustering of replicated microarray data
1 Department of Environmental Health, Center for Genome Information, University of Cincinnati Medical Center, 3223 Eden Avenue ML 56, Cincinnati, OH 45267-0056, USA and 2 Department of Microbiology, Box 358070, University of Washington, Seattle, WA 98195, USA
Received on July 14, 2003; revised on November 4, 2003; accepted on November 5, 2003
Advance Access Publication February 10, 2004
Motivation: Identifying patterns of co-expression in microarray data by cluster analysis has been a productive approach to uncovering molecular mechanisms underlying biological processes under investigation. Using experimental replicates can generally improve the precision of the cluster analysis by reducing the experimental variability of measurements. In such situations, Bayesian mixtures allow for an efficient use of information by precisely modeling between-replicates variability.
Results: We developed different variants of Bayesian mixture based clustering procedures for clustering gene expression data with experimental replicates. In this approach, the statistical distribution of microarray data is described by a Bayesian mixture model. Clusters of co-expressed genes are created from the posterior distribution of clusterings, which is estimated by a Gibbs sampler. We define infinite and finite Bayesian mixture models with different between-replicates variance structures and investigate their utility by analyzing synthetic and the real-world datasets. Results of our analyses demonstrate that (1) improvements in precision achieved by performing only two experimental replicates can be dramatic when the between-replicates variability is high, (2) precise modeling of intra-gene variability is important for accurate identification of co-expressed genes and (3) the infinite mixture model with the ‘elliptical’ between-replicates variance structure performed overall better than any other method tested. We also introduce a heuristic modification to the Gibbs sampler based on the ‘reverse annealing’ principle. This modification effectively overcomes the tendency of the Gibbs sampler to converge to different modes of the posterior distribution when started from different initial positions. Finally, we demonstrate that the Bayesian infinite mixture model with ‘elliptical’ variance structure is capable of identifying the underlying structure of the data without knowing the ‘correct’ number of clusters.
Availability: The MS WindowsTM based program named Gaussian Infinite Mixture Modeling (GIMM) implementing the Gibbs sampler and corresponding C++ code are available at http://homepages.uc.edu/~medvedm/GIMM.htm
Supplemental information: http://expression.microslu.washington.edu/expression/kayee/medvedovic2003/medvedovic_bioinf2003.html
Contact: Mario.Medvedovic{at}uc.edu
* To whom correspondence should be addressed.
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  F. Achcar, J.-M. Camadro, and D. Mestivier AutoClass@IJM: a powerful tool for Bayesian classification of heterogeneous data in biology Nucleic Acids Res., July 1, 2009; 37(suppl_2): W63 - W67. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 C.-T. Li, Y. Yuan, and R. Wilson An unsupervised conditional random fields approach for clustering gene expression time series Bioinformatics, November 1, 2008; 24(21): 2467 - 2473. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 D. B. Dunson and J.-H. Park Kernel stick-breaking processes Biometrika, June 1, 2008; 95(2): 307 - 323. [Abstract] [PDF]
|
|
|
|
|
|
L. Brehelin, O. Gascuel, and O. Martin Using repeated measurements to validate hierarchical gene clusters Bioinformatics, March 1, 2008; 24(5): 682 - 688. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
A. Joshi, Y. Van de Peer, and T. Michoel Analysis of a Gibbs sampler method for model-based clustering of gene expression data Bioinformatics, January 15, 2008; 24(2): 176 - 183. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 W. D. Hardie, T. R. Korfhagen, M. A. Sartor, A. Prestridge, M. Medvedovic, T. D. Le Cras, M. Ikegami, S. C. Wesselkamper, C. Davidson, M. Dietsch, et al. Genomic Profile of Matrix and Vasculature Remodeling in TGF-{alpha} Induced Pulmonary Fibrosis Am. J. Respir. Cell Mol. Biol., September 1, 2007; 37(3): 309 - 321. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
D. Zhu, Y. Li, and H. Li Multivariate correlation estimator for inferring functional relationships from replicated genome-wide data Bioinformatics, September 1, 2007; 23(17): 2298 - 2305. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
Z. Xiang, Z. S. Qin, and Y. He CRCView: a web server for analyzing and visualizing microarray gene expression data using model-based clustering Bioinformatics, July 15, 2007; 23(14): 1843 - 1845. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
Y. Lu, X. He, and S. Zhong Cross-species microarray analysis with the OSCAR system suggests an INSR->Pax6->NQO1 neuro-protective pathway in aging and Alzheimer's disease Nucleic Acids Res., July 13, 2007; 35(suppl_2): W105 - W114. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
A. Thalamuthu, I. Mukhopadhyay, X. Zheng, and G. C. Tseng Evaluation and comparison of gene clustering methods in microarray analysis Bioinformatics, October 1, 2006; 22(19): 2405 - 2412. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
Z. S. Qin Clustering microarray gene expression data using weighted Chinese restaurant process Bioinformatics, August 15, 2006; 22(16): 1988 - 1997. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
X. Liu, S. Sivaganesan, K. Y. Yeung, J. Guo, R. E. Bumgarner, and M. Medvedovic Context-specific infinite mixtures for clustering gene expression profiles across diverse microarray dataset Bioinformatics, July 15, 2006; 22(14): 1737 - 1744. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 M. Rattray, X. Liu, G. Sanguinetti, M. Milo, and N. D. Lawrence Propagating uncertainty in microarray data analysis Brief Bioinform, March 1, 2006; 7(1): 37 - 47.
|
|
|
|
|
|
D. de Ridder, F. J. T. Staal, J. J. M. van Dongen, and M. J. T. Reinders Maximum significance clustering of oligonucleotide microarrays Bioinformatics, February 1, 2006; 22(3): 326 - 331. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
M. A. Sartor, A. M. Zorn, J. A. Schwanekamp, D. Halbleib, S. Karyala, M. L. Howell, G. E. Dean, M. Medvedovic, and C. R. Tomlinson A new method to remove hybridization bias for interspecies comparison of global gene expression profiles uncovers an association between mRNA sequence divergence and differential gene expression in Xenopus Nucleic Acids Res., January 5, 2006; 34(1): 185 - 200. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 S. Kann, C. Estes, J. F. Reichard, M.-y. Huang, M. A. Sartor, S. Schwemberger, Y. Chen, T. P. Dalton, H. G. Shertzer, Y. Xia, et al. Butylhydroquinone Protects Cells Genetically Deficient in Glutathione Biosynthesis from Arsenite-Induced Apoptosis Without Significantly Changing Their Prooxidant Status Toxicol. Sci., October 1, 2005; 87(2): 365 - 384. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 A. Reverter, Y. H. Wang, K. A. Byrne, S. H. Tan, G. S. Harper, and S. A. Lehnert Joint analysis of multiple cDNA microarray studies via multivariate mixed models applied to genetic improvement of beef cattle J Anim Sci, December 1, 2004; 82(12): 3430 - 3439. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 K. K. Lin, D. Chudova, G. W. Hatfield, P. Smyth, and B. Andersen Identification of hair cycle-associated genes from time-course gene expression profile data by using replicate variance PNAS, November 9, 2004; 101(45): 15955 - 15960. [Abstract] [Full Text] [PDF]
|
|