Validation and refinement of gene-regulatory pathways on a network of physical interactions
- Chen-Hsiang Yeang†1,
- H Craig Mak†2,
- Scott McCuine2,
- Christopher Workman2,
- Tommi Jaakkola3 and
- Trey Ideker2Email author
© Yeang et al.; licensee BioMed Central Ltd. 2005
Received: 9 March 2005
Accepted: 3 June 2005
Published: 1 July 2005
As genome-scale measurements lead to increasingly complex models of gene regulation, systematic approaches are needed to validate and refine these models. Towards this goal, we describe an automated procedure for prioritizing genetic perturbations in order to discriminate optimally between alternative models of a gene-regulatory network. Using this procedure, we evaluate 38 candidate regulatory networks in yeast and perform four high-priority gene knockout experiments. The refined networks support previously unknown regulatory mechanisms downstream of SOK2 and SWI4.
Recent advances in genomics and computational biology are enabling construction of large-scale models of gene-regulatory networks. High-throughput technologies such as automated sequencing , gene-expression arrays , chromatin immunoprecipitation , and yeast two-hybrid assays , each probe different aspects of the gene-regulatory system through genome-wide datasets. These data have spawned a variety of methods to infer the structure of gene-regulatory networks or to study their high-level properties, as recently reviewed .
Regulatory network models generated thus far in Escherichia coli and budding yeast (Saccharomyces cerevisiae) have been most often validated against functional databases or previous literature [6, 7]. In contrast, only a few studies have attempted to validate or refine models systematically [8–11]. However, if we are to accurately model large gene networks in complex organisms, including fly, worm, mouse, and human, automated procedures will be essential for analyzing the network, choosing the best new experiments to test the model, conducting the experiments, and integrating the resulting data.
The problem of choosing the best experiments to estimate a model, termed 'experimental design' or 'active learning', has been a significant area of research in statistics and machine learning [12–14]. Automating the experimental design process can greatly accelerate data collection and model building, leading to substantial savings in time, materials, and human effort. For these reasons, many industries such as electronic circuit fabrication and airplane manufacturing incorporate experimental design as an integral step in the design process [15, 16]. A promising application of experimental design for biological systems was presented by King et al. , who integrated computational modeling and experimental design to reconstruct a small, well studied metabolic pathway. Whether automated experimental design can be useful in a large and poorly characterized biological system with noisy data remains an open question.
We recently reported a procedure for inferring gene-regulatory network models by integrating gene-expression profiles with high-throughput measurements of protein interactions . Here we extend this procedure to incorporate automated design of new experiments. First, we use the previously described modeling procedure to generate a library of models corresponding to different gene-regulatory systems in yeast. Many of these models contain transcriptional interactions for which the regulatory effects (inducer versus repressor) are ambiguous and cannot be determined from publicly available expression profiles. Next, to address these ambiguities we implement a score function that ranks possible genetic perturbation experiments on the basis of their projected information content over the models. We perform four of the highest-ranking perturbations experimentally and integrate the data back into the model. The new data support two out of three novel regulatory pathways predicted to mediate expression changes downstream of the yeast transcriptional regulator SWI4.
Summary of physical regulatory models
We applied a previously described network-modeling procedure  to integrate three complementary sources of gene-regulatory information in yeast: 5,558 promoter-binding interactions for 106 transcription factors measured using chromatin immunoprecipitation followed by microarray chip hybridization (ChIP-chip) ; the set of all 15,116 pairwise protein-protein interactions recorded in the Database of Interacting Proteins as of April 2004 ; and a panel of mRNA expression profiles for 273 individual gene-deletion experiments . Software for performing the network-modeling procedure is available as a plug-in to the Cytoscape package [21, 22] on our supplementary website .
For each gene-deletion experiment, the modeling procedure identified the most probable paths of protein-protein and promoter-binding interactions that connect the deleted gene (the perturbation) to genes that were differentially expressed in response to the deletion (the effects of perturbation). Thus, a path represented one possible physical explanation by which a deleted gene regulates a second gene downstream. From the expression data, each interaction on a path was annotated with its probable direction of information flow and its probable regulatory effect as an inducer or repressor.
In total, the modeling process generated 4,836 paths, each explaining expression changes for a particular gene in one or more knockout experiments. Of the 965 interactions covered by paths, 194 had regulatory effects that were uniquely determined by the data, while regulatory effects of the remaining 771 interactions were ambiguous. For example, Figure 1b includes ambiguous interaction paths through SWI4, SOK2, and MSN4, explaining the observation that many genes for which the promoters are bound by Msn4p are upregulated in a swi4Δ knockout. This observation can be explained by several alternative annotations: one scenario is that SWI4 activates SOK2 and SOK2 represses MSN4 (Figure 1b), whereas another is that SWI4 represses SOK2 and SOK2 activates MSN4 (Figure 1c). These regulatory annotations could be uniquely determined by measuring the expression changes of genes downstream of MSN4 in the model in response to a sok2Δ deletion and an msn4Δ deletion (see below).
Internal validation for 21 of the 38 inferred models
Number of genes
Number of variants
Validated literature pathway
Enriched MIPS functions
Kss1/Fus3-Ste12 (mating response and filamentous growth)
Cell fate (1.48 × 10-7); metabolism (0.0067)
Sok2-Msn4 (PKA pathway)
Tup1-Hhf1 (histone regulation)
Protein synthesis (7.13 × 10-8)
Tup1/Ssn6-Nrg1 (glucose metabolism)
Transport (1.05 x 10.5); metabolism (5.41 × 10-4)
Tup1/Ssn6-α2/Mcm1 (mating response)
Cell fate (1.12 × 10-5);
Rpd3-Abf1 (histone modification)
Swi4-Ndd1-Ace2 (cell cycle)
Cell cycle (0.035)
Slt2-Rlm1/Swi4 (PKC pathway)
Med2-Gal4/Gcn4 (general transcription)
Cmd1-Cna1-Skn7 (calcium signaling)
Cell defense (6.33 × 10-6)
Metabolism (1.49 × 10-6);energy (0.04)
Cell defense (9.62 × 10-5)
Yap1-Cad1 (metal response)
Med2-Srb6-Gal4 (general transcription)
Med2-Gal11-Gal4 (general transcription)
Med2-Srb5-Gal4 (general transcription)
Ste12-Mcm1 (mating response)
Cell fate (4.55 × 10-8); homeostasis (0.0012); cell communication (0.0345)
Top-ranking knock-out experiments proposed for model discrimination
Regulator for meiosis and PKA pathway
Protein kinase of cell cycle
Stress response regulator
Regulator of glucose dependent genes
Regulator of cell cycle
Regulator of cell cycle
Protein kinase of cell wall integrity pathway
Regulator of stress response
Regulator of cellular respiration
Among the highest-priority experiments, Model 1 (Figure 1b) was the most often targeted, containing three of the top 10 highest-scoring genetic perturbations: sok2Δ, yap6Δ, and msn4Δ. A fourth perturbation to Model 1, hap4Δ, was also highly ranked (rank 34). Therefore, Model 1 was chosen for further experimentation.
Knockout strains corresponding to the high-ranking perturbations sok2Δ, yap6Δ, hap4Δ, and msn4Δ were grown in quadruplicate under conditions identical to those for the initial 273 knockouts by Hughes et al. . Gene-expression profiles were obtained for each knockout culture versus wild type using yeast genome microarrays. We sought to test the three regulatory cascades leading from SWI4 to SOK2 to either MSN4, HAP4, or YAP6 (Figure 1b). To verify these cascades independently of the model, we analyzed the expression patterns of gene sets known to be directly regulated by MSN4, HAP4, or YAP6 (obtained from the Proteome BioKnowledge Library ; see Additional data file 1). To normalize between our microarray procedures and those of Hughes et al., we also repeated the original swi4Δ expression profile, and filtered the above sets to select only those genes with expression changes that were reproducible (that is, same direction of change) between the Hughes et al. swi4Δ profile and our new profile. Expression changes were reproducible for 28 of 42 Msn4p-regulated genes, 11 of 29 Hap4p-regulated genes, and 64 of 119 Yap6p-regulated genes. Expression similarity among the genes in each filtered set was captured formally in a measure called 'coherence'; details of the computation of expression coherence and the selection of the gene sets are described further in Materials and methods and .
Results were qualitatively similar for the HAP4 pathway (Figure 3b). The gene set downstream of HAP4 was upregulated in the swi4Δ (p ≤ 10-2) and sok2Δ (p ≤ 9 × 10-4) knockouts but downregulated in hap4Δ (p ≤ 10-4). These results suggest that swi4Δ, sok2Δ, and hap4Δ deletions affect the set of genes immediately downstream of HAP4, supporting the SWI4-SOK2-HAP4 regulatory pathway hypothesis. In contrast to the MSN4 and HAP4 pathways, the gene set downstream of YAP6 had insignificant responses to all follow-up knockout experiments (Figure 3c). Thus, the existence of the SWI4-SOK2-YAP6 regulatory pathway was not supported by our validation experiments.
Automated model refinement
We used our modeling procedure to construct a new physical network model using the original 273 knockout gene-expression experiments of Hughes et al. combined with the new sok2Δ, hap4Δ, msn4Δ, and yap6Δ profiles. Overall, 60 protein-DNA interactions were disambiguated by our data: 50 interactions were resolved as definite inducers or repressors, whereas ten interactions were removed from the model because the expression of downstream genes did not change as a result of the knockout. In the updated Model 1, MSN4 and HAP4 were unambiguously annotated as inducers of downstream genes, SOK2 was annotated as a repressor of MSN4 and HAP4, and SWI4 was annotated as an inducer of SOK2 (Figure 3e). These results agree with our previous manually derived annotations (see 'Model validation' above).
We quantified the efficiency of our information-based approach by comparing it to two other methods of prioritizing gene knockout experiments: prioritizing hubs and prioritizing genes randomly. First, we generated a 'reference' model by fixing each ambiguous interaction in Models 1-37 to be an inducer or repressor. Assignments were chosen arbitrarily from the set of annotations that were consistent with the original knockout data. Next, we used each method (information, hub, or random) to iteratively 'learn' these assignments. In each iteration, we selected the highest-priority knockout experiment, simulated the resulting expression changes (up/down) using the reference model, updated the inferred model, and recorded the number of ambiguous interactions that were resolved. This iterative learning procedure was repeated 100 times.
We have used global expression profiles to validate models of transcriptional regulation inferred from protein-protein interactions, genome-wide location analysis, and expression data. A previously described network inference algorithm  identifies probable paths of physical interactions connecting a gene knockout to genes that are differentially expressed as a result of that knockout. The proposed validation strategy uses information gain as a criterion for choosing optimal knockouts to profile using microarray experiments. This strategy agrees with intuition, in that optimal knockouts typically target intermediate genes along the pathways under consideration. If an intermediate gene knockout fails to affect downstream genes in a pathway, that pathway is removed from the model.
The validated pathways point to a combination of previously documented and novel findings. First, in agreement with previous literature, we confirm that MSN4 and HAP4 are inducers [27, 28] and that SOK2 is a repressor . For instance, SOK2 is known to act downstream of protein kinase A (PKA) to repress genes involved in stress response, glycogen storage, and pseudohyphal growth . However, although SOK2 is thought to control these pathways via a transcriptional cascade, the components of this cascade have remained unclear. Here, we provide evidence for a model in which SOK2 acts as a negative regulator upstream of MSN4 and HAP4. Interestingly, MSN4 has been shown to activate stress-response genes , and HAP4 has been shown to activate genes involved in energy conservation and oxidative carbohydrate metabolism . Thus, we have identified a candidate model for the transcriptional cascade downstream of PKA signaling that mediates stress response. This model includes two novel regulatory pathways from SWI4 to SOK2 to MSN4 and from SWI4 to SOK2 to HAP4. The validation experiments do not support the third predicted pathway from SWI4 to SOK2 to YAP6.
In model simulations, choosing new gene knockout experiments with an information-theoretic approach significantly outperformed both random and hub-based selection. It also outperformed the observed experimental results: approximately 280 interactions were disambiguated after four simulated knockouts (Figure 4), whereas only 60 interactions were resolved due to the four actual knockouts sok2Δ, hap4Δ, msn4Δ, and yap6Δ. This difference in performance stems from key differences between the simulated and actual scenarios. In simulation, the four experiments are performed independently and iteratively, selecting the absolute highest-ranking knockout each time. In the actual study, four high-ranking experiments (but not the highest) are chosen to interrogate and maximally resolve a single pathway model, resulting in experiments that are highly co-dependent and performed simultaneously without intervening rounds of inference and experimental design. In addition, the simulation assumes that all interactions in the model are correct, along with one of the initial sets of inducer/repressor annotations. It therefore isolates the process of learning regulatory role annotations, whereas the actual procedure also serves to distinguish interactions as true versus false positives. Nevertheless, the simulation provides a useful comparison of experimental design methods relative to each other.
An important limitation of the single-gene knockout approach is that single perturbations do not identify pathway intermediates for which loss of function can be compensated by another gene. Furthermore, our approach may not identify regulatory pathways in which several transcription factors independently activate gene expression. Applying knockouts in combination may prove fruitful in these cases. For instance, approximately 4,000 double knockouts have been reported in yeast that lead to synthetic lethality: that is, a lethal phenotype that is not observed in either of the single knockouts individually . These interactions suggest regulatory relationships which could be incorporated into future work.
Scientific discovery is an iterative process of building models to explain experimental observations and validating models with new experiments . Experimental design is the essential link between these two aspects. Here we have explored a framework for modeling transcriptional networks in which experimental design and validation are central features. This framework is based on computational analysis and expression microarrays, both of which are amenable to automation, suggesting a high-throughput strategy for mapping gene-regulatory pathways.
Materials and methods
Model building and inference
Physical mechanisms of transcriptional regulation were modeled using an approach described previously . Briefly, we postulated that the regulatory effects of deleting a gene are propagated along paths of physical interactions (protein-protein and protein-DNA). We formalized the properties of these paths and interactions using a factor graph  and found the most probable set of paths using the max-product algorithm . The resulting set of paths was partitioned into independent network models, also as described previously . The raw data used in the modeling procedure included 5,558 promoter-binding interactions (at p-value < 0.001) for 106 transcription factors , the set of all 15,166 pairwise protein-protein interactions recorded in the Database of Interacting Proteins as of April 2004 , and mRNA expression profiles for 273 individual gene deletion experiments . Expression changes with a p-value < 0.02 were considered significant.
We calculated the expected information gain for each of the 4,756 possible non-lethal single-gene deletion experiments that were not included in the set of 273 deletions used to generate our network models. Intuitively, information gain measures (the logarithm of) the number of ambiguous annotations in the model that are likely to be determined after generating a yeast-genome expression profile in response to a particular gene deletion under consideration. Each gene-deletion experiment predicts a distinct expression profile given a particular configuration of model annotations. Experiments with high information gain are those for which the predicted expression profiles are highly variable over the set of possible annotations. In these cases, only one (or at most a few) of the predicted profiles will match the true observed profile, efficiently constraining the space of possible model annotations.
The information gain discussed above arises from the expected value of information calculations in statistical decision theory . Here we describe the score more directly in terms of reduction of model entropy. The entropy of a set of ambiguous model annotations is given by:
The expected information gain is the difference between the entropies before and after a hypothetical experiment:
where Y e denotes the vector of predicted expression changes for each gene in the model under experiment e. The conditional entropy H(M|Y e ) requires us to consider all possible models and corresponding outcomes resulting from experiment e. Direct enumeration of all values of M and Y e is impractical; instead, we make several simplifying approximations as described at .
Expression profiling experiments were based on the wild-type diploid BY4743 and homozygous gene knockout strains derived from this parent  (Invitrogen), with cultures grown identically to those of Hughes et al. . Labeled cDNA from each gene knockout strain was co-hybridized versus wild type cDNA in quadruplicate two-color microarray hybridizations. Total RNA was isolated by hot acid phenol extraction, purified to mRNA (Ambion PolyAPure kits), and labeled with Cy3 or Cy5 by direct incorporation (Amersham CyScribe First-Strand cDNA Labeling Kit). DNA microarrays were spotted from the Yeast Genome Oligo Set v1.1 (Qiagen) on Corning UltraGAPS slides using an OmniGrid 100 robot (Genomic Solutions). Lyophilized Cy3- and Cy5-labeled samples were resuspended in 50 μl buffer (5× SSC, 0.1% SDS, 1× Denhardt's solution, 25% formamide) and co-hybridized at 42°C beneath a coverslip for 15 h. Arrays were imaged at 10 μm resolution using a ScanArray Lite instrument (PerkinElmer). Raw quantitated background intensities were smoothed using a 7 × 7 median filter, separately for the Cy3 and Cy5 channels, and data were corrected for cyanine-dye dependent bias using a Qspline normalization . The VERA/SAM package  was used to assign a log-likelihood statistic λ with each gene, indicating its significance of differential expression in each experiment. Microarray expression data are deposited in the ArrayExpress database  under accession numbers A-MEXP-217 (Arrays) and E-MEXP-351 (Experiments).
The expression coherence of a set of genes measures whether the expression levels of these genes behave similarly in a particular experiment. Each gene i in gene-deletion experiment e has an expression ratio r ie (versus wild type) and associated p-value p ie of differential expression. First, we filter out insignificant expression changes with a p-value > 0.5. Then, we use the inverse Gaussian cumulative distribution function, Φ-1, to convert each remaining p-value into a z-score [36, 37]:
z ie = Φ-1 (1 - p ie )
Next, we compute a 'signed z-score' by multiplying z by +1 if the expression level is increasing and by -1 if it is decreasing. The average signed z-score for a gene subset of size N is computed as:
Gene sets with expression changes that are significant and in the same direction result in large Z-values. A distribution of Z values obtained from random gene sets of size N was used to determine a p-value for each expression coherence score.
Additional data files
Additonal data is available with the online version of this paper. Additional data file 1 contains Tables S1-S4 and wiring diagram illustrations for Models 0-44. Table S1 gives the internal validation for 17 out of 24 restricted network models; Table S2 lists the correlations between swi4δ and gcn4δ data and Rosetta and the new experiments; Table S3 gives the restricted subsets used to evaluate the reproducibility; and Table S4 gives the gene sets for external validation.
We are grateful to Owen Ozier, Ryan Kelley, and Rowan Christmas for their valuable assistance with model visualization, and to Julia Zeitlinger for commenting on the manuscript. C.M., C.W., and S.M. were supported by NIGMS grant GM070743-01 and NSF grant CCF-0425926. T.I. was supported by a David and Lucille Packard Fellowship award. C.Y. and T.J. were supported in part by NIH grant(s) GM68762 and GM69676.
- Hood LE, Hunkapiller MW, Smith LM: Automated DNA sequencing and analysis of the human genome. Genomics. 1987, 1: 201-212. 10.1016/0888-7543(87)90046-2.PubMedView ArticleGoogle Scholar
- Lockhart D, Winzeler E: Genomics, gene expression and DNA arrays. Nature. 2000, 405: 827-836. 10.1038/35015701.PubMedView ArticleGoogle Scholar
- Lee TI, Rinaldi NJ, Robert F, Odom DT, Bar-Joseph Z, Gerber GK, Hannett NM, Harbison CT, Thompson CM, Simon I, et al: Transcriptional regulatory networks in Saccharomyces cerevisiae. Science. 2002, 298: 799-804. 10.1126/science.1075090.PubMedView ArticleGoogle Scholar
- Uetz P, Giot L, Cagney G, Mansfield TA, Judson RS, Knight JR, Lockshon D, Narayan V, Srinivasan M, Pochart P, et al: A comprehensive analysis of protein-protein interactions in Saccharomyces cerevisiae. Nature. 2000, 403: 623-627. 10.1038/35001009.PubMedView ArticleGoogle Scholar
- de Jong H: Modeling and simulation of genetic regulatory systems: a literature review. J Comput Biol. 2002, 9: 67-103. 10.1089/10665270252833208.PubMedView ArticleGoogle Scholar
- Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95: 14863-14868. 10.1073/pnas.95.25.14863.PubMedPubMed CentralView ArticleGoogle Scholar
- Segal E, Shapira M, Regev A, Pe'er D, Botstein D, Koller D, Friedman N: Module networks: identifying regulatory modules and their condition-specific regulators from gene expression data. Nat Genet. 2003, 34: 166-176.PubMedView ArticleGoogle Scholar
- Akutsu T, Kuhara S, Maruyama O, Miyano S: A system for identifying genetic networks in gene expression patterns produced by gene disruptions and overexpressions. Genome Inform Ser Workshop. 1998, 9: 151-160.Google Scholar
- Wagner A: How to reconstruct a large genetic network from n gene perturbations in fewer than n(2) easy steps. Bioinformatics. 2001, 17: 1183-1197. 10.1093/bioinformatics/17.12.1183.PubMedView ArticleGoogle Scholar
- Ideker T, Thorsson V, Karp RM: Discovery of regulatory interactions through perturbation: inference and experimental design. Pac Symp Biocomput. 2000, 305-316.Google Scholar
- Ideker T, Thorsson V, Ranish JA, Christmas R, Buhler J, Eng JK, Bumgarner R, Goodlett DR, Aebersold R, Hood L: Integrated genomic and proteomic analysis of a systematically perturbed metabolic network. Science. 2001, 292: 929-934. 10.1126/science.292.5518.929.PubMedView ArticleGoogle Scholar
- Raiffa H, Shlaifer R: Applied Statistical Decision Theory. 1962, Cambridge, MA: MIT PressGoogle Scholar
- Fedorov FF: Theory of Optimal Experimental Design. 1972, New York: Academic PressGoogle Scholar
- Tong S, Koller D: Active learning for parameter estimation in Bayesian networks. Proc 13th Conf Neural Information Processing. 2000, Tübingen: Neural Information Processing Systems, 647-563. [http://www.nips.cc]Google Scholar
- Abadir MS, Ferguson J, Kirkland TE: Logic design verification via test generation. IEEE Trans Computer-Aided Design. 1988, 7: 138-148. 10.1109/43.3141.View ArticleGoogle Scholar
- Rea C, Settle MA: An automated test approach for US Air Force fighter engines. IEEE Aerospace Electron Syst Mag. 1996, 11: 24-28. 10.1109/62.486733.View ArticleGoogle Scholar
- King RD, Whelan KE, Jones FM, Reiser PG, Bryant CH, Muggleton SH, Kell DB, Oliver SG: Functional genomic hypothesis generation and experimentation by a robot scientist. Nature. 2004, 427: 247-252. 10.1038/nature02236.PubMedView ArticleGoogle Scholar
- Yeang CH, Ideker T, Jaakkola T: Physical network models. J Comput Biol. 2004, 11: 243-262. 10.1089/1066527041410382.PubMedView ArticleGoogle Scholar
- Deane CM, Salwinski L, Xenarios I, Eisenberg D: Protein interactions: two methods for assessment of the reliability of high throughput observations. Mol Cell Proteomics. 2002, 1: 349-356. 10.1074/mcp.M100037-MCP200.PubMedView ArticleGoogle Scholar
- 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-126. 10.1016/S0092-8674(00)00015-5.PubMedView ArticleGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13: 2498-2504. 10.1101/gr.1239303.PubMedPubMed CentralView ArticleGoogle Scholar
- Cytoscape. [http://www.cytoscape.org]
- Cell Circuits Pathway Database. [http://www.cellcircuits.org/Yeang2005]
- Natarajan K, Meyer MR, Jackson BM, Slade D, Roberts C, Hinnebusch AG, Marton MJ: Transcriptional profiling shows that Gcn4p is a master regulator of gene expression during amino acid starvation in yeast. Mol Cell Biol. 2001, 21: 4347-4368. 10.1128/MCB.21.13.4347-4368.2001.PubMedPubMed CentralView ArticleGoogle Scholar
- Munich Information Center for Protein Sequences. [http://mips.gsf.de]
- Proteome BioKnowledge Library. [http://www.proteome.com]
- Blom J, De Mattos MJ, Grivell LA: Redirection of the respiro-fermentative flux distribution in Saccharomyces cerevisiae by overexpression of the transcription factor Hap4p. Appl Environ Microbiol. 2000, 66: 1970-1973. 10.1128/AEM.66.5.1970-1973.2000.PubMedPubMed CentralView ArticleGoogle Scholar
- Smith A, Ward MP, Garrett S: Yeast PKA represses Msn2p/Msn4p-dependent gene expression to regulate growth, stress response and glycogen accumulation. EMBO J. 1998, 17: 3556-3564. 10.1093/emboj/17.13.3556.PubMedPubMed CentralView ArticleGoogle Scholar
- Ward MP, Gimeno CJ, Fink GR, Garrett S: SOK2 may regulate cyclic AMP-dependent protein kinase-stimulated growth and pseudohyphal development by repressing transcription. Mol Cell Biol. 1995, 15: 6854-6863.PubMedPubMed CentralView ArticleGoogle Scholar
- Tong AH, Lesage G, Bader GD, Ding H, Xu H, Xin X, Young J, Berriz GF, Brost RL, Chang M, et al: Global mapping of the yeast genetic interaction network. Science. 2004, 303: 808-813. 10.1126/science.1091317.PubMedView ArticleGoogle Scholar
- Popper K: The Logic of Scientific Discovery. 1959, New York: Basic BooksGoogle Scholar
- Kschischang F, Frey B, Loeliger H: Factor graphs and the sum-product algorithm. IEEE Trans Inform Theory. 2001, 47: 498-519. 10.1109/18.910572.View ArticleGoogle Scholar
- Winzeler EA, Shoemaker DD, Astromoff A, Liang H, Anderson K, Andre B, Bangham R, Benito R, Boeke JD, Bussey H, et al: Functional characterization of the S. cerevisiae genome by gene deletion and parallel analysis. Science. 1999, 285: 901-906. 10.1126/science.285.5429.901.PubMedView ArticleGoogle Scholar
- Workman C, Jensen LJ, Jarmer H, Berka R, Gautier L, Nielser HB, Saxild HH, Nielsen C, Brunak S, Knudsen S: A new non-linear normalization method for reducing variability in DNA microarray experiments. Genome Biol. 2002, 3: research 0048.1-0048.16. 10.1186/gb-2002-3-9-research0048.View ArticleGoogle Scholar
- ArrayExpress Gene Expression Database. [http://www.ebi.ac.uk/arrayexpress]
- Ideker T, Thorsson V, Siegel A, Hood L: Testing for differentially-expressed genes by maximum likelihood analysis of microarray data. J Comput Biol. 2000, 7: 805-817. 10.1089/10665270050514945.PubMedView ArticleGoogle Scholar
- Ideker T, Ozier O, Schwikowski B, Siegel AF: Discovering regulatory and signalling circuits in molecular interaction networks. Bioinformatics. 2002, 18 (Suppl 1): S233-S240.PubMedView ArticleGoogle Scholar
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.