Abstract
Background
Biclustering of gene expression data searches for local patterns of gene expression. A bicluster (or a twoway cluster) is defined as a set of genes whose expression profiles are mutually similar within a subset of experimental conditions/samples. Although several biclustering algorithms have been studied, few are based on rigorous statistical models.
Results
We developed a Bayesian biclustering model (BBC), and implemented a Gibbs sampling procedure for its statistical inference. We showed that Bayesian biclustering model can correctly identify multiple clusters of gene expression data. Using simulated data both from the model and with realistic characters, we demonstrated the BBC algorithm outperforms other methods in both robustness and accuracy. We also showed that the model is stable for two normalization methods, the interquartile range normalization and the smallest quartile range normalization. Applying the BBC algorithm to the yeast expression data, we observed that majority of the biclusters we found are supported by significant biological evidences, such as enrichments of gene functions and transcription factor binding sites in the corresponding promoter sequences.
Conclusions
The BBC algorithm is shown to be a robust modelbased biclustering method that can discover biologically significant genecondition clusters in microarray data. The BBC model can easily handle missing data via Monte Carlo imputation and has the potential to be extended to integrated study of gene transcription networks.
Background
Clustering gene expression data has been an important problem in computational biology. While traditional clustering methods, such as hierarchical and Kmeans clustering, have been shown useful in analyzing microarray data, they have some limitations. First, a gene or an experimental condition can be assigned to only one cluster. Second, all genes and conditions have to be assigned to clusters. However, biologically a gene or a sample could participate in multiple biological pathways, and a cellular process is generally active only under a subset of genes or experimental conditions. A biclustering scheme that produces gene and condition/sample clusters simultaneously can model the situation where a gene or a condition is involved in several biological functions. Furthermore, a biclustering model can avoid those “noise” genes that are not active in any experimental condition.
Biclustering of microarray data was first introduced by Cheng and Church [1]. They defined a residual score to search for submatrices as biclusters. This is a heuristic method and can not model the cases where two biclusters overlap with each other. Segal et al. [2] proposed a modified version of oneway clustering using a Bayesian model in which genes can belong to multiple clusters or none of the clusters. But it can not simultaneously cluster conditions/samples. Tseng and Wong developed a tight clustering algorithm [3]. It allows some of the genes not to be clustered, but does not select conditions. Bergmann et al [4] introduced the iterative signature algorithm (ISA), which searches bicluster modules iteratively based on two predetermined thresholds. ISA can identify multiple biclusters, but is highly sensitive to the threshold values and tends to select a strong bicluster many times. The plaid model [5] introduces a statistical model assuming that the expression value in a bicluster is the sum of the main effect, the gene effect, the condition effect, and the noise term, i.e.:
where noise ε_{ij} ~ N(0, σ^{2}). It further assumes that the expression values of two overlapping biclusters are the sum of the two module effects. The plaid model uses a greedy search strategy, so errors can accumulate easily. Also in multiple clusters case, the clusters identified by the algorithm tend to overlap to a great extent. Tanay et al. [6] proposed a SAMBA biclustering scheme using bipartite graphs containing both conditions and genes. BenDor et al. [7] attempted to identify orderpreserving sub matrices (OPSMs). Murali and Kasif [8] discretized gene expression data into several symbols and searched for conservative symbol motifs (xMOTIFs). A survey of different biclustering methods can be found in [9].
We here propose a Bayesian biclustering (BBC) model. For a single bicluster, we assume the same model as in the plaid model [5], as described in equation (1). But for multiple clusters, we constrain the overlapping of biclusters to only one direction (i.e., either gene or condition direction). Besides, we use a more flexible error model, which allows the error term of each cluster to have to a different variance. To make the Bayesian inference of biclusters, we implemented an efficient Gibbs sampling algorithm with all effect parameters (except the error variances) integrated out. We compared the performance of the BBC algorithm for several different types of simulated datasets with that of the plaid model [5], the ISA [4], the method of Cheng and Church [1], the SAMBA method [6] and the OPSMs [7]. Finally, we applied the BBC algorithm to the yeast expression dataset and identified many biologically significant biclusters.
Results and discussion
Simulation results
Bayesian biclustering in various simulated scenarios
We simulated a dataset with 400 genes and 50 samples. The background data is i.i.d. from N(0, 0.5). Two clusters of 100 genes and 15 conditions are simulated according to the BBC model with main effects, gene effects, condition effects and error terms as μ_{1} ~ N(5, 0.5), μ_{2} ~ N(7, 0.5), α_{i1}, α_{i2} ~ N(0, 0.5), β_{j1}, β_{j2} ~ N(0, 0.5) and ε_{ij1} ~ N(0, 0.5), ε_{ij2} ~ N(0, 0.7).
We considered three scenarios for datasets with two clusters: the two clusters have some common conditions but distinct genes (Figure 1(a)); the two clusters have some common genes but distinct conditions ( Figure 1(d)); and two clusters have both common genes and conditions (Figure 1(h)), in which case an additive model is assumed for the overlapping part. The results from using a nonoverlapping gene version of the BBC model are shown in Figures 1(b)(c),(e)(g),(i)(k). In all cases the BBC model identified the genes and conditions of the simulated clusters correctly, but grouped them slightly differently because of our model constraints.
Figure 1. Simulated data with two biclusters and the results of the BBC analysis. Bayesian biclustering for simulated datasets. (a) A dataset with two nonoverlapping clusters. (b)(c) The two clusters found by the Bayesian biclustering model from (a). (d) A dataset with two clusters with common genes. (e)(g) The three clusters found by the Bayesian biclustering model from (d). (h) A dataset with two clusters with both common samples and common genes. (i)(k) The three clusters found by the Bayesian biclustering model from (h).
Comparison of biclustering algorithms on data simulated from statistical models
We compared six biclustering methods: the BBC method, the plaid model, ISA, SAMBA, OPSMs, and Cheng and Church's biclustering (CC). We considered both the single cluster case and the multiple clusters case using simulated data from the plaid model. A single cluster dataset is shown in Figure 2(a). The 400 × 50 background noise matrix is simulated according to i.i.d. normal N(0,0.5). We superimposed a cluster of size 100 × 20 according to the plaid model with μ_{1} ~ N(5, 0.5) and α_{i1}, β_{j1} ~ N(0, 0.5). The multiple cluster case is shown in Figure 2(b). The background is the same as above. Two clusters of size 100 × 15 are also simulated according to the plaid model with μ_{1} ~ N(5,0.5), μ_{2} ~ N(7, 0.5) and α_{i1}, β_{ji}, α_{i2}, β_{j2} ~ N(0, 0.5). An additive model is used for the overlapping part of the two clusters.
Figure 2. Datasets simulated according to the plaid model Datasets for comparison. (a) A dataset with one single cluster (b) A dataset with two clusters, of which both genes and samples overlap.
Since each method searches for biclusters with different structures, comparing biclustering results is not very straightforward. In order to carry out a comprehensive comparison among various biclustering results for simulated datasets, we use the following four characteristics: sensitivity, specificity, overlapping rate, and number of clusters. Since we know which genecondition combination belongs to the true biclusters, we use the standard definition for sensitivity and specificity, both of which are values between 0 and 1. A higher sensitivity suggests that more “true” members of the clusters have been identified by the algorithm, while a higher specificity suggests that more background data points are excluded from the clusters. The overlapping rate is defined as
Thus, if there is no overlap between the identified clusters, the overlapping rate is 0. On the other hand, if the identified clusters greatly overlap with each other, the overlapping rate is close to 1.
We used the BicAT software package [10] for ISA, CC, and OPSMs. Different gene and condition thresholds are used for the ISA. We carefully chose a set of thresholds with good performance and then slightly changed the thresholds to test the stability of the ISA. We used default settings for CC's model. The plaid algorithm was implemented using the Plaid package [5]. The Expander package [11] is used for SAMBA biclustering. The results are shown in Table 1, where the lefthandside value in each entry is for the single cluster case, and the righthandside value is for the twoclusters case.
Table 1. Biclustering results of different methods for simulated data using the plaid model
It can be seen that the ISA method is very sensitive to the choice of the thresholds. The performance of ISA also degrades in the case of multiple overlapping clusters. The SAMBA method and the OPSMs method correctly identified almost all background noises, but tends to exclude some meaningful patterns. The CC method includes too much background data in clusters. The plaid model performs well in the single cluster case. But it identifies too many overlapped clusters in the multiple clusters case. Our BBC method performs well in both cases, even though the data generation model for the overlapping part in the second case does not satisfy the BBC model assumption.
Comparison of biclustering algorithms on data simulated with biological characteristics
People are mostly interested in how different biclustering methods perform for real microarray datasets. We next carry out a comparison using simulated microarray datasets with realistic characteristics. As shown in Figure 3, the dataset has 100 genes and 50 conditions. The concentration of three transcription factors (TF) are changed across conditions. We assume that gene 21 to gene 40 are transcribed when both TF1 and TF2 are active, with TF1 serving as an activator and TF2 as an inhibitor. We also assume that gene 41 to gene 60 are transcribed when both TF1 and TF3 are active, where TF1 serves as an inhibitor, TF3 an activator. The expression values are simulated using the biochemistry model presented in [12]. If gene i is regulated by M_{I} inhibitors and M_{A} activators, then
Figure 3. The Simulated dataset with realistic characters
where G_{i} represents the abundance of the mRNA of gene i, s(I_{1}, I_{2} …, I_{MI}, A_{1}, A_{2}, … A_{MA}) is a rate law representing mRNA synthesis, I_{1}, I_{2} …, I_{MI} are inhibitor concentrations, A_{1}, A_{2}…, A_{MA} are activator concentrations, and b(G_{i}) is mRNA breakdown rate. In [12] mRNA synthesis rate is modelled as
where V_{basal} is the basal transcription rate, constants Ki_{j} and Ka_{k} are concentrations at which the effect of the inhibitor or activator is half of its saturating value. The exponents n_{i} and n_{k} regulate the sigmoidicity of the transcription rate curve. We set n_{i} = n_{k} = 1.5, and randomly simulated Ki_{j} and Ka_{k} for the dataset.
We added real noise from the well known Leukemia expression dataset [13]. We first obtained a noise data matrix using all scattered (“noise”) genes excluded by the tight clustering algorithm [3]. Then, we chose 100 rows and 50 columns at random. We also scaled the noise to control signal to noise ratio (SNR). Both the good data quality case (SNR=10) and the bad data quality case (SNR=4) are considered. We simulated 10 datasets for both cases and the average value of each characteristics is shown in Table 2, where the lefthandside value in each entry is for SNR=10 case, and the righthandside value is for SNR=4 case. we chose the threshold values for the ISA model in the same way as in the previous simulation study.
Table 2. Biclustering results of different methods for simulated data with realistic characteristics
The BBC model performed the best among these methods. Again the ISA method was sensitive to thresholds. It also had some false positives. The OPSMs method missed most of the significant patterns. The SAMBA method found some small and tight biclusters of genes and conditions, but also excluded many significant patterns. CC's method misidentified many noisy data points as biclusters. The plaid model recognized almost all significant patterns, but its specificity was low. Interestingly, the plaid model gave better results for low SNR case, which was due to the fewer number of clusters found by the plaid model.
Effects of normalization for Bayesian biclustering model
Data normalization is an important step for microarray analysis. Although some clustering methods such as ISA incorporate the normalization step in their procedures, most clustering methods work on normalized microarray data. The BBC model belongs to the latter. Since the normalization procedure greatly changes the microarray data, different normalization procedures may lead to very different clustering results.
We conducted a study on how normalization methods affect the biclustering results. Five normalization procedures including column standardization (CSN), row standardization (RSN), quantile normalization on gene level (QNGL), the interquartile range normalization (IQRN), and the smallest quartile range normalization (SQRN) were considered. In CSN (or RSN), each column (or row) is recentered and rescaled, so that the sample mean of each column (or row) becomes 0, and the sample variance becomes 1. These are quite crude methods, but are still used in many clustering applications. QNGL used the same technique with quantile normalization [14], which was developed for normalizing oligonucleotide arrays on probe levels. In QNGL, we apply the quantile normalization on gene level for the simulated array data.
IQRN and SQRN are two new methods we propose here. They are inspired by CSN, but are more robust to outliers. In IQRN, one first sorts the data in each column, trims off α/2% of the data from each tail, and computes the α%trimmed mean and standard deviation. Then, all data in that column is standardized by subtracting the trimmed mean and being divided by the trimmed standard deviation. This normalization method can reduce the artificial normalization effect caused by outliers. In SQRN, instead of using the middle (100α)% of the data, one first finds for each column the shortest interval that contains a certain percentage (e.g., 50%) of the data. Then the data of that column is standardized by the sample mean and variance of the data inside the shortest quartile range. If distributions of the data in each column are symmetric and unimodal, then SQRN is equivalent to IQRN. But SQRN gives better results for skewed distributions. We applied the above five normalization methods on the same simulated dataset as in Figure 2(b) before applying our BBC procedure.
As shown in Table 3, both IQRN and SQRN performed very well, whereas the other three methods affected the clustering results significantly. Thus, in our yeast data analysis reported in the following section, we used IQRN before applying the BBC procedure.
Table 3. Comparison of normalization methods for Bayesian Biclustering Model
Bayesian biclustering for yeast datasets
We analyzed the same yeast expression data as in [15] using the BBC procedure. This dataset was derived by combining the environmental stress data of [16] and the cell cycle data of [17]. The combined dataset contains 6108 genes and 250 conditions (or samples). We applied the 90% IQRN procedure across all conditions. Since the dataset contains many missing data, we imputed them along with our BBC iterations. The BBC algorithm was asked to search for K biclusters, with K ranging from 30 to 65. We observed that the BIC [18] achieved the optimal value with K = 57. Out of 6108 genes, 6021 were included in one of the clusters, and all conditions were included in at least one cluster.
We analyzed the clustering results from three aspects. First, we identified the significant categories of experimental conditions for each cluster. More precisely, we classified the 250 experimental conditions into 22 categories according to the biological nature of each experiment. Some examples of categories are heat shock stress, amino acid starvation, and α factor synchronization. Then we searched for significant enrichment for each category. Second, we did functional enrichment test of genes in each cluster using functional information from the MIPS database. Third, we searched the promoter sequences (up to 800 bps upstream) of genes in each bicluster for the enrichment of transcription factor binding sites (TFBS). We applied the TFBS scores of 51 known yeast transcription factors used in [15], which measures how likely a promoter sequence contains a TFBS and was computed using ScanACE [19]. A cutoff value of 0.5 was used to make a presence/absence call for a TFBS. The presence frequency of each TFBS in all genes in the dataset was used as null hypothesis. For all the three types of enrichment analysis, we used the criterion of P value< 0.05 after the Bonferroni correction.
Out of the 57 clusters, 36 have significant gene functions enrichments, 26 have significant TFBS enrichments, 51 have significant experimental condition categories enrichments, 22 have all three types of enrichments and 54 have at least one type of enrichment. We named a few of these biclusters and listed them in Table 4. The enriched gene functions and the significant experimental conditions showed strong correlations. For example, in the cell cycle cluster, the enriched gene function terms are cell cycle, DNA processing and cell division, and the significant experimental conditions are all cellcycle experiments including cln3 and clb2 experiments, α factor, cdc15, cdc28 and elutriation synchronization. Another example is stress response and protein folding cluster, the enriched gene function terms are stress response and protein folding, and the significant conditions are heat shock, diamide and osmolarity stress experiments. The biclustering results are also supported by the TFBS information. In many clusters, the enriched TFBSs correspond to TFs known to be involved in the pathways and biological functions that were found significant from gene functional enrichment analysis. For instance, in the ribosome protein bicluster, the enriched TFBS RAP1 is known to be involved in ribosome protein transcription [20]. The ubiquitin cluster is enriched with TFBS of RPN4, which was shown to be part of the ubiquitin fusion degradation pathway [21]. The nitrogen, sulfur and selenium cluster shows significant overrepresentation of binding sites for TFs CBF1, GCN4, MET31, and MET4. CBF1 is known to induce sulfate assimilation pathway along with MET4 [22], GCN4 is an activator involved with protein biosynthesis, and MET31 is a known transcriptional regulator of sulfur amino acid metabolism [23]. The G1 phase cluster is enriched with MBP1 and SWI4 binding sites. MBP1 and SWI4 are known to act together to regulate late G1specific transcription of targets and genes for DNA synthesis [24]. The oxidative stress cluster is enriched by CAD1 and YAP1 binding sites, where YAP1 activates the transcription of antioxidant genes in response to oxidative stress [25]. The glycolysis regulation cluster is enriched by TFBS GCR1, which is known to be involved in glycolysis [26].
Table 4. Bayesian Biclustering results for yeast expression data
Conclusions
We have presented a rigorous hierarchical Bayes model for clustering microarray data in both the gene and the experimental condition directions. We used Gibbs sampling and Bayesian information criterion to identify biclusters as well as the total number of clusters. Using simulated datasets, we showed that the BBC algorithm outperformed other clustering methods especially when multiple clusters were present. Moreover, the BBC method performed the best for simulated data based on biochemistry models with realistic noise background. We also discussed the impact of normalization procedures on the clustering results, and found that both the interquartile range normalization and the smallestquartile range normalization are robust for our BBC model. When applied to a wellknown yeast microarray dataset, the BBC procedure discovered many biologically significant clusters, from which significant enrichments of gene functions, associated experimental conditions, and related TFBS enrichment were found.
Unlike many other biclustering methods, the BBC is completely modelbased and does not need to fine tune any threshold parameters. Because it is a full Bayesian model, the BBC can handle missing data extremely easily, and can also incorporate likelihoodbased criterion, such as AIC, BIC, maximum likelihood, Bayes factors, etc., for model evaluations and comparisons. In addition, the BBC model has the potential to be extended by incorporating other types of data, such as the promoter sequence information into the model.
Methods
Bayesian biclustering model
Consider a microarray dataset with N genes and P conditions (or samples), in which the expression value of the i^{th} gene and j^{th} condition is denoted as y_{ij}, i = 1, 2, • • •, N, j = 1, 2, • • •, P. We assume that
where K is the total number of clusters (unknown), μ_{k} is the main effect of cluster k, and α_{ik} and β_{jk} are the effects of gene i and condition j, respectively, in cluster k, ε_{ijk} is the noise term for cluster k, and e_{ij} models the data points that do not belong to any cluster. Here δ_{ik} and κ_{jk} are binary variables: δ_{ik} = 1 indicates that row (gene) i belongs to cluster k, and δ_{ik} = 0 otherwise; similarly, κ_{jk} = 1 indicates that condition (column) j is in cluster k, and κ_{jk} = 0 otherwise.
When multiple biclusters are allowed, the original plaid model usually finds biclusters greatly overlapping with each other. This effect is quite artificial and is likely due to the nonidentifiability problem caused by the additive assumption made for overlapping clusters. We solve this problem by allowing biclusters to overlap only in one direction, either the gene or the condition direction, but not both. This results in two versions of the BBC model: nonoverlapping gene biclustering and nonoverlapping condition biclustering. In nonoverlapping condition biclustering, a condition can be in one or none of the clusters, but a gene can be assigned to multiple clusters. Mathematically, this constraint can be written as . In nonoverlapping gene biclustering, a condition can be assigned into multiple clusters, while a gene can only be in no more than one cluster. This corresponds to . Note that in either of these two versions, different biclusters do not overlap. Without loss of generality, we focus our discussions on the nonoverlapping gene biclustering in this paper. Thus the priors of the indicators κ and δ are set so that a condition can be in multiple clusters and a gene be in no more than one cluster, i.e.,
k_{i}_{j} ~ Bernoulli(q_{k})
where p_{k} and q_{k} are set to be constant. We tested different values for p_{k} and q_{k}, and found out that different values do not affect the results much. We used q_{k} =0.1 and
for the yeast dataset.
We assume a priori that
The hyperpriors for the are set to be inverse Gamma distributed.
In our model, an observation y_{ij} can belong to either one or none of the biclusters. Thus, we can rewrite the probability distribution of y_{ij} conditional on the cluster indicators. If y_{ij} belongs to cluster k,(k = 1,2,… ,K), then
If y_{ij} belongs to none of the clusters, then
With Gaussian zeromean priors on the effects parameters, we get the marginal distribution of the y_{ij} conditional on the indicators as:
where Σ is the covariance matrix of Y, and Y = (Y_{0}, Y_{1}, Y_{2}, …, Y_{K})^{T} with Y_{k} = {y_{ij} : δ_{ik}κ_{j}_{k} = 1}, k ≥ 1, being the vector composed of the data points belonging to cluster k, and Y_{0} being the vector of data points belonging to no cluster. More specifically, Σ is a sparse matrix of the form
where Σ_{k} = cov(Y_{k}, Y_{k}) is the covariance matrix of all data points belonging to cluster k whose entries are:
Gibbs sampling for Bayesian biclustering
In order to make inference from the BBC model, we implement a Gibbs sampling method [27] to draw samples from the posterior distribution of the indicator variables, which can be derived by combining equation (8) with a prior distribution on the δ and κ. Initializing from a set of randomly assigned values of δ's and κ's, we sample the column (condition) indicators κ by calculating the following logprobability ratio:
Since each data point belongs to no more than one cluster, we can therefore divide data points into two sets given the current parameters except κ_{jk}. The first set contains data points not in cluster k, i.e., V_{1} = {y_{il} : δ_{ik} = 0 or κ_{lk} =0,l ≠ j}. The second set contains data points that are or can in cluster k, i.e., V_{2} = {y_{il} : δ_{ik} = 1,κ_{lk} = 1,l ≠ j}U{y_{ij} : δ_{ik} = 1}.
Two data points are independent if they belong to different clusters, therefore we can write the joint likelihood of Y as a product of the joint likelihood for data in V_{1} and V_{2}, respectively. As a consequence, the logposterior probability ratio can be simplified as
In order to calculate the likelihood term in the above ratio, we need to take the inverse and determinant of the covariance matrices for the vector V_{2} in both cases. The dimensions of these covariance matrices are huge in practice (in the order of thousands), so a brute force calculation would be expensive. Since the covariance matrices have the special structures as shown in equation (9), we can simplify the likelihood ratio term. The final simplified form only involves multiplications and additions of matrices with dimension I_{k} × I_{k}, where I_{k} is the number of genes in cluster k given current parameters.
Similarly we can obtain the log posterior probability ratio for gene indicators δ_{ik},
Since our model requires that , then the gene indicators δ_{ik} are correlated. We thus need to calculate the logposterior probability ratio for every δ_{ik}, k = 1, 2, • • • K, and sample them jointly.
We also sample the effect parameters based on the indicators δ and κ. The gene effects α_{ik} and condition effects β_{jk} can serve as scores for genes and conditions in a cluster. Moreover the BBC model is very convenient and coherent in handling missing data in microarray datasets: just treat them as additional unknown variables and iteratively impute them in the Gibbs sampling iterations. Suppose at step t, we have sampled both the indicator variables and effects parameters, we can impute the missing data, say y_{ij}, by sampling from distribution
In the above procedure, we preset the value K for the total number of clusters. However, the information of K is not available in general. In practice, we search biclusters for a number of K's and select the best K based on the Beyesian information criterion (BIC) [18].
An executable program for the BBC algorithm is available at http://www.people.fas.harvard.edu/~junliu/BBC webcite
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
JG and JSL designed the models and simulation studies together. JG implemented the method and analyzed the data. Both authors contributed to the writing of the manuscript.
Acknowledgements
This work was supported by grants from the National Institute of Health (R01GM078990) and the National Science Foundation (DMS0706989). We thank the four referees and members of the Liu lab for many helpful comments and discussions.
This article has been published as part of BMC Genomics Volume 9 Supplement 1, 2008: The 2007 International Conference on Bioinformatics & Computational Biology (BIOCOMP'07). The full contents of the supplement are available online at http://www.biomedcentral.com/14712164/9?issue=S1.
References

Cheng Y, Church GM: Biclustering of Expression Data.
Proc Int Conf Intell Syst Mol Biol 2000, 93103. PubMed Abstract

Segal E, Battle A, Koller D: Decomposing gene expression into cellular processes. In In Proceedings of the Pacific Symposium on Biocomputing; USA. Edited by Klein E. World Scientific; 2003:89100.

Tseng G, Wong W: Tight clustering: a resamplingbased approach for identifying stable and tight patterns in Data.
Biometrics 2005, 61:1016. PubMed Abstract  Publisher Full Text

Bergmann S, Ihmels J, Barkai N: Iterative signature algorithm for the analysis of largescale gene expression data.

Tanay A, Sharan R, Shamir R: Discovering statistically significant biclusters in gene expression data.
Bioinformatics 2002, 18:S136S144. PubMed Abstract  Publisher Full Text

BenDor A, Chor B, Karp R, Yakhini Z: Discovering local structure in gene expression data: the orderpreserving sumatrix problem. In Proceedings of the sixth international coference on computational biology. Washington, DC, USA, ACM; 2002:89100.

Murali TM, Kasif S: Exacting convserved gene expression motifs from gene exression data. In In Proceedings of the Pacific Symposium on Biocomputing; USA. Edited by Edited by Klein E. World Scientific; 2003:7788.

Madeira S, Oliveira A: Biclustering algorithms for biological data analysis: a survey.
IEEE Trans. on Comp. Biol. and Bioinfo 2004, 1:2445. Publisher Full Text

Barkow S, Bleuler S, Prelic A, Zimmermann E, Pand Zitzler: BicAT: biclustering analysis toolbox.
Bioinformatics 2006, 22:12821283. PubMed Abstract  Publisher Full Text

Sharan R, MaronKatz A, Shamir R: CLICK and EXPANDER: a system for clustering and visualizing gene expression data.
Bioinformatics 2003, 19:17871799. PubMed Abstract  Publisher Full Text

Mendes P, Sha W, Ye K: Artifical gene networks for objective comparison of analysis algorithms.
Bioinformatics 2003, 19:ii122ii129. PubMed Abstract  Publisher Full Text

Golub , et al.: Comprenhensive identification of cell cycleregulated genes of the yeast Saccharomyces cerevisiase by microarray hybridization.
Mol. Biol. Cell 1998, 9:32733297. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Bolstad BM, Irizarry RA, Astrand M, Speed TP: A Comparison of Normalization Methods for High Density Oligonucleotide Array Data Based on Bias and Variance.
Bioinformatics 2003, 19(2):185193. PubMed Abstract  Publisher Full Text

Beer M, Tavazoie S: Predicting gene expression from sequence.
Cell 2004, 117:185198. PubMed Abstract  Publisher Full Text

Gasch A, et al.: Genomic expression programs in the response of yeast cells to environmental changes.
Mol. Biol. Cell 2000, 11:42414257. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Spellman PT, et al.: Comprenhensive identification of cell cycleregulated genes of the yeast Saccharomyces cerevisiase by microarray hybridization.
Mol. Biol. Cell 1998, 9:32733297. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Schwarz G: Estimating the dimension of a model.
Annals of statistics 1978, 6:461464. Publisher Full Text

Hughes J, Estep P, Tavazoie S, Church G: Compuational identifcation of cisregulatory elements associated with groups of functionally related genes in Saccharomyces cerevisiae.
J. Mol. Biol 2000, 296:12051214. PubMed Abstract  Publisher Full Text

Klein C, Struhl K: Protein kinase A mediates growthregulated expression of yeast ribosomal protein genes by modulating RAP1 transcriptional activity.
Mol. Cell. Biol 1994, 14:19208. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Mannhaupt G, et al.: Rpn4p acts as a transcription factor by binding to PACE, a nonamer box found upstream of 26S proteasomal and other genes in yeast.
FEBS Lett 1999, 2734. PubMed Abstract  Publisher Full Text

Thomas D, Jacquemin I, SurdinKerjan Y: MET4, a leucine zipper protein, and centromerebinding factor 1 are both required for transcriptional activation of sulfur metabolism in Saccharomyces cerevisiae.
Mol. Cell. Biol 1992, 12:171927. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Thomas D, SurdinKerjan Y: Metabolism of sulfur amino acids in Saccharomyces cerevisiae.
Microbiol. Mol. Biol. Rev 1997, 61:503532. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Nasmyth K, Dirick L: The role of SWI4 and SWI6 in the activity of G1 cyclins in yeast.
Cell 1991, 66:9951013. PubMed Abstract  Publisher Full Text

Toone W, Jones N: AP1 transcription factors in yeast.
Curr Opin Genet Dev 1999, 9:5561. PubMed Abstract  Publisher Full Text

Clifton D, Weinstock S, Fraenkel D: Glycolysis mutants in Saccharomyces cerevisiae.
Genetics 1978, 88:111. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Liu JS: Monte Carlo Strategies in Scientific Computing. In Genetics. New York: SpringerVerlag; 2001.