Determination of strongly overlapping signaling activity from microarray data
© Bidaut et al; licensee BioMed Central Ltd. 2006
Received: 22 September 2005
Accepted: 28 February 2006
Published: 28 February 2006
As numerous diseases involve errors in signal transduction, modern therapeutics often target proteins involved in cellular signaling. Interpretation of the activity of signaling pathways during disease development or therapeutic intervention would assist in drug development, design of therapy, and target identification. Microarrays provide a global measure of cellular response, however linking these responses to signaling pathways requires an analytic approach tuned to the underlying biology. An ongoing issue in pattern recognition in microarrays has been how to determine the number of patterns (or clusters) to use for data interpretation, and this is a critical issue as measures of statistical significance in gene ontology or pathways rely on proper separation of genes into groups.
Here we introduce a method relying on gene annotation coupled to decompositional analysis of global gene expression data that allows us to estimate specific activity on strongly coupled signaling pathways and, in some cases, activity of specific signaling proteins. We demonstrate the technique using the Rosetta yeast deletion mutant data set, decompositional analysis by Bayesian Decomposition, and annotation analysis using ClutrFree. We determined from measurements of gene persistence in patterns across multiple potential dimensionalities that 15 basis vectors provides the correct dimensionality for interpreting the data. Using gene ontology and data on gene regulation in the Saccharomyces Genome Database, we identified the transcriptional signatures of several cellular processes in yeast, including cell wall creation, ribosomal disruption, chemical blocking of protein synthesis, and, criticially, individual signatures of the strongly coupled mating and filamentation pathways.
This works demonstrates that microarray data can provide downstream indicators of pathway activity either through use of gene ontology or transcription factor databases. This can be used to investigate the specificity and success of targeted therapeutics as well as to elucidate signaling activity in normal and disease processes.
Many diseases develop because of errors in signaling, and newer therapeutics specifically target proteins involved in cellular signaling [1, 2]. However, these therapies are not always effective , and the reason for failure, whether inherent poor interaction or complex cellular response, is unknown. In order to understand the development of disease and drug resistance in these cases, the recovery of the process that led to the specific cellular malfunction must be identified. Such errors generally involve the cellular signaling networks that control cell growth, differentiation, apoptosis, and motility [4, 5]. Because of the extreme underlying biological complexity of these pathways, diseases that involve errors in signaling processes arise from a myriad of different cellular malfunctions, for example in cancers [6, 7] and diabetes [8, 9]. It is from this complex background that functional genomics attempts to glean insight to improve our understanding of diseases.
One of the major uses of microarrays has been elucidation of gene expression in cancer, often focused on refining cancer identification using computational and statistical approaches [10–12]. In addition, the discovery of biomarkers in the form of differential levels of production of mRNA has been a focus in a number of studies [13–15]. The fact that determination of the mRNA levels of a single gene is easier than using an entire array has driven the shift to the use of arrays to generate potential biomarkers, so that the expression levels of these individual genes can be screened for in a more economical way (see, for example, ). For diabetes, microarrays have been used to elucidate gene expression in both type I and type II diseases, and customized chips targeting genes of interest have been developed .
Many tools for statistical inference, pattern recognition, and data mining have been developed for microarray data analysis. Statistical tests include SAM , VERAandSAM , ANOVA techniques [20, 21], Bayesian approaches [22, 23], and rank tests . Pattern recognition and data mining techniques comprise both unsupervised techniques, such as hierarchical clustering , singular value decomposition , multidimensional scaling , Bayesian mixture models [28, 29], and other clustering methods [30–34], and supervised techniques, such as support vector machines  and artificial neural networks , (for a review see ).
While these techniques are useful, they have certain limitations as regards more advanced uses in the elucidation of mechanisms operating in diseased tissues. New therapeutics specifically target proteins involved in cellular signaling [1–3, 38–40]. As noted above, these therapies are not always effective, and a method to understand the reason for their ineffectiveness is highly desirable. If the failure modes for the targeted therapeutics are understood, new therapeutics can be designed or combination therapies undertaken. In addition, to design new therapies that work alone or in combination with other therapies, an understanding of signaling networks is required. Microarray measurements can provide insight into these issues.
Unfortunately, the recovery of pathway information from transcriptional data requires complex analysis, since signaling protein activity is not generally linked to the mRNA expression levels of genes encoding the signaling proteins themselves , nor are protein levels tightly coupled to transcript levels even in yeast [42, 43]. This makes it impossible to directly link an increase in mRNA expression of the gene encoding a signaling protein, such as the therapeutic target, with activity of the protein and therefore of the signaling pathway. Instead, an analysis must treat changes in mRNA levels as downstream indicators of activity.
An important issue to resolve in order to correctly interpret patterns in microarray data is the underlying dimensionality of the data, since statistical analysis of genes in groups relies on correct separation. The dimensionality provides an estimate of the number of patterns required to explain the variation in the data not related to noise, which is equivalent to the number of basis vectors required mathematically to describe the data or the number of principal components required to span the data.
Dimensionality of the data
The freedom to move between branches leads also to some loss of consistency in the annotations as one moves down a branch. This contrasts with the behavior of basis vectors obtained for less than 15 patterns where biological functions split logically as the number of patterns increases. We also observed a reduction in the number of genes related to each basis vector for 16 patterns in comparison to 15. Also, the standard deviation across samples of the obtained vectors is significantly higher for 16 patterns (1.7 × 10-3) than for 15 patterns (7.9 × l0-4) indicating that the Markov chain sampling is not as tightly constrained by the probability distribution. The behavior observed occurs because of the potential to overfit the data with 16 basis vectors allowing the algorithm to find multiple configurations to explain the variation in the data.
Identifying patterns and functions
Bayesian Decomposition retrieves the two linked matrices: the P matrix (pattern matrix) groups mutants that share cellular functions, which can be deduced from the genes linked to each pattern contained in the A matrix. Each mutant (a column of the P matrix) can belong to multiple patterns, which models the fact that each mutant will have many cellular functions active. Each gene (a row of the A matrix) can be assigned to multiple patterns, reflecting the fact that evolution has led to genes being involved in multiple cellular processes. Interpretation of the results involves identifying cellular processes from the genes that are significantly expressed in a pattern (i.e., within a column of A).
We proceed by using the dimensionality estimate of 15 patterns and exploring for each pattern the genes associated with that pattern. These genes are interpreted using the MIPS ontology for yeast  in order to predict the cellular processes associated with a pattern. In addition, for patterns that can be linked to signaling pathways, we discuss the use of data on genes regulated by specific transcription factors and validate the results by analysis of specific key deletion mutants. For each pattern that shows enhancement of ontological terms we provide the terms, the enhancement (as defined in the methods section), and the p value for a hypergeometric test on the term.
We summarize the results in terms of patterns previously identified in other studies using this data set, then we present the new results isolating signatures for activity of the mating and filametation pathways.
Patterns identified in previous studies
Examination of pattern 1 shows expression of the overall common minimal processes necessary for survival, with 386 annotated genes associated with this pattern at a 3σ level. Measure of enhancement, e, of cellular functions, reveals two highly represented functional groups: 1) groups related to protein synthesis and 2) groups related to DNA synthesis. Group 1 includes genes enhanced in Protein Targeting, Sorting and Translocation (Term 14.04, e = 1.84, p = 0.0022), Protein Synthesis (Term 12, e = 1.55, p = 0.024), and Ribosome Biogenesis (Term 12.01, e = 1.60, p = 0.059). Group 2 includes DNA Processing (Term 10.01, e = 1.32, p = 0.097), DNA Recombination and DNA Repair (Term 10.01.05, e = 1.31, p = 0.16), and DNA Synthesis (Term 10.01.03, e = 1.58, p = 0.16). The p-values for the ontology terms remain high, due to the large number of genes associated with this pattern.
This pattern, which essentially includes genes necessary for viability, contains all the mutants of the dataset, although the Ssn6Δ mutant shows a lower level for this pattern than other mutants. As the Ssn6Δ mutant exhibits substantially greater overall expression than any other mutant (including the Tup1Δ mutant with the second highest level), this may reflect the high association of the Ssn6Δ mutant seen in almost all patterns, which will have some gene overlap with this pattern.
Pattern 5 contains 172 annotated genes. Highly enhanced ontologies include transport-related functions: Transported Compounds (Term 20.01, e = 2.02, p < 10-4), C-compound and Carbohydrate Transport (Term 20.01.03, e = 2.60, p < 3 × 10-4), and Cellular Transport, Transport Facilitation and Cellular Routes (Term 20, e = 1.62, p < 6 × 10-4), in addition to other transport terms at p < 10-3. The pattern contains the two deletion mutants, Ssn6Δ and Tup1Δ, and represents the strong response seen in the original study . Ssn6p and Tup1p form a system of transcriptional repression that appears to be highly conserved in eukaryotes . In yeast, the complex acts as a global transcriptional repressor over a large number of genes (more that 150), coordinating several cellular systems, including haploid specific genes, glucose repressible genes, and oxygen utilization genes . Turning off this repression leads to a large overall increase in gene expression (the overall expression in these two mutants is many fold higher than in other mutants).
Pattern 7 is related to the lack of cell wall functions (Cell Wall, Term 42.01, e = 0.0), as 28 cell wall genes (of 32 total) are absent from this pattern, while the other four genes have multiple annotations suggesting they have roles unrelated to Cell Wall function. Enhancement is present for Protein Modification (Term 14.07, e = 4.7, p < 10-4) and Fermentation (Term 02.16, e = 4.5, p < 0.01). This pattern contains the mutants Gas1Δ and Fks1Δ, which impair cell-wall synthesis, as well as the mutant YER083cΔ, annotated as disrupting the cell wall in the original study . The pattern contains other mutants disrupting ergosterol biosynthesis as well, including Erg2Δ, She4Δ, as well as YER044cΔ. In addition, the pattern includes yeast treated with the drugs that are known to disrupt the cell wall, such as Tunicamycin and Glucosamine.
Pattern 11 is related to ribosomal function, with enhancements in terms for Ribosome Biogenesis (Term 12.01, e = 6.32, p < 4 × 10-4) and Protein Synthesis (Term 12, e = 3,89, p < 6 × 10-3). The pattern contains 8 mutants related to ribosomal proteins, Rpll2aΔ, Rpl27aΔ, Rpl34aΔ, Rpl6bΔ, Rp18aΔ, Rps24aΔ, Rps24aΔ (haploid), and Rps27bΔ, as well as some mutants with deleted ORFs of unknown function, YOR078wΔ, YMR269wΔ, and YHR034cΔ, proposed to be involved in ribosomal functions .
Patterns related to cellular signalling pathways
The most enhanced gene ontology terms in patterns 13 and 15. Each term is presented together with a measure of how overrepresented it is compared to a random draw of the same number of genes. These were also confirmed to be significant by hypergeometric tests.
Pheromone response, mating-type determination
Transposable elements, viral and plasmid
Pheromone response, mating-type determination
Transmembrane signal transduction
G-protein mediated signal transduction
Transmembrane signal transduction
Chemoperception and response
Cellular sensing and response
Interaction with cellular environment
Chemoperception and response
Cellular sensing and response
Enzyme mediated signal transduction
Development of asco- basidio- or zygospore
Interaction with cellular environment
Enzyme mediated signal transduction
Protein kinase cascades
G protein mediated signal transduction
Budding, cell polarity, filamentation
The genes most strongly associated with patterns 13 and 15 in order of strength of association. For pattern 13, it is noted whether the genes are known to be regulated in the mating process, and whether the gene is known to be directly regulated by Stel2p. For pattern 15, the gene function is shown. All data is from the Saccharomyces Genome Database [73, 74].
Transposable element gene
Transposable element gene
Verified ORF, Prm5
Transposable element gene
Transposable element gene
Removed from SGD
Transposable element gene
Transposable element gene
Transposable element gene
In order to validate that the patterns were actually measuring activity of the mating and filamentation pathways, we explored deletion mutants related to these pathways [61, 62]. The mating response in S. cerevisiae is mediated via a MAPK signaling cascade initiated by binding to the Ste2p or Ste3p membrane receptors (Figure 3). The signal is transduced through Ste11p, Ste7p, and Fus3p with Ste5p serving as a scaffolding protein. The signal activates the Ste12p transcription factor, leading to transcription of mating response genes. In addition, the signal is transduced to the MAPK cascade from the membrane receptor by a G protein complex or through the Ste20p protein. Pattern 13 shows near zero signal for the deletion mutants Ste11Δ, Ste7Δ Fus3Δ, Ste12Δ, Ste5Δ, and Ste2Δ, while showing signal for deletion mutants of Ste20Δ and Tec1Δ. This is exactly as expected, with the membrane receptor, all signaling proteins in the cascade, and the transcription factor necessary to generate the transcriptional response related to the mating signal (note that the Ste3Δ mutant is not in the data set). Ste20p is not necessary to raise the mating response, since the G-protein complex can trigger activation of Ste11p directly. For pattern 15, the response is very similar. The signal is near zero for the deletion mutants Ste2Δ, Ste11Δ, Ste7Δ, and Ste12Δ. The Fus3Δ mutant shows a signal for pattern 15, as appropriate, while the Fus3Δ, Kss1Δ double deletion mutant does not. In addition, the Tec1Δ mutant shows no signal for pattern 15, indicating that Tec1p is required for filamentation . Finally, the signal is greatly reduced for the Ste20Δ deletion mutant in pattern 15 relative to pattern 13, which agrees with previous work suggesting that the filamentation pathway is more dependent on Ste20p signaling than is the mating pathway .
Microarrays and GeneChips™ have become the tools of choice for the investigation of genome-wide transcription in most biological systems. The resulting data comprises noisy estimates of transcription levels for roughly 6,000 genes in yeast to more than 20,000 genes in typical mammalian studies. Numerous statistical and data mining methods have been applied to this data in order to identify individual genes showing differential expression, to identify patterns related to physiological states, and to identify groups of genes comprising biological processes. These studies generally have not focused on the estimation of cellular signaling from the data, despite the prevalence of cellular signaling in many diseases.
As noted in the introduction, the recovery of signaling pathway information from transcription data requires complex analysis, since protein levels do not correlate well with mRNA levels and signaling protein activity is not generally linked to the mRNA expression levels of genes encoding the signaling proteins themselves. As such, changes in mRNA levels are limited to being a downstream indicator of activity. If a complete model for the transcription of genes, including all known transcriptional regulators and biological processes regulating transcription, was available, the inference of activity would be straightforward. Unfortunately, the network models and even gene annotations are still far from being complete. In addition, the growing evidence supporting the important role of non-coding RNAs in regulation of gene expression (including antisense transcripts and micro-RNAs, see for example [63–65]) further undermines the potential of using mRNA species as markers for proteins and their activities .
In order to overcome this incompleteness, we have created the method described here. We couple identification of transcriptional signatures with our Bayesian Decomposition algorithm to a consistency analysis for gene assignment to patterns determined by comparison of different dimensionalities using ClutrFree. This allows the identification of the correct dimensionality to be applied to subsequent ontology and transcription factor analyses. ClutrFree is also used to determine the ontological terms enhanced within each pattern and to obtain a list of genes tied to this pattern, which can then be linked to specific transcription factors. In this way, the biological processes associated with conditions can be identified, and inferences can be made on the activity of specific transcription factors. This then allows inference on the activity of signaling pathways, which cannot be obtained with methods previously applied to microarray data. Overall, the method requires many separate steps, each modeling an aspect of the biological system, in order to make proper inferences on signaling from the data.
In the application to the Compendium data presented here, our analysis was able to extract the common features for a set of mutants that eliminated related pathways. As in previous studies, the global transcriptional repressor complex Ssn6-Tupl has been isolated in a single group. In addition, patterns for cell-wall synthesis, ribosomal function, and the global functions necessary for continued viability of yeast were isolated. In contrast to previous analyses of this data, two pathways related to the MAPK cascade were isolated, one related to mating and the other to filamentation. Once the correct dimensionality was determined, Bayesian Decomposition was able to identify transcriptional signatures unique for each pathway. The assignment was validated by an investigation of the deletion mutants known to adversely affect these pathways.
Microarray studies have been widespread in biological and medical research, often focusing on identification of genes significantly correlated with various disease states. However, many diseases arise from disruptions in cellular signaling, and in these cases gene expression only provides a downstream indicator of signaling activity. This greatly complicates the analysis. The new approach introduced here recovered signatures allowing us to make validated inferences on strongly overlapping signaling pathways.
The results demonstrate that for Saccharomyces cerevisiae, the mating and filamentation pathways can be distinguished from transcriptional signatures determined from analysis of microarray data, despite the intrinsic high noise, confounding transcriptional activity, and tightly coupled nature of the pathways. The next step will be to apply these methods to more complex signaling networks in worms, flies, and mammals.
The Rosetta deletion mutant data set
The Rosetta Compendium comprises 300 conditions, including 276 deletion mutants, 11 tetracycline regulated genes, and 13 drug treatments, in S. cerevisiae growing in rich medium . The data were generated from a two color cDNA microarray hybridization assay , and transcriptional profiles were measured both with technical replication and biological replication (151 mutants). In parallel with the 300 experiments, 63 controls of wild-type S. cerevisiae were grown in identical conditions and compared against each other, permitting creation of a gene-specific error model. The data was downloaded from Rosetta Inpharmatics.
The data was filtered to retain only conditions characterized by at least a variation of 3 fold in a minimum of 2 genes. Then all genes that did not vary by 3-fold in at least 2 conditions were also removed, leading to a data matrix comprising 764 genes and 228 conditions.
The data used by Bayesian Decomposition included both the mean log ratio for each data point and the uncertainty in this measurement determined by the Rosetta error model. Since Bayesian Decomposition as applied here requires positivity , the log ratios were converted to ratios and the uncertainties propagated to uncertainties on the ratios. Although analysis of residuals suggested that seven dimensions fit the data , the analysis presented here suggests that this is due to overestimation of uncertainty in the data. Bayesian Decomposition is not highly sensitive to minor misestimations of noise however, so that this should not be a problem for this analysis.
Analysis with Bayesian Decomposition
Bayesian Decomposition (BD) has been applied to multiple types of data: in vivo spectroscopic data , medical imaging , microarray data from single cell organisms [44, 45], mammalian model organisms , humans , and on phylogenomic sequence data . A detailed description of the algorithm  and a review of applications  have been published.
Briefly, BD models the microarray data, comprising a matrix of estimates of the ratio of expression between the experimental condition and a control, as the result of the multiplication of two matrices describing behaviors across conditions (the P or pattern matrix) and the distribution of genes within those behaviors (the A or amplitude matrix). Naturally, the data, D, includes noise, so that the full relationship is defined by
where D ij is the estimated ratio for gene i in mutant j, A ik is the strength of gene i in pattern k, P kj is the strength of mutant j in pattern k, and ε ij is the noise for gene i in mutant j estimated by the Rosetta error model. BD estimates A and P by a Markov chain Monte Carlo (MCMC) simulation. For the fixed noise estimate, ε, A and P are inferred from the marginal probability distribution
p (D | A,P) p (A,P) 
where p(A,P) is the prior probability and p (D | A,P) is given by the likelihood such that
The prior is used here to require positivity and to minimize structure in the estimates of A and P .
The analysis with BD is similar mathematically to an analysis with singular value decomposition (SVD) or with principal component analysis (PCA), since all methods estimate two matrices that together reconstruct the data. In both PCA and SVD, orthogonality conditions force each row of P to be linearly independent, deriving P from either the data matrix (SVD) or the covariance matrix (PCA) using deterministic algorithms. Since patterns of expression related to biological processes will not generally be independent, BD uses the MCMC approach to avoid orthogonality. The resulting rows of P are usually easier to relate to biological processes than those from SVD or PCA.
BD was run at each posited dimensionality, K (as in equation 1), between 3 and 25, generating a mean estimate and an uncertainty (i.e., standard deviation from samples of the posterior probability distribution) for each element of A and P. The dimensionality is equivalent to the number of patterns, as the patterns can be viewed as nonorthogonal basis vectors for D.
Persistence and dimensionality
We assessed the dimensionality of the data using measurements of persistence. The persistence was measured for analyses from 3 – 25 patterns and the dimension chosen where a significant drop occurred in an otherwise slow monotonic decline, which was expected due to the branching nature of the tree. Figure 2 shows the significant drop between 15 and 16 dimensions, so 15 patterns were chosen for further analysis.
Ontology and function
To assign ontological terms to the genes contained in basis vectors, we annotated our data using the gene ontologies from the Comprehensive Yeast Genome Database (CYGD) hosted at the Munich Information center for Protein Sequences (MIPS) [51, 72]. The analysis here utilizes the Functional Catalog (FunCat) format that describes each gene using a hierarchical ontological model.
Similar to persistence, we defined enhancement as a measure of the over-representation, or under representation, of a gene function in a subset of the data . It is the ratio of the frequency of occurrence of genes annotated by a particular ontological term in the pattern to the frequency of occurrence of the same term in the whole dataset,
with g p being the number of genes annotated with the term t in pattern p, n p the total number of genes in the pattern p, G the number of genes annotated by the term t in the data, and N being the total number of genes in the dataset. In addition, we apply a hypergeometric test to estimate a p-value for each term. Function was then determined by inspection of enhanced ontological terms.
Transcription factor analysis
The genes were also analyzed for the patterns determined to be related to signaling pathways by exploration of the ten genes most strongly associated with the pattern. Each gene was analyzed using the Saccharomyces Genome Database  to determine whether it was known to be associated with mating or filamentation processes and to determine if it was directly regulated by the Ste12p transcription factor.
We thank the National Institutes of Health, National Library of Medicine (LM008309 to mfo) and National Cancer Institute (CCCG CA06927 supporting mfo), and the Pennsylvania Department of Health (grant to mfo). JMC and KS acknowledge the support from the Réseau National des Génopoles (RNG).
- Mauro MJ, Druker BJ: STI571: targeting BCR-ABL as therapy for CML. Oncologist 2001, 6: 233–8.View ArticlePubMed
- Repka T, Chiorean EG, Gay J, Herwig KE, Kohl VK, Yee D, Miller JS: Trastuzumab and interleukin-2 in HER2-positive metastatic breast cancer: a pilot study. Clin Cancer Res 2003, 9: 2440–6.PubMed
- von Mehren M: Recent advances in the management of gastrointestinal stromal tumors. Curr Oncol Rep 2003, 5: 288–94.View ArticlePubMed
- Jacks T, Weinberg RA: Taking the study of cancer cell survival to a new dimension. Cell 2002, 111: 923–5.View ArticlePubMed
- Kolch W: Meaningful relationships: the regulation of the Ras/Raf/MEK/ERK pathway by protein interactions. Biochem J 2000, 351(Pt 2):289–305.PubMed CentralView ArticlePubMed
- Cooper GM: Elements of Human Cancer. Boston: Jones and Bartlett Publishers; 1992.
- Macdonald F, Ford CHJ: Molecular Biology of Cancer. Oxford: BIOS Scientific Publishers, Ltd; 1997.
- Zdychova J, Komers R: Emerging role of Akt kinase/protein kinase B signaling in pathophysiology of diabetes and its complications. Physiol Res 2005, 54: 1–16.PubMed
- Leng Y, Karlsson HK, Zierath JR: Insulin signaling defects in type 2 diabetes. Rev Endocr Metab Disord 2004, 5: 111–7.View ArticlePubMed
- Alizadeh AA, Eisen MB, Davis RE, Ma C, Lossos IS, Rosenwald A, Boldrick JC, Sabet H, Tran T, Yu X, et al.: Distinct types of diffuse large B-cell lymphoma identified by gene expression profiling. Nature 2000, 403: 503–11.View ArticlePubMed
- Golub TR, Slonim DK, Tamayo P, Huard C, Gaasenbeek M, Mesirov JP, Coller H, Loh ML, Downing JR, Caligiuri MA, et al.: Molecular classification of cancer: class discovery and class prediction by gene expression monitoring. Science 1999, 286: 531–7.View ArticlePubMed
- Zhang H, Yu CY, Singer B, Xiong M: Recursive partitioning for tumor classification with gene expression microarray data. Proc Natl Acad Sci USA 2001, 98: 6730–5.PubMed CentralView ArticlePubMed
- Williams NS, Gaynor RB, Scoggin S, Verma U, Gokaslan T, Simmang C, Fleming J, Tavana D, Frenkel E, Becerra C: Identification and validation of genes involved in the pathogenesis of colorectal cancer using cDNA microarrays and RNA interference. Clin Cancer Res 2003, 9: 931–46.PubMed
- Kikuchi T, Daigo Y, Katagiri T, Tsunoda T, Okada K, Kakiuchi S, Zembutsu H, Furukawa Y, Kawamura M, Kobayashi K, et al.: Expression profiles of non- small cell lung cancers on cDNA microarrays: identification of genes for prediction of lymph-node metastasis and sensitivity to anti-cancer drugs. Oncogene 2003, 22: 2192–205.View ArticlePubMed
- Carr KM, Bittner M, Trent JM: Gene-expression profiling in human cutaneous melanoma. Oncogene 2003, 22: 3076–80.View ArticlePubMed
- Frolov A, Chahwan S, Ochs M, Arnoletti JP, Pan ZZ, Favorova O, Fletcher J, von Mehren M, Eisenberg B, Godwin AK: Response markers and the molecular mechanisms of action of Gleevec in gastrointestinal stromal tumors. Mol Cancer Ther 2003, 2: 699–709.PubMed
- Scearce LM, Brestelli JE, McWeeney SK, Lee CS, Mazzarelli J, Pinney DF, Pizarro A, Stoeckert CJ Jr, Clifton SW, Permutt MA, et al.: Functional genomics of the endocrine pancreas: the pancreas clone set and PancChip, new resources for diabetes research. Diabetes 2002, 51: 1997–2004.View ArticlePubMed
- Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A 2001, 98: 5116–21.PubMed CentralView ArticlePubMed
- Ideker T, Thorsson V, Siegel AF, Hood LE: Testing for differentially-expressed genes by maximum-likelihood analysis of microarray data. J Comput Biol 2000, 7: 805–17.View ArticlePubMed
- Kerr MK, Martin M, Churchill GA: Analysis of variance for gene expression microarray data. J Comput Biol 2000, 7: 819–37.View ArticlePubMed
- Kerr MK, Afshari CA, Bennett L, Bushel P, Martinez J, Walker NJ, Churchill GA: Statistical analysis of a gene expression microarray experiment with replication. Statistica Sinica 2002, 12: 203–218.
- Newton MA, Kendziorski CM, Richmond CS, Blattner FR, Tsui KW: On differential variability of expression ratios: improving statistical inference about gene expression changes from microarray data. J Comput Biol 2001, 8: 37–52.View ArticlePubMed
- Parmigiani G, Garrett E, Anbazhagan R, Gabrielson E: A statistical framework for expression-based molecular classification in cancer. Journal of the Royal Statistical Society, B 2002, 64: 717–736.View Article
- Troyanskaya OG, Garber ME, Brown PO, Botstein D, Altman RB: Nonparametric methods for identifying differentially expressed genes in microarray data. Bioinformatics 2002, 18: 1454–61.View ArticlePubMed
- Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A 1998, 95: 14863–8.PubMed CentralView ArticlePubMed
- Alter O, Brown PO, Botstein D: Singular value decomposition for genome-wide expression data processing and modeling. Proc Natl Acad Sci U S A 2000, 97: 10101–6.PubMed CentralView ArticlePubMed
- Khan J, Simon R, Bittner M, Chen Y, Leighton SB, Pohida T, Smith PD, Jiang Y, Gooden GC, Trent JM, et al.: Gene expression profiling of alveolar rhabdomyosarcoma with cDNA microarrays. Cancer Res 1998, 58: 5009–13.PubMed
- Medvedovic M, Sivaganesan S: Bayesian infinite mixture model based clustering of gene expression profiles. Bioinformatics 2002, 18: 1194–206.View ArticlePubMed
- Medvedovic M, Yeung KY, Bumgarner RE: Bayesian mixture model based clustering of replicated microarray data. Bioinformatics 2004.
- Gasch AP, Eisen MB: Exploring the conditional coregulation of yeast gene expression through fuzzy k-means clustering. Genome Biol 2002, 3: RESEARCH0059.PubMed CentralView ArticlePubMed
- Getz G, Levine E, Domany E: Coupled two-way clustering analysis of gene microarray data. Proc Natl Acad Sci U S A 2000, 97: 12079–84.PubMed CentralView ArticlePubMed
- Ben-Dor A, Shamir R, Yakhini Z: Clustering gene expression patterns. J Comput Biol 1999, 6: 281–97.View ArticlePubMed
- Heyer LJ, Kruglyak S, Yooseph S: Exploring expression data: identification and analysis of coexpressed genes. Genome Res 1999, 9: 1106–15.PubMed CentralView ArticlePubMed
- Lukashin AV, Fuchs R: Analysis of temporal gene expression profiles: clustering by simulated annealing and determining the optimal number of clusters. Bioinformatics 2001, 17: 405–14.View ArticlePubMed
- Brown MP, Grundy WN, Lin D, Cristianini N, Sugnet CW, Furey TS, Ares M Jr, Haussler D: Knowledge-based analysis of microarray gene expression data by using support vector machines. Proc Natl Acad Sci U S A 2000, 97: 262–7.PubMed CentralView ArticlePubMed
- Khan J, Wei JS, Ringner M, Saal LH, Ladanyi M, Westermann F, Berthold F, Schwab M, Antonescu CR, Peterson C, et al.: Classification and diagnostic prediction of cancers using gene expression profiling and artificial neural networks. Nat Med 2001, 7: 673–9.PubMed CentralView ArticlePubMed
- Ochs MF, Godwin AK: Microarrays in cancer: research and applications. Biotechniques 2003, 34: S4-S15.
- Kato-Stankiewicz J, Hakimi I, Zhi G, Zhang J, Serebriiskii I, Guo L, Edamatsu H, Koide H, Menon S, Eckl R, et al.: Inhibitors of Ras/Raf-1 interaction identified by two-hybrid screening revert Ras-dependent transformation phenotypes in human cancer cells. Proc Natl Acad Sci U S A 2002, 99: 14398–403.PubMed CentralView ArticlePubMed
- Strumberg D, Seeber S: Raf kinase inhibitors in oncology. Onkologie 2005, 28: 101–7.View ArticlePubMed
- Heim M, Scharifi M, Zisowsky J, Jaehde U, Voliotis D, Seeber S, Strumberg D: The Raf kinase inhibitor BAY 43–9006 reduces cellular uptake of platinum compounds and cytotoxicity in human colorectal carcinoma cell lines. Anticancer Drugs 2005, 16: 129–36.View ArticlePubMed
- Chen G, Gharib TG, Huang CC, Taylor JM, Misek DE, Kardia SL, Giordano TJ, lannettoni MD, Orringer MB, Hanash SM, et al.: Discordant protein and mRNA expression in lung adenocarcinomas. Mol Cell Proteomics 2002, 1: 304–13.View ArticlePubMed
- Gygi SP, Rochon Y, Franza BR, Aebersold R: Correlation between protein and mRNA abundance in yeast. Mol Cell Biol 1999, 19: 1720–30.PubMed CentralView ArticlePubMed
- Griffin TJ, Gygi SP, Ideker T, Rist B, Eng J, Hood L, Aebersold R: Complementary profiling of gene expression at the transcriptome and proteome levels in Saccharomyces cerevisiae. Mol Cell Proteomics 2002, 1: 323–33.View ArticlePubMed
- Bidaut G, Moloshok TD, Grant JD, Manion FJ, Ochs MF: Bayesian Decomposition analysis of gene expression in yeast deletion mutants. In Methods of Microarray Data Analysis II. Edited by: Johnson K, Lin S. Boston: Kluwer Academic; 2002:105–122.View Article
- Moloshok TD, Klevecz RR, Grant JD, Manion FJ, Speier WFt, Ochs MF: Application of Bayesian Decomposition for analysing microarray data. Bioinformatics 2002, 18: 566–75.View ArticlePubMed
- Ochs MF: Bayesian Decomposition. In The Analysis of Gene Expression Data: Methods and Software. Edited by: Parmigiani G, Garrett E, Irizarry R, Zeger S. New York: Springer Verlag; 2003.
- Moloshok TD, Datta D, Kossenkov AV, Ochs MF: Bayesian Decomposition classification of the Project Normal data set. In Methods of Microarray Data Analysis III. Edited by: Johnson KF, LIn SM. Boston: Kluwer Academic; 2003:211–232.
- Ochs MF, Moloshok TD, Bidaut G, Toby G: Bayesian Decomposition: Analyzing microarray data within a biological context. Annals of the New York Academy of Sciences 2004, 1020: 212–226.View ArticlePubMed
- Bidaut G, Ochs MF: ClutrFree: cluster tree visualization and interpretation. Bioinformatics 2004, 20: 2869–71.View ArticlePubMed
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al.: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 2000, 25: 25–9.PubMed CentralView ArticlePubMed
- Mewes HW, Amid C, Arnold R, Frishman D, Guldener U, Mannhaupt G, Munsterkotter M, Pagel P, Strack N, Stumpflen V, et al.: MIPS: analysis and annotation of proteins from whole genomes. Nucleic Acids Res 2004, 32: D41–4.PubMed CentralView ArticlePubMed
- Matys V, Fricke E, Geffers R, Gossling E, Haubrock M, Hehl R, Hornischer K, Karas D, Kel AE, Kel-Margoulis OV, et al.: TRANSFAC: transcriptional regulation, from patterns to profiles. Nucleic Acids Res 2003, 31: 374–8.PubMed CentralView ArticlePubMed
- Hughes TR, Marton MJ, Jones AR, Roberts CJ, Stoughton R, Armour CD, Bennett HA, Coffey E, Dai H, He YD, et al.: Functional discovery via a compendium of expression profiles. Cell 2000, 102: 109–26.View ArticlePubMed
- Kim PM, Tidor B: Subsystem identification through dimensionality reduction of large-scale gene expression data. Genome Res 2003, 13: 1706–18.PubMed CentralView ArticlePubMed
- Smith RL, Johnson AD: Turning genes off by Ssn6-Tupl: a conserved system of transcriptional repression in eukaryotes. Trends Biochem Sci 2000, 25: 325–30.View ArticlePubMed
- Keleher CA, Redd MJ, Schultz J, Carlson M, Johnson AD: Ssn6-Tupl is a general repressor of transcription in yeast. Cell 1992, 68: 709–19.View ArticlePubMed
- Kusari AB, Molina DM, Sabbagh W Jr, Lau CS, Bardwell L: A conserved protein interaction network involving the yeast MAP kinases Fus3 and Kss1. J Cell Biol 2004, 164: 267–77.PubMed CentralView ArticlePubMed
- Madhani HD, Fink GR: The riddle of MAP kinase signaling specificity. Trends Genet 1998, 14: 151–5.View ArticlePubMed
- Schwartz MA, Madhani HD: Principles of MAP kinase signaling specificity in Saccharomyces cerevisiae. Annu Rev Genet 2004, 38: 725–48.View ArticlePubMed
- Mosch HU, Fink GR: Dissection of filamentous growth by transposon mutagenesis in Saccharomyces cerevisiae. Genetics 1997, 145: 671–84.PubMed CentralPubMed
- Morillon A, Springer M, Lesage P: Activation of the Kss1 invasive-filamentous growth pathway induces Ty1 transcription and retrotransposition in Saccharomyces cerevisiae. Mol Cell Biol 2000, 20: 5766–76.PubMed CentralView ArticlePubMed
- Posas F, Takekawa M, Saito H: Signal transduction by MAP kinase cascades in budding yeast. Curr Opin Microbiol 1998, 1: 175–82.View ArticlePubMed
- Katayama S, Tomaru Y, Kasukawa T, Waki K, Nakanishi M, Nakamura M, Nishida H, Yap CC, Suzuki M, Kawai J, et al.: Antisense transcription in the mammalian transcriptome. Science 2005, 309: 1564–6.View ArticlePubMed
- Claverie JM: Fewer genes, more noncoding RNA. Science 2005, 309: 1529–30.View ArticlePubMed
- Carninci P, Kasukawa T, Katayama S, Gough J, Frith MC, Maeda N, Oyama R, Ravasi T, Lenhard B, Wells C, et al.: The transcriptional landscape of the mammalian genome. Science 2005, 309: 1559–63.View ArticlePubMed
- McManus MT, Sharp PA: Gene silencing in mammals by small interfering RNAs. Nat Rev Genet 2002, 3: 737–47.View ArticlePubMed
- Schena M, Shalon D, Davis RW, Brown PO: Quantitative monitoring of gene expression patterns with a complementary DNA microarray. Science 1995, 270: 467–70.View ArticlePubMed
- Ochs MF, Stoyanova RS, Arias-Mendoza F, Brown TR: A new method for spectral decomposition using a bilinear Bayesian approach. J Magn Reson 1999, 137: 161–76.View ArticlePubMed
- Ochs MF, Stoyanova RS, Brown TR, Rooney WD, Springer CS Jr: A Bayesian Markov chain Monte Carlo solution of the bilinear problem. In Bayesian Inference and Maximum Entropy Methods in Science and Engineering: 19th International Workshop. Edited by: Rychert JT, Erickson GJ, Smith CR. Melville: American Institute of Physics; 2001:274–284.
- Kossenkov A, Bidaut G, Ochs MF: Genes associated with prognosis in adenocarcinoma across studies at multiple institutions. In Methods of Microarray Data Analysis IV. Edited by: Johnson K, Lin S. Boston: Kluwer Academic; 2005:239.View Article
- Bidaut G, Suhre K, Claverie JM, Ochs MF: Bayesian decomposition analysis of bacterial phylogenomic profiles. Am J Pharmacogenomics 2005, 5: 63–70.View ArticlePubMed
- Mewes HW, Heumann K, Kaps A, Mayer K, Pfeiffer F, Stocker S, Frishman D: MIPS: a database for genomes and protein sequences. Nucleic Acids Res 1999, 27: 44–8.PubMed CentralView ArticlePubMed
- Christie KR, Weng S, Balakrishnan R, Costanzo MC, Dolinski K, Dwight SS, Engel SR, Feierbach B, Fisk DG, Hirschman JE, et al.: Saccharomyces Genome Database (SGD) provides tools to identify and analyze sequences from Saccharomyces cerevisiae and related sequences from other organisms. Nucleic Acids Res 2004, 32: D311–4.PubMed CentralView ArticlePubMed
- Cherry JM, Ball C, Weng S, Juvik G, Schmidt R, Adler C, Dunn B, Dwight S, Riles L, Mortimer RK, et al.: Genetic and physical maps of Saccharomyces cerevisiae. Nature 1997, 387: 67–73.PubMed CentralView ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.