Summary
Packaging of DNA into chromatin has a profound impact on gene expression. To understand how changes in chromatin influence transcription, we analyzed 165 mutants of chromatin machinery components in Saccharomyces cerevisiae. mRNA expression patterns change in 80% of mutants, always with specific effects, even for loss of widespread histone marks. The data is assembled into a network of chromatin interaction pathways. The network is function-based, has a branched, interconnected topology and lacks strict one-to-one relationships between complexes. Chromatin pathways are not separate entities for different gene sets, but share many components. The study evaluates which interactions are important for which genes and predicts new interactions, for example between Paf1C and Set3C, as well as a role for Mediator in subtelomeric silencing. The results indicate the presence of gene-dependent effects that go beyond context-dependent binding of chromatin factors and provide a framework for understanding how specificity is achieved through regulating chromatin.
Introduction
In eukaryotes, DNA is densely packed into a higher order structure called chromatin. This has a profound impact on processes that work on DNA, such as replication or gene expression (Campos and Reinberg, 2009; Narlikar et al., 2002). Cells therefore contain various protein complexes that regulate chromatin structure. The basic unit of chromatin is a nucleosome, formed by 147 base pairs of DNA wrapped around an octamer of two copies of the histones H2A, H2B, H3 and H4 (Richmond and Davey, 2003). Chromatin regulators include nucleosome remodelers, histone chaperones and histone modifying complexes (Narlikar et al., 2002). Nucleosome remodeling complexes slide or evict nucleosomes and are also involved in deposition of histones and their variants. Histone modifier complexes covalently modify histone tails with different marks. Besides influencing nucleosome turn-over and altering physical properties such as chromatin condensation, specific modifications also serve as recognition sites for other proteins. These effectors further regulate chromatin structure or facilitate the process of gene expression itself (Campos and Reinberg, 2009).
An elegant model has been put forward to explain the consequences of chromatin modifications (Strahl and Allis, 2000). In the histone code hypothesis, different combinations of modifications form a code that is read by other proteins to influence downstream events. Although the location of many histone marks correlate with particular expression states, stringent evidence for causal relationships is often missing (Rando and Chang, 2009). Furthermore, the discovery that the same histone modification may be bound by different effectors, each mediating different downstream events, also calls into question the existence of a strictly rigid code (Berger, 2007). The consequences of histone modifications are currently being explained by context-dependent binding of effector complexes (Campos and Reinberg, 2009; Lee et al., 2010). The nature of this context is only starting to be investigated. As with the histone code hypothesis itself, proposals about context-dependent binding are based mainly on studies of individual genes. One purpose of this study is therefore to determine to what extent either context-dependency or a code applies to different chromatin interactions when assayed across an entire genome.
A related question is how different chromatin interactions work together. The general architecture of chromatin interaction pathways and how this may vary for different genes, is as yet unexplored. To understand the effects of different chromatin states, the focus of many studies is on the binding of effector proteins. While this is crucial for understanding mechanism, it can result in ignoring the question of whether a binding event has further downstream consequences, for example on gene expression. A second aim of this study is therefore to investigate interactions as manifested by their downstream consequences on gene expression.
Genome-wide expression analysis has previously been applied to study the role of many individual regulators. The use of different microarray platforms, different genetic backgrounds and different growth conditions in these previous studies, confounds proper comparative analyses. Here we analyze the interplay between gene expression and chromatin by determining expression profiles for perturbing the majority of chromatin regulatory machineries in Saccharomyces cerevisiae under identical conditions. This was achieved by DNA microarray expression-profiling 165 yeast strains, each bearing a mutation in a different chromatin factor. The results show a remarkable degree of specificity, also for mutants resulting in loss of widespread histone marks. The data is analyzed at three levels of complexity: analysis of individual profiles to determine cellular roles, analysis of protein complexes to examine subunit relationships and analysis of relationships between complexes to investigate the architecture of interaction pathways. The result is a first function-based network of chromatin interactions. The network reveals that individual chromatin regulators are almost all functionally connected to others and form pathways that branch and interconnect at different levels. The study shows how elements of the histone code and context-dependent binding of chromatin are superimposed to form chromatin interaction pathways. Removal of individual chromatin factors has much more specific and restricted effects on gene expression than is predicted by location. This suggests the presence of additional gene-dependent mechanisms that go beyond context-dependent binding to achieve specificity. The network and underlying data therefore provide a framework for investigating how globally acting chromatin regulators facilitate specific responses.
Results and Discussion
DNA microarray expression profiles were generated for viable deletion mutants of chromatin regulators in S. cerevisiae. Besides factors that directly regulate chromatin, such as nucleosome remodelers (e.g. SWI/SNF, RSC, INO80), histone chaperones/nucleosome assembly factors (e.g. FACT, CAF-1, Rtt106) and histone modifiers (e.g. Set1/COMPASS, Rpd3L/S, NuA4, NuA3, HAT-B), putative chromatin factors were included, as well as coregulators such as Tup1-Ssn6, Ccr4-Not, the Paf1 complex (Paf1C) and SAGA. The result is coverage of over 30 complexes with 174 mutants (Tables S1 and S2). Each mutant was analyzed four times, from two independent cultures on dual-channel microarrays using a batch of wild type (wt) RNA as common reference. To counter other sources of variation, sets of mutants were grown alongside additional wt cultures and these were all profiled in parallel using automated, robotic procedures (Experimental Procedures). Aneuploidy, incorrect deletions and spurious mutations were identified in 13% of strains. These were remade and reprofiled, resulting in data for 165 mutants that passed all quality controls (Supplemental Experimental Procedures, Table S2).
The majority of chromatin regulators have specific effects on gene expression
Throughout the study a p-value of 0.05 in combination with a fold-change (FC) of 1.7 is applied as a threshold for calling a change in mRNA expression significant. The reported FC and p-value is based on the average of the four replicates versus the average of 200 wt profiles. Based on the variation observed within the wt cultures, a threshold can be applied to distinguish mutants with an expression-profile different from wt (Experimental Procedures). By these criteria, 80% of the mutants have altered gene expression. Strikingly, only 29% of kinase deletions and only 15% of phosphatase deletions exhibit expression-profiles different from wt (van Wageningen et al., 2010). This underscores the importance of chromatin complexes in regulation of gene expression.
A second general outcome of this study is that inactivation of chromatin factors result in quite specific effects on gene expression. A two-dimensional hierarchical cluster diagram summarizes all the changes (Fig. 1). In the 80% of mutants with an altered profile, on average 116 genes (2%) show significant changes (p<0.05, FC>1.7). The largest effect is for ssn6Δ with 973 differentially expressed genes (16%). The possibility that data normalization had masked global effects was ruled out by applying a spiked-in external RNA control strategy (van de Peppel et al., 2003), (Supplemental Experimental Procedures). Focus on viable deletion mutants results in coverage of an estimated 78% of known chromatin regulators (Table S1). Although inviable mutants likely have broader effects, the specificity observed here is nevertheless noteworthy since several of these factors are responsible for laying down general marks such as H3K4 (Set1), H3K36 (Set2) and H3K79 (Dot1) methylation. Genome-wide location maps of these histone marks show widespread location (Pokholok et al., 2005). The effects of their removal result in much more specific effects, with changes in expression of 55 genes in set1Δ, 72 genes in set2Δ and 2 genes in dot1Δ (Fig. 1). Disconnect between location and effect has been noted in individual cases before (Venkatasubrahmanyam et al., 2007; Rando and Chang, 2009). The specific nature of all the perturbation signatures (Fig. 1) suggests that this disconnect is a general feature of chromatin regulation. Lack of an effect does not necessarily imply that a factor is not operating on such a gene. Rather, the specific nature of the effects indicates that the function of the factor is important for only a specific subset of genes under this condition.
Fig. 1. Perturbation of chromatin regulators leads to specific effects on gene expression.
Unsupervised hierarchical cluster diagram (cosine correlation) of all the mutants with profiles that differ from wt (Experimental Procedures) and all genes that change significantly (p<0.05, FC>1.7) in any mutant. The dendrograms indicate relationships between genes (top) and between mutants (right). The red part (right) depicts the relationship between Paf1C, Rad6/Bre and Set1C. Mutants are color-coded according to protein complex. Fold-change is indicated by the color-scale, with yellow for up-, blue for downregulation and black for no change, versus the average wt. dot1Δ is not included because the number of genes changing classifies it as similar to wt.
Assigning cellular roles to chromatin regulators
The specificity of the perturbation signatures (Fig. 1) can be exploited in several ways. In a first systematic analysis, cellular roles were assigned to chromatin factors by determining which functional categories of Gene Ontology (GO) are enriched in the signatures (Fig. S1, summarized in Fig. 2A). Many previously reported regulatory roles are confirmed. Regulation of middle-sporulation genes by Sum1/Rfm1/Hst1 (Xie et al., 1999) is reflected by enrichment for the GO Slim term ‘sporulation’ in these profiles (Fig. 2A). Signatures of the Rpd3L complex are enriched for ‘meiosis’ and signatures of the Sir complex are enriched for ‘conjugation’ (Fig. 2A), all as expected (Mallory and Strich, 2003; Rine and Herskowitz, 1987). New roles are also indicated. For example, eight complexes have two or more signatures enriched for GO terms associated with mating (Fig. S1, Fig. 2A “conjugation”). Reduced mating efficiency has been shown for four of these complexes: Sir, Mediator, Rad6/Bre1 and SAGA (Holstege et al., 1998; Rine and Herskowitz, 1987; Huang et al., 1997; Roberts and Winston, 1996). Deletion of two subunits in each of the remaining four complexes (Rpd3L, Isw2C, RSC, SWI/SNF) were tested. All show a significant reduction in mating (Fig. 2B) indicating involvement in this process. Other new roles are predicted (Fig. 2A, Fig. S1) illustrating how the data can be employed to link chromatin factors to specific cellular roles.
Fig. 2. Specific cellular roles for chromatin machinery components.
(A) Enrichment for Gene Ontology Slim terms (top) depicted for each mutant (right), by plotting the significance (hypergeometric test). Yellow indicates enrichment in upregulated genes, blue indicates enrichment in downregulated genes, with a threshold of p<0.01. Fig. S1 depicts all Gene Ontology terms. (B) Histogram showing the mean mating efficiency for mutants (bottom), two from each complex with no previously reported mating defect, with ste7Δ as control. Error bars reflect standard deviations of duplicates. Significance was calculated by two-sided T-test. *p<0.05, **p<0.01, ***p<0.005.
Structure-function analyses of protein complexes that regulate chromatin
The specificity of the signatures also facilitates investigation of functional relationships between protein complex subunits. This second systematic analysis is inspired by the structure of the dendrogram that summarizes similarities between signatures (Fig. 1). The dendrogram is automatically generated upon hierarchical clustering, with no prior knowledge of protein complex composition. Each signature is color-coded according to protein complex (Fig. 1, right), revealing that complex composition is a primary determinant of the dendrogram structure. This indicates that subunits derived from the same protein complex often have similar perturbation signatures. More detailed analysis reveals that this holds especially for complexes that consist of four or less different subunits (Fig. 3, Fig. S2 for all 32 complexes). Throughout the manuscript, for each set of mutants grouped together, the results for all genes with a significant change (p<0.05, FC>1.7) in any single mutant are depicted, rather than selecting for similarly behaving genes. As an example, the HIR complex is a histone chaperone involved in replication-independent deposition of histones throughout the cell-cycle (Green et al., 2005). HIR contributes to setting up a remodeler-resistant, repressive chromatin structure (Prochasson et al., 2005). Deletion of the different HIR subunits results in upregulation of the same set of 9–10 genes (Fig. 3C). Besides exemplifying data consistency, the HIR profiles show that each component is equally important. This agrees with the finding that all HIR subunits are required for interaction with the histone chaperone Asf1, essential for the role in nucleosome assembly (Green et al., 2005). Only a small group of genes show strong mRNA changes in the HIR signatures (Fig. 3C), likely reflecting the function of redundant mechanisms of histone deposition (Kaufman et al., 1998; Formosa et al., 2002). Nearly identical subunit signatures are also observed for the HDAC, SAS and Sir complexes (Fig. 3A, B, D), as well as for the Sum1/Rfm1/Hst1, CAF-1 and Ku complexes (Fig. S2). Although signatures derived from the same complex are generally similar, there are also numerous interesting exceptions, especially in complexes consisting of more than four subunits. This reveals the presence of auxiliary, peripheral and shared subunits, as well as the submodular organization of larger complexes (Fig. 3E–H).
Fig. 3. Structure-function analysis of chromatin regulatory complexes.
For each complex (AH), all genes with significant changes (p<0.05, FC>1.7) in any single mutant are shown, yellow indicating up-, blue indicating downregulation and black indicating no change, relative to the average wt and scaled as in Fig. 1. If a gene is not depicted, then no change in expression is observed outside the thresholds p<0.05, FC>1.7. Fig. S2 depicts all 32 complexes. (A) HDAC (B) SAS (C) HIR (D) Sir (E) SWR1, without Bdf1 as different isolates of bdf1Δ gave different results. The apparent lack of “downregulation” of the deleted gene in swc7Δ and swc3Δ (right hand box) is due to very low expression levels. All such mutants were PCR-verified for carrying the deletion (Supplemental Experimental Procedures). (F) Isw1a and Isw1b (G) Rpd3S and Rpd3L. Ume6 is not shown due to aneuploidy. (H) SAGA. The dendrogram was generated by hierarchical clustering. The grey line indicates that no probe for SUS1 is present on the array. (I) A network of SAGA subunit relationships generated using an edge-weighted, spring-embedded network algorithm (Shannon et al., 2003), revealing connections between submodules.
The presence of auxiliary subunits is exemplified by the SWR1 complex, involved in deposition of the histone variant H2A.Z (Kobor et al., 2004; Krogan et al., 2003; Mizuguchi et al., 2004). In contrast to the other subunits, deletion of SWC7 yields a profile similar to wt (Fig. 3E). This agrees with recent biochemical analyses showing that Swc7 is not required for complex integrity or H2A.Z replacement in vitro (Wu et al., 2009). Swc3 has also been classified as an auxiliary subunit, with little contribution to H2A.Z deposition (Wu et al., 2005, 2009). Deletion of SWC3 results in the same hallmark signature observed for the other components (Fig. 3E), indicating a role not captured by the in vitro analyses. Fig. S2 contains further examples of auxiliary function as well as examples of peripheral roles. The signature of a peripheral subunit is a subset of the effects observed in the others, indicating requirement in only one aspect of the function of a complex (Fig. 3G).
Other deviating subunit signatures result from subunits shared between different complexes. A shared subunit may functionally contribute to both complexes. This is illustrated by ISW1, a component common to the chromatin remodeling complexes Isw1a and Isw1b (Vary et al., 2003). ISW1 deletion yields a signature that largely combines the effects observed upon deletion of the two complexes that it is part of (Fig. 3F). Similarly, Rpd3S and Rpd3L are two related histone deacetylase complexes (Carrozza et al., 2005b). Deletion of the shared components RPD3 or SIN3, results in a signature that is mostly a combination of the profiles from the complex-specific subunits (Fig. 3G). A practical outcome is the identification of model genes suitable for distinguishing between the function of the two complexes. Deletions of the Rpd3L associated proteins, Ash1 and Cti6 (Carrozza et al., 2005a; Puig et al., 2004), show only a subset of the changes observed in other Rpd3L deletion mutants (Fig. 3G), exemplifying peripheral roles, required for particular genes only.
Larger complexes can often be subdivided into distinct submodules, also amenable to analysis by expression profiling. Analysis of the coregulator SAGA indicates four submodules (Fig. 3H and I). SAGA has previously been subdivided into a TATA-binding module, an acetylation module and a structural integrity module (Fig. 3I) (Sterner et al., 1999). The different signatures indicate that the recently characterized deubiquitination components form a functionally distinct submodule, with Sgf73 and Sgf29 connecting this module to the rest of the complex (Fig. 3I) (Rodriguez-Navarro, 2009; Bonnet et al., 2010). As with the other examples (Fig. 3, Fig. S2), this illustrates various ways in which the signatures can provide information on the function of individual subunits within complexes.
A correlation network of chromatin machinery interactions
The data was next analyzed to uncover global properties, focusing on how different complexes work together. The structure of the hierarchical cluster dendrogram is determined by similarities between subunits, but also by similarities between different complexes as a whole (Fig. 1). This is illustrated by the red part of the dendrogram, indicating relationships between Paf1C, Rad6/Bre1 and Set1C. To systematically capture such relationships, a correlation network was derived from the expression data (Fig. 4). In this network, each node represents a mutated component, color-coded according to protein complex. Edges connecting two nodes are drawn if the correlation between the profiles is higher than a specific threshold. Distances between nodes are dependent on the similarity of all connections. The threshold was set to capture those relationships resulting in identical, similar, or broadly overlapping signatures (Suppl. Experimental Procedures).
Fig. 4. A correlation network of perturbation signatures reveals chromatin machinery interactions.
Nodes represent mutated components, color-coded according to complex composition. Edges between nodes are drawn if the correlation between the two signatures is more than 0.4, capturing all identical, similar or broadly overlapping profiles. Fig S3A depicts a heatmap with all correlations. The network was generated with an edge-weighted, spring-embedded algorithm, so that the distances between nodes represent the similarity of connections. Interactions between several complexes are marked by the Roman numerals I (SAGA and SWI/SNF), II (Rad6/Bre1 and Set1C), III (Tup1 and HDAC) and IV (Paf1C and Set3C).
Three properties indicate network reliability. First, as expected, subunits of complexes are often grouped together. Second, the chromatin perturbation network is scale-free, distinct from a network that has a random topology (Fig. S3B) (Barabasi and Oltvai, 2004). Third, the network recapitulates many previously established interactions. For example, SAGA and SWI/SNF are tightly connected (indicated by I in Fig. 4), reflecting that SWI/SNF preferentially remodels nucleosomes acetylated by SAGA (Chandy et al., 2006). Similarly, there are many connections between Rad6/Bre1 and Set1C (Fig. 4, II). This agrees with the requirement for Rad6/Bre1 to ubiquitinate H2BK123 prior to H3K4 methylation by Set1C (Dover et al., 2002; Sun and Allis, 2002). The network indicates that like SAGA-SWI/SNF, this is also a general interaction, with functional consequences on the many genes shared between the corresponding signatures. The network also evaluates interactions between coregulators and chromatin complexes. Tup1-Ssn6 is a global coregulator (Smith and Johnson, 2000). Various repressive mechanisms have been proposed for Tup1 (Zhang and Reese, 2004b). The network indicates that repression through recruitment of the Hda1 deacetylase complex (HDAC) is most prevalent (Fig. 4, III). Evidence for another, less general mechanism, is presented below. The network reflects many individual interactions, facilitating evaluation of those that have previously been characterized on individual genes only.
The network also shows how different chromatin regulators are functionally connected and indicates that almost all regulators are linked to others. This shows that chromatin interactions do not form separate pathways that occur in isolation on different gene sets, but are interconnected. This agrees with the many genetic interactions found between chromatin regulators (Collins et al., 2007). Strictly linear, one-to-one pathway relationships between different regulators are sparse. Perturbation of linear signaling pathways, such as the mating pheromone or High Osmolarity Glycerol (HOG) MAPK cascades, do yield one-to-one relationships. In such pathways, essentially identical mRNA changes are observed upon perturbation of any component (Fig. S4). Here, only the interaction between the Set2 H3K36 methyltransferase and the Rpd3S histone deacetylase complexes (Carrozza et al., 2005b; Joshi and Struhl, 2005; Keogh et al., 2005) shows such a linear, one-to-one relationship (Fig. 5A). All other pathway relationships between chromatin regulators branch (Fig. 4). In other words and as is exemplified below (Fig. 5B and Fig. 6), the signature of one pathway component only partially overlaps or forms a subset of the signature of another pathway component. The varying degrees of overlap in the perturbation signatures are likely caused by the participation of complexes in different interactions, important for different genes.
Fig. 5. The Leo1 subunit of Paf1C contributes to Set3C recruitment.
(A) Identical signatures of Set2 and Rpd3S, indicating a linear pathway. (B) Hierarchical cluster diagram of expression signatures of Bur2, Paf1C, Rad6/Bre1C, COMPASS and Set3C. The genes are colored and scaled as in Fig. 1. (C) and (D) ChIP of Hos2-TAP in wt, paf1Δ, leo1Δ and cdc73Δ. PCR was carried out on RPS13 (C) and PYK1 (D). Numbered primer locations are shown schematically (top). The left panel shows PCR products and the right panel shows quantification of replicates, with error bars for standard deviation. Signals are normalized to input. (E) and (G) ChIP from wt and set3Δ with anti-H3 or anti-acetyl H4 on RPS13 (E) and PYK1 (G). Results are normalized to total H3 signal. (F) and (H) ChIP from wt or leo1Δ with anti-H3, anti-acetyl H4, anti-H3K4me3 or anti-H3K4me2 on RPS13 (F) and PYK1 (H). (I) Previously described (black arrows) and proposed (blue arrow) pathways connecting Bur1/2 and Set3C. The pathways branch at Paf1C and regulate Set3C in two ways: via dimethylation of H3K4 by COMPASS (black arrows) and via the proposed H3K4 methylation-independent mechanism (blue arrow).
Fig. 6. Specific subsets of genes are sensitive to particular interactions.
Previously reported interactions not represented in the network of global correlations, but showing significant overlap on specific subsets of genes. All genes with significant changes (p<0.05, FC>1.7) in any single mutant are shown, except for B. Genes are colored and scaled as in Fig. 1. (A) Tup1, Isw2C and HDAC. (B) RNR3 expression and overlap between Isw2C and Tup1. Here all genes are depicted with significant upregulation in tup1Δ and in any Isw2C signature. (C) SAS and SWR1C. Genes marked by an asterisk (*) are subtelomerically located (<50 kb from chromosome end).
The Leo1 subunit of Paf1C contributes to a distinct branch of Set3C recruitment
To investigate branching and interconnectivity, we first focused on the connections observed between the Paf1C, Rad6/Bre1C, COMPASS/Set1C and Set3C complexes (Fig. 4II and IV). A single chromatin interaction pathway can be assembled from previous studies. Paf1C associates with RNA polymerase II and regulates ubiquitination of H2BK123 by Rad6/Bre1C (Ng et al., 2003; Wood et al., 2003; Robzyk et al., 2000). These processes are dependent on the Bur1/2 cyclin dependent kinase complex (Wood et al., 2005). Subsequent to H2B ubiquitination, COMPASS/Set1 methylates H3K4 (Krogan et al., 2003, 2002; Dover et al., 2002). It has recently been shown that dimethylated H3K4 is bound by Set3C (Kim and Buratowski, 2009). These findings can be summarized in a linear pathway (Fig. 5I, black lines). Many aspects of this pathway are reflected in the network (Fig. 4, II and IV) and in the underlying perturbation signatures (Fig. 5B). BUR2 deletion clusters with Paf1C (Fig. 5B), highlighting its role in recruitment of Paf1C and regulation of Rad6/Bre1 (Wood et al., 2005). Going further down the pathway, the signatures of downstream complexes (e.g. Set1C, Set3C) are subsets of upstream ones (e.g. Paf1C, Rad6/Bre1, Fig 5B). Interestingly, one Paf1C subunit, rtf1Δ, clusters closer to the Rad6/Bre1complex than to members of its own complex, indicating that the sole role of Rtf1 is to promote Rad6/Bre1C function (Ng et al., 2003). Since Set3C recognizes H3K4 dimethylation by Set1C (Kim and Buratowski, 2009), and is furthest downstream, it is surprising that Set3C signatures are more similar to two Paf1C subunits (leo1Δ and cdc73Δ) than to Set1C (Fig. 5B). This is also reflected in the network (Fig. 4, IV). Similarity between leo1Δ and Set3C is even more surprising given that loss of the Paf1C subunits Rtf1, Ctr9, Paf1, all result in defective H3K4 dimethylation, required for Set3C recruitment, while loss of Leo1 has no effect on this mark (Xiao et al., 2005). Together, this suggests that the Leo1 subunit of Paf1C may also contribute to Set3C function in a manner distinct from the pathway through Set1C and Rad6/Bre1 (Fig. 5I, blue arrow).
A functional link between Paf1C and Set3C, separate from the previously established pathway, was tested by chromatin immunoprecipitation (ChIP, Fig. 5C–H). Loss of Leo1 or Cdc73 results in reduced recruitment of Set3C (Fig 5C, D). Set3C is a histone deacetylase and acetylation levels are increased in leo1Δ (Fig. 5F, H). To ensure that the reduced Set3C presence observed in leo1Δ is not simply a reflection of the role of Paf1C in recruitment of Set1C, H3K4 dimethylation patterns were determined. As reported previously (Xiao et al., 2005), leo1Δ has no effect on H3K4 dimethylation (Fig. 5F, H). This agrees with a contribution towards recruiting Set3C that is distinct from the Rad6/Bre1-Set1C pathway. PAF1 deletion leads to a larger decrease in Set3C recruitment (Fig 5C, D), suggesting possible loss of both modes of Set3C recruitment in paf1Δ. The data therefore support a model in which the pathways branch at the level of Paf1C. Paf1C contributes to Set3C recruitment in two ways (Fig. 5I), through the previously described pathway via H3K4 methylation and in a separate manner dependent on Leo1 (Fig. 5I, blue arrow). The latter was predicted by analysis of the chromatin perturbation signatures (Fig. 4, 5B) and perhaps reflects the role of Leo1 in RNA binding (Dermody and Buratowski, 2010). Besides illustrating branching and interconnectivity of a chromatin interaction pathway, this also exemplifies ways in which the signatures can be applied to investigate interaction pathways further.
Specific subsets of genes are sensitive to particular interactions
The threshold applied to generate the correlation network (Fig. 4) only identifies interactions with broad overlaps in the respective signatures. The branched network topology reveals that different combinations of factors are important for different gene sets. Functionally relevant interactions will not necessarily result in strongly overlapping signatures. We therefore also investigated previously reported interactions not represented in the network of global correlations. Besides functioning through the Hda1 deacetylase complex (Fig 4, III), the corepressor Tup1 has also been shown to interact with the nucleosome remodeler Isw2C on RNR3 (Zhang and Reese, 2004a, 2004b). Although the overall correlation is lower than the threshold applied to generate the network, tup1Δ and isw2Δ signatures do overlap significantly (Fig. 6A, p=9*10−7, hypergeometric test). In agreement with previous observations (Zhang and Reese, 2004a, 2004b), upregulation of RNR3 in isw2Δ is very weak (1.3-fold, p<1*10−6). The analysis indicates that the interaction between Tup1 and Isw2C, previously characterized on RNR3, has functional implications for a specific subset of genes. Compared to RNR3, several other genes show a much stronger response in isw2Δ (Fig. 6B). Besides providing a further illustration of branching (Fig. 6A), a practical outcome is therefore the ability to pinpoint model genes that are most sensitive to a particular interaction.
Regulation in the context of telomeric location
Many of the individual interactions captured in the global network represent histone code “writer” – “reader” relationships (Fig. 4). The results only loosely fit the histone code hypothesis however. There is a large amount of branching between different regulators. This is due to the partial nature of the overlaps between signatures of different complexes. Combinatorial effects are part of the original histone code theory (Strahl and Allis, 2000), but would only agree with the branched network topology if different histone marks also exhibit partially overlapping location. Analyses of histone marks across the yeast genome have not yet resulted in the discovery of gene set specific patterns (Pokholok et al., 2005; Rando and Chang, 2009). We therefore focused on recent proposals suggesting that binding of chromatin factors to specific marks is dependent on context (Berger, 2007; Lee et al., 2010; Campos and Reinberg, 2009).
The nature of such contexts is not well understood and the data offers the possibility to explore these further. One context may be formed by location close to telomeres. Telomere-proximal genes are silenced through the telomere position effect, mediated by the Sir proteins (Perrod and Gasser, 2003). Boundary factors, such as H4K16 acetylation and H2A.Z, function synergistically to prevent Sir spreading on a specific subset of genes located at subtelomeric regions (Shia et al., 2006). Overlapping function, but only in the specific context of subtelomeric location, is indeed recapitulated in the expression signatures of the SAS H4K16 acetyltransferase complex and SWR1-C (Fig. 6C).
To investigate this context further, all signatures were tested for significant enrichment of subtelomerically positioned genes (Suppl. Experimental Procedures). For subtelomerically enriched signatures, the average expression in a 1500 bp sliding window from the chromosome end is depicted (Fig. 7A). The contributions of SAS, SWR1C and Rpd3L (Ehrentraut et al., 2010), as well as of Sir proteins, are as expected. In addition to established factors, several others have subtelomerically enriched signatures (Fig. 7A), including the coregulator Mediator. Members of different Mediator submodules (Guglielmi et al., 2004) affect subtelomeric genes differently. Disruption of the Tail (med2Δ, med3Δ, med15Δ) results in upregulation of subtelomeric genes (p = 7*10−4, p = 6*10−5 and p = 2*10−15 respectively, chi-squared test), while disruption of the Middle (med1Δ, med5Δ, med9Δ) results in downregulation (Fig. 7A, p = 2*10−20, p = 5*10−13 and p = 2*10−10, respectively). Although possible explanations include indirect effects, a more direct mechanism is supported by binding of Med5 to telomeres (Esnault et al., 2008). To investigate these observations, the genome-wide location of Sir3 was assessed in two Mediator deletion mutants, med15Δ from the Tail and med1Δ from the Middle (Fig. 7B and Fig. S5). Loss of MED15 causes decreased Sir3 occupancy at several subtelomeric ends (Fig. 7B and Fig. S5). In contrast, loss of MED1 results in increased Sir3 occupancy, coincident with increased silencing (Fig. 7A). Sir mRNA and protein levels are not affected in the mutants (data not shown). In both mutants, Sir3 occupancy is changed, but no general spreading of Sir3 is observed. This suggests that Mediator does not form a boundary element similar to H4K16 acetylation, but instead influences the magnitude of Sir protein binding. The effect of Mediator mutants varies at individual telomeres, in agreement with varying degrees of Sir protein repression at different telomeres (Pryde and Louis, 1999). The changes on subtelomeric Sir3 occupancy are widespread, not just restricted to the vicinity of genes with changed expression in Mediator mutants (Fig. S5). It is therefore unlikely that the changes in Sir occupancy are indirectly reflecting the coregulatory role of Mediator in gene expression. These results agree with a report that Mediator binds directly to H4K16 deacetylated nucleosomes, competing with Sir for nucleosome binding (C. Gustafsson, under review). Together this indicates that Mediator has different roles in different genomic locations.
Fig. 7. Analysis of subtelomeric expression in chromatin mutants indicates that Mediator is a regulator of subtelomeric silencing.
(A) For subtelomerically enriched signatures (p<0.05, chi-squared test), the average expression in a 1500 bp sliding window (steps of 250 bp) is plotted according to the distance from the chromosome end. (B) Sir3 occupancy measured by ChIP-on-Chip in wt, med1Δ or med15Δ at the right arms of chromosome 2 and 3. The ORFs are positioned above the X-axis for Watson genes and below the axis for Crick genes. Sir3 occupancy for all chromosome ends is shown in Fig. S5.
Chromatin pathway topology and specificity
This study is the first comprehensive analysis of chromatin interaction pathways and the first whereby the functional consequences on gene expression are systematically analyzed. The majority of known chromatin regulators have been included. An important outcome is that chromatin interaction pathways show an interconnected and branched topology (Fig. 4). With the exception of Set2-Rpd3S (Fig. 5A), the functional consequences of perturbing different complexes do not support linear, one-to-one relationships between interacting factors. Instead, interacting partners show varying degrees of overlap in their perturbation signatures. Chromatin interaction pathways are therefore not distinct, separate structures for different groups of genes, but share many components. The partial overlaps between signatures show that not all genes are dependent on all components of a pathway and this study identifies the important interactions for many different gene sets.
This study underscores the complexity of regulation through chromatin. It is important to note that histone code “writer” – “reader” relationships are reflected in the network (Fig. 4), but not in an identical, rigid manner across all genes. While the locations of several marks are widespread (Pokholok et al., 2005), their removal results in effects that are much more specific. Such apparent discrepancies from a strict code are starting to be explained by context-dependent binding of chromatin factors (Berger, 2007; Lee et al., 2010; Campos and Reinberg, 2009). In these proposals, different effectors are capable of binding the same mark, but only do so dependent on the context. The result is different functional outcomes mediated by the same histone modification (Lee et al., 2010). Interestingly, the genome-wide binding patterns of effector proteins (Venters et al., 2011) are still much less specific than the perturbation signatures. This indicates that context-dependent binding of effectors is not sufficient to explain the specificity of regulation through chromatin. Other factors that may contribute to specificity include intrinsic properties of genes such as local nucleosome density and dynamics, the influence of neighbouring transcription units and the spatial location of genes. Functional redundancy may also contribute to specificity. Redundancy is not complete for the majority of factors, since 80% do show defects in gene expression. Redundancy may be partial however, with loss of one mechanism compensated by another, but only on a particular set of genes. This is supported by widespread negative synthetic genetic interactions between chromatin factors (Collins et al., 2007). The specificity observed here bodes well for the development of specific therapies based on drugs that target epigenetic factors. The data is useful for individual analyses of the 165 components included and is made available in a variety of formats. The topology of chromatin interaction pathways is revealed for the first time and the study provides a framework for determining which mechanisms are most important for achieving the specific effects observed for regulation through chromatin.
Experimental Procedures
All procedures are described in detail in the Supplemental Experimental Procedures.
Expression-profiling and deletion strains
Each mutant strain, BY4742 (Table S2), was profiled four times from two independently inoculated cultures and harvested in early mid-log phase in synthetic complete medium with 2% glucose. Sets of mutants were grown alongside wt cultures and processed in parallel. Dual-channel 70-mer oligonucleotide arrays were employed with a common reference wt RNA. All steps after RNA isolation were automated using robotic liquid handlers. These procedures were first optimized for accuracy (correct fold-change) and precision (reproducible result), using spiked-in RNA calibration (van Bakel and Holstege, 2004). After quality control, normalization and dye-bias correction (Margaritis et al., 2009), statistical analysis was performed for each mutant versus the collection of 200 wt cultures. The reported fold-change is an average of the 4 replicate mutant profiles versus the average of all wts. 58 genes showed stochastic changes in wt profiles (wt variable genes) and were excluded from all analyses. Incorrect strains from the collection (13%) as indicated by aneuploidy (6%), incorrect deletion (5%) or additional spurious mutation affecting the profile (2%) were remade and reprofiled (Table S2). Less than 1% of the wt profiles had more than three genes changing compared to the average wt as determined by the same criteria as for the mutants (p<0.05, FC>1.7) and after exclusion of the wt variable genes. This threshold was therefore applied to determine whether a mutant had a profile different from wt.
Mating
The tester strain was CSHL3-10 (trp1-239, ura3-52, mat a). Cells were allowed to mate on a filter for 45 min. Serial dilutions were plated on –TRP and –TRP –HIS –URA plates. After 2 days of growth, colonies were counted. The percentage of diploids versus haploids plus diploids was determined. Experiments were normalized to the total mating efficiency of that experiment and calculated relative to wt.
Chromatin Immunoprecipitation (Fig 5)
Chromatin immunoprecipitations were carried out as described (Kim and Buratowski, 2009), see Supplemental experimental procedures. The sequences of oligonucleotides used are listed in Table S3.
ChIP-on-chip (Fig. 6)
Genome-wide location analyses were performed with an adapted linear amplification method (van Bakel et al., 2008). Biotin labeled samples were hybridized to Affymetrix 1.0R S. cerevisiae microarrays, consisting of over 3.2 million probes. An adapted version of the Model-based Analysis of Tiling-arrays (MAT) algorithm (Johnson et al., 2006) was used to detect enriched regions (Schulze et al., 2009). The Sir3 signal was normalized against mock controls to eliminate background enrichment. MATscores were calculated for each probe using a 300 bp sliding window.
Supplementary Material
Highlights.
Loss of chromatin regulators result in highly specific effects on mRNA expression
These effects are much more specific than the location of regulators or of marks
The perturbation signatures identify a network of interactions between regulators
Chromatin interaction pathways are highly branched and interconnected
Acknowledgments
We thank V. Geli and T. van Welsen for strains, H.Th.M. Timmers, M. Vermeulen, F. van Leeuwen and J.A. Lenstra for critical reading. Supported by the Netherlands Bioinformatics Centre (NBIC) and the Netherlands Organization of Scientific Research (NWO), grants 016108607, 81702015, 05071057, 91106009, 021002035 (TLL), 86307007 (PK), 70057407 (JJB).
Footnotes
Database submission
ArrayExpress acc. nr. E-TABM-1074 (mutants) and E-TABM-773/E-TABM-984 (200 wt replicates); GEO acc. nr. GSE25909. The data is also available in flat-file and TreeView format from http://www.holstegelab.nl/publications/chromatin_regulators. Profiles of different culture media and mating types are included to facilitate comparison with other datasets.
References
- van Bakel H, Holstege FCP. In control: systematic assessment of microarray performance. EMBO Rep. 2004;5:964–969. doi: 10.1038/sj.embor.7400253. [DOI] [PMC free article] [PubMed] [Google Scholar]
- van Bakel H, van Werven FJ, Radonjic M, Brok MO, van Leenen D, Holstege FCP, Timmers HTM. Improved genome-wide localization by ChIP-chip using double-round T7 RNA polymerase-based amplification. Nucleic Acids Res. 2008;36:e21. doi: 10.1093/nar/gkm1144. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Barabasi AL, Oltvai ZN. Network biology: understanding the cell’s functional organization. Nat Rev Genet. 2004;5:101–13. doi: 10.1038/nrg1272. [DOI] [PubMed] [Google Scholar]
- Berger SL. The complex language of chromatin regulation during transcription. Nature. 2007;447:407–412. doi: 10.1038/nature05915. [DOI] [PubMed] [Google Scholar]
- Bonnet J, Wang Y, Spedale G, Atkinson RA, Romier C, Hamiche A, Pijnappel WWMP, Timmers HTM, Tora L, Devys D, et al. The structural plasticity of SCA7 domains defines their differential nucleosome-binding properties. EMBO Rep. 2010;11:612–618. doi: 10.1038/embor.2010.98. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Campos EI, Reinberg D. Histones: annotating chromatin. Annu Rev Genet. 2009;43:559–99. doi: 10.1146/annurev.genet.032608.103928. [DOI] [PubMed] [Google Scholar]
- Carrozza MJ, Florens L, Swanson SK, Shia W, Anderson S, Yates J, Washburn MP, Workman JL. Stable incorporation of sequence specific repressors Ash1 and Ume6 into the Rpd3L complex. Biochim Biophys Acta. 2005a;1731:77–87. doi: 10.1016/j.bbaexp.2005.09.005. discussion 75–76. [DOI] [PubMed] [Google Scholar]
- Carrozza MJ, Li B, Florens L, Suganuma T, Swanson SK, Lee KK, Shia W, Anderson S, Yates J, Washburn MP, et al. Histone H3 methylation by Set2 directs deacetylation of coding regions by Rpd3S to suppress spurious intragenic transcription. Cell. 2005b;123:581–592. doi: 10.1016/j.cell.2005.10.023. [DOI] [PubMed] [Google Scholar]
- Chandy M, Gutiérrez JL, Prochasson P, Workman JL. SWI/SNF displaces SAGA-acetylated nucleosomes. Eukaryotic Cell. 2006;5:1738–1747. doi: 10.1128/EC.00165-06. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Collins SR, Miller KM, Maas NL, Roguev A, Fillingham J, Chu CS, Schuldiner M, Gebbia M, Recht J, Shales M, et al. Functional dissection of protein complexes involved in yeast chromosome biology using a genetic interaction map. Nature. 2007;446:806–10. doi: 10.1038/nature05649. [DOI] [PubMed] [Google Scholar]
- Dermody JL, Buratowski S. Leo1 subunit of the yeast paf1 complex binds RNA and contributes to complex recruitment. J Biol Chem. 2010;285:33671–33679. doi: 10.1074/jbc.M110.140764. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Dover J, Schneider J, Tawiah-Boateng MA, Wood A, Dean K, Johnston M, Shilatifard A. Methylation of histone H3 by COMPASS requires ubiquitination of histone H2B by Rad6. J Biol Chem. 2002;277:28368–28371. doi: 10.1074/jbc.C200348200. [DOI] [PubMed] [Google Scholar]
- Ehrentraut S, Weber JM, Dybowski JN, Hoffmann D, Ehrenhofer-Murray AE. Rpd3-dependent boundary formation at telomeres by removal of Sir2 substrate. Proc Natl Acad Sci US A. 2010;107:5522–5527. doi: 10.1073/pnas.0909169107. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Esnault C, Ghavi-Helm Y, Brun S, Soutourina J, Van Berkum N, Boschiero C, Holstege F, Werner M. Mediator-dependent recruitment of TFIIH modules in preinitiation complex. Mol Cell. 2008;31:337–346. doi: 10.1016/j.molcel.2008.06.021. [DOI] [PubMed] [Google Scholar]
- Formosa T, Ruone S, Adams MD, Olsen AE, Eriksson P, Yu Y, Rhoades AR, Kaufman PD, Stillman DJ. Defects in SPT16 or POB3 (yFACT) in Saccharomyces cerevisiae cause dependence on the Hir/Hpc pathway: polymerase passage may degrade chromatin structure. Genetics. 2002;162:1557–1571. doi: 10.1093/genetics/162.4.1557. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Green EM, Antczak AJ, Bailey AO, Franco AA, Wu KJ, Yates JR, Kaufman PD. Replication-independent histone deposition by the HIR complex and Asf1. Curr Biol. 2005;15:2044–2049. doi: 10.1016/j.cub.2005.10.053. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Guglielmi B, van Berkum NL, Klapholz B, Bijma T, Boube M, Boschiero C, Bourbon H, Holstege FCP, Werner M. A high resolution protein interaction map of the yeast Mediator complex. Nucleic Acids Res. 2004;32:5379–5391. doi: 10.1093/nar/gkh878. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Holstege FC, Jennings EG, Wyrick JJ, Lee TI, Hengartner CJ, Green MR, Golub TR, Lander ES, Young RA. Dissecting the regulatory circuitry of a eukaryotic genome. Cell. 1998;95:717–728. doi: 10.1016/s0092-8674(00)81641-4. [DOI] [PubMed] [Google Scholar]
- Huang H, Kahana A, Gottschling DE, Prakash L, Liebman SW. The ubiquitin-conjugating enzyme Rad6 (Ubc2) is required for silencing in Saccharomyces cerevisiae. Mol Cell Biol. 1997;17:6693–6699. doi: 10.1128/mcb.17.11.6693. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Johnson WE, Li W, Meyer CA, Gottardo R, Carroll JS, Brown M, Liu XS. Model-based analysis of tiling-arrays for ChIP-chip. Proc Natl Acad Sci US A. 2006;103:12457–12462. doi: 10.1073/pnas.0601180103. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Joshi AA, Struhl K. Eaf3 chromodomain interaction with methylated H3-K36 links histone deacetylation to Pol II elongation. Mol Cell. 2005;20:971–978. doi: 10.1016/j.molcel.2005.11.021. [DOI] [PubMed] [Google Scholar]
- Kaufman PD, Cohen JL, Osley MA. Hir proteins are required for position-dependent gene silencing in Saccharomyces cerevisiae in the absence of chromatin assembly factor I. Mol Cell Biol. 1998;18:4793–4806. doi: 10.1128/mcb.18.8.4793. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Keogh M, Kurdistani SK, Morris SA, Ahn SH, Podolny V, Collins SR, Schuldiner M, Chin K, Punna T, Thompson NJ, et al. Cotranscriptional set2 methylation of histone H3 lysine 36 recruits a repressive Rpd3 complex. Cell. 2005;123:593–605. doi: 10.1016/j.cell.2005.10.025. [DOI] [PubMed] [Google Scholar]
- Kim T, Buratowski S. Dimethylation of H3K4 by Set1 recruits the Set3 histone deacetylase complex to 5′ transcribed regions. Cell. 2009;137:259–72. doi: 10.1016/j.cell.2009.02.045. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Kobor MS, Venkatasubrahmanyam S, Meneghini MD, Gin JW, Jennings JL, Link AJ, Madhani HD, Rine J. A protein complex containing the conserved Swi2/Snf2-related ATPase Swr1p deposits histone variant H2A into euchromatin. PLoS Biol. 2004;2:E131. doi: 10.1371/journal.pbio.0020131. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Krogan NJ, Dover J, Wood A, Schneider J, Heidt J, Boateng MA, Dean K, Ryan OW, Golshani A, Johnston M, et al. The Paf1 complex is required for histone H3 methylation by COMPASS and Dot1p: linking transcriptional elongation to histone methylation. Mol Cell. 2003;11:721–9. doi: 10.1016/s1097-2765(03)00091-1. [DOI] [PubMed] [Google Scholar]
- Krogan NJ, Dover J, Khorrami S, Greenblatt JF, Schneider J, Johnston M, Shilatifard A. COMPASS, a histone H3 (Lysine 4) methyltransferase required for telomeric silencing of gene expression. J Biol Chem. 2002;277:10753–10755. doi: 10.1074/jbc.C200023200. [DOI] [PubMed] [Google Scholar]
- Krogan NJ, Keogh M, Datta N, Sawa C, Ryan OW, Ding H, Haw RA, Pootoolal J, Tong A, Canadien V, et al. A Snf2 family ATPase complex required for recruitment of the histone H2A variant Htz1. Mol Cell. 2003;12:1565–1576. doi: 10.1016/s1097-2765(03)00497-0. [DOI] [PubMed] [Google Scholar]
- Lee J, Smith E, Shilatifard A. The language of histone crosstalk. Cell. 2010;142:682–685. doi: 10.1016/j.cell.2010.08.011. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Mallory MJ, Strich R. Ume1p represses meiotic gene transcription in Saccharomyces cerevisiae through interaction with the histone deacetylase Rpd3p. J Biol Chem. 2003;278:44727–44734. doi: 10.1074/jbc.M308632200. [DOI] [PubMed] [Google Scholar]
- Margaritis T, Lijnzaad P, van Leenen D, Bouwmeester D, Kemmeren P, van Hooff SR, Holstege FCP. Adaptable gene-specific dye bias correction for two-channel DNA microarrays. Mol Syst Biol. 2009;5:266. doi: 10.1038/msb.2009.21. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Mizuguchi G, Shen X, Landry J, Wu W, Sen S, Wu C. ATP-driven exchange of histone H2AZ variant catalyzed by SWR1 chromatin remodeling complex. Science. 2004;303:343–348. doi: 10.1126/science.1090701. [DOI] [PubMed] [Google Scholar]
- Narlikar GJ, Fan HY, Kingston RE. Cooperation between complexes that regulate chromatin structure and transcription. Cell. 2002;108:475–87. doi: 10.1016/s0092-8674(02)00654-2. [DOI] [PubMed] [Google Scholar]
- Ng HH, Dole S, Struhl K. The Rtf1 component of the Paf1 transcriptional elongation complex is required for ubiquitination of histone H2B. J Biol Chem. 2003;278:33625–8. doi: 10.1074/jbc.C300270200. [DOI] [PubMed] [Google Scholar]
- van de Peppel J, Kemmeren P, van Bakel H, Radonjic M, van Leenen D, Holstege FCP. Monitoring global messenger RNA changes in externally controlled microarray experiments. EMBO Rep. 2003;4:387–393. doi: 10.1038/sj.embor.embor798. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Perrod S, Gasser SM. Long-range silencing and position effects at telomeres and centromeres: parallels and differences. Cell Mol Life Sci. 2003;60:2303–2318. doi: 10.1007/s00018-003-3246-x. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Pokholok DK, Harbison CT, Levine S, Cole M, Hannett NM, Lee TI, Bell GW, Walker K, Rolfe PA, Herbolsheimer E, et al. Genome-wide map of nucleosome acetylation and methylation in yeast. Cell. 2005;122:517–27. doi: 10.1016/j.cell.2005.06.026. [DOI] [PubMed] [Google Scholar]
- Prochasson P, Florens L, Swanson SK, Washburn MP, Workman JL. The HIR corepressor complex binds to nucleosomes generating a distinct protein/DNA complex resistant to remodeling by SWI/SNF. Genes Dev. 2005;19:2534–2539. doi: 10.1101/gad.1341105. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Pryde FE, Louis EJ. Limitations of silencing at native yeast telomeres. EMBO J. 1999;18:2538–2550. doi: 10.1093/emboj/18.9.2538. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Puig S, Lau M, Thiele DJ. Cti6 is an Rpd3-Sin3 histone deacetylase-associated protein required for growth under iron-limiting conditions in Saccharomyces cerevisiae. J Biol Chem. 2004;279:30298–30306. doi: 10.1074/jbc.M313463200. [DOI] [PubMed] [Google Scholar]
- Rando OJ, Chang HY. Genome-wide views of chromatin structure. Annu Rev Biochem. 2009;78:245–71. doi: 10.1146/annurev.biochem.78.071107.134639. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Richmond TJ, Davey CA. The structure of DNA in the nucleosome core. Nature. 2003;423:145–150. doi: 10.1038/nature01595. [DOI] [PubMed] [Google Scholar]
- Rine J, Herskowitz I. Four genes responsible for a position effect on expression from HML and HMR in Saccharomyces cerevisiae. Genetics. 1987;116:9–22. doi: 10.1093/genetics/116.1.9. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Roberts SM, Winston F. SPT20/ADA5 encodes a novel protein functionally related to the TATA-binding protein and important for transcription in Saccharomyces cerevisiae. Mol Cell Biol. 1996;16:3206–3213. doi: 10.1128/mcb.16.6.3206. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Robzyk K, Recht J, Osley MA. Rad6-dependent ubiquitination of histone H2B in yeast. Science. 2000;287:501–504. doi: 10.1126/science.287.5452.501. [DOI] [PubMed] [Google Scholar]
- Rodriguez-Navarro S. Insights into SAGA function during gene expression. EMBO Rep. 2009;10:843–50. doi: 10.1038/embor.2009.168. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Schulze JM, Jackson J, Nakanishi S, Gardner JM, Hentrich T, Haug J, Johnston M, Jaspersen SL, Kobor MS, Shilatifard A. Linking cell cycle to histone modifications: SBF and H2B monoubiquitination machinery and cell-cycle regulation of H3K79 dimethylation. Mol Cell. 2009;35:626–641. doi: 10.1016/j.molcel.2009.07.017. [DOI] [PMC free article] [PubMed] [Google 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. doi: 10.1101/gr.1239303. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Shia WJ, Li B, Workman JL. SAS-mediated acetylation of histone H4 Lys 16 is required for H2A. incorporation at subtelomeric regions in Saccharomyces cerevisiae. Genes Dev. 2006;20:2507–12. doi: 10.1101/gad.1439206. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Smith RL, Johnson AD. Turning genes off by Ssn6-Tup1: a conserved system of transcriptional repression in eukaryotes. Trends Biochem Sci. 2000;25:325–330. doi: 10.1016/s0968-0004(00)01592-9. [DOI] [PubMed] [Google Scholar]
- Sterner DE, Grant PA, Roberts SM, Duggan LJ, Belotserkovskaya R, Pacella LA, Winston F, Workman JL, Berger SL. Functional organization of the yeast SAGA complex: distinct components involved in structural integrity, nucleosome acetylation, and TATA-binding protein interaction. Mol Cell Biol. 1999;19:86–98. doi: 10.1128/mcb.19.1.86. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Strahl BD, Allis CD. The language of covalent histone modifications. Nature. 2000;403:41–45. doi: 10.1038/47412. [DOI] [PubMed] [Google Scholar]
- Sun Z, Allis CD. Ubiquitination of histone H2B regulates H3 methylation and gene silencing in yeast. Nature. 2002;418:104–108. doi: 10.1038/nature00883. [DOI] [PubMed] [Google Scholar]
- Vary JC, Gangaraju VK, Qin J, Landel CC, Kooperberg C, Bartholomew B, Tsukiyama T. Yeast Isw1p forms two separable complexes in vivo. Mol Cell Biol. 2003;23:80–91. doi: 10.1128/MCB.23.1.80-91.2003. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Venkatasubrahmanyam S, Hwang WW, Meneghini MD, Tong AHY, Madhani HD. Genome-wide, as opposed to local, antisilencing is mediated redundantly by the euchromatic factors Set1 and H2A. Proc Natl Acad Sci US A. 2007;104:16609–16614. doi: 10.1073/pnas.0700914104. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Venters BJ, Wachi S, Mavrich TN, Andersen BE, Jena P, Sinnamon AJ, Jain P, Rolleri NS, Jiang C, Hemeryck-Walsh C. A Comprehensive Genomic Binding Map of Gene and Chromatin Regulatory Proteins in Saccharomyces. Molecular Cell. 2011;41:480–492. doi: 10.1016/j.molcel.2011.01.015. [DOI] [PMC free article] [PubMed] [Google Scholar]
- van Wageningen S, Kemmeren P, Lijnzaad P, Margaritis T, Benschop JJ, de Castro IJ, van Leenen D, Groot Koerkamp MJA, Ko CW, Miles AJ, et al. Functional overlap and regulatory links shape genetic interactions between signaling pathways. Cell. 2010;143:991–1004. doi: 10.1016/j.cell.2010.11.021. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Wood A, Schneider J, Dover J, Johnston M, Shilatifard A. The Paf1 complex is essential for histone monoubiquitination by the Rad6-Bre1 complex, which signals for histone methylation by COMPASS and Dot1p. J Biol Chem. 2003;278:34739–42. doi: 10.1074/jbc.C300269200. [DOI] [PubMed] [Google Scholar]
- Wood A, Schneider J, Dover J, Johnston M, Shilatifard A. The Bur1/Bur2 complex is required for histone H2B monoubiquitination by Rad6/Bre1 and histone methylation by COMPASS. Mol Cell. 2005;20:589–599. doi: 10.1016/j.molcel.2005.09.010. [DOI] [PubMed] [Google Scholar]
- Wu W, Alami S, Luk E, Wu C, Sen S, Mizuguchi G, Wei D, Wu C. Swc2 is a widely conserved H2AZ-binding module essential for ATP-dependent histone exchange. Nat Struct Mol Biol. 2005;12:1064–1071. doi: 10.1038/nsmb1023. [DOI] [PubMed] [Google Scholar]
- Wu W, Wu C, Ladurner A, Mizuguchi G, Wei D, Xiao H, Luk E, Ranjan A, Wu C. N terminus of Swr1 binds to histone H2AZ and provides a platform for subunit assembly in the chromatin remodeling complex. J Biol Chem. 2009;284:6200–6207. doi: 10.1074/jbc.M808830200. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Xiao T, Kao CF, Krogan NJ, Sun ZW, Greenblatt JF, Osley MA, Strahl BD. Histone H2B ubiquitylation is associated with elongating RNA polymerase II. Mol Cell Biol. 2005;25:637–51. doi: 10.1128/MCB.25.2.637-651.2005. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Xie J, Pierce M, Gailus-Durner V, Wagner M, Winter E, Vershon AK. Sum1 and Hst1 repress middle sporulation-specific gene expression during mitosis in Saccharomyces cerevisiae. EMBO J. 1999;18:6448–6454. doi: 10.1093/emboj/18.22.6448. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Zhang Z, Reese JC. Redundant mechanisms are used by Ssn6-Tup1 in repressing chromosomal gene transcription in Saccharomyces cerevisiae. J Biol Chem. 2004a;279:39240–39250. doi: 10.1074/jbc.M407159200. [DOI] [PubMed] [Google Scholar]
- Zhang Z, Reese JC. Ssn6-Tup1 requires the ISW2 complex to position nucleosomes in Saccharomyces cerevisiae. EMBO J. 2004b;23:2246–2257. doi: 10.1038/sj.emboj.7600227. [DOI] [PMC free article] [PubMed] [Google Scholar]
Associated Data
This section collects any data citations, data availability statements, or supplementary materials included in this article.