
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
1 Biological Engineering Division and Center for Environmental Health Sciences and 2 Whitehead Institute for Biomedical Research, Massachusetts Institute of Technology, Cambridge, MA
Requests for reprints: Leona D. Samson, Biological Engineering Division and Center for Environmental Health Sciences, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139. Phone: (617) 258-7813; Fax: (617) 253-8099. E-mail: lsamson{at}mit.edu
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
Here we have systematically tested a set of 1615 haploid S. cerevisiae deletion strains to identify those deletion strains that display altered viable recovery after exposure to four classical DNA damaging agents, chosen because they are known to induce different kinds of DNA lesions (1). This approach has been termed genomic phenotyping because it systematically interrogates the genome of an organism to identify individual gene products that affect a particular phenotype, in this case the sensitive or resistant phenotype of S. cerevisiae strains on exposure to mutagens that produce alkylation damage, oxidative damage, or radiation damage. We have analyzed the contribution of one third of the yeast genome (1,615 S. cerevisiae gene products) to recovery after mutagen exposure in the framework of 12,232 protein-protein and protein-DNA interactions comprising the known yeast interactome [derived from both low- and high-throughput protein-protein and protein-DNA interaction studies (713)]. Our results reveal that many unexpected and hitherto uncharacterized pathways influence the recovery of eukaryotic cells after damaging agent exposure.
| Results and Discussion |
|---|
|
|
|---|
80 deletion strains, plus 6 control strains, namely the BY4741 parent (in triplicate), its
MMS-sensitive mag1
derivative, and its MMS-, 4NQO-, and UV-sensitive rev1
and rad14
derivatives
(1416). Although a t-BuOOH-sensitive control strain was not included among the controls, numerous
t-BuOOH-sensitive strains were detected in the screen (see below). Following a 60-h incubation, recovery was recorded by digital imaging of
colony growth, and each strain was scored for sensitivity or resistance (compared to wild type) using image analysis and a simple algorithm. In
essence, after adjusting for varying growth rates in the absence of exposure, treated strains displaying <67% of wild-type growth in at
least two of the three replicates were scored as sensitive, and strains displaying >150% of wild-type growth in at least two replicates
were scored as resistant; these parameters were chosen to minimize false positives while optimizing detection of true positives and were based on
extensive reconstruction experiments (Fig. 2 ). Moreover, all of the sensitivity and resistance calls made by the simple algorithm were
confirmed by visual inspection of the digitally imaged colonies.
|
|
Generation of a Genomic Phenotyping Database
To facilitate independent visual confirmation of the
computational assignment of sensitive and resistant strains, we compiled the imaged data generated from this study (i.e., >90,000 digital
images of the outgrowth from each spotted strain) into a genomic phenotyping database that is publicly accessible at http://GenomicPhenotyping.mit.edu (System Requirements: PC running
Internet Explorer). In this database, the imaged colonies for each spotted strain are arranged from left to right with increasing dose of agent,
generating visual killing curves for MMS, t-BuOOH, 4NQO, and UV. Examples are shown in Fig.
1C for the parent (BY4741), three controls (mag1
, rev1
, rad14
), a strain scored as
sensitive to MMS, t-BuOOH, and 4NQO (yml007w), and a strain scored as MMS resistant (yel033w). Individual web pages were
compiled for each of the 1615 deletion strains tested, displaying their triplicate visual killing curves for all four agents aligned with triplicates
for the parental strain.
The genomic phenotyping database was used to visually confirm all the phenotype assignments, and these are categorized in Table 1. We identified a total of 416 MMS-sensitive, 67 t-BuOOH-sensitive, 149 4NQO-sensitive, 44 UV-sensitive, 23 MMS-resistant, 6 t-BuOOH-resistant, 39 4NQO-resistant, and zero UV-resistant strains. The data were compared to those reported in other screens (17, 18) and, among the limited number of strains and treatments in common, 71% of the assignments were in agreement; considering differences in exposure conditions, the agreement is very good. It is important to point out that for the other screens (17, 18), a different range of mutagens was used and only single-exposure doses tested. In addition, a third study claimed to identify more than 100 genes that affect recovery after exposure to four damaging agents (nonoverlapping with those used here), but in their report only provided the identities of 20 of those genes (19). Here we have tested a wide range of mutagen doses and all of our results, including the primary data, are publicly available.
|
|
|
One of the largest categories of agent-sensitive strains is that comprised of strains deficient in various aspects of protein metabolism, namely protein degradation, processing, synthesis and targeting, as well as amino acid metabolism. Furthermore, processes involved in trafficking proteins from the endoplasmic reticulum (ER) to the Golgi for secretion and to the vacuole for degradation, appear to profoundly influence cellular recovery, especially after MMS exposure. At least part of the protein-related responses involves the unfolded protein response (UPR) signaled by the Ire1 endoplasmic reticulum-membrane-bound protein kinase (21), but it also involves protein degradation via both the proteosome; and the vacuole. From these phenotypic data, and from previous analyses of global transcriptional responses to MMS exposure (36), we propose that eukaryotic cells mount a concerted response to the induction of protein damage on exposure to these so-called DNA damaging agents and that in the absence of these responses, cells are more likely to die.
Several gene deletions affecting chromatin assembly, structure, and silencing were classified as MMS sensitive (msi1, spt8, htz1, yor023c, hho1, san1, ydr363w, yhr154w, hst3, hfi1, Ycl010c) and 4NQO sensitive (msi1, spt3, spt8, ydr363w, yhr154w). These gene products affect chromatin by influencing nucleosome assembly and by modulating histone acetylation (22, 23). Indeed, the Spt3, Spt8, Hfi1, and Ycl010c proteins comprise part of the Spt-Ada-Gcn5-histone acetyltransferase transcriptional activation complex (SAGAA), and it seems that SAGAA disruption produces a damage-sensitive phenotype. Proteins involved in histone acetylation may affect damage sensitivity by influencing the expression of genes needed for recovery. However, it is also possible that overall changes in acetylation of S. cerevisiae histones could both affect the reactivity of DNA with damaging agents and affect access for DNA repair proteins, analogous to affecting access for transcriptional factors (24, 25). Alternatively, in wild-type cells, it may be important to specifically modify chromatin structure adjacent to sites of DNA damage to provide targeted access for DNA repair proteins. Chromatin structural proteins can also influence agent toxicity, because deletion of the histone proteins, Hho1 and Htz1, generates MMS sensitivity. Again, we infer that such structural changes might influence the transcription of genes important for recovery, influence the reactivity of DNA with MMS, or influence the access of DNA repair proteins to damaged DNA. Whatever the mechanisms, it is clear that chromatin structure has a strong influence on damage recovery.
Some of the most sensitive mutants identified in this screen turn out to be deficient in the synthesis of ergosterol, a constituent of the cell membrane thought to contribute to membrane permeability and fluidity (erg3, 4, 5, 6, and 24). Increased permeability can readily explain increased sensitivity to chemical damaging agents, but we were surprised to find that erg6 cells are also modestly UV sensitive. It is possible that, in addition to affecting permeability, ergosterol deficiency diminishes membrane-mediated signal transduction events needed for damage recovery. Indeed, UV light initiates membrane receptor-mediated signaling in mammalian cells (2629), and erg6 mutants display altered receptor-mediated signaling in S. cerevisiae (30).
Manual Interactome Mapping of Genomic
Phenotyping Data
Navigating lists of genes affecting particular phenotypes becomes unwieldy at the genomic level, in particular when trying to
integrate such information with other genomic data sets. Nevertheless, we set out to analyze our genomic phenotyping data in the framework of 12,232
protein-protein and protein-DNA interactions comprising the known yeast interactome. Initial integration of our data with the yeast interactome was
facilitated using the Curagen protein-protein interaction database (http://www.curagen.com). Manual searching of interaction maps for the proteins the absence of which produced a
phenotype, revealed a number of interesting multiprotein interaction hubs, one of which is illustrated in Fig.
4
. From amongst the top 40 MMS-sensitive strains, we identified three gene products
that lie in the same interaction network, namely the Ymr032w/Cyk2 protein involved in cytokinesis, and two proteins of hitherto unknown function
(Yhl006c and Ydr078c). Also included in this network are the Hda1 histone deacetylase, the Hpr5/Sgs2 DNA helicase, plus another protein of unknown
function, Ylr423c, which were not represented in the 1615 strains tested. Because of their association with other proteins that contribute to MMS
resistance, we tested Hda1-, Hpr5/Sgs2-, and Ylr423c-deficient strains for their MMS sensitivity. All three were very sensitive (Fig. 4, data not shown), and so collectively, these six proteins contribute to a hitherto unknown damage recovery
pathway, now ripe for further investigation.
|
The highest scoring subnetworks are shown in Fig. 5B. Subnetworks (1)(4) are associated with MMS sensitivity and subnetworks (5) and (6) are associated with sensitivity to any one of the other three agents. These six large subnetworks comprised of between 18 and 29 proteins emerged from the Cytoscape analysis with P values ranging from 0.0013 to 0.014. Almost all of the proteins in each network contribute to damage recovery (indicated by the green notation). Importantly, the group of interacting proteins described in Fig. 4 was identified as a component of subnetwork (2) in Fig. 5B. In addition to the six large networks shown, several other smaller protein networks required for recovery after damage were also identified (shown at http://GenomicPhenotyping.mit.edu). Interestingly, searches for subnetworks of connected proteins the gene deletions of which lead to damage resistance, did not reveal any significant clusters within the interaction network.
Like the information in Table 1, Fig. 5B illustrates the engagement of many different cellular processes to aid the recovery of S. cerevisiae from damage; the processes embraced by each subnetwork are indicated in Table 1. Subnetwork (1) in Fig. 5B contains a group of DNA damage response proteins (Mec1, Rad17, Swi6, and Cln2) that are indirectly linked to the SAGAA chromatin remodeling complex proteins (Spt8, Spt7, Hfi1, and Ycl010c), via two proteins involved in transcriptional regulation (Yap1 and Med2). That a group of DNA damage response proteins should be physically linked to transcription factors, one of which is known to activate stress response genes (Yap1) (32), as well as to the SAGAA chromatin remodeling complex, supports the idea that the SAGAA complex may provide damage resistance by participating in transcriptional activation of genes whose products aid recovery. Subnetwork (2) is predominated by gene products involved in the degradation of proteins (Vsp9, Aap1, Pep3, Vsp41, Vma22, Vma6, Pep12, and Snf7), but in addition, as mentioned, contains the interacting hub of proteins found by brute force screening and shown in Fig. 4. Subnetwork (2) also contains several proteins of unknown function, and proteins involved in chromatin silencing (Spp1, Swd3, and Bre2), RNA processing (Npl3, Mud13, and Snu7), and amino acid metabolism (Asn2 and Lys14). Subnetworks (3) and (4) contain a large number of proteins involved in chromosome segregation and the associated changes in cytoskeleton required to achieve this; these include Rif2, Cin8, Tpd3, Bem1, Spo11, Mck1, Gbp2, Spo16, Rec104 for subnetwork (3) and Ssp1, Tem1, Hof1, Ste50, Nip100, and Arp1 for subnetwork (4). In addition, subnetworks (3) and (4) contain a number of proteins involved in protein synthesis and protein trafficking, namely Sse1, Spo7, His7, Rpl5, Ard1, Ssf1 for subnetwork (3) and Spc2, Tlg2 and Rsm27 for subnetwork (4). The absence of virtually any one of the proteins in each subnetwork renders cells less able to recover from MMS-induced damage.
Fig. 5B also shows two subnetworks, (5) and (6), for proteins identified as contributing to recovery from t-BuOOH, 4NQO, or UV exposure (phenotypic data for all three agents were combined). Both networks contain a number of proteins involved in protein synthesis and protein trafficking (Yke2, Ylr119w, Rpl18b, Rps17a, Ard1, Rpl8A, Rpl7A, Tef4, Sse1, Ilv1, Mrp4, Tlg2, Rsm27). Subnetwork (5), like subnetwork (1), contains a hub of proteins surrounding the Mec3 DNA damage response protein. However, in this case, the Mec3 hub does not link to the SAGAA histone acetylation complex, but rather links to a group of proteins involved in telomere maintenance and chromosome segregation (Rap1, Rif1, Cin8, and Top1). We infer that connection of the Mec3 DNA damage response hub of proteins to a different part of the interactome reflects the different pathways that are engaged for recovery from t-BuOOH, 4NQO, or UV versus recovery from MMS.
It is important to point out that the 23 gray nodes included in the protein subnetworks shown in Figs. 4 and 5B represent proteins which were not initially tested, simply because they were not represented in the set of 1615 deletion strains screened in this study. We therefore determined whether the inclusion of an untested protein in a subnetwork heavily populated by proteins known to contribute to a damage resistance phenotype, predicts whether that protein also contributes to damage resistance. We therefore tested the mutagen sensitivity of 16 of 23 proteins represented by gray nodes (the other 7 were encoded by essential genes and could not be tested). Eleven of the 16 scored as mutagen sensitive, representing a significant enrichment (with a P value of <0.0035) compared to the overall frequency of mutagen sensitivity observed in this study (535/1615). Thus, combining genomic phenotyping with interactome mapping can provide a predictive component to an analysis with incomplete genome coverage. Such predictive capability is important because it will likely be some years until entire mammalian genomes are represented in libraries of mutant mammalian cell lines, which will presumably be generated by systematic application of RNAi technology (33, 34).
Concluding
Remarks
Global analyses of biological systems are steadily unveiling new roles for old proteins and uncovering hitherto unsuspected
connections between diverse cellular pathways. Indeed, the connectedness of so many proteins in the yeast interactome, albeit dynamic and certain to
change depending on environmental cues, has fundamentally changed the way one must think about both the proximal and distal effects of eliminating
proteins from the network. The present study combines genomic phenotyping with a newly developed computational method for merging genomics databases
(Cytoscape) and demonstrates how systematic phenotypic assays may be directly linked to underlying molecular mechanisms. We have uncovered
an extremely diverse set of mechanisms for cellular recovery on mutagen exposure, mechanisms that may now be viewed as influencing the accumulation
of mutations in cell populations. Our next goal is to accelerate throughput for the biological screen to generate genomic phenotyping profiles for
very large numbers of agents, exposures, and growth conditions, for both yeast and mammalian cells. Hence, we will be able to address the dynamics of
how cells deploy interconnected protein hubs embedded in their interactome, viewing the biological system as a
whole.
| Materials and Methods |
|---|
|
|
|---|
Growth Conditions for Liquid Cultures
Ninety-six-well master plates containing individual deletion strains were thawed
and a 96-pin tool was used to transfer them to a sub-master plate (Corning, Corning, NY) containing 150 µl of YPD supplemented with G418 at 200
µg/ml. There were at least six empty wells on each master plate. The wild-type controls were spotted into empty wells in the plate along with
three damage-sensitive control strains, mag1
, rad14
, and rev1
. Strains in sub-master plates were
grown for 2 days at 30°C. Settled yeast cells in the sub-master plate were mixed using the 96-pin tool and replicated to another 96-well plate
containing 150 µl of YPD supplemented with G418 at 200 µg/ml and grown for 36 h at 30°C.
Growth Conditions for Spotted
Strains
Settled cells in each position of the 96-well plate were resuspended with 60 µl bursts of forced air from a Hydra liquid handling
apparatus (Robbins Scientific, Sunnyvale, CA), which on average provided a suspension of 107 cells/ml. Thirty microliters of cell
suspensions were drawn into each syringe of the Hydra and ninety-six 1-µl samples were spotted simultaneously onto an agar-containing plate. MMS,
t-BuOOH, and 4NQO were purchased from Aldrich, St. Louis, MO. UV radiation (254 nm) was supplied from a UV Stratalinker 2400 (Stratagene, La
Jolla, CA). Plates containing up to 96 strains were tested using the following conditions: no treatment, 0.01% MMS, 0.02% MMS,
0.025% MMS, 0.03% MMS, 0.25 mM t-BuOOH, 0.50 mM t-BuOOH, 0.75 mM t-BuOOH, 1
mM t-BuOOH, 0.1 µg/ml 4NQO, 0.2 µg/ml 4NQO, 0.3 µg/ml 4NQO, 0.4 µg/ml 4NQO, 25 J/m2 UV, 50 J/m2
UV, 75 J/m2 UV, and 100 J/m2 UV. Chemicals were added to cooled agar on plate preparation. UV treatments were performed after
strains had been spotted and dried on the plates. Strains were grown for 60 h at 30°C and imaged using a Gel Doc 1000 from BioRad (Hercules, CA)
running Quantity One software. Images were analyzed using ScanAlyze software to quantitate the pixel intensity of each spotted colony. All screens
were performed in triplicate using fresh liquid cultures.
Database Construction
At website http://GenomicPhenotyping.mit.edu, we have compiled images from 1275
plates and batch processed them using the software program U Lead Smart Saver Pro, version 3.0. This program divided each plate into 96 separate
plate pieces, each containing an image of a single strain; these were recompiled into strain-specific rows in a hierarchical website containing
>1700 web pages using in-house-prepared visual basic scripts via Microsoft applications. Table 1,
hierarchical clustering, and Cytoscape-generated figures, with links to relevant gene lists, can be found at this address. In addition, lists of gene
deletion strains tested and agent-specific phenotype assignments are found here, along with basic and advanced search methods that allow for
efficient mining of the database.
Interactome Mapping
A composite interaction network of 4,232 proteins was constructed from
12,232 previously characterized molecular interactions in yeast, including 5,003 two-hybrid interactions catalogued in BIND (11), 6,925 additional protein-protein interactions determined by coimmunoprecipitation studies (12), and 304 protein-DNA interactions recorded in the TRANSFAC database (13).
To evaluate phenotypic data on the filtered interaction network, we used a network scoring and search procedure as described by Ideker et
al. (31). Briefly, each protein was assigned a score of +1.25, -0.25, or +0.25
representing sensitivity, nonsensitivity, or missing data, respectively. Subnetwork scores zA were computed using the
sum
(where k is the
number of proteins in the subnetwork and zi is the score of subnetwork protein i) and were calibrated against a
"background" distribution of random sets of k proteins drawn independently of the network. Subnetworks with significantly high
scores were identified using a heuristic search algorithm based on simulated annealing. We ran the algorithm with parameters [N =
108, Ti = 1, Tf = 0.01, M = 20, dmin = 70] and required that the
identified subnetworks contain less than 30 proteins.
| Acknowledgements |
|---|
|
|
|---|
| Notes |
|---|
|
|
|---|
2 A Baltimore Whitehead Fellow. ![]()
3 Ellison American Cancer Society Research Professor. ![]()
Received October 29, 2002; accepted November 4, 2002.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
W. Ma, M. A. Resnick, and D. A. Gordenin Apn1 and Apn2 endonucleases prevent accumulation of repair-associated DNA breaks in budding yeast as revealed by direct chromosomal analysis Nucleic Acids Res., April 1, 2008; 36(6): 1836 - 1846. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. S. Plosky, E. G. Frank, D. A. Berry, G. P. Vennall, J. P. McDonald, and R. Woodgate Eukaryotic Y-family polymerases bypass a 3-methyl-2'-deoxyadenosine analog in vitro and methyl methanesulfonate-induced DNA damage in vivo Nucleic Acids Res., April 1, 2008; 36(7): 2152 - 2162. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y.-Q. Zhang and R. Rao Global Disruption of Cell Cycle Progression and Nutrient Response by the Antifungal Agent Amiodarone J. Biol. Chem., December 28, 2007; 282(52): 37844 - 37853. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. C. Reis and J. L. Campbell Contribution of Trf4/5 and the Nuclear Exosome to Genome Stability Through Regulation of Histone mRNA Levels in Saccharomyces cerevisiae Genetics, March 1, 2007; 175(3): 993 - 1010. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Komili and F. P. Roth Genetic interaction screens advance in reverse Genes & Dev., January 15, 2007; 21(2): 137 - 142. [Full Text] [PDF] |
||||
![]() |
H. C. Mak, M. Daly, B. Gruebel, and T. Ideker CellCircuits: a database of protein network models Nucleic Acids Res., January 12, 2007; 35(suppl_1): D538 - D545. [Abstract] [Full Text] [PDF] |
||||
![]() |
D.J. Tweats, A.D. Scott, C. Westmoreland, and P.L. Carmichael Determination of genetic toxicity and potential carcinogenicity in vitro--challenges post the Seventh Amendment to the European Cosmetics Directive Mutagenesis, January 1, 2007; 22(1): 5 - 13. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. M. Roberts Mining literature for systems biology Brief Bioinform, December 1, 2006; 7(4): 399 - 406. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. T. Workman, H. C. Mak, S. McCuine, J.-B. Tagne, M. Agarwal, O. Ozier, T. J. Begley, L. D. Samson, and T. Ideker A systems approach to mapping DNA damage response pathways. Science, May 19, 2006; 312(5776): 1054 - 1059. [Abstract] [Full Text] [PDF] |
||||
![]() |
R.-Y. Huang, M. Eddy, M. Vujcic, and D. Kowalski Genome-Wide Screen Identifies Genes Whose Inactivation Confer Resistance to Cisplatin in Saccharomyces cerevisiae Cancer Res., July 1, 2005; 65(13): 5890 - 5897. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Deng, J. A. Brown, D. You, and J. M. Brown Multiple Endonucleases Function to Repair Covalent Topoisomerase I Complexes in Saccharomyces cerevisiae Genetics, June 1, 2005; 170(2): 591 - 600. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. W. Xie, F. Jin, H. Hwang, S. Hwang, V. Anand, M. C. Duncan, and J. Huang Insights into TOR function and rapamycin response: Chemical genomic profiling by using a high-density cell array method PNAS, May 17, 2005; 102(20): 7215 - 7220. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Strange The end of "naive reductionism": rise of systems biology or renaissance of physiology? Am J Physiol Cell Physiol, May 1, 2005; 288(5): C968 - C974. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. C. Mallory, G. Crudden, B. L. Johnson, C. Mo, C. A. Pierson, M. Bard, and R. J. Craven Dap1p, a Heme-Binding Protein That Regulates the Cytochrome P450 Protein Erg11p/Cyp51p in Saccharomyces cerevisiae Mol. Cell. Biol., March 1, 2005; 25(5): 1669 - 1679. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. C. Game, M. S. Williamson, and C. Baccari X-Ray Survival Characteristics and Genetic Analysis for Nine Saccharomyces Deletion Mutants That Show Altered Radiation Sensitivity Genetics, January 1, 2005; 169(1): 51 - 63. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. R. Said, T. J. Begley, A. V. Oppenheim, D. A. Lauffenburger, and L. D. Samson Global network analysis of phenotypic effects: Protein networks and toxicity modulation in Saccharomyces cerevisiae PNAS, December 28, 2004; 101(52): 18006 - 18011. [Abstract] [Full Text] [PDF] |
||||
![]() |
Session 5: Signalling Networks in Tumorigenesis Toxicol Pathol, October 1, 2004; 32(6): 737 - 738. [PDF] |
||||
![]() |
H. I. Wu, J. A. Brown, M. J. Dorie, L. Lazzeroni, and J. M. Brown Genome-Wide Identification of Genes Conferring Resistance to the Anticancer Agents Cisplatin, Oxaliplatin, and Mitomycin C Cancer Res., June 1, 2004; 64(11): 3940 - 3948. [Abstract] [Full Text] [PDF] |