Nucleic Acids Research Advance Access originally published online on March 9, 2009
Nucleic Acids Research 2009 37(7):e54; doi:10.1093/nar/gkp140
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Nucleic Acids Research, 2009, Vol. 37, No. 7 e54
Published by Oxford University Press 2009
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/2.0/uk/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
Methods Online |
Construction and application of a protein and genetic interaction network (yeast interactome)
1Laboratory of Molecular Genetics, National Institute of Environmental Health Sciences and 2Life Sciences Division, U.S. Army Research Office, Research Triangle Park, NC 27709, USA
*To whom correspondence should be addressed. Tel: +1 919 541 4792; Fax: +1 919 541 7613; Email: copelan1{at}niehs.nih.gov
Received October 28, 2008. Revised February 18, 2009. Accepted February 19, 2009.
| ABSTRACT |
|---|
|
|
|---|
Cytoscape is a bioinformatic data analysis and visualization platform that is well-suited to the analysis of gene expression data. To facilitate the analysis of yeast microarray data using Cytoscape, we constructed an interaction network (interactome) using the curated interaction data available from the Saccharomyces Genome Database (www.yeastgenome.org) and the database of yeast transcription factors at YEASTRACT (www.yeastract.com). These data were formatted and imported into Cytoscape using semi-automated methods, including Linux-based scripts, that simplified the process while minimizing the introduction of processing errors. The methods described for the construction of this yeast interactome are generally applicable to the construction of any interactome. Using Cytoscape, we illustrate the use of this interactome through the analysis of expression data from a recent yeast diauxic shift experiment. We also report and briefly describe the complex associations among transcription factors that result in the regulation of thousands of genes through coordinated changes in expression of dozens of transcription factors. These cells are thus able to sensitively regulate cellular metabolism in response to changes in genetic or environmental conditions through relatively small changes in the expression of large numbers of genes, affecting the entire yeast metabolome.
| INTRODUCTION |
|---|
|
|
|---|
Cytoscape (www.cytoscape.org) is an open source bioinformatics software platform originally intended for, but not limited to, the analysis of molecular interaction data associated with changes in gene expression and other data (1). Cytoscape's core distribution provides a basic set of features for data integration and visualization, with additional features available as plugins. Additionally, the visual display properties are highly customizable, including the use of annotation files that allow additional information to be visually represented in a more meaningful manner (Figure 1).
|
Several years ago, we determined that disruption of the POS5 gene in Saccharomyces cerevisiae results in a 50-fold increase in the reversion of a frameshift deletion in mtDNA and demonstrated that POS5 encodes a NAD(H) kinase, the sole source of NADP+ and NADPH in the mitochondrion of S. cerevisiae (2). In a recent follow-up study, we used a yeast microarray to evaluate the changes in gene expression in S. cerevisiae due to genetic and environmental factors associated with oxidative stress (3). To facilitate those analyses, we created a high-quality yeast interaction network (interactome) suitable for use in Cytoscape, illustrated through the analysis of data from a recent diauxic shift experiment in wild-type yeast cells that serves as an in-house reference source of yeast expression data. Analyses of these data additionally revealed that transcription factors and their target genes form highly complex, interconnected networks affecting all aspects of cellular metabolism in S. cerevisiae.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Strain
Saccharomyces cerevisiae strain YPH925 (ade2-101 cyh2 his3-
200 kar1-
15 leu2-
1 lys2-801 trp1-
63 ura3-52) (4) was employed for this work, and for convenience is referred to as being wild-type.
Bioinformatic platforms
Microarray fluorescence data (3) were imported into Rosetta Resolver (Rosetta Biosoftware, Seattle, WA, USA) for the estimation of random error by application of an error model that calculated the confidence limits (P-values) for the expression values. The Agilent GeneSpring Analysis Platform (Agilent Technologies, Palo Alto, CA, USA) was used for LOWESS data normalization. To aid the visualization and the analysis of the expression data, a S. cerevisiae protein–protein and protein–DNA interactome was constructed as described in the following section. The microarray expression data were then mapped onto this interactome using Cytoscape.
Genes associated with the highest scoring subnetworks were identified using the Cytoscape jActiveModules plugin. Active Modules are connected subnetworks within the interaction network whose genes show significant coordinated changes in mRNA expression over particular experimental conditions (1). The algorithm iteratively reduces network complexity by pinpointing regions whose states are perturbed by the conditions of interest, while removing false-positive interactions and interactions not involved in the perturbation response. Genes present in each of the five top-scoring networks in wild-type cells shifted to growth in glycerol were identified using the jActiveModules algorithm. Since many of the genes in each of these five subnetworks (186–187 genes each) were present in two or more of these subnetworks, for simplicity these groups of genes were combined, resulting in a single jActiveModules gene list.
Construction of a high-quality yeast interactome
The interactome described in this article was constructed in April 2008 using semi-automated methods to format the interaction data from the Saccharomyces Genome Database (SGD: www.yeastgenome.org) and the transcription factor data from YEASTRACT (5) (yeastract.com) into a form suitable for use with Cytoscape, as described in detail in the Supplementary Material (Supplementary file: Constructing a Yeast Interactome). Briefly, the interactions.tab file downloaded from SGD was processed by deleting unnecessary text (e.g. Bait and Hit) and columns and reclassifying the various interaction types as either pi (physical interactions, e.g. protein–protein) or gi (genetic interactions, e.g. synthetic lethality). This step dramatically simplified the visual display of the interaction types (edges) between the various nodes (genes; proteins), while allowing us to assign weights to each of the edges, reflecting the numbers of interactions documented between nodes. These interaction weights provided a measure of the number of times that two genes were found to interact with one another in a specific manner, from among the data curated at SGD from various sources.
Next, an interaction file containing a list of S. cerevisiae transcription factors and their documented target genes, RegulationTwoColumnTable_Documented_2008410_1839_1043605408.tsv, was downloaded from YEASTRACT. In preparation for use in Cytoscape, the yeast common gene names provided by YEASTRACT were converted to their systematic names, as found at SGD (e.g. POS5 was converted to YPL188W). Next, all letters in the gene names appearing in lowercase were converted to uppercase, again a requirement for Cytoscape. This list of genes was then processed through the Batch Download tool at SGD to identify rogue genes (e.g. MAL63 is not present in the in the systematic sequence of the SGD reference strain S288c; or, two or more genes sharing the same common name). Next, a column of interactions weights (all equal to 1) was appended to the transcription factor interaction file, for compatibility with the weighted SGD interaction data.
The SGD and YEASTRACT plain-text, tab-delimited interaction files were then concatenated as a single file (Supplementary file: pp_gi_tf.tab) and imported into Cytoscape using the import tool located under the File menu. At this time various annotation files, including gene-expression data and lists of stress response and mitochondrial genes, were also imported into Cytoscape. Lastly, the visual display properties of the nodes and edges were defined using the Cytoscape VizMapper tool.
The computer used for this work employed an Intel Pentium 4 CPU operating at 3.0 GHz, 1.5 GB of RAM, and the Microsoft Windows XP Professional Version 2002 Service Pack 2 operating system. To facilitate the steps summarized above associated with manipulating and formatting the raw interaction data files, simple perl and awk scripts were employed using Cygwin (http://www.cygwin.com/), a Linux-like environment for Windows [GNU bash shell, version 3.2.33(18)-release (i686-pc-cygwin]. On Macintosh and Linux-based operating systems, the awk and perl programming languages can be implemented directly in a command shell. Cytoscape is available for any major computer platform, including the Windows, Macintosh and Linux operating systems. All of the tools and source data described in this article are freely available from the indicated sources, while the yeast interactome described in this article is provided as Supplementary Material (Supplementary file: pp_gi_tf.cys) a Cytoscape session file with 6,188 nodes (genes/proteins) and 109,179 edges (interactions), that also includes the sample microarray expression data from our yeast diauxic shift experiment, plus the VizMapper visual display settings. For use with older versions of Cytoscape (or for use with other platforms), data files separately containing the diauxic shift expression data, the interactome (Supplemental file: pp_gi_tf.tab), node annotation files (lists of genes associated with the mitochondrion, the response to stress or transcription factors; common gene names; SGD gene descriptions), as well as the VizMapper visual display properties (vizmap.props) file, and lastly an Excel look-up table that can be used to convert common yeast names to their systematic gene name are included in the supplementary interactome files: (Supplementary interactome files: Mitochondrial_Gene_Names.txt; Stress_Response_Genes.txt; Transcription_Factor_Gene_List.txt; Common_Gene_Names.txt; SGD_Gene_Descriptions.txt; vizmap.props; Common_to_Systematic_Name_Lookup_Table.xls; WT_YPG_Shift_Expression_Data.pvals).
| RESULTS |
|---|
|
|
|---|
Changes in expression occurring during a glycerol-induced diauxic shift
In a related study (3), we examined and compared changes in gene expression in cells containing a deletion of the POS5 gene using Cytoscape (1), a bioinformatic data analysis and visualization tool. Here, we describe the construction of the robust yeast interactome used in those analyses. To better describe the application and versatility of this interactome, we describe the results from a parallel study from our laboratory that examined gene expression in a wild-type yeast strain grown to the mid-logarithmic phase of growth, then shifted to growth in glycerol for 2 h.
Similar to previous reports (6–12), we observed profound changes in gene expression following the switch from a fermentable carbon source (glucose) to the nonfermentable carbon source (glycerol), as summarized in Figure 2 and Table 1. Specifically, 3777 of the 6256 genes on the Agilent yeast chip (60.4%) were found to be differentially expressed at a significance level of <10–4 (at this level of significance,
0.6 false-positives were expected). To reduce this list of genes to a more meaningful and manageable dataset, we used the Cytoscape jActiveModules plugin to identify genes showing coordinated, significant changes in expression. The five top-scoring Active Modules subnetworks (jAM5-1 through jAM5-5) each contained 186 or 187 genes, as indicated in Figure 1. Since these subnetworks partially overlapped, these genes were combined into a single list which was further simplified by selecting genes that were greater than 5-fold up- or down-regulated, plus their associated transcription factors (Table 1).
|
|
Examining these genes, we noted the up-regulation of genes associated with mitochondrial function, gluconeogenesis, the TCA cycle, the β-oxidation of fatty acids, transport (including the uptake of amino acids), cell wall stability, copper and iron utilization (both required as prosthetic groups in the cytochromes in the electron transport chain) and glycerol and lactate utilization. Conversely, we observed down-regulation of genes associated with the accumulation of iron in the mitochondrion (required for heme and cytochrome biosyntheses), ribosomal subunit biosynthesis, and cellular growth—a response to glucose starvation.
The increased mitochondrial activity associated with oxidative phosphorylation, required for respiratory growth on nonfermentable carbon sources including ethanol and glycerol, summarized by Maris et al. (8), is accompanied by increased production of reactive oxygen species. In response, we found that numerous genes associated with the response to oxidative stress were up-regulated in our diauxic-shifted cells (Figure 3; Table 2), most notably CTA1 (catalase A, present in the peroxisomal and mitochondrial matrices; 40.2-fold up-regulated), HSP12 (plasma membrane heat-shock protein;18.4-fold), CTT1 (cytosolic catalase T; 7.2-fold), PRX1 (mitochondrial peroxiredoxin; 5.5-fold), MCR1 (mitochondrial NADH-cytochrome b5 reductase; 5.0-fold) and GPX1 (phospholipid hydroperoxide glutathione peroxidase; 4.6-fold). The superoxide dismutases encoded by SOD1 (cytosol; mitochondrial intermembrane space) and SOD2 (mitochondrial matrix) were modestly up-regulated (2.1 and 1.9-fold, respectively), indicating that the burden of the response to increased reactive oxygen species in this strain under these conditions is shared by the other antioxidant defense mechanisms (Cta1, Ctt1, etc.). Interestingly, many genes associated with the response to oxidative stress were down-regulated (Figure 3; Table 2), most notably GPX2 (cytoplasmic phospholipid hydroperoxide glutathione peroxidase; 5.7-fold), TRR1 (cytoplasmic thioredoxin reductase; 4.0-fold) and GSH1 (glutathione biosynthesis; 3.0-fold).
|
|
The numbers of genes differentially expressed in these cells, the modest changes in expression of SOD1 and SOD2, and the down-regulation of oxidative stress-related genes GPX2, TRR1 and GSH1 suggests that the regulation of gene expression in these cells is rather complex. Examining our yeast interactome, it is readily apparent that the regulation of gene expression in yeast is extraordinarily complex, with most genes simultaneously regulated by two or more transcription factors, as indicated in Figures 1–3
| DISCUSSION |
|---|
|
|
|---|
Our laboratory has a long-standing interest in exploring mitochondrial function and maintenance, including the stability and replication of the mitochondrial genome, and mutations and naturally occurring nucleotide polymorphisms associated with mitochondrial disease (13–15). Among the tools that we employ for these studies is the model organism, S. cerevisiae (2,3,13,16).
Energy demands in the facultative anaerobe S. cerevisiae are met under different physiological states when cells are grown on fermentable carbon sources such as glucose versus non-fermentable carbon sources such as ethanol or glycerol (8,17,18). Fermentation and glycolysis supply the cell with energy through the breakdown of glucose and other simple fermentable sugars; however, when the glucose concentration drops below
0.2%, the cells stop growing for a few hours as they undergo the diauxic shift, accompanied by transcriptional and translational changes including the mitochondrial biosynthesis. The cells then resume slower growth for a few generations by oxidative phosphorylation using the ethanol, glycerol and other byproducts accumulated during the fermentive stage of growth, before entering the stationary phase of growth.
Oxidative phosphorylation, which is dependent on mitochondrial activity and oxygen metabolism, provides the most efficient means of energy production (19,20). However, a deleterious consequence of oxidative phosphorylation is the production of reactive oxygen species in the mitochondrion—including the superoxide anion, hydrogen peroxide (H2O2) and the hydroxyl radical—that must be detoxified to minimize damage to nucleic acids, proteins, carbohydrates and lipids (21). The superoxide anion radical is reduced to H2O2 by superoxide dismutases (SOD), and further reduced to water by the antioxidant glutathione (GSH) and the enzymatic activity of catalases and peroxidases (21,22).
Cells switched to growth dependent on mitochondrial activity—the utilization of glycerol via oxidative phosphorylation—experience elevated levels of reactive oxygen species, evidenced by the increased nuclear mutational rates in cells grown in YPG versus YPD (2), a 28-fold increase in oxidative damage to mitochondrial proteins (2), and the changes in gene expression observed in this study, facilitated using Cytoscape and our yeast interactome. The genes that were most significantly affected by the switch from growth on glucose to growth on glycerol were identified using our interactome and the Cytoscape jActiveModules plugin (Figure 2; Table 1). As shown in Figure 2, transcription factors directly associated with the regulation of these genes include the stress response transcription factors Yap1 (required for oxidative stress tolerance), Rpn4 (stimulates expression of proteasome genes), Cad1 (multiple stress responses), Arr1 (resistance to arsenic compounds), Cat8 (derepression of genes following the diauxic shift), Sok2 (signal transduction), Stp2 (external amino acid permease) and Met4 (sulfur amino acid pathway). The involvement of several of these transcription factors in regulating the cellular response to the diauxic shift (e.g. Yap1) is not surprising given the increases in oxidative phosphorylation and reactive oxygen species as a consequence of increased mitochondrial function. Transcription factor Arr1 (Yap8), normally associated with the transcription of genes involved in resistance to arsenic compounds, appeared to directly coregulate the expression of 39 of the 95 genes shown in Figure 2, including up-regulated genes associated with the β-oxidation of fatty acids, carbohydrate metabolism and the TCA cycle, and the response to diauxic shift, suggesting a substantial role for this transcription factor in diauxic-shifted cells.
The complex regulatory nature of gene expression in S. cerevisiae is readily apparent from the interaction data displayed in Figures 1–3![]()
. It appears that transcription factors in S. cerevisiae form a highly interconnected self-regulatory subnetwork, while additionally regulating at least 5734 additional genes (our interactome: data not shown), indicating that substantial redundancy exists among the regulation and utilization of metabolic pathways. These cells may thus be able to respond quickly to changes in their external (e.g. adverse growth conditions) or internal (e.g. nonlethal mutations) environment by adjusting the regulation of their cellular metabolism via modest changes in gene expression involving hundreds or thousands of genes.
The transcription factors present in our yeast interactome appear to regulate most, but not all of the genes present in this interactome. Subtracting the 168 transcription factors and their regulated genes from the yeast interactome reveals that 286 of the 6188 genes present in the interactome (4.6%) are currently not associated with any transcription factor (data not shown). These genes, ranging from 3.1-fold up-regulated (TAM41) to 6.2-fold down-regulated (PAM18), passed through the SGD GO Slim Mapper, are variously associated with unknown biological processes (58 of 286 genes; 20.3%), transport (53/286; 18.5%), transcription (16/286; 5.6%), the cell cycle (11/286 genes; 3.8%), signal transduction (10/286; 3.5%) and amino acid metabolism (3/286 genes; 1.0%). Sixty-four of these genes (64/286; 22.4%) are annotated by SGD as being associated with the mitochondrion.
There is increasing interest regarding the application of bioinformatics and systems biology to the study of organisms and their regulatory mechanisms and metabolic profiles (23–28). The data provided in this study suggest that most genes are regulated in a highly complex manner by more than one transcription factor, and that bioinformatic tools such as Cytoscape—in conjunction with a robust interactome—may provide a useful framework for additional avenues of investigation. For example, by noting the transcription factors associated with specific groups of genes that are differentially expressed, the effect of deleting these transcription factors may be determined, at least partly. Finally, by applying methods similar to those used in the construction of the interactome described in this article, additional types of interaction data—for example those associated with protein kinases and their targets—can be readily incorporated.
| SUPPLEMENTARY DATA |
|---|
|
|
|---|
Supplementary Data are available at NAR Online.
| FUNDING |
|---|
|
|
|---|
National Institutes of Health, National Institute of Environmental Health Sciences (Z01-ES065078 to W.C.C.); US Army Research Office (W991NF-04-D-0001 and DAAG55-98-D-0002 to G.R.S.); National Research Council Research Associateship Award. Funding for open access charge: Intramural Research Program, National Institutes of Health, National Institute of Environmental Health Sciences.
Conflict of interest statement. None declared.
| ACKNOWLEDGEMENTS |
|---|
The authors would like to thank Astrid Haugen (NIEHS) and Leroy Worth (NIEHS) for critical reading of the article, and Astrid Haugen for discussions regarding the use of Cytoscape.
| REFERENCES |
|---|
|
|
|---|
- 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.
[Abstract/Free Full Text] - Strand MK, Stuart GR, Longley MJ, Graziewicz MA, Dominick OC, Copeland WC. POS5 gene of Saccharomyces cerevisiae encodes a mitochondrial NADH kinase required for stability of mitochondrial DNA. Eukaryot. Cell (2003) 2:809–820.
[Abstract/Free Full Text] - Stuart GA, Humble MM, Strand MK, Copeland WC. Transcriptional response to mitochondrial NADH kinase deficiency in Saccharomyces cerevisiae. Mitochondrion (2009) In press.
- Spencer F, Hugerat Y, Simchen G, Hurko O, Connelly C, Hieter P. Yeast kar1 mutants provide an effective method for YAC transfer to new hosts. Genomics (1994) 22:118–126.[CrossRef][Web of Science][Medline]
- Teixeira MC, Monteiro P, Jain P, Tenreiro S, Fernandes AR, Mira NP, Alenquer M, Freitas AT, Oliveira AL, Sá-Correia I. The YEASTRACT database: a tool for the analysis of transcription regulatory associations in Saccharomyces cerevisiae. Nucleic Acids Res. (2006) 34:D446–D451.
[Abstract/Free Full Text] - Bonander N, Ferndahl C, Mostad P, Wilks MDB, Chang C, Showe L, Gustafsson L, Larsson C, Bill RM. Transcriptome analysis of a respiratory Saccharomyces cerevisiae strain suggests the expression of its phenotype is glucose insensitive and predominantly controlled by Hap4, Cat8 and Mig1. BMC Genomics (2008) 9:365.[CrossRef][Medline]
- DeRisi JL, Iyer VR, Brown PO. Exploring the metabolic and genetic control of gene expression on a genomic scale. Science (1997) 278:680–686.
[Abstract/Free Full Text] - Maris AF, Assumpção ALK, Bonatto D, Brendel M, Henriques JAP. Diauxic shift-induced stress resistance against hydroperoxides in Saccharomyces cerevisiae is not an adaptive stress response and does not depend on functional mitochondria. Curr. Genet. (2001) 39:137–149.[CrossRef][Web of Science][Medline]
- Kuhn KM, DeRisi JL, Brown PO, Sarnow P. Global and specific translational regulation in the genomic response of Saccharomyces cerevisiae to a rapid transfer from a fermentable to a nonfermentable carbon source. Mol. Cell Biol. (2001) 21:916–927.
[Abstract/Free Full Text] - Young ET, Dombek KM, Tachibana C, Ideker T. Multiple pathways are co-regulated by the protein kinase Snf1 and the transcription factors Adr1 and Cat8. J. Biol. Chem. (2003) 278:26146–26158.
[Abstract/Free Full Text] - Ohlmeier S, Kastaniotis AJ, Hiltunen JK, Bergmann U. The yeast mitochondrial proteome, a study of fermentative and respiratory growth. J. Biol. Chem. (2004) 279:3956–3979.
[Abstract/Free Full Text] - Xie X, Wilkinson HH, Correa A, Lewis ZA, Bell-Pedersen D, Ebbole DJ. Transcriptional response to glucose starvation and functional analysis of a glucose transporter of Neurospora crassa. Fungal Genet. Biol. (2004) 41:1104–1119.[CrossRef][Web of Science][Medline]
- Stuart GR, Santos JH, Strand MK, Van Houten B, Copeland WC. Mitochondrial and nuclear DNA defects in Saccharomyces cerevisiae with mutations in DNA polymerase
associated with progressive external ophthalmoplegia. Hum. Mol. Genet. (2006) 15:363–374.[Abstract/Free Full Text] - Graziewicz MA, Longley MJ, Copeland WC. DNA polymerase
in mitochondrial DNA replication and repair. Chem. Rev. (2006) 106:383–405.[CrossRef][Web of Science][Medline] - Chan SSL, Copeland WC. DNA polymerase gamma and mitochondrial disease: Understanding the consequence of POLG mutations. Biochimica et Biophysica Acta (2008) In press [29 October, 2008, Epub ahead of print].
- Strand MK, Copeland WC. Measuring mtDNA mutation rates in Saccharomyces cerevisiae using the mtArg8 assay. Methods Mol. Biol. 197:151–157.
- Lagunas R. Energy metabolism of Saccharomyces cerevisiae discrepancy between ATP balance and known metabolic functions. Biochim. Biophys. Acta (1976) 440:661–674.[Medline]
- Rosenfeld E, Beauvoit B. Role of the non-respiratory pathways in the utilization of molecular oxygen by Saccharomyces cerevisiae. Yeast (2003) 20:1115–1144.[CrossRef][Web of Science][Medline]
- Lin S-J, Kaeberlein M, Andalis AA, Sturtz LA, Defossez P-A, Culotta VC, Fink GR, Guarente L. Calorie restriction extends Saccharomyces cerevisiae lifespan by increasing respiration. Nature (2002) 418:344–348.[CrossRef][Medline]
- Johnston M, Kim J-H. Glucose as a hormone: receptor-mediated glucose sensing in the yeast Saccharomyces cerevisiae. Biochem. Soc. Trans. (2005) 33:247–252.[CrossRef][Web of Science][Medline]
- Jackson MJ, Papa S, Bolaños J, Bruckdorfer R, Carlsen H, Elliott RM, Flier J, Griffiths HR, Heales S, Holst B, et al. Antioxidants, reactive oxygen and nitrogen species, gene induction and mitochondrial function. Mol. Aspects Med. (2002) 23:209–285.[CrossRef][Medline]
- Grant CM, Perrone G, Dawes IW. Glutathione and catalase provide overlapping defenses for protection against hydrogen peroxide in the yeast Saccharomyces cerevisiae. Biochem. Biophys. Res. Commun. (1998) 253:893–898.[CrossRef][Web of Science][Medline]
- Futcher B. Transcriptional regulatory networks and the yeast cell cycle. Curr. Opin. Cell Biol. (2002) 14:676–683.[CrossRef][Web of Science][Medline]
- 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.
[Abstract/Free Full Text] - Begley TJ, Samson LD. Network responses to DNA damaging agents. DNA Repair (2004) 3:1123–1132.[Medline]
- Tong AHY, 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.
[Abstract/Free Full Text] - Wu W-S, Li W-H, Chen B-S. Computational reconstruction of transcriptional regulatory modules of the yeast cell cycle. BMC Bioinformatics (2006) 7:421.[CrossRef][Medline]
- Boone C, Bussey H, Andrews BJ. Exploring genetic interactions and networks with yeast. Nat. Rev. Genet (2007) 8:437–449.[CrossRef][Medline]
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||


5-fold, plus the associated transcription factors, displayed using Cytoscape, and also shown at lower magnification in the lower left frame in 