Identifying a bicluster, or submatrix of a gene expression dataset wherein the genes express similar behavior over the columns, is useful for discovering novel functional gene interactions. In this article, we introduce a new algorithm for finding biClusters with Linear Patterns (CLiP). Instead of solely maximizing Pearson correlation, we introduce a fitness function that also considers the correlation of complementary genes and conditions. This eliminates the need for a priori determination of the bicluster size. We employ both greedy search and the genetic algorithm in optimization, incorporating resampling for more robust discovery. When applied to both real and simulation datasets, our results show that CLiP is superior to existing methods. In analyzing RNA-seq fly and worm time-course data from modENCODE, we uncover a set of similarly expressed genes suggesting maternal dependence. Supplementary Material is available online (at www.liebertonline.com/cmb).
Get full access to this article
View all access options for this article.
References
1.
Aguilar-RuizJ.2005. Shifting and scaling patterns from gene expression data. Bioinformatics, 21:3840–3845.
2.
BanfieldJ., RafteryA.1993. Model-based Gaussian and non-Gaussian clustering. Biometrics, 49:803–821.
3.
Ben-DorA., ChorB., KarpR.2003. Discovering local structure in gene expression data: the order-preserving submatrix problem. J. Comput. Biol., 10:803–821.
4.
BhattacharyaA., DeR.2009. Bi-correlation clustering algorithm for determining a set of co-regulated genes. Bioinformatics, 25:2795–2801.
5.
BusyginS., JacobsenG., KramerE.2002. Double conjugated clustering applied to leukemia microarry data. Proc. 2nd SIAM ICDM Workshop Clustering High Dimensional Data.
6.
CalifanoA., StolovitzkyG., TuY.2004. Analysis of gene expression microarrays for phenotype classification. Proc. Int. Conf. Intel. Syst. Mol. Biol., 194:1625–1638.
7.
ChengY., ChurchG.2000. Biclustering of expression data. Proc. Int. Cof. Intell. Syst. Mol. Biol., 8:93–103.
8.
FraleyC., RafteryA.2002. Model-based clustering, discriminant analysis and density estimation. J. Am. Stat. Assoc., 97:611–631.
9.
GanX., LiewA., YanH.2005. Biclustering gene expression data based on a high dimensional geometric method. Proc. Int. Conf. Mach. Learn. Cybernet., 3388–3393.
10.
GarX., LiewA., YanH.2008. Discovering biclusters in gene expression data based on high-dimensional linear geometries. BMC Bioinform., 9:209.
11.
GersteinM.B.et al.2010. Integrative analysis of the Caenorhabditis elegans genome by the modENCODE Project. Science, 330:1775–1787.
12.
GetzG., LevineE., DomanyE.2000. Coupled two-way clustering analysis of gene microarray data. Proc. Natl. Acad. Sci. USA, 97:12079–12084.
13.
HartiganJ.1972. Direct clustering of a data matrix. J. Am. Stat. Assoc., 67:123–129.
14.
HartiganJ.1975. Clustering Algorithms. John Wiley and Sons, Inc.: New York.
HuangD., ShermanB., LempickiR.2009a. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res., 37:1–13.
17.
HuangD., ShermanB., LempickiR.2009b. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc., 4:44–57.
18.
IhmelsJ., FriedlanderG., BergmannS.2002. Revealing modular organization in the yeast transcriptional network. Nat. Genet., 31:370.
NepomucenoJ., TroncosoA., Aguilar-RuizJ.2011. Biclustering of gene expression data by correlation-based scatter search. BioData Mining, 4:3.
35.
PrelicA., BleulerS., ZimmermannP.2006. A systematic comparison and evaluation of biclustering methods for gene expression data. Bioinformatics, 22:1122–1129.
ShengQ., MoreauY., MoorB.2003. Biclustering microarray data by Gibbs sampling. Bioinformatics, 19:ii196–ii205.
38.
SingsonA., MercerK., L'HernaultS.1998. The C. elegans spe-9 gene encodes a sperm transmembrane protein that contains egf-like repeats and is required for fertilization. Cell, 93:71–79.
39.
SpellmanP., SherlockG., ZhangM.et al.1998. Comprehensive identification of cell cycle–regulated genes of the yeast Saccharomyces cerevisiae by microarray hybridization. Mol. Biol. Cell, 9:3273–3297.
40.
TamaryoP., SlonimD., MesirovJ.1999. Interpreting patterns of gene expression with self-organizing maps: methods and application to hematopoietic differentiation. Proc. Natl. Acad. Sci. USA, 96:2907–2912.
The modENCODE Consortium et al. Identification of functional elements and regulatory circuits by Drosophila modENCODE. 2010. Science, 330:1787–1797.
45.
TrapnellC.2010. Transcript assembly and quantification by RNA-seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotech. Mol., 28:511–515.
46.
WangH., WangW., YangJ.et al.2002. Clustering by pattern similarity in large data sets. Proc. ACM SIGMOD Int. Conf. Manage. Data, 394–405.
47.
YangJ., WangH., WangW.et al.2003. Enhanced biclustering on expression data. Proc. 3rd IEEE Int. Symp. Bioinform. BioEng., 321–327.
48.
ZhaoH., LiewA., XieX.et al.2008. A new geometric biclustering algorithm based on the Hough transform for analysis of large-scale microarray data. J. Theoret. Biol., 251:264–274.
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
